Changeset 4068 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Feb 12, 2026, 12:32:15 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:
Following r4065, some safeguards forgotten in "io_netcdf.F90" + 'tsurf' is not needed anymore in the "start1D.txt" + few forgotten small updates.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (4 diffs)

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

@@ -265 +265 @@
     allocate(character(8) :: num)
 end if
-allocate(delta_icetable(ngrid),icetable_depth_old(ngrid,nslope),is_co2ice_disappeared(ngrid,nslope))
+allocate(delta_icetable(ngrid),icetable_depth_old(ngrid,nslope),is_co2ice_disappeared(ngrid,nslope),tsoil_ts_old(ngrid,nsoil,nslope,nday))
 is_co2ice_disappeared(:,:) = .false.
 delta_icetable(:) = 0._dp
@@ -290 +290 @@
         ! Conversion to surface ice
         call layering2surfice(layerings_map,h2o_ice,co2_ice,h2oice_depth)
-
     else
         zlag(:,:) = 0._dp
@@ -324 +323 @@
 
         ! Evolve soil temperature
-        allocate(tsoil_ts_old(ngrid,nsoil,nslope,nday))
         call evolve_soil_temp(tsoil_avg,tsurf_avg,tsoil_ts,tsoil_ts_old,h2o_soildensity_avg)
 
@@ -491 +489 @@
 
 ! Write the "start.nc"
-call write_restart(ps4PCM,pa4PCM,preff4PCM,tsurf4PCM,q4PCM,teta4PCM,air_mass4PCM)
+call write_restart(ps4PCM,pa4PCM,preff4PCM,q4PCM,teta4PCM,air_mass4PCM)
 
 ! Update the duration information of the simulation