Changeset 3154 for trunk/LMDZ.MARS/libf

Timestamp:

Dec 10, 2023, 6:32:46 PM (20 months ago)

Author:

llange

Message:

Mars PCM
Small change following previous commit: when computing the virtual surface temperature, if frost is at the surface, qvap_surf = qsat. else; it is qvap in the first layer
LL

File:

• trunk/LMDZ.MARS/libf/phymars/vdif_cd_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/vdif_cd_mod.F90

@@ -140 +140 @@
       tol_iter =  0.01
       mu_h2o = 18e-3
-      
+      tol_frost = 1e-4
       reynolds(:) = 0.
       pz0t = 0.
@@ -162 +162 @@
          DO ig = 1,ngrid
             temp_v(ig) = ph(ig,1)*(1.+pqvap(ig,1)/(mu_h2o/mumean(ig)))/(1.+pqvap(ig,1))
-            call watersat(1,pts(ig),pp(ig,1),qsat(ig))
-            tsurf_v(ig) = pts(ig)*(1.+qsat(ig)/(mu_h2o/mumean(ig)))/(1.+qsat(ig))
+            IF(pqsurf(ig).gt.tol_frost) then
+               call watersat(1,pts(ig),pp(ig,1),qsat(ig))
+               tsurf_v(ig) = pts(ig)*(1.+qsat(ig)/(mu_h2o/mumean(ig)))/(1.+qsat(ig))
+            ELSE
+               tsurf_v(ig) = pts(ig)*(1.+pqvap(ig,1)/(mu_h2o/mumean(ig)))/(1.+pqvap(ig,1))
+            ENDIF
          ENDDO
       ELSE