Changeset 3075 for trunk/LMDZ.MARS/libf

Timestamp:

Oct 6, 2023, 4:48:17 PM (21 months ago)

Author:

jnaar

Message:

Added coefficient in improvedclouds (adapt_ts subroutine) to ensure that the computed number of subtimestep, as a integer subdivision of ptimestep, is properly scaled to the chosen timestep (scaled from standard iphysiq=5)

File:

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

trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F ¶

@@ -738 +738 @@
       real,intent(in) :: potcond ! condensible vapor / sublimable ice(kg/kg)
       real :: alpha, beta ! Coefficients for power law
+      real :: defstep,coef ! Default ptimestep of 7.5 mins (iphysiq=5)
       integer,intent(out) :: zimicro ! number of ptimestep division
 
+c      Default ptimestep : defstep (7.5 mins)
+       defstep=88775.*5./960.
+       coef=ptimestep/defstep
 c      Conservative coefficients :
        alpha=1.81846504e+11
        beta=1.54550140e+00
-       zimicro=ceiling(min(max(alpha*abs(potcond)**beta,5.),7000.))
+       zimicro=ceiling(min(max(coef*alpha*abs(potcond)**beta,5.),7000.))
 
       END SUBROUTINE adapt_imicro