Skip to content

Commit

Permalink
Merge pull request #344 from SpeedyWeather/mk/howto
Browse files Browse the repository at this point in the history
absolute paths docs img
  • Loading branch information
milankl authored Jun 1, 2023
2 parents b8a5ef7 + 0b5b1b2 commit e3df11d
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions docs/src/how_to_run_speedy.md
Original file line number Diff line number Diff line change
Expand Up @@ -66,7 +66,7 @@ The `run!` command will always return the prognostic variables, which, by defaul
plotted for surface relative vorticity with a unicode plot. The resolution of the plot
is not necessarily representative but it lets us have a quick look at the result

![Barotropic vorticity unicode plot](../img/barotropic_vorticity.jpg)
![Barotropic vorticity unicode plot](https://raw.githubusercontent.com/SpeedyWeather/SpeedyWeather.jl/main/docs/img/barotropic_vorticity.jpg)

Woohoo! I can see turbulence! You could pick up where this simulation stopped by simply
doing `run!(simulation,n_days=50)` again. We didn't store any output, which
Expand Down Expand Up @@ -108,7 +108,7 @@ Now we construct a model, but this time a `ShallowWaterModel`
julia> model = ShallowWaterModel(;spectral_grid, orography, initial_conditions);
julia> simulation = initialize!(model);
```
![Galewsky jet unicode plot](../img/galewsky.jpg)
![Galewsky jet unicode plot](https://raw.githubusercontent.com/SpeedyWeather/SpeedyWeather.jl/main/docs/img/galewsky.jpg)

Oh yeah. That looks like the wobbly jet in their paper. Let's run it again for another 6 days
but this time also store [NetCDF output](@ref).
Expand Down Expand Up @@ -141,7 +141,7 @@ julia> pcolormesh(lon,lat,vor')
```
Which looks like

![Galewsky jet pyplot](../img/galewsky_nc_6days.png)
![Galewsky jet pyplot](https://raw.githubusercontent.com/SpeedyWeather/SpeedyWeather.jl/main/docs/img/galewsky_nc_6days.png)

You see that the unicode plot heavily coarse-grains the simulation, well it's unicode after all!
And now the last time step, that means time=12days is
Expand All @@ -150,7 +150,7 @@ julia> vor = ds["vor"][:,:,1,25];
julia> pcolormesh(lon,lat,vor')
```

![Galewsky jet pyplot](../img/galewsky_nc_12days.png)
![Galewsky jet pyplot](https://raw.githubusercontent.com/SpeedyWeather/SpeedyWeather.jl/main/docs/img/galewsky_nc_12days.png)

The jet broke up into many small eddies, but the turbulence is still confined to the northern hemisphere, cool!
How this may change when we add mountains (we had `NoOrography` above!), say Earth's orography, you may ask?
Expand Down Expand Up @@ -179,7 +179,7 @@ Weather is speedy: run 0003 100%|███████████████
```
This time the run got the id "0003", but otherwise we do as before.

![Galewsky jet pyplot](../img/galewsky_nc_12days_mountains.png)
![Galewsky jet pyplot](https://raw.githubusercontent.com/SpeedyWeather/SpeedyWeather.jl/main/docs/img/galewsky_nc_12days_mountains.png)

Interesting! The initial conditions have zero velocity in the southern hemisphere, but still, one can see
some imprint of the orography on vorticity. You can spot the coastline of Antarctica; the Andes and
Expand Down

0 comments on commit e3df11d

Please sign in to comment.