Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F =================================================================== --- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F (revision 2942) +++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F (revision 2943) @@ -1,4 +1,4 @@ SUBROUTINE lect_start_archive(ngrid,nlayer,nqtot, - & date,tsurf,tsoil,albedo,emis,q2, + & date,tsurf,tsoil,inertiesoil,albedo,emis,q2, & t,ucov,vcov,ps,h,phisold_newgrid, & q,qsurf,tauscaling,totcloudfrac,surfith,nid, @@ -49,4 +49,5 @@ REAL,INTENT(OUT) :: tsurf(ngrid,nslope) ! surface temperature REAL,INTENT(OUT) :: tsoil(ngrid,nsoilmx,nslope) ! soil temperature + REAL,INTENT(OUT) :: inertiesoil(ngrid,nsoilmx,nslope) ! soil thermal inertia REAL,INTENT(OUT) :: albedo(ngrid,2,nslope) ! surface albedo REAL,INTENT(OUT) :: emis(ngrid,nslope) ! ground emissivity @@ -218,4 +219,5 @@ integer :: ndims integer, dimension(:), allocatable :: dimids + integer :: islope c======================================================================= @@ -1354,4 +1356,8 @@ & inertiedatS,inertiedat) + + do islope = 1,nslope + inertiesoil(:,:,islope) = inertiedat(:,:) + enddo c----------------------------------------------------------------------- c 6.2.2 Soil temperature Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F =================================================================== --- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (revision 2942) +++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (revision 2943) @@ -437,5 +437,5 @@ write(*,*) 'Reading file START_ARCHIVE' CALL lect_start_archive(ngridmx,llm,nqtot, - & date,tsurf,tsoil,albedo,emis,q2, + & date,tsurf,tsoil,inertiesoil,albedo,emis,q2, & t,ucov,vcov,ps,teta,phisold_newgrid,q,qsurf, & tauscaling,totcloudfrac,surfith,nid,watercap) @@ -1984,6 +1984,6 @@ call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot, & dtphys,hour_ini, - & tsurf,tsoil,albedo,emis,q2,qsurf,tauscaling, - & totcloudfrac,wstar,watercap) + & tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf, + & tauscaling,totcloudfrac,wstar,watercap) c======================================================================= Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F =================================================================== --- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F (revision 2942) +++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F (revision 2943) @@ -66,4 +66,5 @@ REAL,ALLOCATABLE :: tsurf(:,:) ! Surface temperature REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature + REAL,ALLOCATABLE :: inertiesoil(:,:,:) ! Soil thermal inertia (!= inertiedat which is for present day climate) REAL,ALLOCATABLE :: watercap(:,:) ! h2o ice layer REAL :: tauscaling(ngridmx) ! dust conversion factor @@ -83,5 +84,6 @@ REAL,ALLOCATABLE :: tsurfS(:,:) REAL,ALLOCATABLE :: tsoilS(:,:,:) - REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia + REAL,ALLOCATABLE :: inertiesoilS(:,:,:)! Variable Soil Thermal Inertia (obtained from PEM) + REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate) REAL,ALLOCATABLE :: watercapS(:,:) REAL :: tauscalingS(ip1jmp1) @@ -165,4 +167,5 @@ allocate(tsurf(ngridmx,nslope)) allocate(tsoil(ngridmx,nsoilmx,nslope)) + allocate(inertiesoil(ngridmx,nsoilmx,nslope)) allocate(watercap(ngridmx,nslope)) allocate(emis(ngridmx,nslope)) @@ -172,4 +175,5 @@ allocate(tsurfS(ip1jmp1,nslope)) allocate(tsoilS(ip1jmp1,nsoilmx,nslope)) + allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope)) allocate(watercapS(ip1jmp1,nslope)) allocate(emisS(ip1jmp1,nslope)) @@ -264,4 +268,5 @@ c watercap --> watercapS c tsoil --> tsoilS +c inertiesoil --> inertiesoilS c inertiedat --> ithS c emis --> emisS @@ -281,4 +286,6 @@ call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope), & tsoilS(:,:,islope)) + call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiesoil(:,:,islope), + & inertiesoil(:,:,islope)) ! Note: thermal inertia "inertiedat" is in comsoil.h call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis(:,islope), @@ -451,7 +458,9 @@ call write_archive(nid,ntime,'tsoil','Soil temperature', & 'K',-3,tsoilS(:,:,:)) - ! Write soil thermal inertia - call write_archive(nid,ntime,'inertiedat', + call write_archive(nid,ntime,'inertiesoil','Soil TI', + & 'J.s-1/2.m-2.K-1',-3,inertiesoilS(:,:,:)) +! Write soil thermal inertia for current climate + call write_archive(nid,ntime,'inertiedat - present day TI', & 'Soil thermal inertia', & 'J.s-1/2.m-2.K-1',-3,ithS) Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive_SSO.F =================================================================== --- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive_SSO.F (revision 2942) +++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive_SSO.F (revision 2943) @@ -69,4 +69,5 @@ REAL,ALLOCATABLE :: tsurf(:,:) ! Surface temperature REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature + REAL,ALLOCATABLE :: inertiesoil(:,:,:) ! Soil thermal inertia (!= inertiedat which is for present day climate) REAL,ALLOCATABLE :: watercap(:,:) ! h2o ice layer REAL tauscaling(ngridmx) ! dust conversion factor @@ -89,5 +90,6 @@ REAL,ALLOCATABLE :: tsurfS(:,:) REAL,ALLOCATABLE :: tsoilS(:,:,:) - REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia + REAL,ALLOCATABLE :: inertiesoilS(:,:,:)! Variable Soil Thermal Inertia (obtained from PEM) + REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate) REAL,ALLOCATABLE :: watercapS(:,:) REAL tauscalingS(ip1jmp1) @@ -183,4 +185,5 @@ allocate(tsurf(ngridmx,nslope)) allocate(tsoil(ngridmx,nsoilmx,nslope)) + allocate(inertiesoil(ngridmx,nsoilmx,nslope)) allocate(watercap(ngridmx,nslope)) allocate(emis(ngridmx,nslope)) @@ -190,4 +193,5 @@ allocate(tsurfS(ip1jmp1,nslope)) allocate(tsoilS(ip1jmp1,nsoilmx,nslope)) + allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope)) allocate(watercapS(ip1jmp1,nslope)) allocate(emisS(ip1jmp1,nslope)) @@ -281,4 +285,5 @@ c tsurf --> tsurfS c tsoil --> tsoilS +c inertiesoil ---> inertiesoilS c emis --> emisS c q2 --> q2S @@ -296,4 +301,6 @@ call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope), & tsoilS(:,:,islope)) + call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiesoil(:,:,islope), + & inertiesoil(:,:,islope)) ! Note: thermal inertia "inertiedat" is in comsoil.h call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis(:,islope), @@ -490,6 +497,8 @@ call write_archive(nid,ntime,'tsoil','Soil temperature', & 'K',-3,tsoilS(:,:,:)) - ! Write soil thermal inertia + call write_archive(nid,ntime,'inertiesoil','Soil TI', + & 'J.s-1/2.m-2.K-1',-3,inertiesoilS(:,:,:)) +! Write soil thermal inertia for current climate call write_archive(nid,ntime,'inertiedat', & 'Soil thermal inertia', Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F =================================================================== --- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F (revision 2942) +++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/write_archive.F (revision 2943) @@ -81,5 +81,5 @@ ! Creation de la variable si elle n'existait pas - + write (*,*) "=====================" write (*,*) "creation de ",nom @@ -224,5 +224,5 @@ & nom.eq."ps" .or. nom.eq."q2surf" .or. nom.eq."ZMEA" .or. & nom.eq."ZSTD" .or. nom.eq."ZSIG" .or. nom.eq."ZTHE" .or. - & nom.eq."ZGAM" .or. nom.eq."albedodat" .or. + & nom.eq."ZGAM" .or. nom.eq."albedodat" .or. & nom.eq."z0" .or. nom.eq."summit" .or. nom.eq."hmons" & .or. nom.eq."base") then @@ -238,5 +238,4 @@ write (*,*) "=====================" write (*,*) "creation de ",nom - call def_var(nid,nom,titre,unite,3,id,varid,ierr) @@ -276,5 +275,4 @@ write (*,*) "=====================" write (*,*) "creation de ",nom - call def_var(nid,nom,titre,unite,4,id,varid,ierr) @@ -323,5 +321,4 @@ write (*,*) "=====================" write (*,*) "creation de ",nom - call def_var(nid,nom,titre,unite,1,id,varid,ierr)