using ClimaSeaIce
using ClimaSeaIce.SeaIceMomentumEquations
using ClimaSeaIce.Rheologies
using Oceananigans
using Oceananigans.Units
using Oceananigans.Operators
using Printf
using CairoMakieA solid block of ice moving against a triangular coastline in a periodic channel
Lx = 512kilometers
Ly = 256kilometers
Nx = 256
Ny = 256
y_max = Ly / 2
arch = CPU()
𝓋ₐ = 10.0 # m / s
Cᴰ = 1.2e-3 # Atmosphere - sea ice drag coefficient
ρₐ = 1.3 # kg/m³1.32 km domain
grid = RectilinearGrid(arch; size = (Nx, Ny),
x = (-Lx/2, Lx/2),
y = (0, Ly),
halo = (4, 4),
topology = (Periodic, Bounded, Flat))
bottom(x, y) = ifelse(y > y_max, 0,
ifelse(abs(x / Lx) * Nx + y / Ly * Ny > 24, 0, 1))
grid = ImmersedBoundaryGrid(grid, GridFittedBoundary(bottom))
#####
##### Setup atmospheric and oceanic forcing
#####256×256×1 ImmersedBoundaryGrid{Float64, Oceananigans.Grids.Periodic, Oceananigans.Grids.Bounded, Oceananigans.Grids.Flat} on Oceananigans.Architectures.CPU with 4×4×0 halo:
├── immersed_boundary: Oceananigans.ImmersedBoundaries.GridFittedBoundary{Oceananigans.Fields.Field{Oceananigans.Grids.Center, Oceananigans.Grids.Center, Oceananigans.Grids.Center, Nothing, Oceananigans.Grids.RectilinearGrid{Float64, Oceananigans.Grids.Periodic, Oceananigans.Grids.Bounded, Oceananigans.Grids.Flat, Oceananigans.Grids.StaticVerticalDiscretization{Nothing, Nothing, Float64, Float64}, Float64, Float64, OffsetArrays.OffsetVector{Float64, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}, Int64}}, OffsetArrays.OffsetVector{Float64, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}, Int64}}, Oceananigans.Architectures.CPU}, Tuple{Colon, Colon, Colon}, OffsetArrays.OffsetArray{Bool, 3, Array{Bool, 3}}, Bool, Oceananigans.BoundaryConditions.FieldBoundaryConditions{Oceananigans.BoundaryConditions.BoundaryCondition{Oceananigans.BoundaryConditions.Periodic, Nothing}, Oceananigans.BoundaryConditions.BoundaryCondition{Oceananigans.BoundaryConditions.Periodic, Nothing}, Oceananigans.BoundaryConditions.BoundaryCondition{Oceananigans.BoundaryConditions.Flux, Nothing}, Oceananigans.BoundaryConditions.BoundaryCondition{Oceananigans.BoundaryConditions.Flux, Nothing}, Nothing, Nothing, Oceananigans.BoundaryConditions.BoundaryCondition{Oceananigans.BoundaryConditions.Flux, Nothing}}, Nothing, Oceananigans.Fields.FieldBoundaryBuffers{Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing}}}
├── underlying_grid: 256×256×1 RectilinearGrid{Float64, Oceananigans.Grids.Periodic, Oceananigans.Grids.Bounded, Oceananigans.Grids.Flat} on Oceananigans.Architectures.CPU with 4×4×0 halo
├── Periodic x ∈ [-256000.0, 256000.0) regularly spaced with Δx=2000.0
├── Bounded y ∈ [0.0, 256000.0] regularly spaced with Δy=1000.0
└── Flat z Atmosphere - sea ice stress
Ua = XFaceField(grid)
τᵤ = Field(- ρₐ * Cᴰ * Ua^2)
set!(Ua, 𝓋ₐ)
compute!(τᵤ)
Oceananigans.BoundaryConditions.fill_halo_regions!(τᵤ)
τᵥ = 0.0
#####
##### Ocean stress (a zero-velocity ocean with a drag)
#####
τₒ = SemiImplicitStress()
#####
##### Numerical details
#####ClimaSeaIce.SeaIceMomentumEquations.SemiImplicitStress{Oceananigans.Fields.ZeroField{Float64, 3}, Oceananigans.Fields.ZeroField{Float64, 3}, Float64}(ZeroField{Float64}, ZeroField{Float64}, 1026.0, 0.0055)We use an elasto-visco-plastic rheology and WENO seventh order for advection of h and ℵ
dynamics = SeaIceMomentumEquation(grid;
top_momentum_stress = (u=τᵤ, v=τᵥ),
bottom_momentum_stress = τₒ,
rheology = ElastoViscoPlasticRheology(),
solver = SplitExplicitSolver(substeps=150))
u_bcs = FieldBoundaryConditions(top = nothing, bottom = nothing,
north = ValueBoundaryCondition(0),
south = ValueBoundaryCondition(0))
#Define the model!
model = SeaIceModel(grid;
advection = WENO(order=7),
dynamics = dynamics,
boundary_conditions = (; u=u_bcs),
ice_thermodynamics = nothing)SeaIceModel{Oceananigans.Architectures.CPU, Oceananigans.ImmersedBoundaries.ImmersedBoundaryGrid}(time = 0 seconds, iteration = 0)
├── grid: 256×256×1 ImmersedBoundaryGrid{Float64, Oceananigans.Grids.Periodic, Oceananigans.Grids.Bounded, Oceananigans.Grids.Flat} on Oceananigans.Architectures.CPU with 4×4×0 halo
├── ice_thermodynamics: Nothing
├── advection: WENO(order=7)
└── external_heat_fluxes:
├── top: Nothing
└── bottom: 0We start with a concentration of ℵ = 1 everywhere
set!(model, h = 1)
set!(model, ℵ = 1)
#####
##### Setup the simulation
#####run the model for 10 day
simulation = Simulation(model, Δt = 2minutes, stop_time=2days)Simulation of SeaIceModel
├── Next time step: 2 minutes
├── Elapsed wall time: 0 seconds
├── Wall time per iteration: NaN days
├── Stop time: 2 days
├── Stop iteration: Inf
├── Wall time limit: Inf
├── Minimum relative step: 0.0
├── Callbacks: OrderedDict with 3 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)
├── Output writers: OrderedDict with no entries
└── Diagnostics: OrderedDict with no entriesContainer to hold the data
htimeseries = []
ℵtimeseries = []
utimeseries = []
vtimeseries = []
# Callback function to collect the data from the `sim`ulation
function accumulate_timeseries(sim)
h = sim.model.ice_thickness
ℵ = sim.model.ice_concentration
u = sim.model.velocities.u
v = sim.model.velocities.v
push!(htimeseries, deepcopy(Array(interior(h))))
push!(ℵtimeseries, deepcopy(Array(interior(ℵ))))
push!(utimeseries, deepcopy(Array(interior(u))))
push!(vtimeseries, deepcopy(Array(interior(v))))
end
wall_time = Ref(time_ns())
function progress(sim)
h = sim.model.ice_thickness
ℵ = sim.model.ice_concentration
u = sim.model.velocities.u
v = sim.model.velocities.v
hmax = maximum(interior(h))
ℵmin = minimum(interior(ℵ))
umax = maximum(interior(u)), maximum(interior(v))
step_time = 1e-9 * (time_ns() - wall_time[])
@info @sprintf("Time: %s, Iteration %d, Δt %s, max(vel): (%.2e, %.2e), max(trac): %.2f, %.2f, wtime: %s \n",
prettytime(sim.model.clock.time),
sim.model.clock.iteration,
prettytime(sim.Δt),
umax..., hmax, ℵmin, prettytime(step_time))
wall_time[] = time_ns()
end
simulation.callbacks[:progress] = Callback(progress, IterationInterval(10))
simulation.callbacks[:save] = Callback(accumulate_timeseries, IterationInterval(10))
run!(simulation)[ Info: Initializing simulation...
[ Info: Time: 0 seconds, Iteration 0, Δt 2 minutes, max(vel): (0.00e+00, 0.00e+00), max(trac): 1.00, 0.00, wtime: 1.032 minutes
[ Info: ... simulation initialization complete (29.882 seconds)
[ Info: Executing initial time step...
[ Info: ... initial time step complete (29.267 seconds).
[ Info: Time: 20 minutes, Iteration 10, Δt 2 minutes, max(vel): (1.25e-01, 1.36e-02), max(trac): 1.02, 0.00, wtime: 39.162 seconds
[ Info: Time: 40 minutes, Iteration 20, Δt 2 minutes, max(vel): (1.70e-01, 1.61e-02), max(trac): 1.09, 0.00, wtime: 11.046 seconds
[ Info: Time: 1 hour, Iteration 30, Δt 2 minutes, max(vel): (1.81e-01, 2.15e-02), max(trac): 1.17, 0.00, wtime: 10.911 seconds
[ Info: Time: 1.333 hours, Iteration 40, Δt 2 minutes, max(vel): (1.82e-01, 1.92e-02), max(trac): 1.26, 0.00, wtime: 11.122 seconds
[ Info: Time: 1.667 hours, Iteration 50, Δt 2 minutes, max(vel): (1.89e-01, 2.10e-02), max(trac): 1.34, 0.00, wtime: 11.025 seconds
[ Info: Time: 2 hours, Iteration 60, Δt 2 minutes, max(vel): (1.85e-01, 1.89e-02), max(trac): 1.41, 0.00, wtime: 10.861 seconds
[ Info: Time: 2.333 hours, Iteration 70, Δt 2 minutes, max(vel): (1.92e-01, 2.30e-02), max(trac): 1.47, 0.00, wtime: 10.951 seconds
[ Info: Time: 2.667 hours, Iteration 80, Δt 2 minutes, max(vel): (1.85e-01, 1.76e-02), max(trac): 1.53, 0.00, wtime: 11.301 seconds
[ Info: Time: 3 hours, Iteration 90, Δt 2 minutes, max(vel): (1.85e-01, 1.99e-02), max(trac): 1.58, 0.00, wtime: 11.090 seconds
[ Info: Time: 3.333 hours, Iteration 100, Δt 2 minutes, max(vel): (1.83e-01, 1.88e-02), max(trac): 1.63, 0.00, wtime: 10.856 seconds
[ Info: Time: 3.667 hours, Iteration 110, Δt 2 minutes, max(vel): (1.87e-01, 1.71e-02), max(trac): 1.67, 0.00, wtime: 10.807 seconds
[ Info: Time: 4 hours, Iteration 120, Δt 2 minutes, max(vel): (1.83e-01, 1.99e-02), max(trac): 1.71, 0.00, wtime: 10.804 seconds
[ Info: Time: 4.333 hours, Iteration 130, Δt 2 minutes, max(vel): (1.82e-01, 1.85e-02), max(trac): 1.74, 0.00, wtime: 10.964 seconds
[ Info: Time: 4.667 hours, Iteration 140, Δt 2 minutes, max(vel): (1.89e-01, 2.22e-02), max(trac): 1.78, 0.00, wtime: 11.345 seconds
[ Info: Time: 5 hours, Iteration 150, Δt 2 minutes, max(vel): (1.85e-01, 1.96e-02), max(trac): 1.81, 0.00, wtime: 11.124 seconds
[ Info: Time: 5.333 hours, Iteration 160, Δt 2 minutes, max(vel): (1.84e-01, 1.98e-02), max(trac): 1.84, 0.00, wtime: 11.017 seconds
[ Info: Time: 5.667 hours, Iteration 170, Δt 2 minutes, max(vel): (1.91e-01, 1.69e-02), max(trac): 1.87, 0.00, wtime: 11.007 seconds
[ Info: Time: 6 hours, Iteration 180, Δt 2 minutes, max(vel): (1.88e-01, 2.41e-02), max(trac): 1.89, 0.00, wtime: 11.124 seconds
[ Info: Time: 6.333 hours, Iteration 190, Δt 2 minutes, max(vel): (1.84e-01, 1.75e-02), max(trac): 1.92, 0.00, wtime: 10.889 seconds
[ Info: Time: 6.667 hours, Iteration 200, Δt 2 minutes, max(vel): (1.82e-01, 1.97e-02), max(trac): 1.94, 0.00, wtime: 10.847 seconds
[ Info: Time: 7 hours, Iteration 210, Δt 2 minutes, max(vel): (1.83e-01, 1.78e-02), max(trac): 1.97, 0.00, wtime: 10.778 seconds
[ Info: Time: 7.333 hours, Iteration 220, Δt 2 minutes, max(vel): (1.81e-01, 2.07e-02), max(trac): 1.98, 0.00, wtime: 10.839 seconds
[ Info: Time: 7.667 hours, Iteration 230, Δt 2 minutes, max(vel): (1.89e-01, 1.88e-02), max(trac): 2.00, 0.00, wtime: 10.822 seconds
[ Info: Time: 8 hours, Iteration 240, Δt 2 minutes, max(vel): (1.80e-01, 2.08e-02), max(trac): 2.02, 0.00, wtime: 10.832 seconds
[ Info: Time: 8.333 hours, Iteration 250, Δt 2 minutes, max(vel): (1.87e-01, 2.23e-02), max(trac): 2.04, 0.00, wtime: 10.881 seconds
[ Info: Time: 8.667 hours, Iteration 260, Δt 2 minutes, max(vel): (1.82e-01, 2.02e-02), max(trac): 2.05, 0.00, wtime: 10.808 seconds
[ Info: Time: 9 hours, Iteration 270, Δt 2 minutes, max(vel): (1.82e-01, 1.98e-02), max(trac): 2.07, 0.00, wtime: 10.784 seconds
[ Info: Time: 9.333 hours, Iteration 280, Δt 2 minutes, max(vel): (1.80e-01, 1.98e-02), max(trac): 2.08, 0.00, wtime: 10.903 seconds
[ Info: Time: 9.667 hours, Iteration 290, Δt 2 minutes, max(vel): (1.82e-01, 2.01e-02), max(trac): 2.09, 0.00, wtime: 10.925 seconds
[ Info: Time: 10 hours, Iteration 300, Δt 2 minutes, max(vel): (1.84e-01, 2.11e-02), max(trac): 2.10, 0.00, wtime: 10.861 seconds
[ Info: Time: 10.333 hours, Iteration 310, Δt 2 minutes, max(vel): (1.85e-01, 1.95e-02), max(trac): 2.11, 0.00, wtime: 10.932 seconds
[ Info: Time: 10.667 hours, Iteration 320, Δt 2 minutes, max(vel): (1.86e-01, 1.86e-02), max(trac): 2.12, 0.00, wtime: 11.057 seconds
[ Info: Time: 11 hours, Iteration 330, Δt 2 minutes, max(vel): (1.85e-01, 1.99e-02), max(trac): 2.13, 0.00, wtime: 10.987 seconds
[ Info: Time: 11.333 hours, Iteration 340, Δt 2 minutes, max(vel): (1.84e-01, 1.78e-02), max(trac): 2.14, 0.00, wtime: 10.955 seconds
[ Info: Time: 11.667 hours, Iteration 350, Δt 2 minutes, max(vel): (1.85e-01, 1.79e-02), max(trac): 2.14, 0.00, wtime: 10.855 seconds
[ Info: Time: 12 hours, Iteration 360, Δt 2 minutes, max(vel): (1.85e-01, 2.27e-02), max(trac): 2.15, 0.00, wtime: 10.980 seconds
[ Info: Time: 12.333 hours, Iteration 370, Δt 2 minutes, max(vel): (1.82e-01, 2.00e-02), max(trac): 2.16, 0.00, wtime: 11.408 seconds
[ Info: Time: 12.667 hours, Iteration 380, Δt 2 minutes, max(vel): (1.87e-01, 1.92e-02), max(trac): 2.16, 0.00, wtime: 11.486 seconds
[ Info: Time: 13 hours, Iteration 390, Δt 2 minutes, max(vel): (1.87e-01, 1.99e-02), max(trac): 2.17, 0.00, wtime: 11.151 seconds
[ Info: Time: 13.333 hours, Iteration 400, Δt 2 minutes, max(vel): (1.82e-01, 2.40e-02), max(trac): 2.18, 0.00, wtime: 11.089 seconds
[ Info: Time: 13.667 hours, Iteration 410, Δt 2 minutes, max(vel): (1.86e-01, 2.16e-02), max(trac): 2.18, 0.00, wtime: 10.968 seconds
[ Info: Time: 14 hours, Iteration 420, Δt 2 minutes, max(vel): (1.85e-01, 2.14e-02), max(trac): 2.19, 0.00, wtime: 11.243 seconds
[ Info: Time: 14.333 hours, Iteration 430, Δt 2 minutes, max(vel): (1.86e-01, 1.90e-02), max(trac): 2.19, 0.00, wtime: 11.165 seconds
[ Info: Time: 14.667 hours, Iteration 440, Δt 2 minutes, max(vel): (1.87e-01, 2.00e-02), max(trac): 2.19, 0.00, wtime: 11.140 seconds
[ Info: Time: 15 hours, Iteration 450, Δt 2 minutes, max(vel): (1.85e-01, 2.13e-02), max(trac): 2.20, 0.00, wtime: 11.221 seconds
[ Info: Time: 15.333 hours, Iteration 460, Δt 2 minutes, max(vel): (1.82e-01, 2.32e-02), max(trac): 2.20, 0.00, wtime: 11.486 seconds
[ Info: Time: 15.667 hours, Iteration 470, Δt 2 minutes, max(vel): (1.85e-01, 2.06e-02), max(trac): 2.21, 0.00, wtime: 10.954 seconds
[ Info: Time: 16 hours, Iteration 480, Δt 2 minutes, max(vel): (1.84e-01, 1.99e-02), max(trac): 2.21, 0.00, wtime: 10.883 seconds
[ Info: Time: 16.333 hours, Iteration 490, Δt 2 minutes, max(vel): (1.84e-01, 2.00e-02), max(trac): 2.22, 0.00, wtime: 11.181 seconds
[ Info: Time: 16.667 hours, Iteration 500, Δt 2 minutes, max(vel): (1.89e-01, 1.86e-02), max(trac): 2.22, 0.00, wtime: 10.970 seconds
[ Info: Time: 17 hours, Iteration 510, Δt 2 minutes, max(vel): (1.88e-01, 2.03e-02), max(trac): 2.22, 0.00, wtime: 10.880 seconds
[ Info: Time: 17.333 hours, Iteration 520, Δt 2 minutes, max(vel): (1.88e-01, 1.96e-02), max(trac): 2.23, 0.00, wtime: 10.834 seconds
[ Info: Time: 17.667 hours, Iteration 530, Δt 2 minutes, max(vel): (1.86e-01, 1.84e-02), max(trac): 2.23, 0.00, wtime: 10.826 seconds
[ Info: Time: 18 hours, Iteration 540, Δt 2 minutes, max(vel): (1.84e-01, 2.16e-02), max(trac): 2.23, 0.00, wtime: 10.850 seconds
[ Info: Time: 18.333 hours, Iteration 550, Δt 2 minutes, max(vel): (1.84e-01, 1.94e-02), max(trac): 2.23, 0.00, wtime: 10.890 seconds
[ Info: Time: 18.667 hours, Iteration 560, Δt 2 minutes, max(vel): (1.84e-01, 2.40e-02), max(trac): 2.23, 0.00, wtime: 10.771 seconds
[ Info: Time: 19 hours, Iteration 570, Δt 2 minutes, max(vel): (1.84e-01, 1.94e-02), max(trac): 2.24, 0.00, wtime: 10.901 seconds
[ Info: Time: 19.333 hours, Iteration 580, Δt 2 minutes, max(vel): (1.86e-01, 2.07e-02), max(trac): 2.24, 0.00, wtime: 10.878 seconds
[ Info: Time: 19.667 hours, Iteration 590, Δt 2 minutes, max(vel): (1.86e-01, 2.10e-02), max(trac): 2.24, 0.00, wtime: 10.925 seconds
[ Info: Time: 20 hours, Iteration 600, Δt 2 minutes, max(vel): (1.88e-01, 2.31e-02), max(trac): 2.24, 0.00, wtime: 11.020 seconds
[ Info: Time: 20.333 hours, Iteration 610, Δt 2 minutes, max(vel): (1.88e-01, 2.21e-02), max(trac): 2.24, 0.00, wtime: 10.879 seconds
[ Info: Time: 20.667 hours, Iteration 620, Δt 2 minutes, max(vel): (1.83e-01, 2.23e-02), max(trac): 2.24, 0.00, wtime: 10.826 seconds
[ Info: Time: 21 hours, Iteration 630, Δt 2 minutes, max(vel): (1.82e-01, 1.93e-02), max(trac): 2.24, 0.00, wtime: 10.860 seconds
[ Info: Time: 21.333 hours, Iteration 640, Δt 2 minutes, max(vel): (1.82e-01, 2.00e-02), max(trac): 2.24, 0.00, wtime: 10.919 seconds
[ Info: Time: 21.667 hours, Iteration 650, Δt 2 minutes, max(vel): (1.83e-01, 1.91e-02), max(trac): 2.24, 0.00, wtime: 10.960 seconds
[ Info: Time: 22 hours, Iteration 660, Δt 2 minutes, max(vel): (1.85e-01, 2.08e-02), max(trac): 2.24, 0.00, wtime: 10.874 seconds
[ Info: Time: 22.333 hours, Iteration 670, Δt 2 minutes, max(vel): (1.82e-01, 2.29e-02), max(trac): 2.24, 0.00, wtime: 10.825 seconds
[ Info: Time: 22.667 hours, Iteration 680, Δt 2 minutes, max(vel): (1.84e-01, 2.11e-02), max(trac): 2.24, 0.00, wtime: 10.872 seconds
[ Info: Time: 23 hours, Iteration 690, Δt 2 minutes, max(vel): (1.83e-01, 2.21e-02), max(trac): 2.24, 0.00, wtime: 10.875 seconds
[ Info: Time: 23.333 hours, Iteration 700, Δt 2 minutes, max(vel): (1.84e-01, 2.06e-02), max(trac): 2.24, 0.00, wtime: 10.832 seconds
[ Info: Time: 23.667 hours, Iteration 710, Δt 2 minutes, max(vel): (1.87e-01, 1.82e-02), max(trac): 2.24, 0.00, wtime: 10.800 seconds
[ Info: Time: 1 day, Iteration 720, Δt 2 minutes, max(vel): (1.82e-01, 1.88e-02), max(trac): 2.24, 0.00, wtime: 10.879 seconds
[ Info: Time: 1.014 days, Iteration 730, Δt 2 minutes, max(vel): (1.92e-01, 2.29e-02), max(trac): 2.24, 0.00, wtime: 10.949 seconds
[ Info: Time: 1.028 days, Iteration 740, Δt 2 minutes, max(vel): (1.96e-01, 2.19e-02), max(trac): 2.24, 0.00, wtime: 10.883 seconds
[ Info: Time: 1.042 days, Iteration 750, Δt 2 minutes, max(vel): (1.84e-01, 2.30e-02), max(trac): 2.24, 0.00, wtime: 10.824 seconds
[ Info: Time: 1.056 days, Iteration 760, Δt 2 minutes, max(vel): (1.85e-01, 2.41e-02), max(trac): 2.24, 0.00, wtime: 10.823 seconds
[ Info: Time: 1.069 days, Iteration 770, Δt 2 minutes, max(vel): (1.89e-01, 2.29e-02), max(trac): 2.24, 0.00, wtime: 10.895 seconds
[ Info: Time: 1.083 days, Iteration 780, Δt 2 minutes, max(vel): (1.85e-01, 2.01e-02), max(trac): 2.24, 0.00, wtime: 10.843 seconds
[ Info: Time: 1.097 days, Iteration 790, Δt 2 minutes, max(vel): (1.89e-01, 1.94e-02), max(trac): 2.24, 0.00, wtime: 10.776 seconds
[ Info: Time: 1.111 days, Iteration 800, Δt 2 minutes, max(vel): (1.85e-01, 2.44e-02), max(trac): 2.24, 0.00, wtime: 10.768 seconds
[ Info: Time: 1.125 days, Iteration 810, Δt 2 minutes, max(vel): (1.84e-01, 2.14e-02), max(trac): 2.24, 0.00, wtime: 10.764 seconds
[ Info: Time: 1.139 days, Iteration 820, Δt 2 minutes, max(vel): (1.84e-01, 2.26e-02), max(trac): 2.24, 0.00, wtime: 10.855 seconds
[ Info: Time: 1.153 days, Iteration 830, Δt 2 minutes, max(vel): (1.86e-01, 1.95e-02), max(trac): 2.24, 0.00, wtime: 10.902 seconds
[ Info: Time: 1.167 days, Iteration 840, Δt 2 minutes, max(vel): (1.88e-01, 1.95e-02), max(trac): 2.24, 0.00, wtime: 10.923 seconds
[ Info: Time: 1.181 days, Iteration 850, Δt 2 minutes, max(vel): (1.85e-01, 1.99e-02), max(trac): 2.25, 0.00, wtime: 10.958 seconds
[ Info: Time: 1.194 days, Iteration 860, Δt 2 minutes, max(vel): (1.84e-01, 2.24e-02), max(trac): 2.25, 0.00, wtime: 11.030 seconds
[ Info: Time: 1.208 days, Iteration 870, Δt 2 minutes, max(vel): (1.88e-01, 2.25e-02), max(trac): 2.25, 0.00, wtime: 10.881 seconds
[ Info: Time: 1.222 days, Iteration 880, Δt 2 minutes, max(vel): (1.85e-01, 2.13e-02), max(trac): 2.26, 0.00, wtime: 10.894 seconds
[ Info: Time: 1.236 days, Iteration 890, Δt 2 minutes, max(vel): (1.97e-01, 2.20e-02), max(trac): 2.26, 0.00, wtime: 10.901 seconds
[ Info: Time: 1.250 days, Iteration 900, Δt 2 minutes, max(vel): (1.86e-01, 2.05e-02), max(trac): 2.26, 0.00, wtime: 11.080 seconds
[ Info: Time: 1.264 days, Iteration 910, Δt 2 minutes, max(vel): (1.84e-01, 1.92e-02), max(trac): 2.27, 0.00, wtime: 10.983 seconds
[ Info: Time: 1.278 days, Iteration 920, Δt 2 minutes, max(vel): (1.96e-01, 2.19e-02), max(trac): 2.27, 0.00, wtime: 10.928 seconds
[ Info: Time: 1.292 days, Iteration 930, Δt 2 minutes, max(vel): (1.88e-01, 2.02e-02), max(trac): 2.27, 0.00, wtime: 10.914 seconds
[ Info: Time: 1.306 days, Iteration 940, Δt 2 minutes, max(vel): (1.88e-01, 1.90e-02), max(trac): 2.28, 0.00, wtime: 10.866 seconds
[ Info: Time: 1.319 days, Iteration 950, Δt 2 minutes, max(vel): (1.93e-01, 2.39e-02), max(trac): 2.28, 0.00, wtime: 10.851 seconds
[ Info: Time: 1.333 days, Iteration 960, Δt 2 minutes, max(vel): (1.85e-01, 2.34e-02), max(trac): 2.28, 0.00, wtime: 10.862 seconds
[ Info: Time: 1.347 days, Iteration 970, Δt 2 minutes, max(vel): (1.90e-01, 2.13e-02), max(trac): 2.29, 0.00, wtime: 10.851 seconds
[ Info: Time: 1.361 days, Iteration 980, Δt 2 minutes, max(vel): (1.94e-01, 2.27e-02), max(trac): 2.29, 0.00, wtime: 10.898 seconds
[ Info: Time: 1.375 days, Iteration 990, Δt 2 minutes, max(vel): (1.88e-01, 2.25e-02), max(trac): 2.29, 0.00, wtime: 10.868 seconds
[ Info: Time: 1.389 days, Iteration 1000, Δt 2 minutes, max(vel): (1.89e-01, 2.49e-02), max(trac): 2.30, 0.00, wtime: 10.888 seconds
[ Info: Time: 1.403 days, Iteration 1010, Δt 2 minutes, max(vel): (1.89e-01, 2.14e-02), max(trac): 2.30, 0.00, wtime: 11.055 seconds
[ Info: Time: 1.417 days, Iteration 1020, Δt 2 minutes, max(vel): (1.85e-01, 2.00e-02), max(trac): 2.30, 0.00, wtime: 10.982 seconds
[ Info: Time: 1.431 days, Iteration 1030, Δt 2 minutes, max(vel): (1.90e-01, 1.93e-02), max(trac): 2.30, 0.00, wtime: 11.000 seconds
[ Info: Time: 1.444 days, Iteration 1040, Δt 2 minutes, max(vel): (1.91e-01, 2.02e-02), max(trac): 2.31, 0.00, wtime: 11.053 seconds
[ Info: Time: 1.458 days, Iteration 1050, Δt 2 minutes, max(vel): (1.85e-01, 1.98e-02), max(trac): 2.31, 0.00, wtime: 11.214 seconds
[ Info: Time: 1.472 days, Iteration 1060, Δt 2 minutes, max(vel): (1.88e-01, 2.31e-02), max(trac): 2.31, 0.00, wtime: 11.057 seconds
[ Info: Time: 1.486 days, Iteration 1070, Δt 2 minutes, max(vel): (1.86e-01, 2.10e-02), max(trac): 2.32, 0.00, wtime: 11.033 seconds
[ Info: Time: 1.500 days, Iteration 1080, Δt 2 minutes, max(vel): (1.86e-01, 1.91e-02), max(trac): 2.32, 0.00, wtime: 10.908 seconds
[ Info: Time: 1.514 days, Iteration 1090, Δt 2 minutes, max(vel): (1.86e-01, 2.33e-02), max(trac): 2.32, 0.00, wtime: 10.863 seconds
[ Info: Time: 1.528 days, Iteration 1100, Δt 2 minutes, max(vel): (1.87e-01, 2.24e-02), max(trac): 2.32, 0.00, wtime: 10.826 seconds
[ Info: Time: 1.542 days, Iteration 1110, Δt 2 minutes, max(vel): (1.89e-01, 2.17e-02), max(trac): 2.32, 0.00, wtime: 10.880 seconds
[ Info: Time: 1.556 days, Iteration 1120, Δt 2 minutes, max(vel): (1.91e-01, 2.09e-02), max(trac): 2.33, 0.00, wtime: 10.848 seconds
[ Info: Time: 1.569 days, Iteration 1130, Δt 2 minutes, max(vel): (1.86e-01, 2.11e-02), max(trac): 2.33, 0.00, wtime: 10.871 seconds
[ Info: Time: 1.583 days, Iteration 1140, Δt 2 minutes, max(vel): (1.87e-01, 2.31e-02), max(trac): 2.33, 0.00, wtime: 10.813 seconds
[ Info: Time: 1.597 days, Iteration 1150, Δt 2 minutes, max(vel): (1.86e-01, 2.15e-02), max(trac): 2.33, 0.00, wtime: 10.806 seconds
[ Info: Time: 1.611 days, Iteration 1160, Δt 2 minutes, max(vel): (1.85e-01, 2.02e-02), max(trac): 2.33, 0.00, wtime: 10.806 seconds
[ Info: Time: 1.625 days, Iteration 1170, Δt 2 minutes, max(vel): (1.85e-01, 2.14e-02), max(trac): 2.34, 0.00, wtime: 10.962 seconds
[ Info: Time: 1.639 days, Iteration 1180, Δt 2 minutes, max(vel): (1.87e-01, 2.06e-02), max(trac): 2.34, 0.00, wtime: 10.958 seconds
[ Info: Time: 1.653 days, Iteration 1190, Δt 2 minutes, max(vel): (1.86e-01, 2.34e-02), max(trac): 2.34, 0.00, wtime: 10.886 seconds
[ Info: Time: 1.667 days, Iteration 1200, Δt 2 minutes, max(vel): (1.85e-01, 2.58e-02), max(trac): 2.34, 0.00, wtime: 10.931 seconds
[ Info: Time: 1.681 days, Iteration 1210, Δt 2 minutes, max(vel): (1.86e-01, 2.35e-02), max(trac): 2.35, 0.00, wtime: 10.914 seconds
[ Info: Time: 1.694 days, Iteration 1220, Δt 2 minutes, max(vel): (1.84e-01, 2.26e-02), max(trac): 2.35, 0.00, wtime: 10.922 seconds
[ Info: Time: 1.708 days, Iteration 1230, Δt 2 minutes, max(vel): (1.86e-01, 2.18e-02), max(trac): 2.35, 0.00, wtime: 11.063 seconds
[ Info: Time: 1.722 days, Iteration 1240, Δt 2 minutes, max(vel): (1.87e-01, 2.23e-02), max(trac): 2.35, 0.00, wtime: 10.875 seconds
[ Info: Time: 1.736 days, Iteration 1250, Δt 2 minutes, max(vel): (1.84e-01, 2.21e-02), max(trac): 2.35, 0.00, wtime: 10.831 seconds
[ Info: Time: 1.750 days, Iteration 1260, Δt 2 minutes, max(vel): (1.83e-01, 2.26e-02), max(trac): 2.36, 0.00, wtime: 10.919 seconds
[ Info: Time: 1.764 days, Iteration 1270, Δt 2 minutes, max(vel): (1.86e-01, 1.94e-02), max(trac): 2.36, 0.00, wtime: 10.883 seconds
[ Info: Time: 1.778 days, Iteration 1280, Δt 2 minutes, max(vel): (1.85e-01, 1.93e-02), max(trac): 2.36, 0.00, wtime: 10.776 seconds
[ Info: Time: 1.792 days, Iteration 1290, Δt 2 minutes, max(vel): (1.88e-01, 2.31e-02), max(trac): 2.36, 0.00, wtime: 10.832 seconds
[ Info: Time: 1.806 days, Iteration 1300, Δt 2 minutes, max(vel): (1.85e-01, 2.27e-02), max(trac): 2.36, 0.00, wtime: 10.826 seconds
[ Info: Time: 1.819 days, Iteration 1310, Δt 2 minutes, max(vel): (1.87e-01, 2.13e-02), max(trac): 2.36, 0.00, wtime: 10.854 seconds
[ Info: Time: 1.833 days, Iteration 1320, Δt 2 minutes, max(vel): (1.86e-01, 2.11e-02), max(trac): 2.37, 0.00, wtime: 10.897 seconds
[ Info: Time: 1.847 days, Iteration 1330, Δt 2 minutes, max(vel): (1.84e-01, 2.30e-02), max(trac): 2.37, 0.00, wtime: 10.828 seconds
[ Info: Time: 1.861 days, Iteration 1340, Δt 2 minutes, max(vel): (1.84e-01, 2.10e-02), max(trac): 2.37, 0.00, wtime: 10.811 seconds
[ Info: Time: 1.875 days, Iteration 1350, Δt 2 minutes, max(vel): (1.83e-01, 2.25e-02), max(trac): 2.37, 0.00, wtime: 10.843 seconds
[ Info: Time: 1.889 days, Iteration 1360, Δt 2 minutes, max(vel): (1.87e-01, 2.10e-02), max(trac): 2.37, 0.00, wtime: 10.830 seconds
[ Info: Time: 1.903 days, Iteration 1370, Δt 2 minutes, max(vel): (1.88e-01, 2.10e-02), max(trac): 2.37, 0.00, wtime: 10.866 seconds
[ Info: Time: 1.917 days, Iteration 1380, Δt 2 minutes, max(vel): (1.85e-01, 2.31e-02), max(trac): 2.37, 0.00, wtime: 10.828 seconds
[ Info: Time: 1.931 days, Iteration 1390, Δt 2 minutes, max(vel): (1.86e-01, 2.37e-02), max(trac): 2.37, 0.00, wtime: 10.905 seconds
[ Info: Time: 1.944 days, Iteration 1400, Δt 2 minutes, max(vel): (1.87e-01, 2.16e-02), max(trac): 2.37, 0.00, wtime: 10.893 seconds
[ Info: Time: 1.958 days, Iteration 1410, Δt 2 minutes, max(vel): (1.88e-01, 2.16e-02), max(trac): 2.37, 0.00, wtime: 10.839 seconds
[ Info: Time: 1.972 days, Iteration 1420, Δt 2 minutes, max(vel): (1.88e-01, 1.93e-02), max(trac): 2.37, 0.00, wtime: 10.906 seconds
[ Info: Time: 1.986 days, Iteration 1430, Δt 2 minutes, max(vel): (1.88e-01, 1.99e-02), max(trac): 2.38, 0.00, wtime: 11.311 seconds
[ Info: Simulation is stopping after running for 27.194 minutes.
[ Info: Simulation time 2 days equals or exceeds stop time 2 days.
[ Info: Time: 2 days, Iteration 1440, Δt 2 minutes, max(vel): (1.86e-01, 2.15e-02), max(trac): 2.38, 0.00, wtime: 10.885 secondsVisualize!
Nt = length(htimeseries)
iter = Observable(1)
hi = @lift(htimeseries[$iter][:, :, 1])
ℵi = @lift(ℵtimeseries[$iter][:, :, 1])
ui = @lift(utimeseries[$iter][:, :, 1])
vi = @lift(vtimeseries[$iter][:, :, 1])
fig = Figure(size = (1000, 500))
ax = Axis(fig[1, 1], title = "sea ice thickness")
heatmap!(ax, hi, colormap = :magma, colorrange = (0.0, 2.0))
ax = Axis(fig[1, 2], title = "sea ice concentration")
heatmap!(ax, ℵi, colormap = Reverse(:deep), colorrange = (0.0, 1))
ax = Axis(fig[2, 1], title = "zonal velocity")
heatmap!(ax, ui, colorrange = (0, 0.12), colormap = :balance)
ax = Axis(fig[2, 2], title = "meridional velocity")
heatmap!(ax, vi, colorrange = (-0.025, 0.025), colormap = :bwr)
CairoMakie.record(fig, "sea_ice_advected_on_coastline.mp4", 1:Nt, framerate = 8) do i
iter[] = i
@info "doing iter $i"
end[ Info: doing iter 1
[ Info: doing iter 2
[ Info: doing iter 3
[ Info: doing iter 4
[ Info: doing iter 5
[ Info: doing iter 6
[ Info: doing iter 7
[ Info: doing iter 8
[ Info: doing iter 9
[ Info: doing iter 10
[ Info: doing iter 11
[ Info: doing iter 12
[ Info: doing iter 13
[ Info: doing iter 14
[ Info: doing iter 15
[ Info: doing iter 16
[ Info: doing iter 17
[ Info: doing iter 18
[ Info: doing iter 19
[ Info: doing iter 20
[ Info: doing iter 21
[ Info: doing iter 22
[ Info: doing iter 23
[ Info: doing iter 24
[ Info: doing iter 25
[ Info: doing iter 26
[ Info: doing iter 27
[ Info: doing iter 28
[ Info: doing iter 29
[ Info: doing iter 30
[ Info: doing iter 31
[ Info: doing iter 32
[ Info: doing iter 33
[ Info: doing iter 34
[ Info: doing iter 35
[ Info: doing iter 36
[ Info: doing iter 37
[ Info: doing iter 38
[ Info: doing iter 39
[ Info: doing iter 40
[ Info: doing iter 41
[ Info: doing iter 42
[ Info: doing iter 43
[ Info: doing iter 44
[ Info: doing iter 45
[ Info: doing iter 46
[ Info: doing iter 47
[ Info: doing iter 48
[ Info: doing iter 49
[ Info: doing iter 50
[ Info: doing iter 51
[ Info: doing iter 52
[ Info: doing iter 53
[ Info: doing iter 54
[ Info: doing iter 55
[ Info: doing iter 56
[ Info: doing iter 57
[ Info: doing iter 58
[ Info: doing iter 59
[ Info: doing iter 60
[ Info: doing iter 61
[ Info: doing iter 62
[ Info: doing iter 63
[ Info: doing iter 64
[ Info: doing iter 65
[ Info: doing iter 66
[ Info: doing iter 67
[ Info: doing iter 68
[ Info: doing iter 69
[ Info: doing iter 70
[ Info: doing iter 71
[ Info: doing iter 72
[ Info: doing iter 73
[ Info: doing iter 74
[ Info: doing iter 75
[ Info: doing iter 76
[ Info: doing iter 77
[ Info: doing iter 78
[ Info: doing iter 79
[ Info: doing iter 80
[ Info: doing iter 81
[ Info: doing iter 82
[ Info: doing iter 83
[ Info: doing iter 84
[ Info: doing iter 85
[ Info: doing iter 86
[ Info: doing iter 87
[ Info: doing iter 88
[ Info: doing iter 89
[ Info: doing iter 90
[ Info: doing iter 91
[ Info: doing iter 92
[ Info: doing iter 93
[ Info: doing iter 94
[ Info: doing iter 95
[ Info: doing iter 96
[ Info: doing iter 97
[ Info: doing iter 98
[ Info: doing iter 99
[ Info: doing iter 100
[ Info: doing iter 101
[ Info: doing iter 102
[ Info: doing iter 103
[ Info: doing iter 104
[ Info: doing iter 105
[ Info: doing iter 106
[ Info: doing iter 107
[ Info: doing iter 108
[ Info: doing iter 109
[ Info: doing iter 110
[ Info: doing iter 111
[ Info: doing iter 112
[ Info: doing iter 113
[ Info: doing iter 114
[ Info: doing iter 115
[ Info: doing iter 116
[ Info: doing iter 117
[ Info: doing iter 118
[ Info: doing iter 119
[ Info: doing iter 120
[ Info: doing iter 121
[ Info: doing iter 122
[ Info: doing iter 123
[ Info: doing iter 124
[ Info: doing iter 125
[ Info: doing iter 126
[ Info: doing iter 127
[ Info: doing iter 128
[ Info: doing iter 129
[ Info: doing iter 130
[ Info: doing iter 131
[ Info: doing iter 132
[ Info: doing iter 133
[ Info: doing iter 134
[ Info: doing iter 135
[ Info: doing iter 136
[ Info: doing iter 137
[ Info: doing iter 138
[ Info: doing iter 139
[ Info: doing iter 140
[ Info: doing iter 141
[ Info: doing iter 142
[ Info: doing iter 143
[ Info: doing iter 144
[ Info: doing iter 145This page was generated using Literate.jl.