Changeset 2800 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Oct 13, 2022, 5:04:49 PM (2 years ago)

Author:

mvals

Message:

Mars GCM:
watercloud_mod.F: implementation of the HDO ice and vapor fluxes computation in the case of supersat=false.
MV

File:

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

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

@@ -115 +115 @@
       
       !$OMP THREADPRIVATE(firstcall)
-
+! HDO cycle
+      REAL :: alpha(ngrid,nlay) ! fractionation coefficient for HDO
 ! Representation of sub-grid water ice clouds A. Pottier 2013
       REAL :: ztclf(ngrid, nlay)
@@ -586 +587 @@
              pdqcloud(ig,l,igcm_h2o_ice) = 
      &         -pdqcloud(ig,l,igcm_h2o_vap)
+!            !! HDO ice flux has to be modified in consequence
+             IF (hdo) THEN
+!               !! Logically only condensation can happen in this case
+               IF (pdqcloud(ig,l,igcm_h2o_ice) .gt. 0.0) THEN
+                 IF ( zq0(ig,l,igcm_h2o_vap) .gt. qperemin ) THEN
+             alpha(ig,l) = exp(13525./zt(ig,l)**2.-5.59d-2)
+             pdqcloud(ig,l,igcm_hdo_ice) = 
+     &         pdqcloud(ig,l,igcm_h2o_ice)*alpha(ig,l)*
+     &         ( zq0(ig,l,igcm_hdo_vap)
+     &           /zq0(ig,l,igcm_h2o_vap) )
+             pdqcloud(ig,l,igcm_hdo_ice) =
+     &            min(pdqcloud(ig,l,igcm_hdo_ice),
+     &              zq0(ig,l,igcm_hdo_vap)/ptimestep)
+                 ELSE
+             pdqcloud(ig,l,igcm_hdo_ice) = 0.
+                 ENDIF
+               !! sublimation
+               ELSE
+                 IF ( zq0(ig,l,igcm_h2o_ice) .gt. qperemin ) THEN
+             pdqcloud(ig,l,igcm_hdo_ice) = 
+     &         pdqcloud(ig,l,igcm_h2o_ice)*
+     &         ( zq0(ig,l,igcm_hdo_ice)
+     &           /zq0(ig,l,igcm_h2o_ice) )
+             pdqcloud(ig,l,igcm_hdo_ice) =
+     &          max(pdqcloud(ig,l,igcm_hdo_ice),
+     &            -zq0(ig,l,igcm_hdo_ice)/ptimestep)
+                 ELSE
+             pdqcloud(ig,l,igcm_hdo_ice) = 0.
+                 ENDIF
+               ENDIF !IF (pdqcloud(ig,l,igcm_h2o_ice).gt.0.)
+             pdqcloud(ig,l,igcm_hdo_vap) = 
+     &            -pdqcloud(ig,l,igcm_hdo_ice)
+             ENDIF !IF (hdo)
              ! no need to correct ccn_number, updaterad can handle this properly.
           ENDIF