Changeset 3635

Timestamp:

Feb 20, 2025, 10:31:38 AM (5 months ago)

Author:

tbertrand

Message:

Pluto: small fix to correctly read the soil thermal inertia from the starts
TB

Location:

trunk/LMDZ.PLUTO/libf

Files:

• dynphy_lonlat/phypluto/start2archive.F (modified) (4 diffs)
• phypluto/phyetat0_mod.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F

@@ -71 +71 @@
       REAL tsurf(ngridmx)        ! Surface temperature
       REAL,ALLOCATABLE :: tsoil(:,:) ! Soil temperature
+      REAL,ALLOCATABLE :: inertiedat_tmp(:,:) ! Soil temperature
       REAL n2ice(ngridmx)        ! N2 ice layer
       REAL q2(ngridmx,llm+1)
@@ -203 +204 @@
       ! allocate arrays of nsoilmx size
       allocate(tsoil(ngridmx,nsoilmx))
+      allocate(inertiedat_tmp(ngridmx,nsoilmx))
       allocate(tsoilS(ip1jmp1,nsoilmx))
       allocate(ithS(ip1jmp1,nsoilmx))
@@ -234 +236 @@
       CALL phyetat0(.true.,ngridmx,llm,fichnom,0,Lmodif,nsoilmx,nqtot,
      .      day_ini_fi,timefi,
-     .      tsurf,tsoil,emis,q2,qsurf,inertiedat)
+     .      tsurf,tsoil,emis,q2,qsurf,inertiedat_tmp)
 !       change FF 05/2011
 c     .       cloudfrac,totalcloudfrac,hice,
 !       change BC 05/2014
 c     .       rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
-
-
 
 
@@ -338 +338 @@
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil,tsoilS)
       ! Note: thermal inertia "inertiedat" is in comsoil.h
-      call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiedat,ithS)
+      call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiedat_tmp,ithS)
 
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,zmea,zmeaS)

trunk/LMDZ.PLUTO/libf/phypluto/phyetat0_mod.F90

@@ -23 +23 @@
   use write_field_phy, only: Writefield_phy
 !!
-  use comsoil_h, only: nsoilmx
+  use comsoil_h, only: nsoilmx,inertiedat
   use tabfi_mod, only: tabfi
   USE tracer_h, ONLY: noms
@@ -273 +273 @@
     call soil_settings(nid_start,ngrid,nsoil,tsurf,tsoil,indextime)
   endif
+  therm_inertia(:,:) = inertiedat(:,:)
 else
   thermal_inertia_if_no_startfi=400 ! default value
@@ -279 +280 @@
     !AF24
 endif ! of if (startphy_file)
-
 !
 ! close file: