phylmd/CV30_routines/cv30_param.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/cv3_routines.F,v 1.5 2005/07/11 15:20:02 lmdzadmin Exp $
!
!
!
      SUBROUTINE cv3_param(nd,delt)
      use conema3_m
            use cvparam3
      implicit none

!------------------------------------------------------------
! Set parameters for convectL for iflag_con = 3
!------------------------------------------------------------

!
!   ***  PBCRIT IS THE CRITICAL CLOUD DEPTH (MB) BENEATH WHICH THE ***
!   ***      PRECIPITATION EFFICIENCY IS ASSUMED TO BE ZERO     ***
!   ***  PTCRIT IS THE CLOUD DEPTH (MB) ABOVE WHICH THE PRECIP. ***
!   ***            EFFICIENCY IS ASSUMED TO BE UNITY            ***
!   ***  SIGD IS THE FRACTIONAL AREA COVERED BY UNSATURATED DNDRAFT  ***
!   ***  SPFAC IS THE FRACTION OF PRECIPITATION FALLING OUTSIDE ***
!   ***                        OF CLOUD                         ***
!
! [TAU: CHARACTERISTIC TIMESCALE USED TO COMPUTE ALPHA & BETA]
!   ***    ALPHA AND BETA ARE PARAMETERS THAT CONTROL THE RATE OF ***
!   ***                 APPROACH TO QUASI-EQUILIBRIUM           ***
!   ***    (THEIR STANDARD VALUES ARE 1.0 AND 0.96, RESPECTIVELY) ***
!   ***           (BETA MUST BE LESS THAN OR EQUAL TO 1)        ***
!
!   ***    DTCRIT IS THE CRITICAL BUOYANCY (K) USED TO ADJUST THE ***
!   ***                 APPROACH TO QUASI-EQUILIBRIUM           ***
!   ***                     IT MUST BE LESS THAN 0              ***


      integer nd
      real, intent(in):: delt ! timestep (seconds)

! noff: integer limit for convection (nd-noff)
! minorig: First level of convection

! -- limit levels for convection:

      noff    = 1
      minorig = 1
      nl=nd-noff
      nlp=nl+1
      nlm=nl-1

! -- "microphysical" parameters:

      sigd   = 0.01
      spfac  = 0.15
      pbcrit = 150.0
      ptcrit = 500.0
!IM cf. FH     epmax  = 0.993

      omtrain = 45.0 ! used also for snow (no disctinction rain/snow)

! -- misc:

      dtovsh = -0.2 ! dT for overshoot
      dpbase = -40. ! definition cloud base (400m above LCL)
      dttrig = 5.   ! (loose) condition for triggering

! -- rate of approach to quasi-equilibrium:

      dtcrit = -2.0
      tau    = 8000.
      beta   = 1.0 - delt/tau
      alpha  = 1.5E-3 * delt/tau
! increase alpha to compensate W decrease:
      alpha  = alpha*1.5

! -- interface cloud parameterization:

      delta=0.01  ! cld

! -- interface with boundary-layer (gust factor): (sb)

      betad=10.0   ! original value (from convect 4.3)

      return
      end
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/phylmd/cv3_routines.F,v 1.5 2005/07/11 15:20:02 lmdzadmin Exp $
3			!
4	guez	47	!
5			!
6	guez	3	SUBROUTINE cv3_param(nd,delt)
7			use conema3_m
8	guez	47	use cvparam3
9	guez	3	implicit none
10
11	guez	47	!------------------------------------------------------------
12			! Set parameters for convectL for iflag_con = 3
13			!------------------------------------------------------------
14	guez	3
15	guez	47	!
16			! * PBCRIT IS THE CRITICAL CLOUD DEPTH (MB) BENEATH WHICH THE *
17			! * PRECIPITATION EFFICIENCY IS ASSUMED TO BE ZERO *
18			! * PTCRIT IS THE CLOUD DEPTH (MB) ABOVE WHICH THE PRECIP. *
19			! * EFFICIENCY IS ASSUMED TO BE UNITY *
20			! * SIGD IS THE FRACTIONAL AREA COVERED BY UNSATURATED DNDRAFT *
21			! * SPFAC IS THE FRACTION OF PRECIPITATION FALLING OUTSIDE *
22			! * OF CLOUD *
23			!
24			! [TAU: CHARACTERISTIC TIMESCALE USED TO COMPUTE ALPHA & BETA]
25			! * ALPHA AND BETA ARE PARAMETERS THAT CONTROL THE RATE OF *
26			! * APPROACH TO QUASI-EQUILIBRIUM *
27			! * (THEIR STANDARD VALUES ARE 1.0 AND 0.96, RESPECTIVELY) *
28			! * (BETA MUST BE LESS THAN OR EQUAL TO 1) *
29			!
30			! * DTCRIT IS THE CRITICAL BUOYANCY (K) USED TO ADJUST THE *
31			! * APPROACH TO QUASI-EQUILIBRIUM *
32			! * IT MUST BE LESS THAN 0 *
33	guez	3
34
35			integer nd
36	guez	12	real, intent(in):: delt ! timestep (seconds)
37	guez	3
38	guez	47	! noff: integer limit for convection (nd-noff)
39			! minorig: First level of convection
40	guez	3
41	guez	47	! -- limit levels for convection:
42	guez	3
43			noff = 1
44			minorig = 1
45			nl=nd-noff
46			nlp=nl+1
47			nlm=nl-1
48
49	guez	47	! -- "microphysical" parameters:
50	guez	3
51			sigd = 0.01
52			spfac = 0.15
53			pbcrit = 150.0
54			ptcrit = 500.0
55	guez	47	!IM cf. FH epmax = 0.993
56	guez	3
57			omtrain = 45.0 ! used also for snow (no disctinction rain/snow)
58
59	guez	47	! -- misc:
60	guez	3
61			dtovsh = -0.2 ! dT for overshoot
62			dpbase = -40. ! definition cloud base (400m above LCL)
63	guez	47	dttrig = 5. ! (loose) condition for triggering
64	guez	3
65	guez	47	! -- rate of approach to quasi-equilibrium:
66	guez	3
67			dtcrit = -2.0
68			tau = 8000.
69			beta = 1.0 - delt/tau
70			alpha = 1.5E-3 * delt/tau
71	guez	47	! increase alpha to compensate W decrease:
72	guez	3	alpha = alpha*1.5
73
74	guez	47	! -- interface cloud parameterization:
75	guez	3
76			delta=0.01 ! cld
77
78	guez	47	! -- interface with boundary-layer (gust factor): (sb)
79	guez	3
80			betad=10.0 ! original value (from convect 4.3)
81
82			return
83			end