| 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 = T |
|---|
| 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 | ### Thresholds related parameters |
|---|
| 59 | ### ----------------------------- |
|---|
| 60 | # Aerosol spherical mode total number min. threshold |
|---|
| 61 | m0as_min = 1e-10 |
|---|
| 62 | # Aerosol spherical mode min. caracteristic radius threshold |
|---|
| 63 | rcs_min = 1e-9 |
|---|
| 64 | # Aerosol fractal mode total number min. threshold |
|---|
| 65 | m0af_min = 1e-10 |
|---|
| 66 | # Aerosol fractal mode min. caracteristic radius threshold (updated to monomer radius at runtime if needed) |
|---|
| 67 | rcf_min = 1e-9 |
|---|
| 68 | # cloud drop total number min. threshold |
|---|
| 69 | m0n_min = 1e-10 |
|---|
| 70 | |
|---|
| 71 | |
|---|
| 72 | # =================== # |
|---|
| 73 | # alpha_k cofficients # |
|---|
| 74 | # =================== # |
|---|
| 75 | # alpha_X sections contain the parameters of the inter-moments relation function for the mode X (spherical (s) or fractal (f)) |
|---|
| 76 | [alpha_s] |
|---|
| 77 | a = 1.5044478E-02, -2.0948806E-01, -1.5824302E+02, 1.1597818E-01, 9.9502283E-02, -1.1223796E-01 |
|---|
| 78 | b = -2.8622255E-01, 7.7089599E+00, -1.7000626E+02, 2.6012143E+00, 5.5138784E-01, 9.2024747E-01 |
|---|
| 79 | c = -3.0205020E-02, -3.5510239E+01, -2.0306468E+02, -1.3605159E+01, -4.1653422E+00, -4.2571698E+00 |
|---|
| 80 | |
|---|
| 81 | [alpha_f] |
|---|
| 82 | a = 1.5044478E-02, -2.0948806E-01, -1.5824302E+02, 1.1597818E-01, 9.9502283E-02, -1.1223796E-01 |
|---|
| 83 | b = -2.8622255E-01, 7.7089599E+00, -1.7000626E+02, 2.6012143E+00, 5.5138784E-01, 9.2024747E-01 |
|---|
| 84 | c = -3.0205020E-02, -3.5510239E+01, -2.0306468E+02, -1.3605159E+01, -4.1653422E+00, -4.2571698E+00 |
|---|
| 85 | |
|---|
| 86 | # ================= # |
|---|
| 87 | # b^T_k cofficients # |
|---|
| 88 | # ================= # |
|---|
| 89 | # This section gathers the values of all the btk coefficient used in the coagulation equations for the free-molecular regime |
|---|
| 90 | [btks] |
|---|
| 91 | bt0 = 0.73d0, 0.73d0, 0.75d0, 0.99d0, 0.00d0 |
|---|
| 92 | bt3 = 0.97d0, 0.97d0, 0.00d0, 0.99d0, 0.99d0 |
|---|
