Changeset 3711

Timestamp:

Apr 3, 2025, 10:38:43 AM (3 months ago)

Author:

evos

Message:

adding the growth of the lag layer
changing the stopping critiria to take into account the top startum layer
removing the case where if we have 0 ice the PEM stops.

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

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

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

@@ -632 +632 @@
 updating h2o ice recomp tendency to reduce the ice table loss rate due increase of depth and a decrease in the humidity around the ice.
 
-
+== 3/4/2025 == EV
+adding the growth of the lag layer
+changing the stopping critiria to take into account the top startum layer
+removing the case where if we have 0 ice the PEM stops.
+
+
+

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

@@ -1022 +1022 @@
     if (soil_pem) deallocate(icetable_depth_old,tsoil_PEM_timeseries_old)
     deallocate(vmr_co2_PEM_phys)
+    write(*,*) "> Updating the H2O sub-surface ice depth"
+    do ig = 1,ngrid
+        do islope = 1,nslope
+           if (icetable_depth(ig,islope) > 0) then
+               icetable_depth(ig,islope)=icetable_depth(ig,islope) + d_h2oice(ig,islope)*0.01 !!! 0.01 is for 1 percent dust, needs to be updated !!!
+            endif
+        enddo
+    enddo
+
 
 !------------------------
@@ -1028 +1037 @@
 !------------------------
     write(*,*) "> Checking the stopping criteria"
+    h2o_ice=sum(stratif(ig,islope)%top%h2oice_volfrac*stratif(ig,islope)%top%thickness*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.))
+    co2_ice=sum(stratif(ig,islope)%top%h2oice_volfrac*stratif(ig,islope)%top%thickness*subslope_dist(ig,:)/cos(pi*def_slope_mean(:)/180.))
     call stopping_crit_h2o_ice(cell_area,h2oice_ini_surf,is_h2oice_sublim_ini,h2o_ice,stopPEM,ngrid)
     call stopping_crit_co2(cell_area,co2ice_sublim_surf_ini,is_co2ice_sublim_ini,co2_ice,stopPEM,ngrid,ps_avg_global_ini,ps_avg_global_new,nslope)

trunk/LMDZ.COMMON/libf/evolution/stopping_crit_mod.F90 ¶

@@ -69 +69 @@
 endif
 
-if (abs(h2oice_ini_surf) < 1.e-5) stopPEM = 0
+!if (abs(h2oice_ini_surf) < 1.e-5) stopPEM = 0
 
 END SUBROUTINE stopping_crit_h2o_ice
@@ -133 +133 @@
 endif
 
-if (abs(co2ice_sublim_surf_ini) < 1.e-5) stopPEM = 0
+!if (abs(co2ice_sublim_surf_ini) < 1.e-5) stopPEM = 0
 
 if (ps_avg_global < ps_avg_global_ini*(1. - ps_criterion)) then