-
Notifications
You must be signed in to change notification settings - Fork 29
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
Rossby-Haurwitz wave initial conditions +docs #591
Conversation
Sadly, I don't have a reference for the cut-off amplitude and don't exactly remember why I added it. I think it might be that I added it because at first I had the Rossby-Haurwitz wave given by three spherical harmonics coefficients for |
Well it acts now as an example even if it's just you who came up with it ;) I hope that's okay. If you want to make any changes feel free to create a pull request, otherwise you find the Rossby-Haurwitz wave test case now in the documentation: |
Tobias Ehlert @rivwgk had implemented the Rossby-Haurwitz wave as initial conditions for SpeedyWeather, see Williamson et al. 1992, i.e. initial conditions for vorticity are (eq. 145)
With our
set!
function this is now some easy three lines of code that take care of the transform between grid and spectral space for youSo only differences to the mathematical version from above is that we use degrees (
sind
,cosd
) as coordinates, andσ
is our vertical coordinate that needs to be specified for compatibility with our 3D primitive equations models (yes, it's exactly the same syntax inPrimitiveDryModel
andPrimitiveWetModel
!). Looks like@rivwgk can you provide a reference for the cut-off amplitude? It's not discussed in Williamson et al., 1992. You call it$$\sigma$$ but we use that for the vertical coordinate, so it's
c
here. Anyway, you can do that as easy asno need for the loop you wrote, very similar to numpy if you ask me ;)
This pull request
RossbyHaurwitzWave <: AbstractInitialConditions
to have these initial conditions predefined@clessig you wanted to be kept updated ☝🏼