Changeset 3337 for trunk/LMDZ.MARS

Timestamp:

May 21, 2024, 5:59:03 PM (6 months ago)

Author:

jnaar

Message:

Threshold in qsurf(ice) in the computation of the subtimestep to avoid calling it when ice is not present on the surface (was already the case, but the condition was 0 and it could thus happen with very very low quantities, which is useless). JN

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/vdifc_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4637 +4637 @@
 == 21/05/2024 == JBC
 Addition of the paleoclimate variables in the change the number of slopes by "newtart.F" + some simplifications of the way it is done.
+
+== 21/05/2024 == JN
+In vdifc_mod, adding a threshold on the subtimestep of water vapor
+sublimation/condensation to avoid computation when surface ice is too low.
+Also nsubtimestep is now stored in diagfi at the end of the subtimestep for
+clarity of the code.

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

@@ -1016 +1016 @@
 c           nsubtimestep(ig)=1 !for debug
            subtimestep = ptimestep/nsubtimestep(ig)
-                          call write_output('subtimestep',
-     &                'vdifc substimestep length','s',subtimestep)
             ztsrf(ig)=ptsrf(ig,islope)  !  +pdtsrf(ig)*subtimestep
             zq_tmp_vap(ig,:,:) =zq(ig,:,:)
@@ -1414 +1412 @@
      &                     "Ground ice latent heat flux",
      &                     "W.m-2",surf_h2o_lh(:,iflat))
+         call write_output("nsubtimestep",
+     &                     "Number of subtimesteps in vdifc",
+     &                     "int",nsubtimestep(:))
 
          call write_output('zdqsdif_ssi_frost_tot',
@@ -1498 +1499 @@
 
       do i=1,naersize
-        if ((qsurf(i).eq.0).and.
+        if ((qsurf(i).lt.0.005).and.
      &      (.not.watercaptag(i))) then
           ntsub(i) = 1