Changeset 3985 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Dec 8, 2025, 2:19:07 PM (7 weeks ago)

Author:

jbclement

Message:

PEM:
Addition of a module "phys_constants" to read and store physical parameter of the planet properly, i.e. without going through the module "comcstfi_h" and/or "comconst_mod".
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (modified) (3 diffs)

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

@@ -23 +23 @@
 use callkeys_mod,               only: startphy_file
 use glaciers_mod,               only: rho_co2ice, rho_h2oice
-use comcstfi_h,                 only: r, mugaz, pi
+use phys_constants,             only: pi, r, mugaz
 use surfdat_h,                  only: watercaptag, perennial_co2ice, qsurf
 use metamorphism,               only: frost4PCM, iPCM_h2ofrost, iPCM_co2frost
-use tracers,                    only: mmol
+use tracers,                    only: mmol, iPCM_qh2o
 
 implicit none
@@ -296 +296 @@
             endif !found
 
-            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)*mmol(igcm_h2o_vap)/(mugaz*r)
+            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)*mmol(iPCM_qh2o)/(mugaz*r)
         enddo ! islope
         write(*,*) 'PEMETAT0: TSOIL done'
@@ -469 +469 @@
 
 ! First raw initialization
-            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)*mmol%h2o/(mugaz*r)
+            watersoil_avg(:,nsoil_PCM + 1:nsoil_PEM,islope) = exp(beta_clap_h2o/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope) + alpha_clap_h2o)/tsoil_PEM(:,nsoil_PCM + 1:nsoil_PEM,islope)**mmol(iPCM_qh2o)/(mugaz*r)
         enddo !islope
         write(*,*) 'PEMETAT0: TSOIL done'