/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_param.f
ViewVC logotype

Contents of /trunk/Sources/phylmd/CV30_routines/cv30_param.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 201 - (show annotations)
Mon Jun 6 17:42:15 2016 UTC (7 years, 10 months ago) by guez
File size: 1587 byte(s)
Removed intermediary objects of cv_thermo_m, access suphec_m
directly. Procedure cv_thermo disappeared, all objects are named
constants.

In cv_driver and below, limited extents of arrays to what is needed.

lv, cpn and th in cv30_compress were set at level nl + 1 but lv1, cpn1
and th1 are not defined at this level. This did not lead to an error
because values at nl + 1 were not used.

Removed test on ok_sync in phystokenc because it is not read at run
time. Printing min and max of output NetCDF variables is heavy and
archaic.

Used histwrite_phy in phytrac.

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