Ignore:
Timestamp:
Oct 17, 2023, 9:40:14 AM (15 months ago)
Author:
slebonnois
Message:

BdeBatz? : Cleans microphysics and makes few corrections for physics

Location:
trunk/LMDZ.TITAN/libf/muphytitan/config
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.TITAN/libf/muphytitan/config/config.ne15.cfg

    r1897 r3090  
    1919#    2 - SF interactions
    2020#    4 - FF interactions.
    21 # (for example: 5 = 4+1 --> SS and FF coagulation only)
     21# (for example : 5 = 4+1 --> SS and FF coagulation only)
    2222haze_coag_interactions = 7
    2323# Enable/disable Haze sedimentation process
    2424haze_sedimentation     = T
    2525# Disable Fiadero correction for sedimentation process
    26 no_fiadero             = F
     26no_fiadero             = T
    2727# Fiadero correction minimum ratio threshold
    2828fiadero_min_ratio      = 0.1
     
    3030fiadero_max_ratio      = 10.
    3131# Force settling velocity to M0
    32 wsed_m0 = T
     32wsed_m0                = T
    3333# Force settling velocity to M3
    34 wsed_m3 = F
     34wsed_m3                = F
    3535# Enable/disable clouds sedimentation process
    3636# (automatically set to F if clouds microphysics is not enabled)
     
    4141# Condensible species configuration file
    4242# (not needed if clouds microphysics is not enabled)
    43 specie_cfg             = ../datagcm/microphysics/mp2m_species.cfg
     43specie_cfg             = ../../datagcm/microphysics/mp2m_species.cfg
    4444
    4545# Enable/disable spherical mode transfert probability
     
    4747# Path of the spherical mode transfert probability look-up tables file
    4848# (optional if 'transfert_probability' is False)
    49 ps2s_file             = ../datagcm/microphysics/mmp_ps2s_rm50_ne15.nc
     49ps2s_file             = ../../datagcm/microphysics/mmp_ps2s_rm50_ne15.nc
    5050
    5151# Electric charging coagulation correction
     
    5454# Path of the electric charging correction factor.
    5555# (optional if 'electric_charging' is False)
    56 mq_file               = ../datagcm/microphysics/mmp_qmean_rm50_ne15.nc
     56mq_file               = ../../datagcm/microphysics/mmp_qmean_rm50_ne15.nc
    5757
    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
     58### Thresholds related parameters
     59### -----------------------------
     60# Aerosol spherical mode total number min. threshold
     61m0as_min = 1e-10
     62# Aerosol spherical mode min. caracteristic radius threshold
     63rcs_min  = 1e-9
     64# Aerosol fractal mode total number min. threshold
     65m0af_min = 1e-10
     66# Aerosol fractal mode min. caracteristic radius threshold (updated to monomer radius at runtime if needed)
     67rcf_min  = 1e-9
     68# cloud drop total number min. threshold
     69m0n_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))
    6176[alpha_s]
    6277a =   1.5044478E-02,  -2.0948806E-01,  -1.5824302E+02,   1.1597818E-01,   9.9502283E-02,  -1.1223796E-01
    6378b =  -2.8622255E-01,   7.7089599E+00,  -1.7000626E+02,   2.6012143E+00,   5.5138784E-01,   9.2024747E-01
    6479c =  -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
     80
    7381[alpha_f]
    7482a =   1.5044478E-02,  -2.0948806E-01,  -1.5824302E+02,   1.1597818E-01,   9.9502283E-02,  -1.1223796E-01
    7583b =  -2.8622255E-01,   7.7089599E+00,  -1.7000626E+02,   2.6012143E+00,   5.5138784E-01,   9.2024747E-01
    7684c =  -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
    8585
    8686# ================= #
    8787# b^T_k cofficients #
    8888# ================= #
    89 # This section gathers the values of all the btk coefficient used in the coagulation
    90 # equations for the free-molecular regime.
     89# This section gathers the values of all the btk coefficient used in the coagulation equations for the free-molecular regime
    9190[btks]
    9291bt0 = 0.73d0, 0.73d0, 0.75d0, 0.99d0, 0.00d0
    9392bt3 = 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 
  • trunk/LMDZ.TITAN/libf/muphytitan/config/mp2m_species.cfg

    r1793 r3090  
    1 #----------------------------------------------------------------------------------------
    2 # mp_scpecies.cfg
     1# ----------------------------------------------------------------------------------------
     2# mp2m_scpecies.cfg
    33# Thermodynamic properties of chemical species used in the cloud microphysics processes.
    44# Each condensible specie that should be treated by the model must be entirely described
    55# here.
    66#
    7 # Air specie is always mandatory. It is saved in a special variable in the model and is
    8 # never used as cloud condensible specie.
    9 #-----------------------------------------------------------------------------------------
     7# ----------------------------------------------------------------------------------------
    108# PARAMETER | DESCRIPTION
    119# ----------|-----------------------------------------------------------------------------
    1210#  name     | Specie name, should be the same as the section name
    13 #  mas      | Molecular weight
    14 #  vol      | Molecular volume
    15 #  ray      | Molecular radius
    16 #  masmol   | Molar mass
    17 #  rho      | Density
    18 #  tc       | Critical temperature
    19 #  tb       | Boiling temperature
    20 #  pc       | Critical pressure (in bar !!!)
     11#  mas      | Molecular weight      [kg]
     12#  vol      | Molecular volume      [m3]
     13#  ray      | Molecular radius      [m]
     14#  masmol   | Molar mass            [kg.mol-1]
     15#  rho      | Density               [kg.m-3]
     16#  tc       | Critical temperature  [K]
     17#  tb       | Boiling temperature   [K]
     18#  pc       | Critical pressure     [bar]
     19#  tx_prod  | Production rate       [kg.m-2.s-1] (actually this is not used by the model)
    2120#  w        | Acentric factor
    22 #  a_sat    | Coefficient A of Psat equation (from Reid et al. 1986).
    23 #  b_sat    | Coefficient B of Psat equation (from Reid et al. 1986).
    24 #  c_sat    | Coefficient C of Psat equation (from Reid et al. 1986).
    25 #  d_sat    | Coefficient D of Psat equation (from Reid et al. 1986).
     21#  a_sat    | Coefficient A of Psat equation (from Reid et al. 1986)
     22#  b_sat    | Coefficient B of Psat equation (from Reid et al. 1986)
     23#  c_sat    | Coefficient C of Psat equation (from Reid et al. 1986)
     24#  d_sat    | Coefficient D of Psat equation (from Reid et al. 1986)
    2625#  mteta    | Wettability (free parameter in most cases, from 0 to 1)
    27 #  tx_prod  | Production rate (actually this is not used by the model)
    28 #-----------------------------------------------------------------------------------------
     26# -----------------------------------------------------------------------------------------
    2927
    30 ### List of actual species to be used in the model:
    31 ### WARNING : the list of species specified here should be ordered:
    32 ###    In the model, ice tracers as well as condensible gazs species must have the same
    33 ###    index.
    34 used_species = "CH4", "C2H6", "C2H2"
    35 
    36 # AIR properties (MANDATORY !!!)
    37 ################################
    38 [air]
    39 name    = "air"
    40 mas     = 4.650e-26
    41 vol     = 5.750e-29
    42 ray     = 1.750e-10
    43 masmol  = 28.e-3
    44 rho     = 808.6
    45 tc      = 126.2
    46 tb      = 77.4
    47 pc      = 33.9
    48 w       = 3.9e-2
    49 a_sat   = -6.09676
    50 b_sat   = 1.13670
    51 c_sat   = -1.04072
    52 d_sat   = -1.93306
    53 mteta   = 0.
    54 tx_prod = 0.
     28### List of actual species to be used in the model :
     29### WARNING : The list of species specified here should be ordered :
     30###           In the model, ice tracers as well as condensible gazs species must
     31###           have the same index.
     32used_species = CH4, C2H6, C2H2, HCN
    5533
    5634# CH4 properties (useful for Titan :)
    5735#####################################
    5836[CH4]
    59 name    = "CH4"
     37name    = CH4
    6038mas     = 2.6578e-26
    6139vol     = 6.252e-29
    6240ray     = 2.000e-10
    63 masmol  = 16.e-3
     41masmol  = 16.04e-3
    6442rho     = 425.
    65 tc      = 190.4
    66 tb      = 111.6
    67 pc      = 46.0
     43tc      = 190.551
     44tb      = 111.66
     45pc      = 45.992
    6846w       = 1.1e-2
    69 a_sat   = -6.00435
    70 b_sat   = 1.18850
    71 c_sat   = -0.83408
    72 d_sat   = -1.22833
    73 mteta   = 0.92
     47a_sat   = -6.02242
     48b_sat   = 1.26652
     49c_sat   = -0.5707
     50d_sat   = -1.366
     51mteta   = 0.99
    7452tx_prod = 0.
    7553
     
    7755#################
    7856[C2H6]
    79 name    = "C2H6"
     57name    = C2H6
    8058mas     = 4.983e-26
    8159vol     = 9.094e-29
    8260ray     = 2.220e-10
    83 masmol  = 30.e-3
     61masmol  = 30.07e-3
    8462rho     = 544.6
    85 tc      = 305.4
    86 tb      = 184.6
    87 pc      = 48.8
     63tc      = 305.33
     64tb      = 184.55
     65pc      = 48.71
    8866w       = 9.9e-2
    89 a_sat   = -6.34307
    90 b_sat   = 1.01163
    91 c_sat   = -1.19116
    92 d_sat   = -2.03539
    93 mteta   = 0.92
     67a_sat   = -6.47500
     68b_sat   = 1.41071
     69c_sat   = -1.1440
     70d_sat   = -1.8590
     71mteta   = 0.99
    9472tx_prod = 1.2e-12
    9573
     
    9775#################
    9876[C2H2]
    99 name    = "C2H2"
     77name    = C2H2
    10078mas     = 4.319e-26
    10179vol     = 7.020e-29
    10280ray     = 2.015e-10
    103 masmol  = 26.e-3
     81masmol  = 26.04e-3
    10482rho     = 615.
    105 tc      = 308.8
    106 tb      = 188.4
    107 pc      = 61.4
     83tc      = 308.35
     84tb      = 188.40
     85pc      = 61.39
    10886w       = 19.0e-2
    109 a_sat   = -6.90128
    110 b_sat   = 1.26873
    111 c_sat   = -2.09113
    112 d_sat   = -2.75601
    113 mteta   = 0.92
     87a_sat   = -6.87886
     88b_sat   = 1.30164
     89c_sat   = -1.22474
     90d_sat   = -3.59556
     91mteta   = 0.99
    11492tx_prod = 3.2e-13
    11593
     
    11795#################
    11896[HCN]
    119 name    = "HCN"
     97name    = HCN
    12098mas     = 4.484e-26
    12199vol     = 6.498e-29
     
    131109c_sat   = -3.004
    132110d_sat   = 1635.
    133 mteta   = 0.92
     111mteta   = 0.99
    134112tx_prod = 1e-12
    135 
Note: See TracChangeset for help on using the changeset viewer.