Tapio Schneider

Thanks, @mwarusz. Filtering perturbations from reference state helps, but is not the general solution to the problem. Deviations from reference state can be substantial, and then the problem remains that...

The key thin-shell approximation is replacing the Coriolis force `source.ρu -= SVector(0, 0, 2 * _Omega) × state.ρu` by `source.ρu -= f*vertical_unit_vector × state.ρu`, with Coriolis parameter `f=2\Omega sin(latitude)` (`\Omega`...

The point of ClimaParameters is not just to ensure uniform values of planetary constants across all model components, although that too is important (e.g., to conserve energy in an Earth...

Looks good. A standard test for limiters are step waves in tracers advected with constant velocity. If we do not already have that, we shoud add it as a test....

Thanks, @charleskawczynski. The Guba et al. paper you have describes the horizontal limiter have. In the vertical, there are more choices. * Van Leer limiters (currently not implemented AFAIK but...

You do not need to do anything to the horizontal limiters (Guba); please just run the standard tests on them (not sure what's implemented, perhaps shallow water flow with a...

There's a bit more, @charleskawczynski: * A central result here should be a tested default configuration of timestepper (ideally SSP, but it can be ARIS if SSP is too slow)...

I'd like to see a clear description of the minimum functionality that we need (e.g., transformation covariant -> contravariant and inverse, in the horizontal and in 3d). That is, put...

Standard solution: Specify start time tstart (in seconds since epoch, e.g., obtained from DateTime object), end time tend, and timestep dt. Calculate number of integer time steps needed to get...

Really nice and helpful. Thank you!