Changeset 3635 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto

Timestamp:

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

Author:

tbertrand

Message:

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

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/start2archive.F (modified) (4 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)