Changeset 4170 for trunk/LMDZ.COMMON/libf/evolution/layered_deposits.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/layered_deposits.F90 (modified) (5 diffs)

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

@@ -33 +33 @@
 logical(k4), protected, private :: impose_dust_ratio            ! Impose dust-to-ice ratio instead of dust tendency (see 'dust2ice_ratio')
 real(dp),    protected, private :: dust2ice_ratio               ! Dust-to-ice ratio when ice condenses (10% by default)
-real(dp),    protected, private :: d_dust                       ! Tendency of dust [kg.m-2.y-1]
+real(dp),    protected, private :: d_dust                       ! Tendency of dust [kg/m2/y]
 real(dp),    parameter, private :: dry_lag_porosity  = 0.2_dp   ! Porosity of dust lag
 real(dp),    parameter, private :: wet_lag_porosity  = 0.15_dp  ! Porosity of water ice lag
@@ -85 +85 @@
 !
 ! DESCRIPTION
-!     Setter for 'layered_deposits' configuration parameters.
+!     Setter for "layered_deposits" configuration parameters.
 !
 ! AUTHORS & DATE
@@ -749 +749 @@
             h2o_ice(ig,islope) = layerings_map(ig,islope)%top%h_h2oice*rho_h2oice
         else
-            call subsurface_ice_layering(layerings_map(ig,islope),h2oice_depth(ig,islope),h2o_ice(ig,islope),co2_ice(ig,islope))
+            call get_ssice_layering(layerings_map(ig,islope),h2oice_depth(ig,islope),h2o_ice(ig,islope),co2_ice(ig,islope))
         end if
     end do
@@ -1034 +1034 @@
 
 !=======================================================================
-SUBROUTINE subsurface_ice_layering(this,h_toplag,h2o_ice,co2_ice)
-!-----------------------------------------------------------------------
-! NAME
-!     subsurface_ice_layering
+SUBROUTINE get_ssice_layering(this,h_toplag,h2o_ice,co2_ice)
+!-----------------------------------------------------------------------
+! NAME
+!     get_ssice_layering
 !
 ! DESCRIPTION
@@ -1095 +1095 @@
 end if
 
-END SUBROUTINE subsurface_ice_layering
+END SUBROUTINE get_ssice_layering
 !=======================================================================