Changeset 3707

Timestamp:

Apr 2, 2025, 9:23:10 AM (3 months ago)

Author:

evos

Message:

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.

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (1 diff)
• recomp_tend_mod.F90 (modified) (1 diff)

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

@@ -630 +630 @@
 == 2/4/2025 == EV
 updating the h2o and co2 ice values that will be used by the PCM according to the top layer in the stratum
-
+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.
+
+

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

@@ -1014 +1014 @@
 !------------------------
     call recomp_tend_co2(ngrid,nslope,timelen,d_co2ice,d_co2ice_ini,co2_ice,emis,vmr_co2_PCM,vmr_co2_PEM_phys,ps_timeseries,ps_avg_global_old,ps_avg_global_new)
-    !write(*,*) "> Updating the H2O sub-surface ice tendency due to lag layer"
-    !do ig = 1,ngrid
-    !    do islope = 1,nslope
-    !        call recomp_tend_h2o(icetable_depth_old(ig,islope),icetable_depth(ig,islope),tsurf_avg(ig,islope),tsoil_PEM_timeseries_old(ig,:,islope,:),tsoil_PEM_timeseries(ig,:,islope,:),d_h2oice(ig,islope))
-    !    enddo
-    !enddo
+    write(*,*) "> Updating the H2O sub-surface ice tendency due to lag layer"
+    do ig = 1,ngrid
+        do islope = 1,nslope
+            call recomp_tend_h2o(icetable_depth_old(ig,islope),icetable_depth(ig,islope),tsurf_avg(ig,islope),tsoil_PEM_timeseries_old(ig,:,islope,:),tsoil_PEM_timeseries(ig,:,islope,:),d_h2oice(ig,islope))
+        enddo
+    enddo
     if (soil_pem) deallocate(icetable_depth_old,tsoil_PEM_timeseries_old)
     deallocate(vmr_co2_PEM_phys)

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

@@ -106 +106 @@
 
 ! Flux correction (decrease)
-d_h2oice = d_h2oice*R_dec*hum_dec
+d_h2oice = d_h2oice/R_dec/hum_dec
 
 END SUBROUTINE recomp_tend_h2o