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

Timestamp:

May 29, 2024, 11:12:58 PM (19 months ago)

Author:

jnaar

Message:

Hardcoded threshold in vdifc subtimestep now changed to callphys variable "frost_albedo_threshold" for clarity and flexibility. JN

File:

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

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

@@ -1487 +1487 @@
 c et en adaptant avec un critere "dtmax" du chauffage a accomoder
 c dtmax est regle empiriquement (pour l'instant) a 0.5 K
+      use surfdat_h, only: frost_albedo_threshold
 
       integer,intent(in) :: naersize
@@ -1499 +1500 @@
 
       do i=1,naersize
-        if ((qsurf(i).lt.0.005).and.
+c     If not on permanent ice or thin frost layer :
+c     do not use a subtimestep (ntsub=1)
+        if ((qsurf(i).lt.frost_albedo_threshold).and.
      &      (.not.watercaptag(i))) then
           ntsub(i) = 1