Changeset 1738 for LMDZ5/trunk/libf/phylmd/calltherm.F90

Timestamp:

Apr 2, 2013, 11:48:33 AM (11 years ago)

Author:

idelkadi

Message:

Modifications for numerical stability of the boundary layer parameterizations.
Concerning yamada4.F :

• option with asymptotic mixing length l0 imposed (iflag_pbl=8 and 9)
• option with a new temporal scheme (iflag_pbl=10 and 11)
• Correction for very stable PBLs (iflag_pbl=10 and 11) iflag_pbl=8 converges numerically with NPv3.1 iflag_pbl=11 -> now starts with NP from start files created by ce0l

-> the model can run with longer time-steps.

Concerning thermals :

Introduction of an implicit computation of vertical advection in
the environment of thermal plumes in thermcell_dq
impl = 0 : explicit, 1 : implicit, -1 : old version
controled by iflag_thermals =

15, 16 run with impl=-1 : numerical convergence with NPv3
17, 18 run with impl=1 : more stable

15 and 17 correspond to the activation of the stratocumulus "bidouille"

Modified routines (phylmd/):
calltherm.F90 : for managing the various options of thermcell_dq
coef_diff_turb_mod.F90 : yamada4 called for iflag_pbl<= 18 instead of 11
physiq.F : desactivation of the vertical diffusion of TKE for iflag_pbl=10
thermcellV0_main.F90 : calling thermcell_dq with implicit scheme
thermcell_dq.F90 : thermcell_dq with optional explicit/implicit scheme
thermcell_main.F90 : various option for vertical transport by thermals
yamada4.F : Yamada scheme with new options

File:

• LMDZ5/trunk/libf/phylmd/calltherm.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/calltherm.F90

@@ -234 +234 @@
      &      ,tau_thermals,Ale,Alp,lalim_conv,wght_th &
      &      ,zmax0,f0,zw2,fraca)
-          else if (iflag_thermals==15.or.iflag_thermals==16) then
+          else if (iflag_thermals>=15.and.iflag_thermals<=18) then
 
 !            print*,'THERM iflag_thermas_ed=',iflag_thermals_ed
@@ -271 +271 @@
 ! fait bien ce qu'on croit.
 
-       flag_bidouille_stratocu=iflag_thermals<=12.or.iflag_thermals==14.or.iflag_thermals==16
+       flag_bidouille_stratocu=iflag_thermals<=12.or.iflag_thermals==14.or.iflag_thermals==16.or.iflag_thermals==18
 
       if (iflag_thermals<=12) then