Changeset 2316 for trunk/LMDZ.MARS/libf/phymars/vdifc_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/vdifc_mod.F (modified) (4 diffs)

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

@@ -835 +835 @@
             CALL hdo_surfex(ngrid,nlay,nq,ptimestep,
      &                      zt,zq,pqsurf,
-     &                      old_h2o_vap,pdqsdif,h2oflux,
+     &                      old_h2o_vap,pdqsdif,dwatercap_dif,
      &                      hdoflux)
             DO ig=1,ngrid
@@ -879 +879 @@
                 if(watercaptag(ig)) then
 c              dwatercap_dif same sign as pdqsdif (pos. toward ground)
-                   dwatercap_dif(ig) = -pdqsdif(ig,igcm_h2o_ice)
-     &                         - pqsurf(ig,igcm_h2o_ice)/ptimestep
+                   dwatercap_dif(ig) = pdqsdif(ig,igcm_h2o_ice)
+     &                         + pqsurf(ig,igcm_h2o_ice)/ptimestep
                    pdqsdif(ig,igcm_h2o_ice)=
      &                         -pqsurf(ig,igcm_h2o_ice)/ptimestep
@@ -892 +892 @@
                    endif
 
-                endif   
 c                 write(*,*)'flux vers le sol=',pdqsdif(ig,nq)
                   z1(ig)=1./(za(ig,1)+ zb(ig,2)*(1.-zd(ig,2)))
                   zc(ig,1)=(za(ig,1)*zq(ig,1,igcm_h2o_vap)+
-     $            zb(ig,2)*zc(ig,2) + (dwatercap_dif(ig)-
-c     $            zb(ig,2)*zc(ig,2) + (-
-     $            pdqsdif(ig,igcm_h2o_ice)) *ptimestep) *z1(ig)
+     $            zb(ig,2)*zc(ig,2) + 
+     $              (-pdqsdif(ig,igcm_h2o_ice)) *ptimestep) *z1(ig)
                   zq1temp(ig)=zc(ig,1)
-              endif ! if (.not.watercaptag(ig))
+                endif  !if .not.watercaptag(ig) 
+              endif ! if sublim more than surface
+
 c             Starting upward calculations for water :
               zq(ig,1,igcm_h2o_vap)=zq1temp(ig)
@@ -925 +925 @@
 
                if(pqsurf(ig,igcm_h2o_ice)
-     &           +(dwatercap_dif(ig)+pdqsdif(ig,igcm_h2o_ice))*ptimestep
+     &           +pdqsdif(ig,igcm_h2o_ice)*ptimestep
      &           .gt.frost_albedo_threshold) ! if there is still ice, T cannot exceed To
      &           pdtsrf(ig) = min(pdtsrf(ig),(To-ptsrf(ig))/ptimestep) ! ice melt case