Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Deep atmosphere inconsistencies

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.


New member

I just started with WRF. I implemented the deep atmosphere into the ICON-IAP model, so I checked wether the deep atmospehre is implemented into the ARW solver. I found that replacing the functional determinant of spherical coordinates according to

sqrt(g) = r^2*cos(latitude) -> a^2*cos(latitutde),

where a is the mean Earth radius, is sufficient to consistently derive the shallow atmosphere approximation (compare , they do not show this in detail, but go into the same direction). Otherwise you will run into inconsistencies concerning spurious divergences and/or vorticities (numerically small but aviodable). In my view the shallow atmosphere has four implications:

1.) radius and related grid properties (faces, volumes, depending on the model formulation) independant on height
2.) non-traditional Coriolis components neglected
3.) gravity independant on height
4.) metric terms ~ uw/r, vw/r neglected

I can show You detailed derivations if You like.

From a look into the documentation and the code, ARW is a deep solver in terms of aspects 2.) and 4.), but still makes approximations 1.) and 3.).

Would it be appreciated if I sat down and implemented a height-dependant gravity as well as height-dependant grid properties into the code? Since I only just begin with WRF, this would take some time.

Our top priority at present is MPAS development and application. The capability for deep atmosphere simulation is available in MPAS. We have no plan to add this capability to WRF. It is completely your decision whether to proceed or not on this work.