Changeset 2761 for LMDZ5/trunk/libf/phylmd/cv3_routines.F90

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)

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F90 (modified) (5 diffs)

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
 
 ! =====================================================================