Changeset 3724 for trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

Timestamp:

Apr 14, 2025, 4:42:01 PM (2 months ago)

Author:

mmaurice

Message:

Generic PCM

Add radioactive tracers in order to compare tracers mixing between 3D and 1D.

MM

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 ¶

@@ -37 +37 @@
                           igcm_h2o_ice, igcm_h2o_vap, igcm_dustbin, &
                           igcm_co2_ice, nesp, is_chim, is_condensable,constants_epsi_generic
+                        !  igcm_co2_ice, nesp, is_chim, is_condensable,constants_epsi_generic, &
+                        !  half_life, top_prod, bot_prod
       use time_phylmdz_mod, only: ecritphy, iphysiq, nday
       use phyetat0_mod, only: phyetat0,tab_cntrl_mod
@@ -364 +366 @@
       REAL,allocatable,save :: zdqchim(:,:,:) ! Calchim_asis routine
       REAL,allocatable,save :: zdqschim(:,:)  ! Calchim_asis routine
+      REAL zdqradio(ngrid,nlayer,nq)   ! Radioactive decay
 !$OMP THREADPRIVATE(zdqchim,zdqschim)
       real zdqvolc(ngrid,nlayer,nq) ! injection by volcanoes (kg/kg_of_air/s)
@@ -2072 +2075 @@
          qsurf_hist(:,:) = qsurf(:,:)
 
+  ! -----------------------------
+  !   VI.7. Radioactive Tracers
+  ! -----------------------------
+         
+         call radioactive_tracers(ngrid,nlayer,nq,ptimestep,pq,zdqradio)
+         pdq(1:ngrid,1:nlayer,1:nq) = pdq(1:ngrid,1:nlayer,1:nq) + zdqradio(1:ngrid,1:nlayer,1:nq)
+
       endif! end of if 'tracer'