Skip to content
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

Merged
merged 1 commit into from
Oct 7, 2024

Conversation

trontrytel
Copy link
Member

@trontrytel trontrytel commented Oct 2, 2024

TODO:

  • confirm that the values we are getting for cloud liquid and ice are reasonable
  • add unit tests
  • update the documentation plots with new function names
  • add plots with all velocities with v(D) and group velocities for comparison

Copy link

codecov bot commented Oct 3, 2024

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 96.98%. Comparing base (0c3e082) to head (06e1a1e).
Report is 1 commits behind head on main.

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              
Components Coverage Δ
src 98.76% <100.00%> (+0.01%) ⬆️
ext 69.79% <ø> (ø)

@trontrytel trontrytel self-assigned this Oct 5, 2024
@trontrytel trontrytel added the enhancement New feature or request label Oct 5, 2024
Copy link
Contributor

@tapios tapios left a 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!

docs/src/Microphysics1M.md Outdated Show resolved Hide resolved
docs/src/Microphysics1M.md Outdated Show resolved Hide resolved
docs/src/Microphysics1M.md Outdated Show resolved Hide resolved
docs/src/Microphysics1M.md Outdated Show resolved Hide resolved
docs/src/Microphysics1M.md Outdated Show resolved Hide resolved
@@ -187,37 +185,31 @@ function terminal_velocity(
q::FT,
) where {FT}
fall_w = FT(0)
# For now we assume B4 table coefficiantes for snow
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
# For now we assume B4 table coefficiantes for snow
# For now, we assume the B4 table coefficients for snow

@@ -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
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
# TODO - we should do an integral by parts
# TODO - we should do partial integrals

a0c = a0 * χa
mec = me + Δm
aec = ae + Δa
# As a next step could keep ϕ(r) under the integrals
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
# As a next step could keep ϕ(r) under the integrals
# As a next step, we could keep ϕ(r) under the integrals


Returns the mass weighted average terminal velocity assuming a
monodisperse size distribution with prescribed number concentration.
Fall velocity of individual particles is parameterized following
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
Fall velocity of individual particles is parameterized following
The fall velocity of individual particles is parameterized following

# 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.
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
# 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>
@trontrytel trontrytel merged commit f225c60 into main Oct 7, 2024
12 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants