Changeset 3337 for trunk/LMDZ.MARS/libf

Timestamp:

May 21, 2024, 5:59:03 PM (8 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

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (3 diffs)

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