-
Notifications
You must be signed in to change notification settings - Fork 7
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Add cloud condensate sedimentation velocity #462
Conversation
Codecov ReportAll modified and coverable lines are covered by tests ✅
Additional details and impacted files@@ Coverage Diff @@
## main #462 +/- ##
==========================================
+ Coverage 96.95% 96.98% +0.03%
==========================================
Files 41 41
Lines 1543 1559 +16
==========================================
+ Hits 1496 1512 +16
Misses 47 47
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me. Just minor cosmetic comments. Thank you for the thorough documentation. Makes it so much easier to review the code too!
src/Microphysics1M.jl
Outdated
@@ -187,37 +185,31 @@ function terminal_velocity( | |||
q::FT, | |||
) where {FT} | |||
fall_w = FT(0) | |||
# For now we assume B4 table coefficiantes for snow |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
# For now we assume B4 table coefficiantes for snow | |
# For now, we assume the B4 table coefficients for snow |
src/Microphysics1M.jl
Outdated
@@ -187,37 +185,31 @@ function terminal_velocity( | |||
q::FT, | |||
) where {FT} | |||
fall_w = FT(0) | |||
# For now we assume B4 table coefficiantes for snow | |||
# and B2 table coefficients for cloud ice. | |||
# TODO - we should do an integral by parts |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
# TODO - we should do an integral by parts | |
# TODO - we should do partial integrals |
src/Microphysics1M.jl
Outdated
a0c = a0 * χa | ||
mec = me + Δm | ||
aec = ae + Δa | ||
# As a next step could keep ϕ(r) under the integrals |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
# As a next step could keep ϕ(r) under the integrals | |
# As a next step, we could keep ϕ(r) under the integrals |
src/MicrophysicsNonEq.jl
Outdated
|
||
Returns the mass weighted average terminal velocity assuming a | ||
monodisperse size distribution with prescribed number concentration. | ||
Fall velocity of individual particles is parameterized following |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Fall velocity of individual particles is parameterized following | |
The fall velocity of individual particles is parameterized following |
src/MicrophysicsNonEq.jl
Outdated
# that is more suited for cloud droplets. For now I'm just multiplying | ||
# by an arbitrary correction factor. | ||
aiu, bi, ciu = CO.Chen2022_vel_coeffs_B1(vel, ρ) | ||
# The 1-m scheme does not assume any cloud droplet size distribution. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
# The 1-m scheme does not assume any cloud droplet size distribution. | |
# The 1M scheme does not assume any cloud droplet size distribution. |
Co-authored-by: Tapio Schneider <tapio@caltech.edu>
5f6914f
to
06e1a1e
Compare
TODO: