Sa/add aerosol activation #3872
Conversation
sajjadazimi
force-pushed
the
sa/add_aerosol_activation
branch
17 times, most recently
from
3f12ff2 to
b3eb3c5
Compare
sajjadazimi
force-pushed
the
sa/add_aerosol_activation
branch
8 times, most recently
from
8b62879 to
9c0eefc
Compare
…r 2M microphysics
sajjadazimi
force-pushed
the
sa/add_aerosol_activation
branch
from
9c0eefc to
f438d40
Compare
| seasalt_num, | ||
| seasalt_mean_radius, | ||
| sulfate_num, | ||
| specific(Y.c.ρq_tot, Y.c.ρ), |
There was a problem hiding this comment.
For what moisture and precipitation configurations this will be tested?
If we are aiming to swap out the 1M for the 2M scheme, then we should be passing q_tot, q_liq, q_ice, q_rai and q_sno. Precipitation is part of the working fluid, so q_vap = q_tot - q_liq - q_ice - q_rai - q_sno.
| dt, | ||
| ) | ||
| # Aerosol activation using prescribed aerosol (Sea salt and sulfate) | ||
| if !(:prescribed_aerosols_field in propertynames(p.tracers)) |
There was a problem hiding this comment.
Is it possible for prescribed_aerosols_field to not be empty, but not contain any of the seasalt or sulfate aerosol fields you are expecting?
| @. seasalt_num += seasalt_mass / seasalt_particle_mass | ||
| @. seasalt_mean_radius += | ||
| seasalt_mass / seasalt_particle_mass * | ||
| log(seasalt_particle_radius) |
There was a problem hiding this comment.
Or I guess I just dont understand how the mean radius is computed here
| aerosol_params = cmc.aerosol | ||
| T = Tₐ(thermo_params, ts) | ||
| p = Pₐ(thermo_params, ts) | ||
| S = CMTDI.supersaturation_over_liquid(thermo_params, qₜ, qₗ, qᵢ, ρ, T) |
There was a problem hiding this comment.
Here q_l should be equal to q_lcl + q_rai and q_i = q_icl + q_sno
|
Left some small comments:
Thank you! |
Purpose
Add aerosol activation based on Abdulrazzak and Ghan (2000).
To-do
Add a good test
Content