Changeset 2800 for trunk/LMDZ.MARS

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

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/README

@@ -3742 +3742 @@
  - atm_wat_profile = -1 : water vapor and ice are free in the column. This is the default mode of the 1d model which conserves total water tracers (qsurf + watercap + h2o_ice + h2o_vap)
  - atm_wat_profile = 0 : dry atmosphere, where q(h2o_ice) = q(h2o_vap) = 0.
+
+== 13/10/2022 == MV
+watercloud_mod.F: implementation of the HDO ice and vapor fluxes computation in the case of supersat=false.

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