Changeset 2316 for trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F

Timestamp:

May 12, 2020, 11:51:22 AM (5 years ago)

Author:

lrossi

Message:

Mars GCM:
Fixing some errors in vdifc_mod related to variable watercap. This variable was also integrated to the hdo cycle.
Also added watercap output for the 1D model.
LR

File:

• trunk/LMDZ.MARS/libf/phymars/hdo_surfex_mod.F (modified) (5 diffs)

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

@@ -7 +7 @@
       subroutine hdo_surfex(ngrid,nlay,nq,ptimestep,
      &                      zt,zq,pqsurf,
-     &                      old_h2o_vap,pdqsdif,h2oflux,
+     &                      old_h2o_vap,pdqsdif,dwatercap_dif,
      &                      hdoflux)
 
@@ -37 +37 @@
       REAL, INTENT(IN) :: old_h2o_vap(ngrid)     ! traceur d'eau avant
                                            !traitement de l'eau (kg/kg)
+      REAL, INTENT(IN) :: dwatercap_dif(ngrid)  ! trend related to permanent ice
       REAL, INTENT(INOUT) :: pdqsdif(ngrid,nq)    ! tendance towards surface 
                                  !   (kg/kg.s-1)
-      REAL, INTENT(IN) :: h2oflux(ngrid)       ! value of vapour flux of H2O
-                                      ! same sign as pdqsdif
 
 c     Output:
@@ -51 +50 @@
       REAL extrasublim ! sublimation in excess of surface ice
       REAL tmpratio(ngrid)   ! D/H ratio in flux to surf
+      REAL h2oflux(ngrid)       ! value of vapour flux of H2O
+                                      ! same sign as pdqsdif
 
       INTEGER ig,l
@@ -62 +63 @@
 
         DO ig=1,ngrid
+              
+            h2oflux(ig) = pdqsdif(ig,igcm_h2o_ice) + 
+     &          dwatercap_dif(ig)
 
             !! IF Sublimation
@@ -87 +91 @@
      &           .gt.pqsurf(ig,igcm_h2o_ice)) then
 
-C               This is how much we sublimate in excess of
+C               dwatercap_dif is how much we sublimate in excess of
 C               pqsurf for H2O                        
-                  extrasublim=
-     &              ( pqsurf(ig,igcm_h2o_ice) /ptimestep 
-     &                   +h2oflux(ig))
-c               extrasublim has same sign conv. as pdqsdif
- 
 C               hdoflux(ig) is the flux of HDO from atm. to surf.
 c               The D/H of the old ice is supposed to be 5 SMOW
 c               We need D/H of the flux to be 5, so we need
-c               extrasublim* 5 * 2 * 155.76e-6 (=1 SMOW)
+c               dwatercap_dif* 5 * 2 * 155.76e-6 (=1 SMOW)
                     hdoflux(ig)= hdoflux(ig)
-     &                       +(extrasublim*(2.*155.76e-6)*5.)
-!    &                       +extrasublim
-!                 hdoflux(ig) = h2oflux(ig) !test
-              endif
-                endif ! watercap
+     &                   +(dwatercap_dif(ig)*(2.*155.76e-6)*5.)
+                endif
+             endif ! watercap
 
             else ! condensation