/[lmdze]/trunk/phylmd/CV30_routines/cv30_param.f90
ViewVC logotype

Contents of /trunk/phylmd/CV30_routines/cv30_param.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 328 - (show annotations)
Thu Jun 13 14:40:06 2019 UTC (4 years, 11 months ago) by guez
File size: 1600 byte(s)
Change all `.f` suffixes to `.f90`. (The opposite was done in revision
82.)  Because of change of philosopy in GNUmakefile: we already had a
rewritten rule for `.f`, so it does not make the makefile longer to
replace it by a rule for `.f90`. And it spares us options of
makedepf90 and of the compiler. Also we prepare the way for a simpler
`CMakeLists.txt`.

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

  ViewVC Help
Powered by ViewVC 1.1.21