Changeset 2324 for trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F

Timestamp:

May 15, 2020, 3:24:17 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Follow-up of the last commit for the transport of the isotopic ratio: cleaning to help the user:
-initracer.F: hdo must be defined as an isotope (meaning its "father" has to be precised in traceur.def) and must be placed at the end of the list of
tracers in traceur.def.
-/deftank/traceur.def.isotopes: writing of traceur.def in the case of hdo=true has been simplified to match both the 3D and the 1D simulations.
-simpleclouds.F,hdo_surfex_mod.F,physiq_mod.F: thresholds of 1.e-16 have been replaced by the variable qperemin defined in tracer_mod.F.
MV

File:

• trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F (modified) (3 diffs)

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

@@ -11 +11 @@
 
       use tracer_mod, only: igcm_h2o_vap, igcm_h2o_ice,
-     &                      igcm_hdo_vap, igcm_hdo_ice
-
+     &                      igcm_hdo_vap, igcm_hdo_ice,
+     &                      qperemin
       use surfdat_h, only: watercaptag
       use geometry_mod, only: longitude_deg,latitude_deg
@@ -70 +70 @@
             if (h2oflux(ig).le.0.) then
 
-               if (pqsurf(ig,igcm_h2o_ice).gt.1e-16) then
+               if (pqsurf(ig,igcm_h2o_ice).gt.qperemin) then
                 pdqsdif(ig,igcm_hdo_ice) =
      &            pdqsdif(ig,igcm_h2o_ice)*
@@ -111 +111 @@
                endif
 
-               if (old_h2o_vap(ig).gt.1.e-16) then
+               if (old_h2o_vap(ig).gt.qperemin) then
                          pdqsdif(ig,igcm_hdo_ice)=
      &                      alpha_c(ig)*pdqsdif(ig,igcm_h2o_ice)*