Changeset 2801 for trunk/LMDZ.MARS

Timestamp:

Oct 14, 2022, 10:21:12 AM (2 years ago)

Author:

mvals

Message:

Mars GCM:
watercloud_mod.F: follow-up of previous commit.
improvedclouds_mod.F: Lamb et al. 2017 formula put as default in the fractionation by condensation calculation.
MV

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/README

@@ -3745 +3745 @@
 == 13/10/2022 == MV
 watercloud_mod.F: implementation of the HDO ice and vapor fluxes computation in the case of supersat=false.
+
+== 14/10/2022 == MV
+watercloud_mod.F: follow-up of previous commit.
+improvedclouds_mod.F: Lamb et al. 2017 formula set as default in the fractionation by condensation calculation.
+

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

@@ -415 +415 @@
      &          * sqrt(1.+mhdo/mco2) )
              !! Calculation of the fractionnation coefficient at equilibrium
-             alpha(ig,l) = exp(16288./zt(ig,l)**2.-9.34d-2)
-c             alpha = exp(13525./zt(ig,l)**2.-5.59d-2)  !Lamb
+c             alpha(ig,l) = exp(16288./zt(ig,l)**2.-9.34d-2)  ! Merlivat and Nief et al. 1967
+             alpha = exp(13525./zt(ig,l)**2.-5.59d-2)  ! Lamb et al. 2017
              !! Calculation of the 'real' fractionnation coefficient
              alpha_c(ig,l) = (alpha(ig,l)*satu)/

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

@@ -117 +117 @@
 ! HDO cycle
       REAL :: alpha(ngrid,nlay) ! fractionation coefficient for HDO
+      REAL :: zq0(ngrid,nlay,nq) ! Initial mixing ratio: intermediate variable for HDO
 ! Representation of sub-grid water ice clouds A. Pottier 2013
       REAL :: ztclf(ngrid, nlay)
@@ -574 +575 @@
        ! In case one does not want to allow supersatured water when using microphysics.
        ! Not done by default.
-       IF(.not.supersat) THEN     
+       IF(.not.supersat) THEN
+!       !! initial mixing ratios for initial D/H ratio calculation
+        zq0(:,:,:) = pq(:,:,:) + pdq(:,:,:)*ptimestep 
         zt  = pt + (pdt+pdtcloud)*ptimestep
         call watersat(ngrid*nlay,zt,pplay,zqsat)
@@ -592 +595 @@
                IF (pdqcloud(ig,l,igcm_h2o_ice) .gt. 0.0) THEN
                  IF ( zq0(ig,l,igcm_h2o_vap) .gt. qperemin ) THEN
+!            !! Lamb et al. 2017
              alpha(ig,l) = exp(13525./zt(ig,l)**2.-5.59d-2)
              pdqcloud(ig,l,igcm_hdo_ice) =