phylmd/CV30_routines/cv30_param.f

module cv30_param_m

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

  ! Set parameters for Emanuel convection scheme. Includes
  ! microphysical parameters and parameters that control the rate of
  ! approach to quasi-equilibrium.

  USE dimphy, ONLY: klev

  implicit none

  integer, parameter:: minorig = 1 
  ! first level of convection (limit levels for convection)

  integer, parameter:: nl = klev - 1
  ! Limit for convection. The maximum number of levels to which
  ! convection can penetrate, plus 1. We should have:
  ! 6 <= nl <= KLEV - 1
  ! (because of locate in cv30_feed).

  real, parameter:: sigd = 0.01 
  ! fractional area covered by unsaturated downdraft
  ! \sigma_d in Emanuel (1991 928)

  real, parameter:: dtcrit = - 2.
  ! Critical buoyancy (K) used to adjust the approach to
  ! quasi-equilibrium. It must be < 0.

  real beta, alpha
  ! Parameters that control the rate of approach to quasi-equilibrium
  ! (their standard values are 1. and 0.96, respectively) (beta must
  ! be less than or equal to 1).

  private klev

contains

  SUBROUTINE cv30_param

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

    use comconst, only: dtphys

    ! Local:

    real, parameter:: tau = 8000.
    ! characteristic timescale used to compute alpha and beta

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

    beta = 1. - dtphys / tau

    alpha = 1.5E-3 * dtphys / tau * 1.5
    ! factor 1.5: increase alpha to compensate W decrease

  end SUBROUTINE cv30_param

end module cv30_param_m
1	guez	185	module cv30_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	195	! Set parameters for Emanuel convection scheme. Includes
6			! microphysical parameters and parameters that control the rate of
7			! approach to quasi-equilibrium.
8	guez	3
9	guez	187	USE dimphy, ONLY: klev
10
11	guez	69	implicit none
12	guez	3
13	guez	195	integer, parameter:: minorig = 1
14			! first level of convection (limit levels for convection)
15	guez	187
16			integer, parameter:: nl = klev - 1
17			! Limit for convection. The maximum number of levels to which
18	guez	195	! convection can penetrate, plus 1. We should have:
19			! 6 <= nl <= KLEV - 1
20			! (because of locate in cv30_feed).
21	guez	187
22	guez	195	real, parameter:: sigd = 0.01
23	guez	201	! fractional area covered by unsaturated downdraft
24			! \sigma_d in Emanuel (1991 928)
25	guez	3
26	guez	195	real, parameter:: dtcrit = - 2.
27			! Critical buoyancy (K) used to adjust the approach to
28			! quasi-equilibrium. It must be < 0.
29	guez	3
30	guez	69	real beta, alpha
31	guez	195	! Parameters that control the rate of approach to quasi-equilibrium
32			! (their standard values are 1. and 0.96, respectively) (beta must
33			! be less than or equal to 1).
34	guez	3
35	guez	187	private klev
36
37	guez	69	contains
38	guez	3
39	guez	195	SUBROUTINE cv30_param
40	guez	3
41	guez	97	! From LMDZ4/libf/phylmd/cv3_routines.F, version 1.5, 2005/07/11 15:20:02
42	guez	3
43	guez	195	use comconst, only: dtphys
44	guez	3
45	guez	195	! Local:
46	guez	3
47	guez	195	real, parameter:: tau = 8000.
48			! characteristic timescale used to compute alpha and beta
49
50	guez	69	!------------------------------------------------------------
51	guez	3
52	guez	195	beta = 1. - dtphys / tau
53	guez	3
54	guez	195	alpha = 1.5E-3 * dtphys / tau * 1.5
55			! factor 1.5: increase alpha to compensate W decrease
56	guez	69
57	guez	185	end SUBROUTINE cv30_param
58	guez	69
59	guez	185	end module cv30_param_m