-
Notifications
You must be signed in to change notification settings - Fork 7
Add Vertical Advection design document #301
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
base: develop
Are you sure you want to change the base?
Add Vertical Advection design document #301
Conversation
|
|
||
| ### 2.3 Requirement: Compute tendencies for mass, momentum, and tracers | ||
|
|
||
| Separate methods are neededto compute the vertical-advection tendencies |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
| Separate methods are neededto compute the vertical-advection tendencies | |
| Separate methods are needed to compute the vertical-advection tendencies |
| for mass, momentum, and tracers. Mass (`LayerThickness`) and horizontal | ||
| momentum (`NormalVelocity`) are stored as 2-D arrays. `LayerThickness` is | ||
| cell-centered, while `NormalVelocity` is defined on edges. These tendencies | ||
| use a first-order flux formulation that ensures strict conservation and energy |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
do we really use first order? I thought thickness advection was centered (2nd order), but I could be misremembering
| This first-order form ensures exact conservation of total column thickness and | ||
| numerical stability. | ||
|
|
||
| ### 3.3 Momentum (Velocity) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
we should call this the velocity equation, it is not momentum as defined
| update the tracer tendencies. | ||
|
|
||
| The `computeFCTVAdvTend` method implements the flux-corrected transport scheme | ||
| following Zalesak 1979. At each interface, both a high-order flux (chosen via |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Are we using zalesak? isn't that 1st order 2nd order? I guess the principle is the same though, but maybe cite this too for the high order part?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
sorry forgot the citation - https://journals.ametsoc.org/view/journals/mwre/139/9/mwr-d-10-05056.1.xml
This PR adds the design document for the vertical advection module. Compiled here.
Checklist