| 1 | # YAMMS model parameters |
|---|
| 2 | # ====================== |
|---|
| 3 | |
|---|
| 4 | ### Model global parameters |
|---|
| 5 | ###------------------------ |
|---|
| 6 | # Monomer radius (m) |
|---|
| 7 | rm = 50e-9 |
|---|
| 8 | # Fractal dimension |
|---|
| 9 | df = 2.0 |
|---|
| 10 | # Aerosol density (kg.m-3) |
|---|
| 11 | rho_aer = 1000. |
|---|
| 12 | # Enable/disable Haze production process |
|---|
| 13 | haze_production = T |
|---|
| 14 | # Enable/disable Haze coagulation process |
|---|
| 15 | haze_coagulation = T |
|---|
| 16 | # Coagulation interactions, a combination of: |
|---|
| 17 | # 0 - no interactions (same as haze_coagulation == F) |
|---|
| 18 | # 1 - SS interactions |
|---|
| 19 | # 2 - SF interactions |
|---|
| 20 | # 4 - FF interactions. |
|---|
| 21 | # (for example: 5 = 4+1 --> SS and FF coagulation only) |
|---|
| 22 | haze_coag_interactions = 7 |
|---|
| 23 | # Enable/disable Haze sedimentation process |
|---|
| 24 | haze_sedimentation = T |
|---|
| 25 | # Disable Fiadero correction for sedimentation process |
|---|
| 26 | no_fiadero = F |
|---|
| 27 | # Fiadero correction minimum ratio threshold |
|---|
| 28 | fiadero_min_ratio = 0.1 |
|---|
| 29 | # Fiadero correction maximum ratio threshold |
|---|
| 30 | fiadero_max_ratio = 10. |
|---|
| 31 | # Force settling velocity to M0 |
|---|
| 32 | wsed_m0 = T |
|---|
| 33 | # Force settling velocity to M3 |
|---|
| 34 | wsed_m3 = F |
|---|
| 35 | # Enable/disable clouds sedimentation process |
|---|
| 36 | # (automatically set to F if clouds microphysics is not enabled) |
|---|
| 37 | clouds_sedimentation = T |
|---|
| 38 | # Enable/disable clouds nucleation and condensation processes |
|---|
| 39 | # (automatically set to F if clouds microphysics is not enabled) |
|---|
| 40 | clouds_nuc_cond = T |
|---|
| 41 | # Condensible species configuration file |
|---|
| 42 | # (not needed if clouds microphysics is not enabled) |
|---|
| 43 | specie_cfg = ../datagcm/microphysics/mp2m_species.cfg |
|---|
| 44 | |
|---|
| 45 | # Enable/disable spherical mode transfert probability |
|---|
| 46 | transfert_probability = T |
|---|
| 47 | # Path of the spherical mode transfert probability look-up tables file |
|---|
| 48 | # (optional if 'transfert_probability' is False) |
|---|
| 49 | ps2s_file = ../datagcm/microphysics/mmp_ps2s_rm50_ne15.nc |
|---|
| 50 | |
|---|
| 51 | # Electric charging coagulation correction |
|---|
| 52 | # If set to .false. then no correction is assumed |
|---|
| 53 | electric_charging = T |
|---|
| 54 | # Path of the electric charging correction factor. |
|---|
| 55 | # (optional if 'electric_charging' is False) |
|---|
| 56 | mq_file = ../datagcm/microphysics/mmp_qmean_rm50_ne15.nc |
|---|
| 57 | |
|---|
| 58 | # alpha_X sections contain the parameters of the inter-moments relation function for |
|---|
| 59 | # the mode X: either spherical (s) or fractal (f) |
|---|
| 60 | # dndr_X sections contain the parameters of the size-distribution law of the mode X |
|---|
| 61 | [alpha_s] |
|---|
| 62 | a = 1.5044478E-02, -2.0948806E-01, -1.5824302E+02, 1.1597818E-01, 9.9502283E-02, -1.1223796E-01 |
|---|
| 63 | b = -2.8622255E-01, 7.7089599E+00, -1.7000626E+02, 2.6012143E+00, 5.5138784E-01, 9.2024747E-01 |
|---|
| 64 | c = -3.0205020E-02, -3.5510239E+01, -2.0306468E+02, -1.3605159E+01, -4.1653422E+00, -4.2571698E+00 |
|---|
| 65 | [dndr_s] |
|---|
| 66 | rc = 4.58219580180634588E-007 |
|---|
| 67 | a0 = 86144.861255561875 |
|---|
| 68 | c = 0d0 |
|---|
| 69 | a = 2.48333861883769357E-040, 1.46076790655632173E-013, 1.71525517568997062E-009, |
|---|
| 70 | 1.80855172875974993E-019, 1.48212594918347503E-047, 6.87247318898338451E-081 |
|---|
| 71 | b = 59.518212357684796, 15.507500262021228, -5.4179933012448069, |
|---|
| 72 | -9.3500794017892854, -18.207927270524777, -27.248924688740562 |
|---|
| 73 | [alpha_f] |
|---|
| 74 | a = 1.5044478E-02, -2.0948806E-01, -1.5824302E+02, 1.1597818E-01, 9.9502283E-02, -1.1223796E-01 |
|---|
| 75 | b = -2.8622255E-01, 7.7089599E+00, -1.7000626E+02, 2.6012143E+00, 5.5138784E-01, 9.2024747E-01 |
|---|
| 76 | c = -3.0205020E-02, -3.5510239E+01, -2.0306468E+02, -1.3605159E+01, -4.1653422E+00, -4.2571698E+00 |
|---|
| 77 | [dndr_f] |
|---|
| 78 | rc = 4.58219580180634588E-007 |
|---|
| 79 | a0 = 86144.861255561875 |
|---|
| 80 | c = 0d0 |
|---|
| 81 | a = 2.48333861883769357E-040, 1.46076790655632173E-013, 1.71525517568997062E-009, |
|---|
| 82 | 1.80855172875974993E-019, 1.48212594918347503E-047, 6.87247318898338451E-081 |
|---|
| 83 | b = 59.518212357684796, 15.507500262021228, -5.4179933012448069, |
|---|
| 84 | -9.3500794017892854, -18.207927270524777, -27.248924688740562 |
|---|
| 85 | |
|---|
| 86 | # ================= # |
|---|
| 87 | # b^T_k cofficients # |
|---|
| 88 | # ================= # |
|---|
| 89 | # This section gathers the values of all the btk coefficient used in the coagulation |
|---|
| 90 | # equations for the free-molecular regime. |
|---|
| 91 | [btks] |
|---|
| 92 | bt0 = 0.73d0, 0.73d0, 0.75d0, 0.99d0, 0.00d0 |
|---|
| 93 | bt3 = 0.97d0, 0.97d0, 0.00d0, 0.99d0, 0.99d0 |
|---|
| 94 | |
|---|
| 95 | [optics] |
|---|
| 96 | optic_file = /path/to/optics_look_up_table.nc |
|---|
| 97 | |
|---|
| 98 | |
|---|
| 99 | |
|---|
