Changeset 3984 for trunk

Timestamp:

Dec 8, 2025, 11:36:38 AM (6 days ago)

Author:

jbclement

Message:

PEM:
Addition of a module "tracers" to retain properties of atmospheric tracers.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• metamorphism.F90 (added)
• pem.F90 (modified) (4 diffs)
• pemetat0.F90 (modified) (3 diffs)
• tracers.F90 (added)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -817 +817 @@
 - Addition of a module "metamorphism" to compute the PCM frost at the PEM beginning and give it back to the PCM at the PEM end. The frost is considered as the ice given by the PCM "startfi.nc" which is above the yearly minimum. Thereby, metamorphism is performed through this operation.
 - Ice reservoirs representation in the PEM is modernized.
+
+== 08/12/2025 == JBC
+Addition of a module "tracers" to retain properties of atmospheric tracers.

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

@@ -82 +82 @@
 use dimradmars_mod,             only: totcloudfrac, albedo
 use dust_param_mod,             only: tauscaling
-use tracer_mod,                 only: noms, mmol, igcm_h2o_vap ! Tracer names and molar masses
+use tracer_mod,                 only: noms ! Tracer names
 use mod_phys_lmdz_para,         only: is_parallel, is_sequential, is_mpi_root, is_omp_root, is_master
 use planete_h,                  only: year_day
@@ -88 +88 @@
 use comcstfi_h,                 only: mugaz
 use metamorphism,               only: ini_frost_id, set_frost4PCM, iPCM_h2ofrost, iPCM_co2frost
+use tracers,                    only: mmol
 
 #ifndef CPP_1D
@@ -416 +417 @@
 call surfini(ngrid,nslope,qsurf)
 call ini_frost_id(nqtot,noms)
-do nnq = 1,nqtot
-    if (noms(nnq) == "h2o_vap") then
-        igcm_h2o_vap = nnq
-        mmol(igcm_h2o_vap) = 18.
-    endif
-enddo
+call ini_tracers_id(nqtot,noms)
 
 !------------------------
@@ -767 +763 @@
                         tsoil_PEM_timeseries(ig,isoil,islope,t) = tsoil_PEM_timeseries(ig,isoil,islope,t)*tsoil_PEM(ig,isoil,islope)/tsoil_avg_old(ig,isoil)
                         ! Update of watersoil density
-                        watersoil_density_PEM_timeseries(ig,isoil,islope,t) = exp(beta_clap_h2o/tsoil_PEM_timeseries(ig,isoil,islope,t) + alpha_clap_h2o)/tsoil_PEM_timeseries(ig,isoil,islope,t)*mmol(igcm_h2o_vap)/(mugaz*r)
+                        watersoil_density_PEM_timeseries(ig,isoil,islope,t) = exp(beta_clap_h2o/tsoil_PEM_timeseries(ig,isoil,islope,t) + alpha_clap_h2o)/tsoil_PEM_timeseries(ig,isoil,islope,t)*mmol%h2o/(mugaz*r)
                         if (isnan(tsoil_PEM(ig,isoil,islope))) call abort_pem("PEM - Update Tsoil","NaN detected in tsoil_PEM",1)
                     enddo

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

@@ -18 +18 @@
 use constants_marspem_mod,      only: alpha_clap_h2o, beta_clap_h2o, TI_breccia, TI_bedrock
 use soil_thermalproperties_mod, only: update_soil_thermalproperties
-use tracer_mod,                 only: mmol, igcm_h2o_vap ! tracer names and molar masses
 use abort_pem_mod,              only: abort_pem
 use compute_soiltemp_mod,       only: ini_tsoil_pem, compute_tsoil_pem
@@ -27 +26 @@
 use surfdat_h,                  only: watercaptag, perennial_co2ice, qsurf
 use metamorphism,               only: frost4PCM, iPCM_h2ofrost, iPCM_co2frost
+use tracers,                    only: mmol
 
 implicit none
@@ -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(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%h2o/(mugaz*r)
         enddo !islope
         write(*,*) 'PEMETAT0: TSOIL done'