Changeset 2312 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

May 6, 2020, 4:46:33 PM (5 years ago)

Author:

lrossi

Message:

MARS GCM:
Implementation of HDO in the physics
New tracers hdo_vap and hdo_ice are added. Cf. traceur.def.isotopes in deftank for exemple of traceur.def.
All HDO related computations are activated by flag hdo=.true. in callphys.def. (see callphys.def.hdo in deftank for an example)
Additional option hdofrac (true by default) allows for turning on/off fractionation (for tests).
For now, only functional with simplified cloud scheme (so microphys=.false. and activice=.false. also recommended).
Initialisation of start files can be done with option 'inihdo' in newstart.
LR

File:

• trunk/LMDZ.MARS/libf/aeronomars/inichim_newstart.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/aeronomars/inichim_newstart.F90

@@ -87 +87 @@
       igcm_h2o_vap=0
       igcm_h2o_ice=0
+      igcm_hdo_vap=0
+      igcm_hdo_ice=0
       igcm_co2=0
       igcm_co=0
@@ -274 +276 @@
            igcm_h2o_ice = iq
            mmol(igcm_h2o_ice) = 18.
+           count = count + 1
+        end if
+        if (noms(iq) == "hdo_vap") then
+           igcm_hdo_vap = iq
+           mmol(igcm_hdo_vap) = 18.
+           count = count + 1
+        end if
+        if (noms(iq) == "hdo_ice") then
+           igcm_hdo_ice = iq
+           mmol(igcm_hdo_ice) = 18.
            count = count + 1
         end if
@@ -561 +573 @@
                   do n = 1,nspe
                      iq = niq(n)
-                     if (iq /= igcm_h2o_vap .or. flagh2o == 1) then
+                     if (iq /= igcm_h2o_vap .or. iq/= igcm_hdo_vap .or. flagh2o == 1) then
                         pq(i,j,l,iq) = vmrint(n)*mmol(iq)/mmean(i,j,l)
                      end if