def __init__(self, settings, products=None):
env = Parcel(
dt=settings.dt,
mass_of_dry_air=settings.mass_of_dry_air,
p0=settings.p0,
q0=settings.q0,
T0=settings.T0,
w=settings.w,
)
builder = Builder(n_sd=settings.n_sd,
backend=CPU(formulae=settings.formulae))
builder.set_environment(env)
attributes = env.init_attributes(n_in_dv=settings.n_in_dv,
kappa=settings.kappa,
r_dry=settings.r_dry)
attributes = {
**attributes,
**settings.starting_amounts,
"pH": np.zeros(settings.n_sd),
}
builder.add_dynamic(AmbientThermodynamics())
builder.add_dynamic(Condensation())
builder.add_dynamic(
AqueousChemistry(
environment_mole_fractions=settings.ENVIRONMENT_MOLE_FRACTIONS,
system_type=settings.system_type,
n_substep=settings.n_substep,
dry_rho=settings.DRY_RHO,
dry_molar_mass=settings.dry_molar_mass,
))
products = products or (
PySDM_products.AmbientRelativeHumidity(name="RH", unit="%"),
PySDM_products.WaterMixingRatio(
name="ql", radius_range=[1 * si.um, np.inf], unit="g/kg"),
PySDM_products.ParcelDisplacement(name="z"),
PySDM_products.AmbientPressure(name="p"),
PySDM_products.AmbientTemperature(name="T"),
PySDM_products.AmbientDryAirDensity(name="rhod"),
PySDM_products.AmbientWaterVapourMixingRatio(name="qv",
unit="g/kg"),
PySDM_products.Time(name="t"),
*(PySDM_products.AqueousMoleFraction(
comp, unit="ppb", name=f"aq_{comp}_ppb")
for comp in AQUEOUS_COMPOUNDS),
*(PySDM_products.GaseousMoleFraction(
comp, unit="ppb", name=f"gas_{comp}_ppb")
for comp in GASEOUS_COMPOUNDS),
PySDM_products.Acidity(
name="pH_pH_number_weighted",
radius_range=settings.cloud_radius_range,
weighting="number",
attr="pH",
),
PySDM_products.Acidity(
name="pH_pH_volume_weighted",
radius_range=settings.cloud_radius_range,
weighting="volume",
attr="pH",
),
PySDM_products.Acidity(
name="pH_conc_H_number_weighted",
radius_range=settings.cloud_radius_range,
weighting="number",
attr="conc_H",
),
PySDM_products.Acidity(
name="pH_conc_H_volume_weighted",
radius_range=settings.cloud_radius_range,
weighting="volume",
attr="conc_H",
),
PySDM_products.TotalDryMassMixingRatio(
settings.DRY_RHO, name="q_dry", unit="ug/kg"),
PySDM_products.PeakSupersaturation(unit="%", name="S_max"),
PySDM_products.ParticleSpecificConcentration(
radius_range=settings.cloud_radius_range,
name="n_c_mg",
unit="mg^-1"),
PySDM_products.AqueousMassSpectrum(
key="S_VI",
dry_radius_bins_edges=settings.dry_radius_bins_edges,
name="dm_S_VI/dlog_10(dry diameter)",
unit='ug / m^3"',
),
)
particulator = builder.build(attributes=attributes, products=products)
self.settings = settings
super().__init__(particulator=particulator)
def __init__(self, settings, products=None):
env = Parcel(dt=settings.dt,
mass_of_dry_air=settings.mass_of_dry_air,
p0=settings.p0,
q0=settings.q0,
T0=settings.T0,
w=settings.w,
g=settings.g)
builder = Builder(n_sd=settings.n_sd,
backend=CPU,
formulae=settings.formulae)
builder.set_environment(env)
attributes = env.init_attributes(n_in_dv=settings.n_in_dv,
kappa=settings.kappa,
r_dry=settings.r_dry)
attributes = {
**attributes,
**settings.starting_amounts, 'pH': np.zeros(settings.n_sd)
}
builder.add_dynamic(AmbientThermodynamics())
builder.add_dynamic(Condensation(kappa=settings.kappa))
builder.add_dynamic(
AqueousChemistry(settings.ENVIRONMENT_MOLE_FRACTIONS,
system_type=settings.system_type,
n_substep=settings.n_substep,
dry_rho=settings.DRY_RHO,
dry_molar_mass=settings.dry_molar_mass))
products = products or (
PySDM_products.RelativeHumidity(),
PySDM_products.WaterMixingRatio(name='ql',
description_prefix='liquid',
radius_range=[1 * si.um, np.inf]),
PySDM_products.ParcelDisplacement(), PySDM_products.Pressure(),
PySDM_products.Temperature(), PySDM_products.DryAirDensity(),
PySDM_products.WaterVapourMixingRatio(), PySDM_products.Time(), *[
PySDM_products.AqueousMoleFraction(compound)
for compound in AQUEOUS_COMPOUNDS.keys()
], *[
PySDM_products.GaseousMoleFraction(compound)
for compound in GASEOUS_COMPOUNDS.keys()
],
PySDM_products.pH(radius_range=settings.cloud_radius_range,
weighting='number',
attr='pH'),
PySDM_products.pH(radius_range=settings.cloud_radius_range,
weighting='volume',
attr='pH'),
PySDM_products.pH(radius_range=settings.cloud_radius_range,
weighting='number',
attr='conc_H'),
PySDM_products.pH(radius_range=settings.cloud_radius_range,
weighting='volume',
attr='conc_H'),
PySDM_products.TotalDryMassMixingRatio(
settings.DRY_RHO), PySDM_products.PeakSupersaturation(),
PySDM_products.CloudDropletConcentration(
radius_range=settings.cloud_radius_range),
PySDM_products.AqueousMassSpectrum("S_VI",
settings.dry_radius_bins_edges))
self.core = builder.build(attributes=attributes, products=products)
self.settings = settings