Changeset 2516 for trunk/LMDZ.MARS/libf/phymars/watercloud_mod.F

Timestamp:

May 6, 2021, 5:31:05 PM (4 years ago)

Author:

mvals

Message:

Mars GCM:
Correction of the threshold value used to prevent from negative values in watercloud_mod.F: the threshold is set to 1e-16 (=qperemin, an already
existing parameter for the water cycle), instead of 1e-8. This former value was at the origin of a strong and unrealistic drying out of the upper
atmosphere, when the photochemistry was active, during the first part of the martian year (months 3 and 4). (bug reported by Francisco
Gonzalez-Galindo).
MV

File:

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

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

@@ -22 +22 @@
      &                      igcm_dust_mass, igcm_dust_number,
      &                      igcm_ccn_mass, igcm_ccn_number,
-     &                      rho_dust, nuice_sed, nuice_ref
+     &                      rho_dust, nuice_sed, nuice_ref,
+     &                      qperemin
       use dimradmars_mod, only: naerkind
       IMPLICIT NONE
@@ -514 +515 @@
           IF (pq(ig,l,igcm_h2o_ice) + ptimestep*
      &       (pdq(ig,l,igcm_h2o_ice) + pdqcloud(ig,l,igcm_h2o_ice)) 
-     &       .le. 1.e-8) THEN
+     &       .le. qperemin) THEN
            pdqcloud(ig,l,igcm_h2o_ice) = 
      &     - pq(ig,l,igcm_h2o_ice)/ptimestep - pdq(ig,l,igcm_h2o_ice)
@@ -526 +527 @@
           IF (pq(ig,l,igcm_h2o_vap) + ptimestep*
      &       (pdq(ig,l,igcm_h2o_vap) + pdqcloud(ig,l,igcm_h2o_vap)) 
-     &       .le. 1.e-8) THEN
+     &       .le. qperemin) THEN
            pdqcloud(ig,l,igcm_h2o_vap) = 
      &     - pq(ig,l,igcm_h2o_vap)/ptimestep - pdq(ig,l,igcm_h2o_vap)