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

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.

File:

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

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)