Changeset 3451 for trunk/LMDZ.VENUS

Timestamp:

Oct 7, 2024, 3:37:34 PM (7 weeks ago)

Author:

emillour

Message:

Venus PCM:
Add "Age of Air" scheme from Maureen Cohen.
Activated with flag "ok_aoa=y" and requires that there is also
an "aoa" tracer in traceur.def
MC

Location:

trunk/LMDZ.VENUS

Files:

• deftank/field_def_physics.xml (modified) (1 diff)
• deftank/physiq.def (modified) (1 diff)
• libf/phyvenus/age_of_air_mod.F90 (added)
• libf/phyvenus/conf_phys.F90 (modified) (4 diffs)
• libf/phyvenus/ini_histmth.h (modified) (1 diff)
• libf/phyvenus/phyetat0.F90 (modified) (1 diff)
• libf/phyvenus/phyredem.F90 (modified) (2 diffs)
• libf/phyvenus/phys_state_var_mod.F90 (modified) (3 diffs)
• libf/phyvenus/physiq_mod.F (modified) (6 diffs)
• libf/phyvenus/write_histmth.h (modified) (1 diff)

trunk/LMDZ.VENUS/deftank/field_def_physics.xml

@@ -281 +281 @@
                    long_name="Kinetic Energy net flux" 
                    unit="W/m2" />
+           <!---- TRACERS ---->
+            <field id="age"
+                   long_name="Age of air" 
+                   unit="s" />
+            <field id="aoa"
+                   long_name="Age of air tracer" 
+                   unit="kg/kg" />
+           <!---- CHEMISTRY TRACERS ---->
             <field id="co2"
                    long_name="CO2" 

trunk/LMDZ.VENUS/deftank/physiq.def

@@ -41 +41 @@
 # Reinitialization of tracer abundances 
 reinit_trac=n
+#
+# Parameters for Age of air
+##~~~~~~~~~~~~~~~~~~~~~~~~~
+# Compute age of air (requires that there is an "aoa" tracer; default: n)
+ok_aoa=n
+# Reinitialize age of air tracer (default: y)
+reinit_aoa=n
+# Emission altitude (in model layers) of aoa tracer emission (default: 1)
+lev_aoa=1
 #
 # parameters for chemistry and microphysics

trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90

@@ -10 +10 @@
    use print_control_mod, only: lunout
    use IOIPSL, only: getin
+   use age_of_air_mod, only: ok_aoa, reinit_aoa, lev_aoa
 
    implicit none
@@ -318 +319 @@
   call getin('tr_scheme',tr_scheme)
 
+!Config Key = ok_aoa
+!Config Desc = Flag to run age of air subroutine
+!Config Def = .false.
+  ok_aoa=.false.
+  call getin('ok_aoa',ok_aoa)
+
+!Config Key = reinit_aoa
+!Config Desc = Flag to reinitiate age of air array to 0
+!Config Def = .true.
+  reinit_aoa=.true.
+  call getin('reinit_aoa',reinit_aoa)
+  
+!Config Key = lev_aoa
+!Config Desc = Model level where tracer associated with age of air is released
+!Config Def = 1 (surface)
+  lev_aoa=1
+  call getin('lev_aoa',lev_aoa)
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !   PARAMETRES DE LA CHIMIE/NUAGE dans physiq.def
@@ -369 +388 @@
     stop
   endif
+
+
+! Check that chemistry and age of air flages are not both active
+  if ((ok_chem) .and. (ok_aoa)) then
+     write(*,*) "Chemistry and age of air cannot be run together"
+     write(*,*) "ok_chem=",ok_chem," / ok_aoa=",ok_aoa
+     stop
+  endif
+
+
 
 !
@@ -566 +595 @@
   write(lunout,*)' Inclinaison =',R_incl
   write(lunout,*)' tr_scheme = ', tr_scheme
+  write(lunout,*)' ok_aoa = ', ok_aoa
+  write(lunout,*)' reinit_aoa = ', reinit_aoa
+  write(lunout,*)' lev_aoa = ', lev_aoa
   write(lunout,*)' iflag_pbl = ', iflag_pbl
   write(lunout,*)' z0 = ',z0 

trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h

@@ -109 +109 @@
 c    .                "ave(X)", zsto,zout)
 
+       IF ((iflag_trac.eq.1) .and. ok_aoa) THEN
+         CALL histdef(nid_mth, "age", "Age of air", "s",
+     .                nbp_lon,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+       ENDIF
 
 c plusieurs traceurs  !!!outputs in [vmr]

trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90

@@ -351 +351 @@
 ELSE
   ancien_ok=.false.
+ENDIF
+
+IF (startphy_file) THEN
+  CALL get_field("age",age,found)
+  IF (.not.found) THEN
+    PRINT*, "phyetat0: Age of air is missing"
+    PRINT*, "Reinitialising age of air to 0"
+    age(:,:) = 0.
+  ENDIF
+ELSE
+  age(:,:) = 0.
 ENDIF
 

trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90

@@ -15 +15 @@
   use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd, &
                                 east_gwstress, west_gwstress
-
-      implicit none
+  use age_of_air_mod, only: ok_aoa 
+  implicit none
 !======================================================================
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
@@ -109 +109 @@
 endif
 
+if (ok_aoa) then
+   call put_field("age","Age of air",age)
+endif
+
 ! close file
 

trunk/LMDZ.VENUS/libf/phyvenus/phys_state_var_mod.F90

@@ -117 +117 @@
       REAL,save,allocatable :: fder(:) ! Derive de flux (sensible et latente) 
 !$OMP THREADPRIVATE(dlw,fder)
+      REAL, ALLOCATABLE, SAVE :: age(:,:) ! Age of air (s)
+!$OMP THREADPRIVATE(age)
+
 CONTAINS
 
@@ -171 +174 @@
       ALLOCATE(topsw0(klon),toplw0(klon),solsw0(klon),sollw0(klon))
       ALLOCATE(dlw(klon), fder(klon))
+      ALLOCATE(age(klon,klev))
 
       ! allocate turb_mod module variables
@@ -212 +216 @@
       deallocate(topsw0,toplw0,solsw0,sollw0)
       deallocate(dlw, fder)
+      deallocate(age)
 
       ! deallocate turb_mod module variables

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -68 +68 @@
       USE cpdet_phy_mod, only: cpdet, t2tpot
       USE chemparam_mod
+      USE age_of_air_mod, ONLY: ok_aoa, reinit_aoa, i_aoa, init_age
+      USE age_of_air_mod, ONLY: aoa_ini, age_of_air
       USE conc
       USE param_v4_h
@@ -454 +456 @@
       
 c======================================================================
-c PREMIER APPEL SEULEMENT
+c DONE ONLY AT FIRST CALL
 c========================
-      IF (debut) THEN
+      IF (debut) THEN         
          allocate(source(klon,nqmax))
          allocate(prod_tr(klon,klev,nqmax))
@@ -473 +475 @@
 
          IF (if_ebil.ge.1) d_h_vcol_phy=0.
-c
-c appel a la lecture du physiq.def 
-c
+!
+! load flag and parameter values from physiq.def 
+!
          call conf_phys(ok_journe, ok_mensuel,
      .                  ok_instan,
@@ -743 +745 @@
       end if
          
+      if ((iflag_trac.eq.1) .and. ok_aoa) then
+          write(lunout,*) 'Initialising age of air subroutine'
+          allocate(init_age(klon,klev))
+          call aoa_ini(ok_chem, tr_scheme) ! get index of age of air tracer in the tracer array 
+          if (reinit_aoa) then
+            age(:,:) = 0. ! Set initial value of age of air to 0 everywhere
+            tr_seri(:,:,i_aoa) = 1e-30 ! Set initial value of tracer to tiny everywhere
+          endif ! reinitialisation loop
+          init_age(:,:) = age(:,:) !  save initial value of age, either read in from restartphy or set to 0
+      endif ! age of air initialisation
+
 !     initialise cloud parameters (for cl_scheme = 1)
 
@@ -1030 +1043 @@
 
       if (iflag_trac.eq.1) then
+
+        if ( ok_aoa ) then
+            call age_of_air(tr_seri(:,:,i_aoa),klon,klev,
+     I                   itap,pdtphys,init_age,
+     O                   age)
+        end if
+
 !====================================================================
 ! Case 1: pseudo-chemistry with relaxation toward fixed profile
@@ -2148 +2168 @@
       if (iflag_trac == 1) then
 
+         if (ok_aoa) then
+            call send_xios_field("age",age)
+            call send_xios_field("aoa",tr_seri(:,:,i_aoa))
+         endif
+
 ! production and destruction rate, cm-3.s-1
 ! Beware of the context*.xml file !!

trunk/LMDZ.VENUS/libf/phyvenus/write_histmth.h

@@ -48 +48 @@
       call histwrite_phy(nid_mth,.false.,"mmean",itau_w,mmean)
 c     call histwrite_phy(nid_mth,.false.,"rho",itau_w,rho)
+      IF ((iflag_trac.eq.1) .and. ok_aoa) THEN
+        call histwrite_phy(nid_mth,.false.,"age",itau_w,age)
+      ENDIF
 
 c plusieurs traceurs  !!!outputs in [vmr]