Changeset 3172 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Jan 8, 2024, 12:12:09 PM (11 months ago)

Author:

llange

Message:

PEM
Fixing bug when rewritting the startfi.nc for the PCM: fluxgeo, read in the run_PEM.def, was not written in the startfi.

LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (2 diffs)
• pemredem.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -192 +192 @@
 == 05/01/2024 == LL
 Add the possibily to output the soil fields during a PEM run in a dedicated diagsoilpem
+
+== 08/01/2024 == LL
+Fixing bug when rewritting the startfi.nc for the PCM: fluxgeo, read in the run_PEM.def,  was not written in the startfi.
+

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -71 +71 @@
 
 #ifndef CPP_STD
-    use comsoil_h,          only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil, nqsoil, qsoil
+    use comsoil_h,          only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil, flux_geo, nqsoil, qsoil
     use surfdat_h,          only: tsurf, emis, qsurf, watercap, ini_surfdat_h, &
                                   albedodat, zmea, zstd, zsig, zgam, zthe,     &
@@ -984 +984 @@
     call interpol_TI_PEM2PCM(ngrid,nslope,nsoilmx_PEM,nsoilmx,TI_PEM,inertiesoil)
     tsoil = tsoil_phys_PEM_timeseries(:,1:nsoilmx,:,timelen)
+#ifndef CPP_STD
+    flux_geo(:,:) = fluxgeo
+#endif
 endif
 

trunk/LMDZ.COMMON/libf/evolution/pemredem.F90

@@ -19 +19 @@
 
 ! create physics restart file and write time-independent variables
-use comsoil_h_PEM,     only: soil_pem,mlayer_PEM,fluxgeo,inertiedat_PEM
+use comsoil_h_PEM,     only: soil_pem,mlayer_PEM,inertiedat_PEM
 use iostart_PEM,       only: open_restartphy, close_restartphy, put_var, put_field, length
 use mod_grid_phy_lmdz, only: klon_glo
@@ -73 +73 @@
 ! write time-dependent variable to restart file
 use iostart_PEM,   only: open_restartphy, close_restartphy, put_var, put_field
-use comsoil_h_PEM, only: mlayer_PEM,fluxgeo, inertiedat_PEM, soil_pem
+use comsoil_h_PEM, only: mlayer_PEM, inertiedat_PEM, soil_pem
 use time_evol_mod, only: year_bp_ini, convert_years