fix: missing (1+qv) factor in calculation of derivative of rho_d when…

[yoctoyotta1024]

… accounting for water vapour

rho_d(z) = rho(z) / (1 + qv(z))
therefore drho_d/dz = drho/dz / (1+ qv) - dqv/dz / (1 + qv)^2
not(!) drho_d/dz = drho/dz / (1+ qv) - dqv/dz / (1 + qv)

[codecov]

Codecov Report

✅ All modified and coverable lines are covered by tests.
✅ Project coverage is 92.69%. Comparing base (a921acc) to head (f0c5b67).
⚠️ Report is 5 commits behind head on main.

Additional details and impacted files

@@           Coverage Diff           @@
##             main     #614   +/-   ##
=======================================
  Coverage   92.69%   92.69%           
=======================================
  Files          39       39           
  Lines        1697     1697           
=======================================
  Hits         1573     1573           
  Misses        124      124           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:

• ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

[github-actions]

Stale pull request message

[yoctoyotta1024]

@slayoo did this ever get merged? Also in PySDM?

[AgnieszkaZaba]

@yoctoyotta1024
In PySDM I found this: https://github.com/open-atmos/PySDM/blob/e4ddacdbf7fdbb01a3e43ab947a152ac55bdcf92/examples/PySDM_examples/Shipway_and_Hill_2012/settings.py#L131-L133

Maybe we can do similar test to this: https://github.com/open-atmos/PySDM/blob/main/tests/smoke_tests/kinematic_1d/shipway_and_hill_2012/test_initial_condition.py - test_density_profile

[yoctoyotta1024]

@AgnieszkaZaba Ah yes nice spot thank you, I had forgotten about this 😅
Good idea for the test too.