Changeset 2856

Timestamp:

Dec 23, 2022, 2:09:37 PM (3 years ago)

Author:

llange

Message:

PEM
CO2 glacier flows added in a module
Maximum thickness of glacier before flow is not anymore hardcoded, but computed with the value of Tcond and the slope angle.
The method is from A.Grau Galfore (LPG), inspired from Nye et al., 2000
Works for high stress regim, will be adapted to consider also low stress regim (Tco2 > 200K);
LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• co2glaciers_mod.F90 (added)
• pem.F90 (modified) (4 diffs)

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

@@ -85 +85 @@
 
       use pemredem, only:  pemdem1
-
+      use co2glaciers_mod,only: co2glaciers_evol
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! SOIL
       use comsoil_h_PEM, only: soil_pem,ini_comsoil_h_PEM,end_comsoil_h_PEM,nsoilmx_PEM, &
@@ -225 +225 @@
       REAL, dimension(:,:),allocatable  :: tendencies_co2_ice_phys_slope   ! physical point xslope field : Tendency of evolution of perenial co2 ice over a year
       REAL, dimension(:,:),allocatable  :: tendencies_co2_ice_phys_slope_ini ! physical point x slope field x nslope: Tendency of evolution of perenial co2 ice over a year in the GCM
-      REAL, dimension(:,:),allocatable  :: tendencies_h2o_ice_phys_slope   ! physical pointx slope  field : Tendency of evolution of perenial co2 ice over a year
-      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 ()
+      REAL, dimension(:,:),allocatable  :: tendencies_h2o_ice_phys_slope   ! physical pointx slope  field : Tendency of evolution of perenial co2 ice 
       REAL , dimension(:,:), allocatable :: flag_co2flow(:,:)   !(ngrid,nslope)          ! To flag where there is a glacier flow
       REAL , dimension(:), allocatable :: flag_co2flow_mesh(:)  !(ngrid)          ! To flag where there is a glacier flow
@@ -483 +480 @@
      PRINT*,'corresponding criterium = ',def_slope_mean(iflat)
 
-! CO2 max thickness (for glaciers flows)
-       allocate(co2_hmax(nslope))
-       if(nslope.eq.7) then ! ugly way to implement that ...
-!        CF documentation that explain how the values are computed
-         co2_hmax(1) = 1.5
-         co2_hmax(7) = co2_hmax(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.
-     elseif(nslope.eq.1) then
-         co2_hmax(1) = 1000000.
-     else
-         print *, "Nslope must be=to 1 or 7 for now!!"
-     endif
 
      allocate(flag_co2flow(ngrid,nslope))
@@ -963 +944 @@
 !------------------------
 
-      print *, "Co2 glacier flow"
-       DO ig = 1,ngrid
-        DO islope = 1,nslope
-          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(islope) * &
-                  cos(pi*def_slope_mean(islope)/180.)) THEN 
-! Second: determine the flatest slopes possible:
-                IF(islope.gt.iflat) THEN
-                  iaval=islope-1
-                ELSE 
-                 iaval=islope+1
-                ENDIF
-                do while ((iaval.ne.iflat).and.  &
-                    (subslope_dist(ig,iaval).eq.0))
-                  IF(iaval.gt.iflat) THEN
-                     iaval=iaval-1
-                  ELSE 
-                     iaval=iaval+1
-                  ENDIF
-                enddo
-              co2ice_slope(ig,iaval) = co2ice_slope(ig,iaval) + &
-               (co2ice_slope(ig,islope) - rho_co2* co2_hmax(islope) *     &
-               cos(pi*def_slope_mean(islope)/180.)) *             &
-               subslope_dist(ig,islope)/subslope_dist(ig,iaval) * &
-               cos(pi*def_slope_mean(iaval)/180.) /               &
-               cos(pi*def_slope_mean(islope)/180.)                
-
-              co2ice_slope(ig,islope)=rho_co2*co2_hmax(islope) *        &
-               cos(pi*def_slope_mean(islope)/180.) 
-
-              flag_co2flow(ig,islope) = 1.
-              flag_co2flow_mesh(ig) = 1.
-            ENDIF ! co2ice > hmax
-          ENDIF ! iflat
-        ENDDO !islope 
-       ENDDO !ig
+      print *, "Co2 glacier flows"
+
+
+
+    call co2glaciers_evol(timelen,ngrid,nslope,iflat,subslope_dist,def_slope_mean,vmr_co2_pem_phys,ps_phys_timeseries,&
+                         global_ave_press_GCM,global_ave_press_new,co2ice_slope,flag_co2flow,flag_co2flow_mesh)
+
+
+
+
 
 !------------------------