Changeset 2761 for LMDZ5/trunk/libf/phylmd

Timestamp:

Jan 11, 2017, 12:39:50 PM (8 years ago)

Author:

jyg

Message:

Two options introduced concerning deep convection
and wakes:

+ deep convection is allowed only when top
temperature is lower than t_top_max.
(Default: t_top_max=1000.).
+ wake number per unit area may now be
different over ocean and over land:
wdens_ref_o and wdens_ref_l.
(Default is:wdens_ref_o=8.e-12,
wdens_ref_l=8.e-12)

Location:

LMDZ5/trunk/libf/phylmd

Files:

• calwake.F90 (modified) (8 diffs)
• cv3_routines.F90 (modified) (5 diffs)
• cv3param.h (modified) (2 diffs)
• cva_driver.F90 (modified) (2 diffs)
• wake.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/calwake.F90

@@ -26 +26 @@
 
   USE dimphy
+  USE phys_state_var_mod, ONLY: pctsrf
+  USE indice_sol_mod, ONLY: is_oce
+  USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level
   IMPLICIT NONE
   ! ======================================================================
@@ -51 +54 @@
   ! ------
   REAL, DIMENSION(klon, klev),   INTENT (OUT)        :: dt_wake, dq_wake
-  INTEGER, DIMENSION(klon),      INTENT (OUT)        :: wake_k
+  REAL, DIMENSION(klon),         INTENT (OUT)        :: wake_k
   REAL, DIMENSION(klon, klev),   INTENT (OUT)        :: wake_d_deltat_gw
   REAL, DIMENSION(klon),         INTENT (OUT)        :: wake_h
@@ -69 +72 @@
   ! -----------------
   INTEGER                                            :: i, l
+  INTEGER, DIMENSION(klon)                           :: znatsurf    ! 0 if pctsrf(is_oce)>0.1; 1 else.
   REAL                                               :: aire
   REAL, DIMENSION(klon, klev)                        :: p,  pi
@@ -97 +101 @@
 
 
-  ! print *, '-> calwake, wake_s ', wake_s(1)
+  IF (prt_level >= 10) THEN
+    print *, '-> calwake, wake_s input ', wake_s(1)
+  ENDIF
 
   rdcp = 1./3.5
+
+  znatsurf(:) = 0
+  DO i = 1,klon
+    IF (pctsrf(i,is_oce) < 0.1) znatsurf(i) = 1
+  ENDDO
 
 
@@ -124 +135 @@
   END DO
 
+!----------------------------------------------------------------
+!         Initialize tendencies to zero
+!----------------------------------------------------------------
+dtls(:,:) = 0.
+dqls(:,:) = 0.
+d_deltat_gw(:,:) = 0.
+d_deltatw(:,:) = 0.
+d_deltaqw(:,:) = 0.
+d_sigmaw(:) = 0.
+d_wdens(:) = 0.
+!
+
   DO i = 1, klon
     sigd0(i) = sigd(i)
@@ -132 +155 @@
   END DO
 
-  DO i = 1, klon
-    ktopw(i) = wake_k(i)
-  END DO
+!!jyg!  DO i = 1, klon                  
+!!jyg!    ktopw(i) = NINT(wake_k(i))    
+!!jyg!  END DO                          
 
   DO i = 1, klon
@@ -178 +201 @@
   END DO
 
-  CALL wake(p, ph, pi, dtime, &
+  CALL wake(znatsurf, p, ph, pi, dtime, &
     te, qe, omgbe, &
     dtdwn, dqdwn, amdwn, amup, dta, dqa, &
@@ -254 +277 @@
     END IF
   END DO
+!
 
 !jyg<  

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -85 +85 @@
      cv_flag_feed=1
      CALL getin_p('cv_flag_feed',cv_flag_feed)
+     T_top_max = 1000.
+     CALL getin_p('t_top_max',T_top_max)
      dpbase=-40.
      CALL getin_p('dpbase',dpbase)
@@ -119 +121 @@
      CALL getin_p('tlcrit',tlcrit)
 
+    WRITE (*, *) 't_top_max=', t_top_max
     WRITE (*, *) 'dpbase=', dpbase
     WRITE (*, *) 'pbcrit=', pbcrit
@@ -1071 +1074 @@
 END SUBROUTINE icefrac
 
-SUBROUTINE cv3_undilute2(nloc, ncum, nd, icb, icbs, nk, &
+SUBROUTINE cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &
                          tnk, qnk, gznk, hnk, t, q, qs, gz, &
                          p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
@@ -1116 +1119 @@
   REAL, DIMENSION (nloc, nd), INTENT (INOUT)         :: tp, tvp, clw   ! Input for k = 1, icb+1 (computed in cv3_undilute1)
                                                                        ! Output above
+  INTEGER, DIMENSION (nloc), INTENT (INOUT)          :: iflag
 
 !outputs:
@@ -1594 +1598 @@
 ! ori          frac(i)=max(frac(i),0.0)
 ! ori 540  continue
+
+! --------------------------------------------------------------------
+!   Prevent convection when top is too hot
+! --------------------------------------------------------------------
+  DO i = 1,ncum
+    IF (t(i,inb(i)) > T_top_max) iflag(i) = 10
+  ENDDO
 
 ! =====================================================================

LMDZ5/trunk/libf/phylmd/cv3param.h

@@ -20 +20 @@
       real dtovsh, dpbase, dttrig
       real dtcrit, tau, beta, alpha, alpha1
+      real T_top_max
       real tau_stop, noconv_stop
       real wbmax
@@ -32 +33 @@
                       ,dtovsh, dpbase, dttrig &
                       ,dtcrit, tau, beta, alpha, alpha1 &
+                      ,T_top_max &
                       ,tau_stop, noconv_stop &
                       ,wbmax &

LMDZ5/trunk/libf/phylmd/cva_driver.F90

@@ -353 +353 @@
 !         9     No moist convection: cloud base is higher
 !               then the level NL-1.
+!        10     No moist convection: cloud top is too warm.
+!
 
 ! ft:   Array of temperature tendency (K/s) of dimension ND, defined at same
@@ -907 +909 @@
         if (prt_level >= 9) &
              PRINT *, 'cva_driver -> cv3_undilute2'
-      CALL cv3_undilute2(nloc, ncum, nd, icb, icbs, nk, &              !na->nd
+      CALL cv3_undilute2(nloc, ncum, nd, iflag, icb, icbs, nk, &        !na->nd
                          tnk, qnk, gznk, hnk, t, q, qs, gz, &
                          p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &

LMDZ5/trunk/libf/phylmd/wake.F90

@@ -2 +2 @@
 ! $Id$
 
-SUBROUTINE wake(p, ph, pi, dtime, &
+SUBROUTINE wake(znatsurf, p, ph, pi, dtime, &
                 te0, qe0, omgb, &
                 dtdwn, dqdwn, amdwn, amup, dta, dqa, &
@@ -124 +124 @@
   ! --------------------
 
+  INTEGER, DIMENSION (klon),        INTENT(IN)          :: znatsurf
   REAL, DIMENSION (klon, klev),     INTENT(IN)          :: p, pi
   REAL, DIMENSION (klon, klev+1),   INTENT(IN)          :: ph
@@ -167 +168 @@
   LOGICAL, SAVE                                         :: first = .TRUE.
   !$OMP THREADPRIVATE(first)
-  REAL, SAVE                                            :: stark, wdens_ref, coefgw, alpk
+!jyg<
+!!  REAL, SAVE                                            :: stark, wdens_ref, coefgw, alpk
+  REAL, SAVE, DIMENSION(2)                              :: wdens_ref
+  REAL, SAVE                                            :: stark, coefgw, alpk
+!>jyg
   REAL, SAVE                                            :: crep_upper, crep_sol  
   !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, crep_upper, crep_sol)
@@ -319 +324 @@
   alpk=0.25
   CALL getin_p('alpk',alpk)
-  wdens_ref=8.E-12
-  CALL getin_p('wdens_ref',wdens_ref)
+!jyg<
+!!  wdens_ref=8.E-12
+!!  CALL getin_p('wdens_ref',wdens_ref)
+  wdens_ref(1)=8.E-12
+  wdens_ref(2)=8.E-12
+  CALL getin_p('wdens_ref_o',wdens_ref(1))    !wake number per unit area ; ocean
+  CALL getin_p('wdens_ref_l',wdens_ref(2))    !wake number per unit area ; land
+!>jyg
   coefgw=4.
   CALL getin_p('coefgw',coefgw)
@@ -326 +337 @@
   WRITE(*,*) 'stark=', stark
   WRITE(*,*) 'alpk=', alpk
-  WRITE(*,*) 'wdens_ref=', wdens_ref
+!jyg<
+!!  WRITE(*,*) 'wdens_ref=', wdens_ref
+  WRITE(*,*) 'wdens_ref_o=', wdens_ref(1)
+  WRITE(*,*) 'wdens_ref_l=', wdens_ref(2)
+!>jyg
   WRITE(*,*) 'coefgw=', coefgw 
 
@@ -339 +354 @@
   ! Les densites peuvent evoluer si les poches debordent
   ! (voir au tout debut de la boucle sur les substeps)
-  wdens(:) = wdens_ref
+!jyg<
+!!  wdens(:) = wdens_ref
+  DO i = 1,klon
+    wdens(i) = wdens_ref(znatsurf(i)+1)
+  ENDDO 
+!>jyg
 
   ! print*,'stark',stark