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

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

File:

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

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) =