Changeset 4170 for trunk/LMDZ.COMMON/libf/evolution/clim_state_rec.F90

Timestamp:

Apr 3, 2026, 4:34:51 PM (7 days ago)

Author:

jbclement

Message:

PEM:

• Deletion of 'flux_ssice' ('zqdsdif_tot') from the "startfi.nc" as it is not needed.
• Using the yearly average flux for the sublimating subsurface ice instead of the value at last timestep.
• Keeping a clear separation between the subsurface ice flux and the surface ice tendency.
• Making sure that subsurface ice depth is well given to the PCM at the end of the PEM (ice table VS layering).

JBC

File:

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

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

@@ -180 +180 @@
 
 !=======================================================================
-SUBROUTINE write_restartfi(is_h2o_perice,h2o_ice4PCM,co2_ice4PCM,tsurf4PCM,tsoil4PCM,inertiesoil4PCM,albedo4PCM,emissivity4PCM,flux_geo4PCM)
+SUBROUTINE write_restartfi(is_h2o_perice,h2o_ice4PCM,co2_ice4PCM,tsurf4PCM,tsoil4PCM,inertiesoil4PCM,albedo4PCM,emissivity4PCM,flux_geo4PCM,h2oice_depth4PCM)
 !-----------------------------------------------------------------------
 ! NAME
@@ -208 +208 @@
 ! ARGUMENTS
 ! ---------
-real(dp),    dimension(:,:),   intent(in) :: h2o_ice4PCM, co2_ice4PCM, albedo4PCM, emissivity4PCM, tsurf4PCM, flux_geo4PCM
+real(dp),    dimension(:,:),   intent(in) :: h2o_ice4PCM, co2_ice4PCM, albedo4PCM, emissivity4PCM, tsurf4PCM, flux_geo4PCM, h2oice_depth4PCM
 real(dp),    dimension(:,:,:), intent(in) :: tsoil4PCM, inertiesoil4PCM
 logical(k4), dimension(:),     intent(in) :: is_h2o_perice
@@ -240 +240 @@
 
 ! Orbital parameters (controle)
-controle(18) = obliquity ! Obliquity
+controle(18) = obliquity  ! Obliquity
 controle(15) = perihelion ! Perihelion
-controle(16) = aphelion ! Aphelion
-controle(17) = date_peri ! Date of perihelion
+controle(16) = aphelion   ! Aphelion
+controle(17) = date_peri  ! Date of perihelion
 call put_var_nc('controle',controle)
 deallocate(controle)
@@ -259 +259 @@
 call put_var_nc('emis',emissivity4PCM,1)
 call put_var_nc('flux_geo',flux_geo4PCM,1)
-
-! h2oice_depth
-! lag_co2_ice
-! zdqsdif_ssi_tot
-! d_coef
-
+call put_var_nc('h2oice_depth',h2oice_depth4PCM)
+
+! Close
 call close_nc('re'//startfi_name)