Changeset 3154 for trunk/LMDZ.MARS

Timestamp:

Dec 10, 2023, 6:32:46 PM (12 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

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/vdif_cd_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4395 +4395 @@
 Add the possibility to compute Cd; Ch based on the virtual potential temperature to account for water flotability.
 To do so, a boolean "virtual", set to false by default, must be set to true (for now hard coded as future modifications will follow).
+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

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