SpeedyWeather / SpeedyWeather.jl

Play atmospheric modelling like it's LEGO.
https://speedyweather.github.io/SpeedyWeather.jl/dev
MIT License
400 stars 24 forks source link

Correct vertical advection of humidity #556

Open simone-silvestri opened 2 weeks ago

simone-silvestri commented 2 weeks ago

I don't know why, at the moment, humidity is reconstructed differently from temperature and velocities when it comes to vertical advection.

This should not be the case (i.e., temperature and humidity should follow the same transport equation). I suspect this was just a mishap given that humidity is separated by the other variables because it is active only in a wet core. This PR corrects this "problem". I will try some tests to make sure everything is correct before merging

milankl commented 1 week ago

Thanks Simone, this is only relevant for WENO advection though right? Because all others use the current time step anyway? Asking because with #555 I've moved the evaluating of all parameterizations to the previous instead of the current timestep which seems to help massively with numerical stability. It also largely makes all flux limiters we previously had redundant. Just checking in relation to the vertical advection whether it would also be important to evaluate more at the previous rather than the current time step?