Changeset 2259 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Apr 15, 2015, 12:51:04 PM (9 years ago)

Author:

jyg

Message:

Changes in Emanuel's deep convection scheme: the
upper bound of deep convection loops is set at
the first level above 22 km.

Modified files:

physiq.F90,
concvl.F90,
cva_driver.F90,
cv3a_compress.F90,
cv3a_uncompress.F90,
cv3_routines.F90

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -316 +316 @@
   SAVE top
   !$OMP THREADPRIVATE(bas, top)
-
+  !------------------------------------------------------------------
+  ! Upmost level reached by deep convection and related variable (jyg)
+  !
+  INTEGER izero
+  INTEGER k_upper_cv
+  !------------------------------------------------------------------
   !
   !=================================================================================================
@@ -528 +533 @@
   !                            !par la masse/airetot (moyglo_pondaima) et la vraie masse (moyglo_pondmass)
   !
-  ! Variables locales
+  !
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+  ! Local variables
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   !
   REAL rhcl(klon,klev)    ! humiditi relative ciel clair
@@ -570 +578 @@
   REAL zx_t, zx_qs, zdelta, zcor, zlvdcp, zlsdcp
   real zqsat(klon,klev)
+!
   INTEGER i, k, iq, ig, j, nsrf, ll, l, iiq
+!
   REAL t_coup
   PARAMETER (t_coup=234.0)
@@ -885 +895 @@
      igout=klon/2+1/klon
      write(lunout,*) 'DEBUT DE PHYSIQ !!!!!!!!!!!!!!!!!!!!'
+     write(lunout,*) 'igout, rlat, rlon ',igout, rlatd(igout)*180./3.141593, rlond(igout)*180./3.141593
      write(lunout,*) &
           'nlon,klev,nqtot,debut,lafin, jD_cur, jH_cur,pdtphys'
@@ -904 +915 @@
 
   if (first) then 
-
+     
      !CR:nvelles variables convection/poches froides
 
@@ -1333 +1344 @@
 #endif
      END IF
-     !
      !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -2230 +2240 @@
      !
      IF (ok_cvl) THEN ! new driver for convectL
-
+     !
+!jyg<
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+     ! Calculate the upmost level of deep convection loops: k_upper_cv
+     !  (near 22 km)
+   izero = igout
+   DO i = 1,klon
+     IF (abs(rlatd(i)) < 4.e-2 .AND. abs(rlond(i)) < 4.e-2) izero = i
+   ENDDO
+!
+   k_upper_cv = klev
+   DO k = klev,1,-1
+     IF (pphi(izero,k) > 22.e4) k_upper_cv = k
+   ENDDO
+   IF (prt_level .ge. 5) THEN
+     Print *, 'upmost level of deep convection loops: k_upper_cv = ',k_upper_cv
+   ENDIF
+     !
+!>jyg
         IF (type_trac == 'repr') THEN
            nbtr_tmp=ntra
@@ -2239 +2267 @@
         !c          CALL concvl (iflag_con,iflag_clos,
         CALL concvl (iflag_clos, &
-             dtime,paprs,pplay,t_undi,q_undi, &
+             dtime, paprs, pplay, k_upper_cv, t_undi,q_undi, &
              t_wake,q_wake,wake_s, &
              u_seri,v_seri,tr_seri,nbtr_tmp, &