Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

Mod_PHY_CM_dat.f90 @ 5336

Last change on this file since 5336 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics

File size: 9.1 KB

Line
1	module Mod_PHY_CM_dat
2
3	!--------------------------------------------------------------------------+
4	! Sun 30-Jun-2013 MAR \|
5	! module Mod_PHY_CM_dat contains the data of \|
6	! Cloud Microphsics Scheme \|
7	! \|
8	! version 3.p.4.1 created by H. Gallee, Sat 23-Mar-2013 \|
9	! Last Modification by H. Gallee, Sun 30-Jun-2013 \|
10	! \|
11	!--------------------------------------------------------------------------+
12
13
14	use Mod_Real
15
16
17	IMPLICIT NONE
18
19
20	integer, parameter :: npt_CM = 1 ! Nb of OUTPUT Points
21	integer, SAVE :: ipt_CM ! Index of OUTPUT Points
22	integer, SAVE, dimension(npt_CM) :: i0__CM ! x-Coordinate of the OUTPUT Point
23	integer, SAVE, dimension(npt_CM) :: j0__CM ! y-Coordinate of the OUTPUT Point
24	integer, SAVE, dimension(npt_CM) :: k0__CM ! z-Coordinate of the OUTPUT Point
25	integer, SAVE, dimension(npt_CM) :: ikl0CM ! v-Coordinate of the OUTPUT Point
26	data (i0__CM(ipt_CM), ipt_CM=1,npt_CM) / 1/
27	data (j0__CM(ipt_CM), ipt_CM=1,npt_CM) / 1/
28	data (k0__CM(ipt_CM), ipt_CM=1,npt_CM) / 1/
29
30	real(kind=real8), SAVE :: qv_MIN = 3.e-6 ! Minimum Specific Humidity (Ch. Tricot) [kg/kg]
31	!real(kind=real8), SAVE :: qh_MIN = 1.e-18 ! Minimum Hydrometeors Concentration [kg/kg]
32	real(kind=real8), SAVE :: qh_MIN = 1.e-9 ! Minimum Hydrometeors Concentration [kg/kg]
33	real(kind=real8), SAVE :: WatIce = 273.16 ! Computation of Saturation pressure
34	real(kind=real8), SAVE :: ExpWat = 5.138 ! Ref.: Dudhia 1989, JAS
35	real(kind=real8), SAVE :: ExpWa2 = 6827. !
36
37	real(kind=real8), SAVE :: RH_MAX = 1.0 ! Maximum Relative Humidity (1 = 100%) [-]
38	real(kind=real8), SAVE :: RHcrit = 1.0 ! Sursaturation Threshold (1 = 100%) [-]
39	real(kind=real8), SAVE :: CFrMIN = 1.e-6 ! Minimum Cloud Fraction when qw OR qi > 0 [-]
40	real(kind=real8), SAVE :: SSImax = 101.0 ! Maximum Sursaturation % ICE (101 ==> RH= 201%) [%]
41	real(kind=real8), SAVE :: n0___s = 3.0e+6 ! intercept parameter / snow gamma distribution
42	real(kind=real8), SAVE :: n0___r = 8.0e+6 ! intercept parameter / rain gamma distribution
43	real(kind=real8), SAVE :: n0___g = 4.0e+4 ! intercept parameter / graupel gamma distribution
44	! (Lin et al. 1983, JCAM 22, p.1068: 1, 2 and 3)
45	real(kind=real8), SAVE :: qi0_DC = 0.10e-3 ! Ice Crystals Critical Mixing Ratio (tuned Dome C) [kg/kg]
46	! qi0_DC = 0.30e-3 ! Ice Crystals Critical Mixing Ratio (standard) [kg/kg]
47	real(kind=real8), SAVE :: qs__D0 = 2.00e-4 ! Smallest Diameter of Particles in the snow Class [m]
48	! ! Ref.: Levkov et al. 1992, Contr.Atm.Phys.65, p.41, para 1
49	real(kind=real8), SAVE :: T_NuId = -5. ! Deposition & Condensation-Freezing Nucleation Temperature [dgC]
50	real(kind=real8), SAVE :: a_NuId= -0.639 ! Deposition & Condensation-Freezing Nucleation Parameters
51	real(kind=real8), SAVE :: b_NuId = 0.1296 ! Ref.: Meyers et al. 1992 p.713
52	real(kind=real8), SAVE :: T_NuIc = -2. ! Contact Freezing Nucleation Temperature [dgC]
53	real(kind=real8), SAVE :: a_NuIc= -2.80 ! Contact Freezing Nucleation Parameters
54	real(kind=real8), SAVE :: b_NuIc = 0.262 ! Ref.: Meyers et al. 1992 p.713
55
56	real(kind=real8), SAVE :: Di_Hex = 1.1e+4 ! Diameter of Hexagonal Plates
57	! Ref.: Emde & Kalig 1989, Ann.Geophys. 7, p.408 (15): D
58
59	real(kind=real8), SAVE :: TmaxHM = -3. ! [dgC]
60	real(kind=real8), SAVE :: TminHM = -8. ! [dgC]
61	real(kind=real8), SAVE :: wa__HM = 1. ! V. Wind > wa__HM => Nucleation II (Hall-Mossop) may occur[m/s]
62
63	real(kind=real8), SAVE :: TqwFrz =-35.e0 ! T < TqwFrz => Instantaneous Freezing [dgC]
64	! Ref.: Levkov et al. 1992, Contr.Atm.Phys.65, p.39
65	real(kind=real8), SAVE :: qw_VOL = 18.e-15 ! Typical Cloud Droplet Volume [m3] (typ. diam.: 32.5 mim)
66	! OR Cloud Droplet Weight [T]
67
68	real(kind=real8), SAVE :: qisMAX = 0.0010 ! Ice Crystals MAX Concentration before autoconversion [kg/kg]
69	! Ref.: Lin et al. 1983, JCAM 22, p.1070 (21)
70	real(kind=real8), SAVE :: qigMAX = 0.0006 ! Ice Crystals MAX Concentration before autoconversion [kg/kg]
71	! Ref.: Lin et al. 1983, JCAM 22, p.1074 (37)
72	real(kind=real8), SAVE :: qwTURB = 0.27 ! qwTURB=1/3 ln(1/k), k=0.8 (droplets dispersion parameter)
73	! Ref.: Martin et al. 1994, JAS 51, p.1823
74	real(kind=real8), SAVE :: C1_EkM = 0.14e-3 ! Partial Condensation Scheme
75	real(kind=real8), SAVE :: C2_EkM = 9.75e+0 ! Ref.: Ek and Mahrt 1991, An.Geoph. 9, 716--724
76
77	real(kind=real8), SAVE :: qw_MAX = 0.0001 ! Cloud Droplets MAX Concentration before autoconversion [kg/kg]
78	real(kind=real8), SAVE :: qw_MAXL= 0.0020 ! Cloud droplets MAX concentration before autoconversion
79	real(kind=real8), SAVE :: rwCrit = 10.0e-6 ! Droplets Autoconversion: Critical Radius
80	! Ref.: Liou and Ou, 1989
81	real(kind=real8), SAVE :: c_Sund = 1.0e-4 ! Droplets Autoconversion: 1/(characteristic time scale)
82	! Ref.: Sundqvist, 1988
83
84	real(kind=real8), SAVE :: cc1 = 1.200e-04 ! Cloud droplets autoconversion parameter
85	real(kind=real8), SAVE :: cc2 = 1.569e-12 ! Cloud droplets autoconversion parameter
86	real(kind=real8), SAVE :: dd0 = 0.15e0 ! Cloud droplets autoconversion parameter
87	! Ref.: Lin et al. 1983, JCAM 22, p.1076 (50)
88
89
90
91
92	! ======================================================================
93	! Bergeron Process Data (given by Koenig, 1971, J.A.S. 28,p235) ========
94
95	real(kind=real8), SAVE, dimension(31) :: aa1
96
97	data aa1/0.7939e-07 , 0.7841e-06 , 0.3369e-05 , 0.4336e-05 , &
98	& 0.5285e-05 , 0.3728e-05 , 0.1852e-05 , 0.2991e-06 , &
99	& 0.4248e-06 , 0.7434e-06 , 0.1812e-05 , 0.4394e-05 , &
100	& 0.9145e-05 , 0.1725e-06 , 0.3348e-04 , 0.1725e-04 , &
101	& 0.9175e-05 , 0.4412e-05 , 0.2252e-05 , 0.9115e-06 , &
102	& 0.4876e-06 , 0.3473e-06 , 0.4758e-06 , 0.6306e-06 , &
103	& 0.8573e-06 , 0.7868e-06 , 0.7192e-06 , 0.6513e-06 , &
104	& 0.5956e-06 , 0.5333e-06 , 0.4834e-06 /
105
106
107	real(kind=real8), SAVE, dimension(31) :: aa2
108
109	data aa2/0.4006e0, 0.4831e0, 0.5320e0, 0.5307e0, 0.5319e0, &
110	& 0.5249e0, 0.4888e0, 0.3894e0, 0.4047e0, 0.4318e0, &
111	& 0.4771e0, 0.5183e0, 0.5463e0, 0.5651e0, 0.5813e0, &
112	& 0.5655e0, 0.5478e0, 0.5203e0, 0.4906e0, 0.4447e0, &
113	& 0.4126e0, 0.3960e0, 0.4149e0, 0.4320e0, 0.4506e0, &
114	& 0.4483e0, 0.4460e0, 0.4433e0, 0.4413e0, 0.4382e0, &
115	& 0.4361e0/
116
117	! Bergeron Process Data (given by Koenig, 1971, J.A.S. 28,p235) ========
118	! ======================================================================
119
120
121
122	end module Mod_PHY_CM_dat

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format