Changeset 2995

Timestamp:

Jul 17, 2023, 7:04:26 PM (18 months ago)

Author:

llange

Message:

Mars PEM
*Implementation of the H2O glacier flow laws
*The algorithm for glacier flow is now more generic and not specific for co2 ice
*Principle: if ice thickness > ice mass, computed from models (cf note attached in the wiki), then the excess of ice is transfered
LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• abort_pem.F (modified) (1 diff)
• co2glaciers_mod.F90 (deleted)
• conf_pem.F90 (modified) (2 diffs)
• glaciers_mod.F90 (added)
• pem.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/abort_pem.F

@@ -32 +32 @@
       character(len=*), intent(in):: message
 
-      write(lunout,*) 'in abort_gcm'
+      write(lunout,*) 'in abort_pem'
 
 #ifdef CPP_XIOS

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

@@ -25 +25 @@
   USE comsoil_h_pem, only: soil_pem,fluxgeo,water_reservoir_nom,depth_breccia,depth_bedrock,reg_thprop_dependp 
   USE adsorption_mod,only: adsorption_pem
-  USE co2glaciers_mod, only: co2glaciersflow
+  USE glaciers_mod, only: co2glaciersflow,h2oglaciersflow
   use ice_table_mod, only: icetable_equilibrium, icetable_dynamic 
 
@@ -79 +79 @@
   co2glaciersflow = .true.
   CALL getin('co2glaciersflow', co2glaciersflow)
+
+  h2oglaciersflow = .true.
+  CALL getin('h2oglaciersflow', h2oglaciersflow)
 
   reg_thprop_dependp = .false.

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

@@ -93 +93 @@
       use conf_pem_mod,     only: conf_pem
       use pemredem,         only:  pemdem0,pemdem1
-      use co2glaciers_mod,        only: co2glaciers_evol,co2glaciersflow
+      use glaciers_mod,        only: co2glaciers_evol,h2oglaciersflow,co2glaciersflow,h2oglaciersflow
       use criterion_pem_stop_mod, only: criterion_waterice_stop,criterion_co2_stop
       use constants_marspem_mod,  only: alpha_clap_co2,beta_clap_co2, alpha_clap_h2o,beta_clap_h2o, &
@@ -942 +942 @@
        call WRITEDIAGFI(ngrid,'co2ice_slope'//str2,'CO2 ice','kg.m-2',2,qsurf(:,igcm_co2,islope))
        call WRITEDIAGFI(ngrid,'tendencies_co2_ice_slope'//str2,'CO2 ice tend','kg.m-2.year-1',2,tendencies_co2_ice(:,islope))
+       call WRITEDIAGFI(ngrid,'Flow_co2ice_slope'//str2,'CO2 ice flow','Boolean',2,flag_co2flow(:,islope))
      ENDDO