Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

Mod_PHY_CM_dat.f90 @ 5442

Last change on this file since 5442 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

Rev	Line
[2089]	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