The sponge layer : Différence entre versions
(Page créée avec « == in a nutshell == The sponge layer is necessary to damp vertically propagating waves, which would otherwise non-physically be reflected downward from the model top. In p... ») |
|||
| (Une révision intermédiaire par le même utilisateur non affichée) | |||
| Ligne 5 : | Ligne 5 : | ||
== Sponge layer parameters == | == Sponge layer parameters == | ||
| − | The parameters controlling the sponge layer are in input file gcm.def | + | The parameters controlling the sponge layer are in input file gcm.def (or sometimes in vert_L**.def as they are strongly linked to the vertical grid) |
* iflag_top_bound (integer): 0 for no sponge, 1 for a sponge over the last 4 topmost layers, 2 to have a sponge layer extending over topmost layers down to a pressure which is 100 times that of the topmost layer's. | * iflag_top_bound (integer): 0 for no sponge, 1 for a sponge over the last 4 topmost layers, 2 to have a sponge layer extending over topmost layers down to a pressure which is 100 times that of the topmost layer's. | ||
* mode_top_bound (integer): 0 for no relaxation, 1 to relax winds (u,v) to zero, 2 to relax winds (u,v) towards their zonal mean, 3 to relax winds (u,v) and potential temperature towards their zonal mean. | * mode_top_bound (integer): 0 for no relaxation, 1 to relax winds (u,v) to zero, 2 to relax winds (u,v) towards their zonal mean, 3 to relax winds (u,v) and potential temperature towards their zonal mean. | ||
| − | * tau_top_bound (real, Hz): inverse of the relaxation characteristic time scale at the topmost layer; halved at each successive descending layer. | + | * tau_top_bound (real, Hz): inverse of the relaxation characteristic time scale at the topmost layer; which is then halved at each successive descending layer. |
== Application of the sponge == | == Application of the sponge == | ||
| − | The sponge is applied right after dissipation at each dissipation step | + | The sponge is applied right after dissipation at each dissipation step, but only if the flag '''ok_strato=y''' |
| − | + | == Sponge routine == | |
| + | The sponge computations are done in routine '''dyn3d/top_bound.F''' (for the serial dynamical core; equivalently in '''dyn3dmem/top_bound_loc.F''' for the parallel dynamical core) | ||
| + | |||
| + | 09/05/2022 | ||
[[Category:inputs]] | [[Category:inputs]] | ||
Version actuelle en date du 9 mai 2022 à 11:10
in a nutshell
The sponge layer is necessary to damp vertically propagating waves, which would otherwise non-physically be reflected downward from the model top. In practice its application is thus limited to a few of the topmost layers (usually 4), with decreasing impact with decreasing altitudes. The sponge layer is a relaxation of main variables (winds and/or potential temperature) towards a given value or their zonal mean along a given time scale.
Sponge layer parameters
The parameters controlling the sponge layer are in input file gcm.def (or sometimes in vert_L**.def as they are strongly linked to the vertical grid)
- iflag_top_bound (integer): 0 for no sponge, 1 for a sponge over the last 4 topmost layers, 2 to have a sponge layer extending over topmost layers down to a pressure which is 100 times that of the topmost layer's.
- mode_top_bound (integer): 0 for no relaxation, 1 to relax winds (u,v) to zero, 2 to relax winds (u,v) towards their zonal mean, 3 to relax winds (u,v) and potential temperature towards their zonal mean.
- tau_top_bound (real, Hz): inverse of the relaxation characteristic time scale at the topmost layer; which is then halved at each successive descending layer.
Application of the sponge
The sponge is applied right after dissipation at each dissipation step, but only if the flag ok_strato=y
Sponge routine
The sponge computations are done in routine dyn3d/top_bound.F (for the serial dynamical core; equivalently in dyn3dmem/top_bound_loc.F for the parallel dynamical core)
09/05/2022
