Changeset 2312 for trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

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/phymars/callsedim_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/callsedim_mod.F

@@ -17 +17 @@
      &                      rho_dust, rho_q, radius, varian,
      &                      igcm_ccn_mass, igcm_ccn_number,
-     &                      igcm_h2o_ice, nuice_sed, nuice_ref,
+     &                      igcm_h2o_ice, igcm_hdo_ice,
+     &                      nuice_sed, nuice_ref,
      &                      igcm_ccnco2_mass,igcm_ccnco2_number,
      &                      igcm_co2_ice, igcm_stormdust_mass, 
@@ -499 +500 @@
 c -----------------------------------------------------------------
            else if ((iq .eq. iccn_mass) .or. (iq .eq. iccn_number)
-     &       .or. (iq .eq. igcm_h2o_ice)) then
+     &       .or. (iq .eq. igcm_h2o_ice)
+     &       .or. (iq .eq. igcm_hdo_ice)) then
             if (microphys) then
               ! water ice sedimentation