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

Last change on this file since 520 was 520, checked in by tnavarro, 13 years ago

10/02/12 == TN

Major update on watercycle: a smaller integration timestep is now used

in watercloud.F, sedimentation of clouds is done in watercloud instead of
callsedim.F

Temperature-dependant contact parameter in nuclea.F
No dust lifting if CO2 ice in vdif.c
Ice integrated column opacity is written in diagfi from physiq.F, instead

of aeropacity.F. Mandatory if iradia is not 1.

New definition of permanent ice in surfini.F and possibility to have an ice

cap in it in 1d.

Update in deftank: callphys.def.outliers,run.def.1d; added traceur.def.scavenging
File size: 2.2 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!     Effective CO2 gas molecular radius (m)
22      DOUBLE PRECISION, PARAMETER :: molco2 = 2.2d-10
23!     Effective H2O gas molecular radius (m)
24      DOUBLE PRECISION, PARAMETER :: molh2o = 1.2d-10
25!     Surface tension of ice/vapor (N.m)
26      DOUBLE PRECISION, PARAMETER :: sigh2o = 0.12
27!     Activation energy for desorption of
28!       water on a dust-like substrate
29!       (J/molecule)
30      DOUBLE PRECISION, PARAMETER :: desorp = 0.288e-19
31!     Jump frequency of a water molecule (s-1)
32      DOUBLE PRECISION, PARAMETER :: nus = 1.e+13
33!     Estimated activation energy for
34!       surface diffusion of water molecules
35!       (J/molecule)
36      DOUBLE PRECISION, PARAMETER :: surfdif = desorp / 10.
37!     Weight of a water molecule (kg)
38      DOUBLE PRECISION, PARAMETER :: m0 = mh2o / nav
39
40!     Contact parameter ( m=cos(theta) )
41!       (initialized in improvedclouds.F)
42      REAL mteta
43
44!     Volume of a water molecule (m3)
45      DOUBLE PRECISION vo1
46!     Radius used by the microphysical scheme (m)
47      DOUBLE PRECISION rad_cld(nbin_cld)
48
49
50! NB: to keep commons aligned: 
51!     split them in groups (reals, integers and characters)
52
53      COMMON/microphys/rad_cld,vo1
54     
55      COMMON/microphys_2/mteta
56     
57!     EXAMPLE:
58!     COMMON/tracer/radius,rho_q,alpha_lift,alpha_devil,mmol,           &
59!    & varian,r3n_q,rho_dust,rho_ice,nuice_ref,nuice_sed,               &
60!    & ref_r0,dryness
61!-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.