The zonal jet initial conditions from Galewsky, 2004 but for many/multiple zonal jets for SpeedyWeather's shallow-water model.
using SpeedyWeather
using ManyZonalJets
spectral_grid = SpectralGrid(trunc=85, nlev=1)
initial_conditions = ZonalJets(latitude=[45,0,-45])
orography = NoOrography(spectral_grid)
model = ShallowWaterModel(;spectral_grid, orography, initial_conditions)
simulation = initialize!(model)
run!(simulation, period=Day(6))
using ManyZonalJets
exports the ZonalJets
type.
spectral_grid
defines the resolution to be T85 (~165km global).
ZonalJets
initial conditions are then created with jets at 45˚N, 0˚, 45˚S.
NoOrography
is used as in Galewsky, 2004.
Then create a ShallowWaterModel
with those components, initialize!
and
run!
the simulation
for 6 days.
At T127 vorticity evolves like
vor.mp4
with mountains via orography = EarthOrography(spectral_grid)
we have instead
vor.mp4
ManyZonalJets.jl is not registered, so you have to install it manually with
Julia's package manager (open it with ]
)
(@v1.10) pkg> add https://github.com/SpeedyWeather/ManyZonalJets.jl
Or alternatively, using Pkg
, then Pkg.add("https://github.com/SpeedyWeather/ManyZonalJets.jl")
(which is equivalent).