Changeset 3404 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Aug 12, 2024, 8:06:46 PM (17 months ago)

Author:

llange

Message:

Mars PCM
Following previous commit, fixing a bug: more frost than what is a the surface could sublimes (when working with lag layer only). It is corrected, not in a elegant way...
LL

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (1 diff)

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

@@ -1121 +1121 @@
      &                .and.lag_layer)) then
 
-     
-!                     zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
+                    if ((-zdqsdif_tot(ig)*subtimestep)
+     &                .ge.(zqsurf(ig))) then
+                          zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
+                    endif
+
                     if((h2o_ice_depth(ig,islope).gt.0)
      &                  .and. (zqsurf(ig) .lt. tol_frost)