PySDM_examples.Arabas_and_Pawlowska_2011.settings
1from PySDM import Formulae 2from PySDM.physics import si 3from PySDM.initialisation.spectra import Lognormal, Sum 4 5 6class Settings: 7 def __init__(self): 8 self.dt = 0.25 * si.s 9 self.mass_of_dry_air = 1000 * si.kg 10 self.p0 = 1000 * si.hPa 11 self.RH0 = 0.99 * si.dimensionless 12 self.T0 = 280 * si.K 13 self.w = 0.25 * si.m / si.s 14 15 self.initial_air_density = 1.245 * si.kg / si.m**3 16 17 self.output_interval = 4 18 self.output_points = 250 19 self.n_sd = 1024 20 21 self.formulae = Formulae() 22 23 self.cloud_range = (1 * si.um, 25 * si.um) 24 25 self.kappa_sea_salt = 1.28 * si.dimensionless 26 self.kappa_sulphate = 0.61 * si.dimensionless 27 28 sea_salt_spectrum = Sum( 29 ( 30 Lognormal( 31 norm_factor=51.1 / si.cm**3, 32 m_mode=0.10 * si.um, 33 s_geom=1.90, 34 ), 35 Lognormal( 36 norm_factor=2.21 / si.cm**3, 37 m_mode=1.00 * si.um, 38 s_geom=2.00, 39 ), 40 Lognormal( 41 norm_factor=1e-5 / si.cm**3, 42 m_mode=6.00 * si.um, 43 s_geom=3.00, 44 ), 45 ) 46 ) 47 48 sulphate_spectrum = Lognormal( 49 norm_factor=100 / si.cm**3, 50 m_mode=0.08 * si.um, 51 s_geom=1.45, 52 ) 53 54 self.total_aerosol_concentration = ( 55 sea_salt_spectrum.norm_factor + sulphate_spectrum.norm_factor 56 ) 57 58 self.aerosol_modes_by_kappa = { 59 self.kappa_sea_salt: sea_salt_spectrum, 60 self.kappa_sulphate: sulphate_spectrum, 61 } 62 63 self.aerosol_mode_names = ( 64 "sea salt", 65 "sulphate", 66 )
class
Settings:
7class Settings: 8 def __init__(self): 9 self.dt = 0.25 * si.s 10 self.mass_of_dry_air = 1000 * si.kg 11 self.p0 = 1000 * si.hPa 12 self.RH0 = 0.99 * si.dimensionless 13 self.T0 = 280 * si.K 14 self.w = 0.25 * si.m / si.s 15 16 self.initial_air_density = 1.245 * si.kg / si.m**3 17 18 self.output_interval = 4 19 self.output_points = 250 20 self.n_sd = 1024 21 22 self.formulae = Formulae() 23 24 self.cloud_range = (1 * si.um, 25 * si.um) 25 26 self.kappa_sea_salt = 1.28 * si.dimensionless 27 self.kappa_sulphate = 0.61 * si.dimensionless 28 29 sea_salt_spectrum = Sum( 30 ( 31 Lognormal( 32 norm_factor=51.1 / si.cm**3, 33 m_mode=0.10 * si.um, 34 s_geom=1.90, 35 ), 36 Lognormal( 37 norm_factor=2.21 / si.cm**3, 38 m_mode=1.00 * si.um, 39 s_geom=2.00, 40 ), 41 Lognormal( 42 norm_factor=1e-5 / si.cm**3, 43 m_mode=6.00 * si.um, 44 s_geom=3.00, 45 ), 46 ) 47 ) 48 49 sulphate_spectrum = Lognormal( 50 norm_factor=100 / si.cm**3, 51 m_mode=0.08 * si.um, 52 s_geom=1.45, 53 ) 54 55 self.total_aerosol_concentration = ( 56 sea_salt_spectrum.norm_factor + sulphate_spectrum.norm_factor 57 ) 58 59 self.aerosol_modes_by_kappa = { 60 self.kappa_sea_salt: sea_salt_spectrum, 61 self.kappa_sulphate: sulphate_spectrum, 62 } 63 64 self.aerosol_mode_names = ( 65 "sea salt", 66 "sulphate", 67 )