Changeset 3075

Timestamp:

Oct 6, 2023, 4:48:17 PM (14 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)

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/improvedclouds_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/changelog.txt

@@ -4245 +4245 @@
 == 06/10/2023 == JBC
 Transformation of 'writerestart1D' into module.
+
+== 06/10/2023 == JN
+Bugfix : adaptative timestep now is insensitive to ptimestep

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