Single-column ocean simulation forced by JRA55 re-analysis

In this example, we simulate the evolution of an ocean water column forced by an atmosphere derived from the JRA55 re-analysis. The simulated column is located at ocean station Papa (144.9ᵒ W and 50.1ᵒ N)

Install dependencies

First let's make sure we have all required packages installed.

using Pkg
pkg"add Oceananigans, ClimaOcean, CairoMakie"
using ClimaOcean
using Oceananigans
using Oceananigans.Units
using Oceananigans.BuoyancyFormulations: buoyancy_frequency
using Oceananigans.Units: Time
using Printf

Construct the grid

First, we construct a single-column grid with 2 meter spacing located at ocean station Papa.

Ocean station papa location

location_name = "ocean_station_papa"
λ★, φ★ = 35.1, 50.1

grid = RectilinearGrid(size = 200,
                       x = λ★,
                       y = φ★,
                       z = (-400, 0),
                       topology = (Flat, Flat, Bounded))
1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── Flat x = 35.1              
├── Flat y = 50.1              
└── Bounded  z ∈ [-400.0, 0.0] regularly spaced with Δz=2.0

An "ocean simulation"

Next, we use ClimaOcean's ocean_simulation constructor to build a realistic ocean simulation on the single-column grid,

ocean = ocean_simulation(grid; Δt=10minutes, coriolis=FPlane(latitude = φ★))
Simulation of HydrostaticFreeSurfaceModel{CPU, RectilinearGrid}(time = 0 seconds, iteration = 0)
├── Next time step: 10 minutes
├── Elapsed wall time: 0 seconds
├── Wall time per iteration: NaN days
├── Stop time: Inf days
├── Stop iteration: Inf
├── Wall time limit: Inf
├── Callbacks: OrderedDict with 4 entries:
│   ├── stop_time_exceeded => Callback of stop_time_exceeded on IterationInterval(1)
│   ├── stop_iteration_exceeded => Callback of stop_iteration_exceeded on IterationInterval(1)
│   ├── wall_time_limit_exceeded => Callback of wall_time_limit_exceeded on IterationInterval(1)
│   └── nan_checker => Callback of NaNChecker for u on IterationInterval(100)
├── Output writers: OrderedDict with no entries
└── Diagnostics: OrderedDict with no entries

which wraps around the ocean model

ocean.model
HydrostaticFreeSurfaceModel{CPU, RectilinearGrid}(time = 0 seconds, iteration = 0)
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── timestepper: QuasiAdamsBashforth2TimeStepper
├── tracers: (T, S, e)
├── closure: CATKEVerticalDiffusivity{VerticallyImplicitTimeDiscretization}
├── buoyancy: SeawaterBuoyancy with g=9.80665 and BoussinesqEquationOfState{Float64} with ĝ = NegativeZDirection()
├── advection scheme: 
│   ├── momentum: Nothing
│   ├── T: Nothing
│   ├── S: Nothing
│   └── e: Nothing
└── coriolis: Oceananigans.Coriolis.FPlane{Float64}

We set initial conditions from ECCO:

set!(ocean.model, T=ECCOMetadata(:temperature), S=ECCOMetadata(:salinity))

A prescribed atmosphere based on JRA55 re-analysis

We build a PrescribedAtmosphere at the same location as the single-colunm grid which is based on the JRA55 reanalysis.

simulation_days = 31
snapshots_per_day = 8 # corresponding to JRA55's 3-hour frequency
last_time = simulation_days * snapshots_per_day
atmosphere = JRA55PrescribedAtmosphere(1:last_time;
                                       longitude = λ★,
                                       latitude = φ★,
                                       backend = InMemory())
2×2×1×248 PrescribedAtmosphere{Float32} on Oceananigans.Grids.LatitudeLongitudeGrid:
├── times: 248-element StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}, Int64}
├── reference_height: 10.0
└── boundary_layer_height: 600.0

This builds a representation of the atmosphere on the small grid

atmosphere.grid
2×2×1 LatitudeLongitudeGrid{Float32, Oceananigans.Grids.Bounded, Oceananigans.Grids.Bounded, Oceananigans.Grids.Flat} on Oceananigans.Architectures.CPU with 2×2×0 halo and with precomputed metrics
├── longitude: Bounded  λ ∈ [34.5938, 35.7188] variably spaced with min(Δλ)=0.5625, max(Δλ)=0.5625
├── latitude:  Bounded  φ ∈ [49.4227, 50.5459] variably spaced with min(Δφ)=0.561619, max(Δφ)=0.561623
└── z:         Flat z                          

Let's take a look at the atmospheric state

ua = interior(atmosphere.velocities.u, 1, 1, 1, :)
va = interior(atmosphere.velocities.v, 1, 1, 1, :)
Ta = interior(atmosphere.tracers.T, 1, 1, 1, :)
qa = interior(atmosphere.tracers.q, 1, 1, 1, :)
t_days = atmosphere.times / days

using CairoMakie

set_theme!(Theme(linewidth=3, fontsize=24))

fig = Figure(size=(800, 600))
axu = Axis(fig[2, 1], xlabel="Days since Jan 1 1990", ylabel="Atmosphere \n velocity (m s⁻¹)")
axT = Axis(fig[3, 1], xlabel="Days since Jan 1 1990", ylabel="Atmosphere \n temperature (K)")
axq = Axis(fig[4, 1], xlabel="Days since Jan 1 1990", ylabel="Atmosphere \n specific humidity")
Label(fig[1, 1], "Atmospheric state over ocean station Papa", tellwidth=false)

lines!(axu, t_days, ua, label="Zonal velocity")
lines!(axu, t_days, va, label="Meridional velocity")
ylims!(axu, -6, 6)
axislegend(axu, framevisible=false, nbanks=2, position=:lb)

lines!(axT, t_days, Ta)
lines!(axq, t_days, qa)

current_figure()

We continue constructing a simulation. For the fluxes computation we use a SkinTemperature formulation that computes the skin temperature from a balance between internal and external heat fluxes.

radiation = Radiation()
similarity_theory = SimilarityTheoryTurbulentFluxes(grid; surface_temperature_type=SkinTemperature())
coupled_model = OceanSeaIceModel(ocean; atmosphere, radiation)
simulation = Simulation(coupled_model, Δt=ocean.Δt, stop_time=30days)

wall_clock = Ref(time_ns())

function progress(sim)
    msg = "Ocean Station Papa"
    msg *= string(", iter: ", iteration(sim), ", time: ", prettytime(sim))

    elapsed = 1e-9 * (time_ns() - wall_clock[])
    msg *= string(", wall time: ", prettytime(elapsed))
    wall_clock[] = time_ns()

    u, v, w = sim.model.ocean.model.velocities
    msg *= @sprintf(", max|u|: (%.2e, %.2e)", maximum(abs, u), maximum(abs, v))

    T = sim.model.ocean.model.tracers.T
    S = sim.model.ocean.model.tracers.S
    e = sim.model.ocean.model.tracers.e

    τx = first(sim.model.fluxes.total.ocean.momentum.u)
    τy = first(sim.model.fluxes.total.ocean.momentum.v)
    Q = first(sim.model.fluxes.total.ocean.heat)

    u★ = sqrt(sqrt(τx^2 + τy^2))

    Nz = size(T, 3)
    msg *= @sprintf(", u★: %.2f m s⁻¹", u★)
    msg *= @sprintf(", Q: %.2f W m⁻²",  Q)
    msg *= @sprintf(", T₀: %.2f ᵒC", first(interior(T, 1, 1, Nz)))
    msg *= @sprintf(", extrema(T): (%.2f, %.2f) ᵒC", minimum(T), maximum(T))
    msg *= @sprintf(", S₀: %.2f g/kg", first(interior(S, 1, 1, Nz)))
    msg *= @sprintf(", e₀: %.2e m² s⁻²", first(interior(e, 1, 1, Nz)))

    @info msg
end

simulation.callbacks[:progress] = Callback(progress, IterationInterval(100))
Callback of progress on IterationInterval(100)

Build flux outputs

τx = coupled_model.fluxes.total.ocean.momentum.u
τy = coupled_model.fluxes.total.ocean.momentum.v
JT = coupled_model.fluxes.total.ocean.tracers.T
Js = coupled_model.fluxes.total.ocean.tracers.S
E  = coupled_model.fluxes.turbulent.fields.water_vapor
Qc = coupled_model.fluxes.turbulent.fields.sensible_heat
Qv = coupled_model.fluxes.turbulent.fields.latent_heat
ρₒ = coupled_model.fluxes.ocean_reference_density
cₚ = coupled_model.fluxes.ocean_heat_capacity

Q = ρₒ * cₚ * JT
ρτx = ρₒ * τx
ρτy = ρₒ * τy
N² = buoyancy_frequency(ocean.model)
κc = ocean.model.diffusivity_fields.κc

fluxes = (; ρτx, ρτy, E, Js, Qv, Qc)
auxiliary_fields = (; N², κc)
fields = merge(ocean.model.velocities, ocean.model.tracers, auxiliary_fields)
(u = 1×1×200 Field{Oceananigans.Grids.Face, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: Flux, top: Flux, immersed: ZeroFlux
└── data: 1×1×206 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:203) with eltype Float64 with indices 1:1×1:1×-2:203
    └── max=0.0, min=0.0, mean=0.0, v = 1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Face, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: Flux, top: Flux, immersed: ZeroFlux
└── data: 1×1×206 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:203) with eltype Float64 with indices 1:1×1:1×-2:203
    └── max=0.0, min=0.0, mean=0.0, w = ZeroField{Int64}, T = 1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: ZeroFlux, top: Flux, immersed: ZeroFlux
└── data: 1×1×206 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:203) with eltype Float64 with indices 1:1×1:1×-2:203
    └── max=17.4568, min=14.8668, mean=16.0725, S = 1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: ZeroFlux, top: Flux, immersed: ZeroFlux
└── data: 1×1×206 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:203) with eltype Float64 with indices 1:1×1:1×-2:203
    └── max=39.0664, min=38.5484, mean=38.9286, e = 1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: ZeroFlux, top: Flux, immersed: ZeroFlux
└── data: 1×1×206 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:203) with eltype Float64 with indices 1:1×1:1×-2:203
    └── max=0.0, min=0.0, mean=0.0, N² = KernelFunctionOperation at (Center, Center, Face)
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── kernel_function: ∂z_b (generic function with 4 methods)
└── arguments: ("SeawaterBuoyancy with g=9.80665 and BoussinesqEquationOfState{Float64}", "(T=1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU, S=1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU, e=1×1×200 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU)"), κc = 1×1×201 Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Face} on Oceananigans.Grids.RectilinearGrid on Oceananigans.Architectures.CPU
├── grid: 1×1×200 RectilinearGrid{Float64, Oceananigans.Grids.Flat, Oceananigans.Grids.Flat, Oceananigans.Grids.Bounded} on Oceananigans.Architectures.CPU with 0×0×3 halo
├── boundary conditions: FieldBoundaryConditions
│   └── west: Nothing, east: Nothing, south: Nothing, north: Nothing, bottom: Nothing, top: Nothing, immersed: ZeroFlux
└── data: 1×1×207 OffsetArray(::Array{Float64, 3}, 1:1, 1:1, -2:204) with eltype Float64 with indices 1:1×1:1×-2:204
    └── max=4.09156e-5, min=0.0, mean=3.04304e-7)

Slice fields at the surface

outputs = merge(fields, fluxes)

filename = "single_column_omip_$(location_name)"

simulation.output_writers[:jld2] = JLD2OutputWriter(ocean.model, outputs; filename,
                                                    schedule = TimeInterval(3hours),
                                                    overwrite_existing = true)

run!(simulation)

filename *= ".jld2"

u  = FieldTimeSeries(filename, "u")
v  = FieldTimeSeries(filename, "v")
T  = FieldTimeSeries(filename, "T")
S  = FieldTimeSeries(filename, "S")
e  = FieldTimeSeries(filename, "e")
N² = FieldTimeSeries(filename, "N²")
κ  = FieldTimeSeries(filename, "κc")

Qv = FieldTimeSeries(filename, "Qv")
Qc = FieldTimeSeries(filename, "Qc")
Js = FieldTimeSeries(filename, "Js")
Ev = FieldTimeSeries(filename, "E")
ρτx = FieldTimeSeries(filename, "ρτx")
ρτy = FieldTimeSeries(filename, "ρτy")

Nz = size(T, 3)
times = Qc.times

ua  = atmosphere.velocities.u
va  = atmosphere.velocities.v
Ta  = atmosphere.tracers.T
qa  = atmosphere.tracers.q
Qlw = atmosphere.downwelling_radiation.longwave
Qsw = atmosphere.downwelling_radiation.shortwave
Pr  = atmosphere.freshwater_flux.rain
Ps  = atmosphere.freshwater_flux.snow

Nt   = length(times)
uat  = zeros(Nt)
vat  = zeros(Nt)
Tat  = zeros(Nt)
qat  = zeros(Nt)
Qswt = zeros(Nt)
Qlwt = zeros(Nt)
Pt   = zeros(Nt)

for n = 1:Nt
    t = times[n]
    uat[n]  =  ua[1, 1, 1, Time(t)]
    vat[n]  =  va[1, 1, 1, Time(t)]
    Tat[n]  =  Ta[1, 1, 1, Time(t)]
    qat[n]  =  qa[1, 1, 1, Time(t)]
    Qswt[n] = Qsw[1, 1, 1, Time(t)]
    Qlwt[n] = Qlw[1, 1, 1, Time(t)]
    Pt[n]   =  Pr[1, 1, 1, Time(t)] + Ps[1, 1, 1, Time(t)]
end

fig = Figure(size=(1800, 1800))

axτ = Axis(fig[1, 1:3], xlabel="Days since Oct 1 1992", ylabel="Wind stress (N m⁻²)")
axQ = Axis(fig[1, 4:6], xlabel="Days since Oct 1 1992", ylabel="Heat flux (W m⁻²)")
axu = Axis(fig[2, 1:3], xlabel="Days since Oct 1 1992", ylabel="Velocities (m s⁻¹)")
axT = Axis(fig[2, 4:6], xlabel="Days since Oct 1 1992", ylabel="Surface temperature (ᵒC)")
axF = Axis(fig[3, 1:3], xlabel="Days since Oct 1 1992", ylabel="Freshwater volume flux (m s⁻¹)")
axS = Axis(fig[3, 4:6], xlabel="Days since Oct 1 1992", ylabel="Surface salinity (g kg⁻¹)")

axuz = Axis(fig[4:5, 1:2], xlabel="Velocities (m s⁻¹)",                ylabel="z (m)")
axTz = Axis(fig[4:5, 3:4], xlabel="Temperature (ᵒC)",                  ylabel="z (m)")
axSz = Axis(fig[4:5, 5:6], xlabel="Salinity (g kg⁻¹)",                 ylabel="z (m)")
axNz = Axis(fig[6:7, 1:2], xlabel="Buoyancy frequency (s⁻²)",          ylabel="z (m)")
axκz = Axis(fig[6:7, 3:4], xlabel="Eddy diffusivity (m² s⁻¹)",         ylabel="z (m)", xscale=log10)
axez = Axis(fig[6:7, 5:6], xlabel="Turbulent kinetic energy (m² s⁻²)", ylabel="z (m)", xscale=log10)

title = @sprintf("Single-column simulation at %.2f, %.2f", φ★, λ★)
Label(fig[0, 1:6], title)

n = Observable(1)

times = (times .- times[1]) ./days
Nt = length(times)
tn = @lift times[$n]

colors = Makie.wong_colors()

ρₒ = coupled_model.fluxes.ocean_reference_density
τx = interior(ρτx, 1, 1, 1, :) ./ ρₒ
τy = interior(ρτy, 1, 1, 1, :) ./ ρₒ
u★ = @. (τx^2 + τy^2)^(1/4)

lines!(axu, times, interior(u, 1, 1, Nz, :), color=colors[1], label="Zonal")
lines!(axu, times, interior(v, 1, 1, Nz, :), color=colors[2], label="Meridional")
lines!(axu, times, u★, color=colors[3], label="Ocean-side u★")
vlines!(axu, tn, linewidth=4, color=(:black, 0.5))
axislegend(axu)

lines!(axτ, times, interior(ρτx, 1, 1, 1, :), label="Zonal")
lines!(axτ, times, interior(ρτy, 1, 1, 1, :), label="Meridional")
vlines!(axτ, tn, linewidth=4, color=(:black, 0.5))
axislegend(axτ)

lines!(axT, times, Tat[1:Nt] .- 273.15,      color=colors[1], linewidth=2, linestyle=:dash, label="Atmosphere temperature")
lines!(axT, times, interior(T, 1, 1, Nz, :), color=colors[2], linewidth=4, label="Ocean surface temperature")
vlines!(axT, tn, linewidth=4, color=(:black, 0.5))
axislegend(axT)

lines!(axQ, times, interior(Qv, 1, 1, 1, 1:Nt),    color=colors[2], label="Sensible",  linewidth=2)
lines!(axQ, times, interior(Qc, 1, 1, 1, 1:Nt),    color=colors[3], label="Latent",    linewidth=2)
lines!(axQ, times, - interior(Qsw, 1, 1, 1, 1:Nt), color=colors[4], label="Shortwave", linewidth=2)
lines!(axQ, times, - interior(Qlw, 1, 1, 1, 1:Nt), color=colors[5], label="Longwave",  linewidth=2)
vlines!(axQ, tn, linewidth=4, color=(:black, 0.5))
axislegend(axQ)

lines!(axF, times, Pt[1:Nt], label="Prescribed freshwater flux")
lines!(axF, times, - interior(Ev, 1, 1, 1, 1:Nt), label="Evaporation")
vlines!(axF, tn, linewidth=4, color=(:black, 0.5))
axislegend(axF)

lines!(axS, times, interior(S, 1, 1, Nz, :))
vlines!(axS, tn, linewidth=4, color=(:black, 0.5))

zc = znodes(T)
zf = znodes(κ)
un  = @lift interior(u[$n],  1, 1, :)
vn  = @lift interior(v[$n],  1, 1, :)
Tn  = @lift interior(T[$n],  1, 1, :)
Sn  = @lift interior(S[$n],  1, 1, :)
κn  = @lift interior(κ[$n],  1, 1, :)
en  = @lift interior(e[$n],  1, 1, :)
N²n = @lift interior(N²[$n], 1, 1, :)

scatterlines!(axuz, un,  zc, label="u")
scatterlines!(axuz, vn,  zc, label="v")
scatterlines!(axTz, Tn,  zc)
scatterlines!(axSz, Sn,  zc)
scatterlines!(axez, en,  zc)
scatterlines!(axNz, N²n, zf)
scatterlines!(axκz, κn,  zf)

axislegend(axuz)

ulim = max(maximum(abs, u), maximum(abs, v))
xlims!(axuz, -ulim, ulim)

Tmax = maximum(interior(T))
Tmin = minimum(interior(T))
xlims!(axTz, Tmin - 0.1, Tmax + 0.1)

Nmax = maximum(interior(N²))
xlims!(axNz, -Nmax/10, Nmax * 1.05)

κmax = maximum(interior(κ))
xlims!(axκz, 1e-9, κmax * 1.1)

emax = maximum(interior(e))
xlims!(axez, 1e-11, emax * 1.1)

Smax = maximum(interior(S))
Smin = minimum(interior(S))
xlims!(axSz, Smin - 0.2, Smax + 0.2)

record(fig, "single_column_profiles.mp4", 1:Nt, framerate=24) do nn
    @info "Drawing frame $nn of $Nt..."
    n[] = nn
end
[ Info: Initializing simulation...
[ Info: Ocean Station Papa, iter: 0, time: 0 seconds, wall time: 2.749 minutes, max|u|: (0.00e+00, 0.00e+00), u★: 0.00 m s⁻¹, Q: 232.78 W m⁻², T₀: 17.38 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.55 g/kg, e₀: 0.00e+00 m² s⁻²
[ Info:     ... simulation initialization complete (5.448 minutes)
[ Info: Executing initial time step...
┌ Warning: Simulation stopped during initialization.
└ @ Oceananigans.Simulations /central/scratch/esm/slurm-buildkite/climaocean-examples/1401/depot/default/packages/Oceananigans/Qooz1/src/Simulations/run.jl:129
[ Info:     ... initial time step complete (15.279 seconds).
[ Info: Ocean Station Papa, iter: 100, time: 16.667 hours, wall time: 6.007 minutes, max|u|: (1.53e-04, 7.64e-04), u★: 0.00 m s⁻¹, Q: 261.26 W m⁻², T₀: 17.34 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 3.19e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 200, time: 1.389 days, wall time: 641.788 ms, max|u|: (3.01e-04, 1.90e-03), u★: 0.00 m s⁻¹, Q: 181.47 W m⁻², T₀: 17.30 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 1.82e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 300, time: 2.083 days, wall time: 565.751 ms, max|u|: (1.70e-04, 9.27e-04), u★: 0.00 m s⁻¹, Q: 191.75 W m⁻², T₀: 17.27 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 2.46e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 400, time: 2.778 days, wall time: 787.054 ms, max|u|: (1.30e-03, 1.29e-03), u★: 0.00 m s⁻¹, Q: 242.76 W m⁻², T₀: 17.24 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 2.27e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 500, time: 3.472 days, wall time: 617.457 ms, max|u|: (1.34e-03, 6.33e-04), u★: 0.00 m s⁻¹, Q: 254.63 W m⁻², T₀: 17.20 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 3.01e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 600, time: 4.167 days, wall time: 937.528 ms, max|u|: (7.67e-04, 1.43e-03), u★: 0.00 m s⁻¹, Q: 212.42 W m⁻², T₀: 17.15 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 2.49e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 700, time: 4.861 days, wall time: 753.056 ms, max|u|: (5.52e-05, 1.88e-03), u★: 0.00 m s⁻¹, Q: 175.63 W m⁻², T₀: 17.13 ᵒC, extrema(T): (14.87, 17.46) ᵒC, S₀: 38.56 g/kg, e₀: 2.16e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 800, time: 5.556 days, wall time: 896.540 ms, max|u|: (8.14e-04, 1.52e-03), u★: 0.00 m s⁻¹, Q: 166.69 W m⁻², T₀: 17.10 ᵒC, extrema(T): (14.87, 17.39) ᵒC, S₀: 38.57 g/kg, e₀: 2.55e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 900, time: 6.250 days, wall time: 735.445 ms, max|u|: (2.19e-03, 1.44e-03), u★: 0.00 m s⁻¹, Q: 186.30 W m⁻², T₀: 17.07 ᵒC, extrema(T): (14.87, 17.36) ᵒC, S₀: 38.57 g/kg, e₀: 2.44e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1000, time: 6.944 days, wall time: 751.437 ms, max|u|: (1.89e-03, 1.80e-04), u★: 0.00 m s⁻¹, Q: 183.07 W m⁻², T₀: 17.04 ᵒC, extrema(T): (14.87, 17.36) ᵒC, S₀: 38.57 g/kg, e₀: 2.63e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1100, time: 7.639 days, wall time: 854.597 ms, max|u|: (2.64e-03, 8.74e-04), u★: 0.00 m s⁻¹, Q: 196.86 W m⁻², T₀: 17.02 ᵒC, extrema(T): (14.87, 17.28) ᵒC, S₀: 38.57 g/kg, e₀: 2.49e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1200, time: 8.333 days, wall time: 780.477 ms, max|u|: (1.88e-03, 5.63e-04), u★: 0.00 m s⁻¹, Q: 85.95 W m⁻², T₀: 16.99 ᵒC, extrema(T): (14.87, 17.27) ᵒC, S₀: 38.57 g/kg, e₀: 1.70e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1300, time: 9.028 days, wall time: 954.475 ms, max|u|: (1.90e-03, 1.24e-03), u★: 0.00 m s⁻¹, Q: 195.95 W m⁻², T₀: 16.96 ᵒC, extrema(T): (14.87, 17.27) ᵒC, S₀: 38.58 g/kg, e₀: 3.03e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1400, time: 9.722 days, wall time: 914.318 ms, max|u|: (1.28e-03, 2.06e-03), u★: 0.00 m s⁻¹, Q: 184.79 W m⁻², T₀: 16.93 ᵒC, extrema(T): (14.87, 17.21) ᵒC, S₀: 38.58 g/kg, e₀: 2.51e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1500, time: 10.417 days, wall time: 1.266 seconds, max|u|: (3.36e-03, 4.39e-03), u★: 0.00 m s⁻¹, Q: 239.17 W m⁻², T₀: 16.89 ᵒC, extrema(T): (14.87, 17.17) ᵒC, S₀: 38.58 g/kg, e₀: 3.04e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1600, time: 11.111 days, wall time: 958.568 ms, max|u|: (1.69e-03, 3.12e-03), u★: 0.00 m s⁻¹, Q: 212.68 W m⁻², T₀: 16.86 ᵒC, extrema(T): (14.87, 17.17) ᵒC, S₀: 38.58 g/kg, e₀: 2.86e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1700, time: 11.806 days, wall time: 1.076 seconds, max|u|: (8.60e-05, 3.05e-03), u★: 0.00 m s⁻¹, Q: 137.42 W m⁻², T₀: 16.84 ᵒC, extrema(T): (14.87, 17.13) ᵒC, S₀: 38.59 g/kg, e₀: 1.90e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1800, time: 12.500 days, wall time: 1.016 seconds, max|u|: (1.55e-03, 3.15e-03), u★: 0.00 m s⁻¹, Q: 202.47 W m⁻², T₀: 16.82 ᵒC, extrema(T): (14.87, 17.09) ᵒC, S₀: 38.59 g/kg, e₀: 2.57e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 1900, time: 13.194 days, wall time: 1.169 seconds, max|u|: (2.32e-03, 3.97e-03), u★: 0.00 m s⁻¹, Q: 247.62 W m⁻², T₀: 16.78 ᵒC, extrema(T): (14.87, 17.08) ᵒC, S₀: 38.59 g/kg, e₀: 3.61e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2000, time: 13.889 days, wall time: 1.011 seconds, max|u|: (1.25e-03, 3.32e-03), u★: 0.00 m s⁻¹, Q: 360.37 W m⁻², T₀: 16.73 ᵒC, extrema(T): (14.87, 17.04) ᵒC, S₀: 38.59 g/kg, e₀: 3.97e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2100, time: 14.583 days, wall time: 1.100 seconds, max|u|: (2.02e-03, 2.70e-03), u★: 0.00 m s⁻¹, Q: 298.53 W m⁻², T₀: 16.69 ᵒC, extrema(T): (14.87, 16.98) ᵒC, S₀: 38.60 g/kg, e₀: 3.42e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2200, time: 15.278 days, wall time: 953.548 ms, max|u|: (1.97e-03, 1.56e-03), u★: 0.00 m s⁻¹, Q: 218.40 W m⁻², T₀: 16.64 ᵒC, extrema(T): (14.87, 16.98) ᵒC, S₀: 38.60 g/kg, e₀: 3.15e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2300, time: 15.972 days, wall time: 1.095 seconds, max|u|: (2.18e-03, 8.61e-04), u★: 0.00 m s⁻¹, Q: 280.61 W m⁻², T₀: 16.60 ᵒC, extrema(T): (14.87, 16.90) ᵒC, S₀: 38.60 g/kg, e₀: 3.56e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2400, time: 16.667 days, wall time: 1.199 seconds, max|u|: (1.89e-03, 3.03e-04), u★: 0.00 m s⁻¹, Q: 223.01 W m⁻², T₀: 16.57 ᵒC, extrema(T): (14.87, 16.90) ᵒC, S₀: 38.61 g/kg, e₀: 2.93e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2500, time: 17.361 days, wall time: 884.027 ms, max|u|: (1.88e-03, 1.26e-03), u★: 0.00 m s⁻¹, Q: 81.78 W m⁻², T₀: 16.54 ᵒC, extrema(T): (14.87, 16.90) ᵒC, S₀: 38.61 g/kg, e₀: 1.83e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2600, time: 18.056 days, wall time: 1.101 seconds, max|u|: (1.14e-03, 1.59e-03), u★: 0.00 m s⁻¹, Q: 262.60 W m⁻², T₀: 16.51 ᵒC, extrema(T): (14.87, 16.88) ᵒC, S₀: 38.61 g/kg, e₀: 3.74e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2700, time: 18.750 days, wall time: 812.141 ms, max|u|: (2.33e-04, 1.24e-03), u★: 0.00 m s⁻¹, Q: 286.11 W m⁻², T₀: 16.49 ᵒC, extrema(T): (14.87, 16.82) ᵒC, S₀: 38.61 g/kg, e₀: 3.79e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2800, time: 19.444 days, wall time: 918.794 ms, max|u|: (4.81e-04, 1.23e-03), u★: 0.00 m s⁻¹, Q: 163.27 W m⁻², T₀: 16.45 ᵒC, extrema(T): (14.87, 16.82) ᵒC, S₀: 38.61 g/kg, e₀: 2.16e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 2900, time: 20.139 days, wall time: 922.204 ms, max|u|: (9.19e-04, 1.13e-03), u★: 0.00 m s⁻¹, Q: 224.43 W m⁻², T₀: 16.42 ᵒC, extrema(T): (14.87, 16.82) ᵒC, S₀: 38.62 g/kg, e₀: 3.48e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3000, time: 20.833 days, wall time: 949.309 ms, max|u|: (1.63e-03, 9.63e-05), u★: 0.00 m s⁻¹, Q: 264.77 W m⁻², T₀: 16.39 ᵒC, extrema(T): (14.87, 16.77) ᵒC, S₀: 38.62 g/kg, e₀: 3.82e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3100, time: 21.528 days, wall time: 700.949 ms, max|u|: (1.49e-03, 5.81e-04), u★: 0.00 m s⁻¹, Q: 220.15 W m⁻², T₀: 16.36 ᵒC, extrema(T): (14.87, 16.74) ᵒC, S₀: 38.62 g/kg, e₀: 2.92e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3200, time: 22.222 days, wall time: 912.599 ms, max|u|: (9.29e-04, 1.01e-03), u★: 0.00 m s⁻¹, Q: 143.97 W m⁻², T₀: 16.33 ᵒC, extrema(T): (14.87, 16.74) ᵒC, S₀: 38.62 g/kg, e₀: 2.68e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3300, time: 22.917 days, wall time: 895.096 ms, max|u|: (5.07e-04, 1.81e-03), u★: 0.00 m s⁻¹, Q: 226.24 W m⁻², T₀: 16.31 ᵒC, extrema(T): (14.87, 16.74) ᵒC, S₀: 38.62 g/kg, e₀: 3.35e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3400, time: 23.611 days, wall time: 858.535 ms, max|u|: (5.92e-04, 2.12e-03), u★: 0.00 m s⁻¹, Q: 202.04 W m⁻², T₀: 16.28 ᵒC, extrema(T): (14.87, 16.68) ᵒC, S₀: 38.63 g/kg, e₀: 2.55e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3500, time: 24.306 days, wall time: 773.777 ms, max|u|: (1.05e-03, 1.18e-03), u★: 0.00 m s⁻¹, Q: 127.61 W m⁻², T₀: 16.25 ᵒC, extrema(T): (14.87, 16.66) ᵒC, S₀: 38.63 g/kg, e₀: 2.47e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3600, time: 25 days, wall time: 805.292 ms, max|u|: (3.03e-03, 2.64e-03), u★: 0.01 m s⁻¹, Q: 410.07 W m⁻², T₀: 16.21 ᵒC, extrema(T): (14.87, 16.66) ᵒC, S₀: 38.63 g/kg, e₀: 4.56e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3700, time: 25.694 days, wall time: 941.486 ms, max|u|: (9.61e-04, 2.46e-03), u★: 0.00 m s⁻¹, Q: 343.23 W m⁻², T₀: 16.17 ᵒC, extrema(T): (14.87, 16.58) ᵒC, S₀: 38.64 g/kg, e₀: 4.45e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3800, time: 26.389 days, wall time: 742.183 ms, max|u|: (8.11e-04, 1.80e-03), u★: 0.00 m s⁻¹, Q: 163.10 W m⁻², T₀: 16.13 ᵒC, extrema(T): (14.87, 16.58) ᵒC, S₀: 38.64 g/kg, e₀: 2.43e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 3900, time: 27.083 days, wall time: 839.337 ms, max|u|: (1.40e-04, 1.55e-03), u★: 0.00 m s⁻¹, Q: 406.33 W m⁻², T₀: 16.08 ᵒC, extrema(T): (14.87, 16.53) ᵒC, S₀: 38.64 g/kg, e₀: 4.84e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 4000, time: 27.778 days, wall time: 786.815 ms, max|u|: (3.16e-04, 1.98e-03), u★: 0.00 m s⁻¹, Q: 395.20 W m⁻², T₀: 16.04 ᵒC, extrema(T): (14.87, 16.53) ᵒC, S₀: 38.65 g/kg, e₀: 4.72e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 4100, time: 28.472 days, wall time: 905.502 ms, max|u|: (4.29e-04, 1.55e-03), u★: 0.00 m s⁻¹, Q: 227.02 W m⁻², T₀: 16.00 ᵒC, extrema(T): (14.87, 16.49) ᵒC, S₀: 38.65 g/kg, e₀: 2.93e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 4200, time: 29.167 days, wall time: 988.850 ms, max|u|: (1.21e-03, 6.53e-04), u★: 0.00 m s⁻¹, Q: 219.16 W m⁻², T₀: 15.97 ᵒC, extrema(T): (14.87, 16.49) ᵒC, S₀: 38.65 g/kg, e₀: 3.18e-04 m² s⁻²
[ Info: Ocean Station Papa, iter: 4300, time: 29.861 days, wall time: 803.514 ms, max|u|: (1.98e-03, 2.05e-04), u★: 0.00 m s⁻¹, Q: 355.56 W m⁻², T₀: 15.94 ᵒC, extrema(T): (14.87, 16.49) ᵒC, S₀: 38.65 g/kg, e₀: 4.08e-04 m² s⁻²
[ Info: Simulation is stopping after running for 6.660 minutes.
[ Info: Simulation time 30 days equals or exceeds stop time 30 days.
[ Info: Drawing frame 1 of 241...
[ Info: Drawing frame 2 of 241...
[ Info: Drawing frame 3 of 241...
[ Info: Drawing frame 4 of 241...
[ Info: Drawing frame 5 of 241...
[ Info: Drawing frame 6 of 241...
[ Info: Drawing frame 7 of 241...
[ Info: Drawing frame 8 of 241...
[ Info: Drawing frame 9 of 241...
[ Info: Drawing frame 10 of 241...
[ Info: Drawing frame 11 of 241...
[ Info: Drawing frame 12 of 241...
[ Info: Drawing frame 13 of 241...
[ Info: Drawing frame 14 of 241...
[ Info: Drawing frame 15 of 241...
[ Info: Drawing frame 16 of 241...
[ Info: Drawing frame 17 of 241...
[ Info: Drawing frame 18 of 241...
[ Info: Drawing frame 19 of 241...
[ Info: Drawing frame 20 of 241...
[ Info: Drawing frame 21 of 241...
[ Info: Drawing frame 22 of 241...
[ Info: Drawing frame 23 of 241...
[ Info: Drawing frame 24 of 241...
[ Info: Drawing frame 25 of 241...
[ Info: Drawing frame 26 of 241...
[ Info: Drawing frame 27 of 241...
[ Info: Drawing frame 28 of 241...
[ Info: Drawing frame 29 of 241...
[ Info: Drawing frame 30 of 241...
[ Info: Drawing frame 31 of 241...
[ Info: Drawing frame 32 of 241...
[ Info: Drawing frame 33 of 241...
[ Info: Drawing frame 34 of 241...
[ Info: Drawing frame 35 of 241...
[ Info: Drawing frame 36 of 241...
[ Info: Drawing frame 37 of 241...
[ Info: Drawing frame 38 of 241...
[ Info: Drawing frame 39 of 241...
[ Info: Drawing frame 40 of 241...
[ Info: Drawing frame 41 of 241...
[ Info: Drawing frame 42 of 241...
[ Info: Drawing frame 43 of 241...
[ Info: Drawing frame 44 of 241...
[ Info: Drawing frame 45 of 241...
[ Info: Drawing frame 46 of 241...
[ Info: Drawing frame 47 of 241...
[ Info: Drawing frame 48 of 241...
[ Info: Drawing frame 49 of 241...
[ Info: Drawing frame 50 of 241...
[ Info: Drawing frame 51 of 241...
[ Info: Drawing frame 52 of 241...
[ Info: Drawing frame 53 of 241...
[ Info: Drawing frame 54 of 241...
[ Info: Drawing frame 55 of 241...
[ Info: Drawing frame 56 of 241...
[ Info: Drawing frame 57 of 241...
[ Info: Drawing frame 58 of 241...
[ Info: Drawing frame 59 of 241...
[ Info: Drawing frame 60 of 241...
[ Info: Drawing frame 61 of 241...
[ Info: Drawing frame 62 of 241...
[ Info: Drawing frame 63 of 241...
[ Info: Drawing frame 64 of 241...
[ Info: Drawing frame 65 of 241...
[ Info: Drawing frame 66 of 241...
[ Info: Drawing frame 67 of 241...
[ Info: Drawing frame 68 of 241...
[ Info: Drawing frame 69 of 241...
[ Info: Drawing frame 70 of 241...
[ Info: Drawing frame 71 of 241...
[ Info: Drawing frame 72 of 241...
[ Info: Drawing frame 73 of 241...
[ Info: Drawing frame 74 of 241...
[ Info: Drawing frame 75 of 241...
[ Info: Drawing frame 76 of 241...
[ Info: Drawing frame 77 of 241...
[ Info: Drawing frame 78 of 241...
[ Info: Drawing frame 79 of 241...
[ Info: Drawing frame 80 of 241...
[ Info: Drawing frame 81 of 241...
[ Info: Drawing frame 82 of 241...
[ Info: Drawing frame 83 of 241...
[ Info: Drawing frame 84 of 241...
[ Info: Drawing frame 85 of 241...
[ Info: Drawing frame 86 of 241...
[ Info: Drawing frame 87 of 241...
[ Info: Drawing frame 88 of 241...
[ Info: Drawing frame 89 of 241...
[ Info: Drawing frame 90 of 241...
[ Info: Drawing frame 91 of 241...
[ Info: Drawing frame 92 of 241...
[ Info: Drawing frame 93 of 241...
[ Info: Drawing frame 94 of 241...
[ Info: Drawing frame 95 of 241...
[ Info: Drawing frame 96 of 241...
[ Info: Drawing frame 97 of 241...
[ Info: Drawing frame 98 of 241...
[ Info: Drawing frame 99 of 241...
[ Info: Drawing frame 100 of 241...
[ Info: Drawing frame 101 of 241...
[ Info: Drawing frame 102 of 241...
[ Info: Drawing frame 103 of 241...
[ Info: Drawing frame 104 of 241...
[ Info: Drawing frame 105 of 241...
[ Info: Drawing frame 106 of 241...
[ Info: Drawing frame 107 of 241...
[ Info: Drawing frame 108 of 241...
[ Info: Drawing frame 109 of 241...
[ Info: Drawing frame 110 of 241...
[ Info: Drawing frame 111 of 241...
[ Info: Drawing frame 112 of 241...
[ Info: Drawing frame 113 of 241...
[ Info: Drawing frame 114 of 241...
[ Info: Drawing frame 115 of 241...
[ Info: Drawing frame 116 of 241...
[ Info: Drawing frame 117 of 241...
[ Info: Drawing frame 118 of 241...
[ Info: Drawing frame 119 of 241...
[ Info: Drawing frame 120 of 241...
[ Info: Drawing frame 121 of 241...
[ Info: Drawing frame 122 of 241...
[ Info: Drawing frame 123 of 241...
[ Info: Drawing frame 124 of 241...
[ Info: Drawing frame 125 of 241...
[ Info: Drawing frame 126 of 241...
[ Info: Drawing frame 127 of 241...
[ Info: Drawing frame 128 of 241...
[ Info: Drawing frame 129 of 241...
[ Info: Drawing frame 130 of 241...
[ Info: Drawing frame 131 of 241...
[ Info: Drawing frame 132 of 241...
[ Info: Drawing frame 133 of 241...
[ Info: Drawing frame 134 of 241...
[ Info: Drawing frame 135 of 241...
[ Info: Drawing frame 136 of 241...
[ Info: Drawing frame 137 of 241...
[ Info: Drawing frame 138 of 241...
[ Info: Drawing frame 139 of 241...
[ Info: Drawing frame 140 of 241...
[ Info: Drawing frame 141 of 241...
[ Info: Drawing frame 142 of 241...
[ Info: Drawing frame 143 of 241...
[ Info: Drawing frame 144 of 241...
[ Info: Drawing frame 145 of 241...
[ Info: Drawing frame 146 of 241...
[ Info: Drawing frame 147 of 241...
[ Info: Drawing frame 148 of 241...
[ Info: Drawing frame 149 of 241...
[ Info: Drawing frame 150 of 241...
[ Info: Drawing frame 151 of 241...
[ Info: Drawing frame 152 of 241...
[ Info: Drawing frame 153 of 241...
[ Info: Drawing frame 154 of 241...
[ Info: Drawing frame 155 of 241...
[ Info: Drawing frame 156 of 241...
[ Info: Drawing frame 157 of 241...
[ Info: Drawing frame 158 of 241...
[ Info: Drawing frame 159 of 241...
[ Info: Drawing frame 160 of 241...
[ Info: Drawing frame 161 of 241...
[ Info: Drawing frame 162 of 241...
[ Info: Drawing frame 163 of 241...
[ Info: Drawing frame 164 of 241...
[ Info: Drawing frame 165 of 241...
[ Info: Drawing frame 166 of 241...
[ Info: Drawing frame 167 of 241...
[ Info: Drawing frame 168 of 241...
[ Info: Drawing frame 169 of 241...
[ Info: Drawing frame 170 of 241...
[ Info: Drawing frame 171 of 241...
[ Info: Drawing frame 172 of 241...
[ Info: Drawing frame 173 of 241...
[ Info: Drawing frame 174 of 241...
[ Info: Drawing frame 175 of 241...
[ Info: Drawing frame 176 of 241...
[ Info: Drawing frame 177 of 241...
[ Info: Drawing frame 178 of 241...
[ Info: Drawing frame 179 of 241...
[ Info: Drawing frame 180 of 241...
[ Info: Drawing frame 181 of 241...
[ Info: Drawing frame 182 of 241...
[ Info: Drawing frame 183 of 241...
[ Info: Drawing frame 184 of 241...
[ Info: Drawing frame 185 of 241...
[ Info: Drawing frame 186 of 241...
[ Info: Drawing frame 187 of 241...
[ Info: Drawing frame 188 of 241...
[ Info: Drawing frame 189 of 241...
[ Info: Drawing frame 190 of 241...
[ Info: Drawing frame 191 of 241...
[ Info: Drawing frame 192 of 241...
[ Info: Drawing frame 193 of 241...
[ Info: Drawing frame 194 of 241...
[ Info: Drawing frame 195 of 241...
[ Info: Drawing frame 196 of 241...
[ Info: Drawing frame 197 of 241...
[ Info: Drawing frame 198 of 241...
[ Info: Drawing frame 199 of 241...
[ Info: Drawing frame 200 of 241...
[ Info: Drawing frame 201 of 241...
[ Info: Drawing frame 202 of 241...
[ Info: Drawing frame 203 of 241...
[ Info: Drawing frame 204 of 241...
[ Info: Drawing frame 205 of 241...
[ Info: Drawing frame 206 of 241...
[ Info: Drawing frame 207 of 241...
[ Info: Drawing frame 208 of 241...
[ Info: Drawing frame 209 of 241...
[ Info: Drawing frame 210 of 241...
[ Info: Drawing frame 211 of 241...
[ Info: Drawing frame 212 of 241...
[ Info: Drawing frame 213 of 241...
[ Info: Drawing frame 214 of 241...
[ Info: Drawing frame 215 of 241...
[ Info: Drawing frame 216 of 241...
[ Info: Drawing frame 217 of 241...
[ Info: Drawing frame 218 of 241...
[ Info: Drawing frame 219 of 241...
[ Info: Drawing frame 220 of 241...
[ Info: Drawing frame 221 of 241...
[ Info: Drawing frame 222 of 241...
[ Info: Drawing frame 223 of 241...
[ Info: Drawing frame 224 of 241...
[ Info: Drawing frame 225 of 241...
[ Info: Drawing frame 226 of 241...
[ Info: Drawing frame 227 of 241...
[ Info: Drawing frame 228 of 241...
[ Info: Drawing frame 229 of 241...
[ Info: Drawing frame 230 of 241...
[ Info: Drawing frame 231 of 241...
[ Info: Drawing frame 232 of 241...
[ Info: Drawing frame 233 of 241...
[ Info: Drawing frame 234 of 241...
[ Info: Drawing frame 235 of 241...
[ Info: Drawing frame 236 of 241...
[ Info: Drawing frame 237 of 241...
[ Info: Drawing frame 238 of 241...
[ Info: Drawing frame 239 of 241...
[ Info: Drawing frame 240 of 241...
[ Info: Drawing frame 241 of 241...


This page was generated using Literate.jl.