Changeset 2812

Timestamp:

Nov 3, 2022, 2:28:56 PM (2 years ago)

Author:

llange

Message:

PEM
Update CO2 glaciers evolution following a parametrization made by LL & Anna Grau (LPG).
Works only for nslope = 7. Hardcoded now, will be improved in future commits (need to continue theoretical work with Anna).
LL

File:

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

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

@@ -230 +230 @@
       REAL, dimension(:,:),allocatable  :: tendencies_h2o_ice_phys_slope   ! physical point field : Tendency of evolution of perenial co2 ice over a year
 
-      REAL, PARAMETER :: co2_hmax = 10                    ! Maximum height  for CO2 deposit on slopes (m)
+      REAL,SAVE,ALLOCATABLE,DIMENSION(:) ::  co2_hmax                          ! Maximum height  for CO2 deposit on slopes (m)
+
       REAL, PARAMETER :: rho_co2 = 1600           ! CO2 ice density (kg/m^3)
       INTEGER :: iaval                            ! Index of the neighboord slope ()
@@ -460 +461 @@
        PRINT*,'Flat slope for islope = ',iflat
        PRINT*,'corresponding criterium = ',def_slope_mean(iflat)
+! CO2 max thickness (for glaciers flows)
+
+       allocate(co2_hmax(nslope)
+       if(nslope.eq.7) ! ugly way to implement that ...
+!        CF documentation that explain how the values are computed
+         co2_hmax(1) = 1.5
+         co2_hmax(7) = co2_max(1)
+         co2_hmax(2) = 2.4
+         co2_hmax(6) = co2_hmax(2)
+         co2_hmax(3) = 5.6
+         co2_hmax(5) = co2_hmax(3)
+         co2_hmax(4) = 1000000.
+       endif
+
+
+
+
 
        flag_co2flow(:,:) = 0.     
@@ -1018 +1036 @@
           IF(islope.ne.iflat) THEN ! ice can be infinite on flat ground
 ! First: check that CO2 ice must flow (excess of ice on the slope), ice can accumulate on flat ground
-            IF(co2ice_slope(ig,islope).ge.rho_co2*co2_hmax * &
+            IF(co2ice_slope(ig,islope).ge.rho_co2*co2_hmax(islope) * &
                   cos(pi*def_slope_mean(islope)/180.)) THEN 
 
@@ -1041 +1059 @@
 
                 co2ice_slope(ig,iaval) = co2ice_slope(ig,iaval) + &
-               (co2ice_slope(ig,islope) - rho_co2* co2_hmax *     &
+               (co2ice_slope(ig,islope) - rho_co2* co2_hmax(islope) *     &
                cos(pi*def_slope_mean(islope)/180.)) *             &
                subslope_dist(ig,islope)/subslope_dist(ig,iaval) * &
@@ -1047 +1065 @@
                cos(pi*def_slope_mean(islope)/180.)                
 
-                co2ice_slope(ig,islope)=rho_co2*co2_hmax *        &
+                co2ice_slope(ig,islope)=rho_co2*co2_hmax(islope) *        &
                   cos(pi*def_slope_mean(islope)/180.)