Changeset 2944 for trunk/LMDZ.COMMON/libf/evolution/co2glaciers_mod.F90

Timestamp:

Apr 18, 2023, 10:41:16 AM (20 months ago)

Author:

llange

Message:

PEM
Introduce module 'constants_marspem': Constants that are used are now put in this module, and load in each subroutine.
LL

File:

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

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

@@ -64 +64 @@
 
      USE comconst_mod, ONLY: pi
+     use comcstfi_h, only: g
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -87 +88 @@
       REAL,INTENT(OUT) :: hmax(ngrid,nslope) ! Physical grid x Slope field: maximum co2 thickness before flaw [m]
 ! Local
-      REAL, PARAMETER :: g = 3.71 ! surface gravity [m/s^2]
       INTEGER,PARAMETER :: n = 7 ! flow law exponent Nye et al., 2000
       REAL,PARAMETER :: Rg = 8.3145 ! gas constant [J/K/mol]
@@ -226 +226 @@
 !!!
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-implicit none
+  
+use constants_marspem_mod,only : alpha_clap_co2,beta_clap_co2
+
+implicit none 
+
 ! arguments:
 ! ----------
@@ -243 +247 @@
       INTEGER :: ig,it ! for loop
       REAL :: ave ! intermediate to compute average
-      REAL :: alpha_clap, beta_clap ! Clapeyron law for CO2
-      alpha_clap = 23.3494 ! James et al. 1992
-      beta_clap = 3182.48  ! James et al. 1992
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -253 +254 @@
         ave = 0
         DO it = 1,timelen
-           ave = ave + beta_clap/(alpha_clap-log(vmr_co2_PEM(ig,it)*ps_GCM(ig,it)*global_ave_ps_GCM/global_ave_ps_PEM/100))
+           ave = ave + beta_clap_co2/(alpha_clap_co2-log(vmr_co2_PEM(ig,it)*ps_GCM(ig,it)*global_ave_ps_GCM/global_ave_ps_PEM/100))
         ENDDO
         Tcond(ig) = ave/timelen