Changeset 2943 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

Timestamp:

Apr 17, 2023, 12:07:33 PM (19 months ago)

Author:

llange

Message:

Mars PCM
Following r-2942: Fix a bug in newstart when rewriting inertiesoil. Inertiesoil is now also managed in startarchive.
When using startarchive or newstart, inertiesoil is set to inertiedat.
LL

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -66 +66 @@
       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 +84 @@
       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 +167 @@
       allocate(tsurf(ngridmx,nslope))
       allocate(tsoil(ngridmx,nsoilmx,nslope))
+      allocate(inertiesoil(ngridmx,nsoilmx,nslope))
       allocate(watercap(ngridmx,nslope))
       allocate(emis(ngridmx,nslope))
@@ -172 +175 @@
       allocate(tsurfS(ip1jmp1,nslope))
       allocate(tsoilS(ip1jmp1,nsoilmx,nslope))
+      allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
       allocate(watercapS(ip1jmp1,nslope))
       allocate(emisS(ip1jmp1,nslope))
@@ -264 +268 @@
 c watercap --> watercapS
 c tsoil --> tsoilS
+c inertiesoil --> inertiesoilS
 c inertiedat --> ithS
 c emis --> emisS
@@ -281 +286 @@
       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 +458 @@
       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)