Changeset 1738 for LMDZ5/trunk/libf/phylmd/thermcellV0_main.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/thermcellV0_main.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/thermcellV0_main.F90

@@ -519 +519 @@
 !------------------------------------------------------------------
 
-      call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse,  &
+      call thermcell_dq(ngrid,nlay,1,ptimestep,fm0,entr0,masse,  &
      &                    zthl,zdthladj,zta,lev_out)
-      call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse,  &
+      call thermcell_dq(ngrid,nlay,1,ptimestep,fm0,entr0,masse,  &
      &                   po,pdoadj,zoa,lev_out)
 
@@ -561 +561 @@
 
 ! calcul purement conservatif pour le transport de V
-         call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse  &
+         call thermcell_dq(ngrid,nlay,1,ptimestep,fm0,entr0,masse  &
      &    ,zu,pduadj,zua,lev_out)
-         call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse  &
+         call thermcell_dq(ngrid,nlay,1,ptimestep,fm0,entr0,masse  &
      &    ,zv,pdvadj,zva,lev_out)
       endif