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

LMDZ5/trunk/libf/phylmd/thermcell_main.F90

@@ -22 +22 @@
 
       USE dimphy
+      USE ioipsl
       USE comgeomphy , ONLY:rlond,rlatd
       IMPLICIT NONE
@@ -44 +45 @@
 !     4. un detrainement
 !
+! Modif 2013/01/04 (FH hourdin@lmd.jussieu.fr)
+!    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"
+!
 !=======================================================================
+
 
 !-----------------------------------------------------------------------
@@ -79 +90 @@
 
       integer icount
+
+      integer, save :: dvdq=1,dqimpl=-1
+!$OMP THREADPRIVATE(dvdq,dqimpl)
       data icount/0/
       save icount
@@ -247 +261 @@
 
       if (debut)  then
+!        call getin('dvdq',dvdq)
+!        call getin('dqimpl',dqimpl)
+
+         if (iflag_thermals==15.or.iflag_thermals==16) then
+            dvdq=0
+            dqimpl=-1
+         else
+            dvdq=1
+            dqimpl=1
+         endif
+
          fm0=0.
          entr0=0.
@@ -593 +618 @@
 !------------------------------------------------------------------
 
-      call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse,  &
+      call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse,  &
      &                    zthl,zdthladj,zta,lev_out)
-      call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse,  &
+      call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse,  &
      &                   po,pdoadj,zoa,lev_out)
 
@@ -620 +645 @@
 
 !IM 090508  
-      if (1.eq.1) then
-!IM 070508 vers. _dq       
-!     if (1.eq.0) then
-
+      if (dvdq == 0 ) then
 
 ! Calcul du transport de V tenant compte d'echange par gradient
@@ -629 +651 @@
 
          call thermcell_dv2(ngrid,nlay,ptimestep,fm0,entr0,masse  &
-     &    ,fraca,zmax  &
+!    &    ,fraca*dvdq,zmax &
+     &    ,fraca,zmax &
      &    ,zu,zv,pduadj,pdvadj,zua,zva,lev_out)
 
@@ -635 +658 @@
 
 ! calcul purement conservatif pour le transport de V
-         call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse  &
+         call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse  &
      &    ,zu,pduadj,zua,lev_out)
-         call thermcell_dq(ngrid,nlay,ptimestep,fm0,entr0,masse  &
+         call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse  &
      &    ,zv,pdvadj,zva,lev_out)
+
       endif