Changeset 3724 for trunk/LMDZ.GENERIC

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

Location:

trunk/LMDZ.GENERIC

Files:

• deftank/traceur.def.radioactive (added)
• libf/phystd/initracer.F90 (modified) (3 diffs)
• libf/phystd/physiq_mod.F90 (modified) (3 diffs)
• libf/phystd/radioactive_tracers.F90 (added)
• libf/phystd/tracer_h.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/initracer.F90

@@ -129 +129 @@
        IF (.NOT. allocated(constants_metallicity_coeff)) allocate(constants_metallicity_coeff(nq))
        IF (.NOT. allocated(constants_RCPV_generic)) allocate(constants_RCPV_generic(nq))
+       IF (.NOT.ALLOCATED(half_life))         ALLOCATE(half_life(nqtot))
+       IF (.NOT.ALLOCATED(top_prod))         ALLOCATE(top_prod(nqtot))
+       IF (.NOT.ALLOCATED(bot_prod))         ALLOCATE(bot_prod(nqtot))
 
        !! initialization
@@ -141 +144 @@
        is_recomb_qset(:) = 0
        is_recomb_qotf(:) = 0
+       half_life(:)      = 0.
+       top_prod(:)       = 0.
+       bot_prod(:)       = 0.
        
        ! Added by JVO 2017 : these arrays are handled later
@@ -716 +722 @@
             write(*,*)'Parameter value (default) : '// &
               'is_rgcs = ',is_rgcs(iq)
-          endif
+          end if
+             ! option top_prod 
+              if (index(tracline,'top_prod=') .ne. 0) then
+                  read(tracline(index(tracline,'top_prod=') &
+                  +len('top_prod='):),*) top_prod(iq)
+                  write(*,*) ' Parameter value (traceur.def) : top_prod=', &
+                  top_prod(iq)
+              else
+                  write(*,*) ' Parameter value (default)     : top_prod=',  &
+                  top_prod(iq)
+              end if
+              ! option bot_prod 
+              if (index(tracline,'bot_prod=') .ne. 0) then
+                  read(tracline(index(tracline,'bot_prod=') &
+                  +len('bot_prod='):),*) bot_prod(iq)
+                  write(*,*) ' Parameter value (traceur.def) : bot_prod=', &
+                  bot_prod(iq)
+              else
+                  write(*,*) ' Parameter value (default)     : bot_prod=',  &
+                  bot_prod(iq)
+              end if
       end subroutine get_tracdat
 

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'
 

trunk/LMDZ.GENERIC/libf/phystd/tracer_h.F90

@@ -55 +55 @@
 !$OMP THREADPRIVATE(constants_Pref,constants_epsi_generic)
 !$OMP THREADPRIVATE(constants_RLVTT_generic,constants_metallicity_coeff,constants_RCPV_generic)
+
+       real, save, allocatable :: half_life(:) ! half-life (s) for radioactive tracers. If 0, tracer is not decaying.
+       real, save, allocatable :: top_prod(:)  ! top production rate (1/s) for radioactive tracers. Non-dim, so use 1.
+       real, save, allocatable :: bot_prod(:)  ! bottom production rate (1/s) for radioactive tracers. Non-dim, so use 1.
+!$OMP THREADPRIVATE(half_life,top_prod,bot_prod)
+
 
 !$OMP THREADPRIVATE(is_condensable,is_rgcs)   !also added by LT