| 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 | ! bachnar 2016 value :1.97d-10 ! old value = 2.2d-10 |
|---|
| 25 | DOUBLE PRECISION, PARAMETER :: molco2 = 1.97d-10 |
|---|
| 26 | ! Effective H2O gas molecular radius (m) |
|---|
| 27 | DOUBLE PRECISION, PARAMETER :: molh2o = 1.2d-10 |
|---|
| 28 | ! Surface tension of ice/vapor (N.m) |
|---|
| 29 | DOUBLE PRECISION, PARAMETER :: sigh2o = 0.12 |
|---|
| 30 | ! Activation energy for desorption of |
|---|
| 31 | ! water on a dust-like substrate |
|---|
| 32 | ! (J/molecule) |
|---|
| 33 | DOUBLE PRECISION, PARAMETER :: desorp = 0.288e-19 |
|---|
| 34 | ! Jump frequency of a water molecule (s-1) |
|---|
| 35 | DOUBLE PRECISION, PARAMETER :: nus = 1.e+13 |
|---|
| 36 | ! Estimated activation energy for |
|---|
| 37 | ! surface diffusion of water molecules |
|---|
| 38 | ! (J/molecule) |
|---|
| 39 | DOUBLE PRECISION, PARAMETER :: surfdif = desorp / 10. |
|---|
| 40 | ! Weight of a water molecule (kg) |
|---|
| 41 | DOUBLE PRECISION, PARAMETER :: m0 = mh2o / nav |
|---|
| 42 | |
|---|
| 43 | ! Contact parameter ( m=cos(theta) ) |
|---|
| 44 | ! (initialized in improvedclouds.F) |
|---|
| 45 | REAL mteta |
|---|
| 46 | |
|---|
| 47 | ! Volume of a water molecule (m3) |
|---|
| 48 | DOUBLE PRECISION vo1 |
|---|
| 49 | ! Radius used by the microphysical scheme (m) |
|---|
| 50 | DOUBLE PRECISION rad_cld(nbin_cld) |
|---|
| 51 | |
|---|
| 52 | |
|---|
| 53 | |
|---|
| 54 | |
|---|
| 55 | !CO2 part |
|---|
| 56 | ! number of bins for nucleation |
|---|
| 57 | INTEGER, PARAMETER :: nbinco2_cld=100 |
|---|
| 58 | ! Surface tension of ice/vapor (J.m-2) |
|---|
| 59 | DOUBLE PRECISION, PARAMETER :: sigco2 = 0.08 |
|---|
| 60 | ! Activation energy for desorption of |
|---|
| 61 | ! water on a dust-like substrate |
|---|
| 62 | ! (J/molecule) |
|---|
| 63 | DOUBLE PRECISION, PARAMETER ::desorpco2=3.07e-20 |
|---|
| 64 | ! bachnar 2016 value 3.07d-20 |
|---|
| 65 | !old value 3.20e-20 |
|---|
| 66 | ! Jump frequency of a co2 molecule (s-1) |
|---|
| 67 | DOUBLE PRECISION, PARAMETER :: nusco2 = 2.9e+12 |
|---|
| 68 | ! Estimated activation energy for |
|---|
| 69 | ! surface diffusion of co2 molecules |
|---|
| 70 | ! (J/molecule) |
|---|
| 71 | DOUBLE PRECISION, PARAMETER :: surfdifco2 = desorpco2 / 10. |
|---|
| 72 | ! Weight of a co2 molecule (kg) |
|---|
| 73 | DOUBLE PRECISION, PARAMETER :: m0co2 = mco2 / nav |
|---|
| 74 | ! Contact parameter ( m=cos(theta) ) |
|---|
| 75 | ! (initialized in improvedCO2clouds.F) |
|---|
| 76 | ! bachnar 2016 value :0.78 |
|---|
| 77 | !old value 0.952 |
|---|
| 78 | REAL, parameter :: mtetaco2 = 0.952 |
|---|
| 79 | ! Volume of a co2 molecule (m3) |
|---|
| 80 | DOUBLE PRECISION :: vo1co2 |
|---|
| 81 | ! Radius used by the microphysical scheme (m) |
|---|
| 82 | DOUBLE PRECISION :: rad_cldco2(nbinco2_cld) |
|---|
| 83 | REAL, parameter :: threshJA = 1 |
|---|
| 84 | ! COMMON/microphys/vo1co2,rad_cldco2 |
|---|
| 85 | |
|---|
| 86 | ! NB: to keep commons aligned: |
|---|
| 87 | ! split them in groups (reals, integers and characters) |
|---|
| 88 | |
|---|
| 89 | COMMON/microphys/rad_cld,vo1,rad_cldco2,vo1co2 |
|---|
| 90 | COMMON/microphys_2/mteta |
|---|
| 91 | |
|---|
| 92 | ! EXAMPLE: |
|---|
| 93 | ! COMMON/tracer/radius,rho_q,alpha_lift,alpha_devil,mmol, & |
|---|
| 94 | ! & varian,r3n_q,rho_dust,rho_ice,nuice_ref,nuice_sed, & |
|---|
| 95 | ! & ref_r0,dryness |
|---|
| 96 | !----------------------------------------------------------------------- |
|---|
