phylmd/CV30_routines/cv30_param.f

module cv3_param_m

  ! From LMDZ4/libf/phylmd/cvparam3.h, version 1.1.1.1, 2004/05/19 12:53:09

  ! Parameters for Emanuel convection scheme:
  ! - microphysical parameters
  ! - parameters that control the rate of approach to quasi-equilibrium

  implicit none

  integer noff ! integer limit for convection (nd - noff)
  integer minorig ! first level of convection
  integer nl, nlp, nlm
  real sigd ! FRACTIONAL AREA COVERED BY UNSATURATED DNDRAFT 
  real spfac ! FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD 

  real pbcrit
  !  CRITICAL CLOUD DEPTH (MB) BENEATH WHICH THE PRECIPITATION
  ! EFFICIENCY IS ASSUMED TO BE ZERO

  real ptcrit
  ! CLOUD DEPTH (MB) ABOVE WHICH THE PRECIPitation EFFICIENCY IS
  ! ASSUMED TO BE UNITY

  real omtrain
  real dtovsh, dpbase, dttrig

  real dtcrit
  ! CRITICAL BUOYANCY (K) USED TO ADJUST THE APPROACH TO
  ! QUASI-EQUILIBRIUM. IT MUST BE LESS THAN 0.

  real beta, alpha
  ! 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)

  real tau ! CHARACTERISTIC TIMESCALE USED TO COMPUTE ALPHA and BETA
  real delta
  real betad

contains

  SUBROUTINE cv3_param(nd, delt)

    ! From LMDZ4/libf/phylmd/cv3_routines.F, version 1.5, 2005/07/11 15:20:02

    ! Set parameters for Emanuel convection scheme

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

    !------------------------------------------------------------

    ! 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
    ! cf. FH epmax = 0.993

    omtrain = 45.0 ! used also for snow (no distinction 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)

  end SUBROUTINE cv3_param

end module cv3_param_m
1	guez	69	module cv3_param_m
2	guez	3
3	guez	97	! From LMDZ4/libf/phylmd/cvparam3.h, version 1.1.1.1, 2004/05/19 12:53:09
4	guez	3
5	guez	182	! Parameters for Emanuel convection scheme:
6	guez	69	! - microphysical parameters
7			! - parameters that control the rate of approach to quasi-equilibrium
8	guez	3
9	guez	69	implicit none
10	guez	3
11	guez	69	integer noff ! integer limit for convection (nd - noff)
12			integer minorig ! first level of convection
13			integer nl, nlp, nlm
14			real sigd ! FRACTIONAL AREA COVERED BY UNSATURATED DNDRAFT
15			real spfac ! FRACTION OF PRECIPITATION FALLING OUTSIDE OF CLOUD
16	guez	3
17	guez	69	real pbcrit
18			! CRITICAL CLOUD DEPTH (MB) BENEATH WHICH THE PRECIPITATION
19			! EFFICIENCY IS ASSUMED TO BE ZERO
20	guez	3
21	guez	69	real ptcrit
22			! CLOUD DEPTH (MB) ABOVE WHICH THE PRECIPitation EFFICIENCY IS
23			! ASSUMED TO BE UNITY
24	guez	3
25	guez	69	real omtrain
26			real dtovsh, dpbase, dttrig
27	guez	3
28	guez	69	real dtcrit
29			! CRITICAL BUOYANCY (K) USED TO ADJUST THE APPROACH TO
30			! QUASI-EQUILIBRIUM. IT MUST BE LESS THAN 0.
31	guez	3
32	guez	69	real beta, alpha
33			! PARAMETERS THAT CONTROL THE RATE OF APPROACH TO QUASI-EQUILIBRIUM
34			! (THEIR STANDARD VALUES ARE 1.0 AND 0.96, RESPECTIVELY) (BETA MUST
35			! BE LESS THAN OR EQUAL TO 1)
36	guez	3
37	guez	69	real tau ! CHARACTERISTIC TIMESCALE USED TO COMPUTE ALPHA and BETA
38			real delta
39			real betad
40	guez	3
41	guez	69	contains
42	guez	3
43	guez	69	SUBROUTINE cv3_param(nd, delt)
44	guez	3
45	guez	97	! From LMDZ4/libf/phylmd/cv3_routines.F, version 1.5, 2005/07/11 15:20:02
46	guez	3
47	guez	182	! Set parameters for Emanuel convection scheme
48	guez	3
49	guez	69	integer, intent(in):: nd
50			real, intent(in):: delt ! timestep (seconds)
51	guez	3
52	guez	69	!------------------------------------------------------------
53	guez	3
54	guez	69	! Limit levels for convection:
55	guez	3
56	guez	69	noff = 1
57			minorig = 1
58	guez	97	nl = nd - noff
59	guez	69	nlp=nl+1
60			nlm=nl-1
61	guez	3
62	guez	69	! "Microphysical" parameters:
63
64			sigd = 0.01
65			spfac = 0.15
66			pbcrit = 150.0
67			ptcrit = 500.0
68			! cf. FH epmax = 0.993
69
70			omtrain = 45.0 ! used also for snow (no distinction rain/snow)
71
72			! Misc:
73
74			dtovsh = -0.2 ! dT for overshoot
75			dpbase = -40. ! definition cloud base (400m above LCL)
76			dttrig = 5. ! (loose) condition for triggering
77
78			! Rate of approach to quasi-equilibrium:
79
80			dtcrit = -2.0
81			tau = 8000.
82			beta = 1.0 - delt/tau
83			alpha = 1.5E-3 * delt/tau
84			! Increase alpha to compensate W decrease:
85			alpha = alpha*1.5
86
87			! Interface cloud parameterization:
88
89			delta=0.01 ! cld
90
91			! Interface with boundary-layer (gust factor): (sb)
92
93			betad=10.0 ! original value (from convect 4.3)
94
95			end SUBROUTINE cv3_param
96
97			end module cv3_param_m