Changeset 2516

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/watercloud_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -3376 +3376 @@
  - hdo_vap is now treated independantly after h2o_vap (should be fine !)
 
-
+== 06/05/2021 == MV
+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).
+

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)