source: trunk/LMDZ.MARS/libf/phymars/microphys.h @ 1617

Last change on this file since 1617 was 1617, checked in by jaudouard, 8 years ago

Added modifications for CO2 clouds scheme in physiq_mod.F and added several routines and variables for CO2 microphysics. October 2016 Joachim Audouard
File size: 3.4 KB
Line 
1!-----------------------------------------------------------------------
2! INCLUDE 'microphys.h'
3! Parameters and physical constants used by the microphysal scheme;
4!-----------------------------------------------------------------------
5
6!     Number of bins
7      INTEGER, PARAMETER :: nbin_cld = 5
8
9!     Reference temperature, T=273.15 K
10      REAL, PARAMETER :: To = 273.15
11!     Avogadro number
12      DOUBLE PRECISION, PARAMETER :: nav = 6.023d23
13!     Perfect gas constant
14      DOUBLE PRECISION, PARAMETER :: rgp = 8.3143
15!     Boltzman constant
16      DOUBLE PRECISION, PARAMETER :: kbz = 1.381d-23
17!     Molecular weight of H2O (kg.mol-1)
18      DOUBLE PRECISION, PARAMETER :: mh2o = 18.d-3
19!     Molecular weight of CO2 (kg.mol-1)
20      DOUBLE PRECISION, PARAMETER :: mco2 = 44.d-3
21!     Molecular weight of N2 (kg.mol-1)
22      DOUBLE PRECISION, PARAMETER :: mn2 = 28.01d-3
23!     Effective CO2 gas molecular radius (m)
24      DOUBLE PRECISION, PARAMETER :: molco2 = 2.2d-10
25!     Effective H2O gas molecular radius (m)
26      DOUBLE PRECISION, PARAMETER :: molh2o = 1.2d-10
27!     Surface tension of ice/vapor (N.m)
28      DOUBLE PRECISION, PARAMETER :: sigh2o = 0.12
29!     Activation energy for desorption of
30!       water on a dust-like substrate
31!       (J/molecule)
32      DOUBLE PRECISION, PARAMETER :: desorp = 0.288e-19
33!     Jump frequency of a water molecule (s-1)
34      DOUBLE PRECISION, PARAMETER :: nus = 1.e+13
35!     Estimated activation energy for
36!       surface diffusion of water molecules
37!       (J/molecule)
38      DOUBLE PRECISION, PARAMETER :: surfdif = desorp / 10.
39!     Weight of a water molecule (kg)
40      DOUBLE PRECISION, PARAMETER :: m0 = mh2o / nav
41
42!     Contact parameter ( m=cos(theta) )
43!       (initialized in improvedclouds.F)
44      REAL mteta
45
46!     Volume of a water molecule (m3)
47      DOUBLE PRECISION vo1
48!     Radius used by the microphysical scheme (m)
49      DOUBLE PRECISION rad_cld(nbin_cld)
50
51
52
53
54!CO2 part
55!      number of bins for nucleation
56      INTEGER, PARAMETER :: nbinco2_cld=40
57!     Surface tension of ice/vapor (J.m-2)
58      DOUBLE PRECISION, PARAMETER :: sigco2 = 0.08
59!     Activation energy for desorption of
60!       water on a dust-like substrate
61!       (J/molecule)
62      DOUBLE PRECISION, PARAMETER :: desorpco2 = 3.25e-20
63!     Jump frequency of a co2 molecule (s-1)
64      DOUBLE PRECISION, PARAMETER :: nusco2 =  2.9e+12
65!     Estimated activation energy for
66!       surface diffusion of co2 molecules
67!       (J/molecule)
68      DOUBLE PRECISION, PARAMETER :: surfdifco2 = desorpco2 / 10.
69!     Weight of a co2 molecule (kg)
70      DOUBLE PRECISION, PARAMETER :: m0co2 = mco2 / nav
71!     Contact parameter ( m=cos(theta) )
72!       (initialized in improvedCO2clouds.F)
73       REAL, parameter :: mtetaco2 = 0.952
74!     Volume of a co2 molecule (m3)
75       DOUBLE PRECISION :: vo1co2
76!     Radius used by the microphysical scheme (m)
77      DOUBLE PRECISION :: rad_cldco2(nbinco2_cld)
78       REAL, parameter :: threshJA = 1.0
79!     COMMON/microphys/vo1co2,rad_cldco2
80
81! NB: to keep commons aligned: 
82!     split them in groups (reals, integers and characters)
83
84      COMMON/microphys/rad_cld,vo1,rad_cldco2,vo1co2
85                  COMMON/microphys_2/mteta
86     
87!     EXAMPLE:
88!     COMMON/tracer/radius,rho_q,alpha_lift,alpha_devil,mmol,           &
89!    & varian,r3n_q,rho_dust,rho_ice,nuice_ref,nuice_sed,               &
90!    & ref_r0,dryness
91!-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.