Restart Simulation

This simulation is split into two pieces, run one after the other. The purpose of this example is to show how to use the restart functionality.

using Subzero, CairoMakie, GeoInterfaceMakie
using JLD2, Random, Statistics

const FT = Float64
const Δt = 20
const nΔt = 2500
const nfloes = 20
const L = 1e5
const Δgrid = 1e4
const hmean = 2
const concentration = 0.7
const uomax = 2

dirs = [joinpath("restart_sim", "run" * string(i)) for i in 1:2]

2-element Vector{String}:
 "restart_sim/run1"
 "restart_sim/run2"

Create Grid

grid = RegRectilinearGrid(; x0 = 0.0, xf = L, y0 = 0.0, yf = L, Δx = Δgrid, Δy = Δgrid)

RegRectilinearGrid{Float64}
  ⊢x extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m
  ∟y extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m

Create Domain

nboundary = PeriodicBoundary(North; grid)
sboundary = PeriodicBoundary(South; grid)
eboundary = PeriodicBoundary(East; grid)
wboundary = PeriodicBoundary(West; grid)
domain = Domain(; north = nboundary, south = sboundary, east = eboundary,west =  wboundary)

Domain
  ⊢Northern boundary of type PeriodicBoundary{North, Float64}
  ⊢Southern boundary of type PeriodicBoundary{South, Float64}
  ⊢Eastern boundary of type PeriodicBoundary{East, Float64}
  ⊢Western boundary of type PeriodicBoundary{West, Float64}
  ∟0-element TopograpahyElement{Float64} list

Create Ocean

ngrid = Int(L/Δgrid) + 1
ygrid = range(0,L,ngrid)
uoprofile = @. uomax * (1 - abs(1 - 2 * ygrid/L))
uvels_ocean = repeat(
    uoprofile,
    outer = (1, ngrid),
)
ocean = Ocean(;
    u = uvels_ocean',
    grid,
    v = 0,
    temp = 0,
)

Ocean{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.90909 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C

Create Atmos

atmos = Atmos(FT; grid, u = 0.0, v = 0.0, temp = 0.0)

Atmos{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.0 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C

Create Floes

floe_settings = FloeSettings(subfloe_point_generator = SubGridPointsGenerator(grid, 2))
floe_arr = initialize_floe_field(
    FT,
    nfloes,
    [concentration],
    domain,
    hmean,
    0;
    rng = Xoshiro(1),
    floe_settings = floe_settings
)

21-element StructArray(::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, ::Vector{Vector{Float64}}, ::Vector{Vector{Vector{Vector{Float64}}}}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Vector{Float64}}, ::Vector{Vector{Float64}}, ::Vector{Vector{Float64}}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Subzero.Status}, ::Vector{Int64}, ::Vector{Int64}, ::Vector{Vector{Int64}}, ::Vector{Vector{Int64}}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Matrix{Float64}}, ::Vector{Float64}, ::Vector{Matrix{Float64}}, ::Vector{Int64}, ::Vector{Matrix{Float64}}, ::Vector{Matrix{Float64}}, ::Vector{Matrix{Float64}}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}, ::Vector{Float64}) with eltype Floe{Float64}:
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 15631.577920172793), … (4) … , (100000.0, 15631.577920172793)])]), [85711.61922005899, 27312.213227680793], [[[100000.0, 15631.577920172793], [71716.46244051738, 22675.03020194567], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [100000.0, 28356.11795020577], [100000.0, 15631.577920172793]]], 2.0, 4.69398005598528e8, 8.636923303012916e11, 18455.21786648942, 7.819040385550921e19, [76.01607055889471, 126.25405949075277, 90.67870268750109, 116.79692289860256, 130.2542443642489], [14288.380779941013, 16003.757384271776, 12832.328307457361, 9660.899230642946, 6489.470153828533, 3318.0410770141184, 146.6120001997042, -3024.81707661471, -6196.246153429124, -9367.675230243538  …  6226.312199831872, 9373.463013403258, -9509.441868025058, -6362.291054453672, -3215.1402408822855, -67.98942731089967, 3079.1613862604863, 6226.312199831872, -67.98942731089967, 3079.1613862604863], [-11680.635307508, -12107.815720582692, -11318.034409790465, -10528.25309899824, -9738.471788206014, -8948.69047741379, -8158.909166621563, -7369.127855829337, -6579.346545037111, -5789.565234244886  …  3564.730798936382, 3564.730798936382, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 10303.53037567539, 10303.53037567539], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 1, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 73544.51422997881), … (2) … , (100000.0, 73544.51422997881)])]), [97870.36737262971, 80465.95765510257], [[[100000.0, 73544.51422997881], [93611.10211788886, 80173.12477744419], [100000.0, 87680.23395788489], [100000.0, 73544.51422997881]]], 2.0, 4.5155834915867805e7, 8.308673624519676e10, 7522.042123017303, 8.810650693918952e17, [43.9450383721843, 95.65563805510841, 40.39932357270729], [2129.632627370287, 3356.4063526085815, 1145.79501946165, -1064.8163136852816, -3275.427646832213, -5486.038979979145, -4259.265254740851, -3113.5562055044256, -1064.8163136852818, 983.9235781338623, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 0.0, 0.0], [-6921.443425123754, -8194.245934181998, -5900.6920427865325, -3607.138151391068, -1313.5842599956027, 979.9696313998622, -292.8328776583803, 1053.4029020242754, 3460.721712561972, 5868.0405230996685, 7214.276302782324, 5446.509349815955, 1913.1140758248414, -1620.281198166272, -5153.676472157385, -1620.281198166272, 1913.1140758248414], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 2, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(76898.32736748773, 77881.43845066903), … (5) … , (76898.32736748773, 77881.43845066903)])]), [80594.00652552379, 67170.8397752156], [[[76898.32736748773, 77881.43845066903], [71678.17724422461, 74000.83776654628], [69731.37176368426, 60788.69094406669], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [76898.32736748773, 77881.43845066903]]], 2.0, 3.549290796181326e8, 6.53069506497364e11, 12598.756195436905, 3.7710184125653115e19, [131.74832702379143, 135.00889306377988, 103.44974080292391, 100.31882102703666, 113.24112293558707, 136.23309514688106], [-3695.6791580360587, -6305.754219667615, -8915.829281299171, -8658.131933678354, -9150.57196883475, -9643.012003991149, -10135.452039147545, -10627.892074303943, -11120.33210946034, -10862.634761839523  …  10466.726802718476, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663, 7206.461034119871, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663], [10710.598675453424, 8770.29833339205, 6829.997991330674, 8578.88108548797, 5236.898483329135, 1894.9158811703005, -1447.0667209885346, -4789.049323147369, -8131.031925306204, -6382.148831148908  …  2523.747090724229, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 8942.831722487055, 8942.831722487055, 8942.831722487055, 8942.831722487055], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 3, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (3) … , (0.0, 81608.99581128299)])]), [1992.8823393862158, 72479.44017635356], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [3770.057775250367, 67773.38988696749], [0.0, 67022.0962985524], [0.0, 81608.99581128299]]], 2.0, 4.472412766760963e7, 8.229239490840172e10, 9344.536698515838, 1.0370393406838369e18, [25.720422825686196, 108.82356664552104, 146.72621968139975, 78.72979084739306], [-1992.8823393862158, -1225.7064094673674, 174.08558838475847, 1573.8775862368843, 2973.66958408901, 3740.8455140078586, 1777.175435864151, -107.85345176103237, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, 0.0, 0.0, 0.0, 0.0], [9129.555634929435, 7536.934809825998, 4631.033047725278, 1725.131285624558, -1180.7704764761622, -2773.391301579599, -4706.0502893860685, -5081.697083593615, -5457.343877801162, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601, 7361.788681963066, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 4, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 2108.7909233259197), … (3) … , (100000.0, 2108.7909233259197)])]), [85349.20058304284, 12759.257210627215], [[[100000.0, 2108.7909233259197], [68671.97979189342, 11357.646577717289], [71716.46244051738, 22675.03020194567], [100000.0, 15631.577920172793], [100000.0, 2108.7909233259197]]], 2.0, 3.8259072868817246e8, 7.039669407862373e11, 18112.933381781473, 6.5333077688377205e19, [73.55201323481629, 91.39132522232647, 91.072732101752, 103.98392944110529], [14650.799416957161, 16346.224342551182, 13189.963428069854, 10033.702513588523, 6877.441599107194, 3721.1806846258646, 564.919770144535, -2591.3411443367945, -5747.602058818124, -8903.862973299454  …  12883.032463990792, -14909.45383818305, -11435.39305041132, -7961.332262639589, -4487.271474867859, -1013.2106870961288, 2460.8501006756014, 5934.910888447332, -11435.39305041132, -7961.332262639589], [-10650.466287301295, -11151.00032378084, -10219.18907583989, -9287.377827898938, -8355.566579957987, -7423.755332017036, -6491.944084076085, -5560.132836135135, -4628.321588194184, -3696.510340253233  …  1335.7238892772807, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 8148.006038352085, 8148.006038352085], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 5, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (4) … , (0.0, 81608.99581128299)])]), [8962.33309769216, 87293.65732020503], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [18547.875224604184, 80169.53795708867], [17908.287958045414, 100000.0], [0.0, 100000.0], [0.0, 81608.99581128299]]], 2.0, 4.543822952720574e8, 8.360634233005857e11, 17881.494978141724, 6.7991925253237375e19, [154.27957717431383, 76.48680563883177, 127.38631014576067, 91.84730704109374, 90.0], [-8962.33309769216, -8195.157167773312, -6795.365169921186, -5395.57317206906, -3995.7811742169347, -3228.6052442980863, -1859.3397652074657, 659.5643380497515, 3178.468441306969, 5697.372544564186  …  -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654, -7194.566144725792, -4192.097880991502, -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654], [-5684.661508922043, -7277.28233402548, -10183.1840961262, -13089.08585822692, -15994.98762032764, -17587.608445431077, -16469.52434851764, -14412.69412639568, -12355.86390427372, -10299.033682151758  …  7593.773574416933, 7593.773574416933, 7593.773574416933, 7593.773574416933, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 6, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(60505.550527273575, 0.0), … (3) … , (60505.550527273575, 0.0)])]), [78082.51384882341, 3876.119291709452], [[[60505.550527273575, 0.0], [68671.97979189342, 11357.646577717289], [100000.0, 2108.7909233259197], [100000.0, 0.0], [60505.550527273575, 0.0]]], 2.0, 2.5731412177668297e8, 4.734579840690967e11, 22257.594208507348, 4.574450329071635e19, [54.28292924502844, 109.26908398978784, 106.44798676518374, 90.0], [-17576.963321549832, -16544.970781905584, -14510.822720128468, -12476.674658351354, -10442.526596574238, -9410.534056929988, -7715.109131335967, -4610.979091678351, -1506.8490520207365, 1597.2809876368785  …  10342.742225448308, 13611.73454970228, -12540.204044329492, -9271.21172007552, -6002.219395821549, -2733.2270715675772, 535.7652526863941, 3804.7575769403657, -9271.21172007552, -6002.219395821549], [-3876.119291709452, -2440.8522788792225, 388.18523847305397, 3217.2227558253303, 6046.260273177607, 7481.527286007837, 6980.993249528292, 6064.572407146928, 5148.151564765563, 4231.730722384199  …  499.0185518517671, 499.0185518517671, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 5713.760333041468, 5713.760333041468], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 7, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(25482.063933135447, 48800.33213577528), … (8) … , (25482.063933135447, 48800.33213577528)])]), [39324.736398807945, 42187.873742157804], [[[25482.063933135447, 48800.33213577528], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [54356.76582028708, 48907.36542038052], [39877.20185576307, 60461.079867567416], [39281.31096026047, 60492.26994952193], [25482.063933135447, 48800.33213577528]]], 2.0, 7.194006581353397e8, 1.3236972109690251e12, 19610.54645883592, 1.5746305872438765e20, [124.64444697758991, 153.14239702539754, 118.89501290853333, 171.41958894582058, 144.47982878226506, 133.37283990700973, 132.9075986376685, 144.40873395402195, 136.72955286169338], [-13842.672465672498, -13669.254773411065, -13378.448592749266, -13087.642412087465, -12796.836231425666, -12506.030050763868, -12215.223870102067, -11924.417689440268, -11750.999997178835, -10509.063520727177  …  -2314.7503588359277, 938.6346924541394, 4192.019743744207, 7445.404795034274, 10698.78984632434, -5568.135410125995, -2314.7503588359277, 938.6346924541394, 4192.019743744207, 938.6346924541394], [6612.458393617475, 4853.2181303189345, 1903.1273110670156, -1046.9635081849035, -3997.0543274368224, -6947.145146688741, -9897.235965940661, -12847.32678519258, -14606.56704849112, -16556.965794839445  …  9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 13055.981688145439, 13055.981688145439, 13055.981688145439, 13055.981688145439, 16536.629254397754], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 8, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(75500.23553787598, 100000.0), … (3) … , (75500.23553787598, 100000.0)])]), [81812.79465207118, 93051.11247855646], [[[75500.23553787598, 100000.0], [73438.92036209526, 95968.05099128203], [77330.45737587528, 79599.16016136229], [94990.72859285733, 100000.0], [75500.23553787598, 100000.0]]], 2.0, 2.2352717453713274e8, 4.112900011483242e11, 14897.817985644646, 1.8595750879226687e19, [117.07817808302875, 139.54859765859973, 54.25478202082039, 49.118442237551136], [-6312.559114195203, -7343.216702085563, -8373.874289975924, -7965.001680463247, -7196.55373177458, -6428.105783085914, -5659.657834397247, -4891.20988570858, -4482.3372761959035, -3325.3382472102185  …  2402.0298254051086, 5404.742212876663, 8407.454600348217, -6606.107337009555, -3603.394949538, -600.6825620664458, 2402.0298254051086, 5404.742212876663, 8407.454600348217, 11410.166987819772], [6948.887521443539, 4932.913017084553, 2916.9385127255664, 1197.1061984605117, -2035.2003518868946, -5267.506902234301, -8499.813452581708, -11732.120002929114, -13451.95231719417, -12115.407051620952  …  1808.0593819361761, 1808.0593819361761, 1808.0593819361761, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 9, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 25731.58204587066), … (2) … , (0.0, 25731.58204587066)])]), [4111.848557080553, 18450.34292670981], [[[0.0, 25731.58204587066], [12335.54567124166, 21299.167510608822], [0.0, 8320.279223649946], [0.0, 25731.58204587066]]], 2.0, 1.0738896057966119e8, 1.975956874665766e11, 10932.771340534675, 4.3667696262100403e18, [70.23559643234292, 66.22020874088905, 43.54419482676801], [-4111.848557080553, -2448.218922594947, 554.5432114952022, 3557.305345585351, 6560.0674796755, 8223.697114161107, 9441.536313379724, 7331.361446282739, 5221.1865791857535, 3111.011712088769  …  589.2556509887896, -2344.0816041141843, 589.2556509887896, 3522.592906091764, -2344.0816041141843, 589.2556509887896, 3522.592906091764, 6455.930161194738, -2344.0816041141843, 589.2556509887896], [7281.239119160848, 6683.462874014483, 5604.508859024781, 4525.5548440350785, 3446.6008290453765, 2848.8245838990115, 4130.178432313476, 1909.950434629504, -310.27756305446826, -2530.5055607384406  …  -4893.354521021502, -1424.4122919495085, -1424.4122919495085, -1424.4122919495085, 2044.529937122485, 2044.529937122485, 2044.529937122485, 2044.529937122485, 5513.472166194479, 5513.472166194479], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 10, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 ⋮
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 40825.89817681637), … (6) … , (0.0, 40825.89817681637)])]), [12697.167115723327, 31418.32758700577], [[[0.0, 40825.89817681637], [4420.1954736259295, 41993.70052733918], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [26292.921098828716, 22649.172216637835], [17962.641911132076, 23027.3563608404], [0.0, 32207.178831597743], [0.0, 40825.89817681637]]], 2.0, 2.7430885631593776e8, 5.047282956213255e11, 17126.016877750342, 3.914516269592779e19, [104.79922854795807, 133.2997300542375, 154.38840686641623, 107.27913607125417, 127.63412732875173, 155.53003091488236, 117.06934021650007], [-12697.167115723327, -10487.06937891036, -8276.971642097396, -6776.2045279990125, -3897.346428026668, -1018.4883280543233, 1860.3697719180213, 4739.227871890366, 7618.085971862711, 10496.944071835054  …  -7769.382791238454, -4609.36541971995, -1449.3480482014465, 1710.6693233170572, -10929.400162756958, -7769.382791238454, -4609.36541971995, -1449.3480482014465, -10929.400162756958, -7769.382791238454], [9407.570589810603, 9991.471765072005, 10575.372940333407, 9641.18983770266, 7849.185890786803, 6057.181943870945, 4265.177996955089, 2473.174050039232, 681.170103123375, -1110.8338437924817  …  2484.008225459597, 2484.008225459597, 2484.008225459597, 2484.008225459597, 5645.807106413317, 5645.807106413317, 5645.807106413317, 5645.807106413317, 8807.605987367038, 8807.605987367038], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 13, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(18708.048610863716, 0.0), … (4) … , (18708.048610863716, 0.0)])]), [29692.671074503003, 10025.935716776967], [[[18708.048610863716, 0.0], [26292.921098828716, 22649.172216637835], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [40564.17237212511, 0.0], [18708.048610863716, 0.0]]], 2.0, 3.4089674709385103e8, 6.272500146526859e11, 15997.228292747279, 4.233716017246853e19, [71.48506477943155, 158.2814367604344, 133.8258510202125, 104.5093946598428, 71.8982527800788], [-10984.622463639287, -10423.264807987005, -9346.238611876932, -8269.21241576686, -7192.186219656787, -6115.160023546715, -5038.133827436642, -3961.1076313265694, -3399.7499756742873, -2138.3740860482612  …  2996.8710595461866, 6050.302702100963, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -56.560583008589674, 2996.8710595461866], [-10025.935716776967, -8349.666769847336, -5133.5610493842405, -1917.4553289211447, 1298.6503915419507, 4514.756112005046, 7730.861832468142, 10946.967552931237, 12623.236499860868, 14114.10436520113  …  4486.231424657013, 4486.231424657013, 7672.331471773916, 7672.331471773916, 7672.331471773916, 10858.431518890819, 10858.431518890819, 10858.431518890819, 14044.531566007721, 14044.531566007721], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 14, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(69731.37176368426, 60788.69094406669), … (6) … , (69731.37176368426, 60788.69094406669)])]), [71288.1799442355, 45938.89146732441], [[[69731.37176368426, 60788.69094406669], [54356.76582028708, 48907.36542038052], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [69731.37176368426, 60788.69094406669]]], 2.0, 6.692501784369707e8, 1.2314203283240261e12, 17364.42744277742, 1.3688265153115973e20, [120.47160929676905, 123.37635799052929, 126.97662415698366, 124.39463803085083, 115.79395304455706, 133.56558891941265, 155.4212285608975], [-1556.8081805512338, -158.04244405797817, -2754.0620749699347, -5350.081705881891, -7946.101336793848, -10542.120967705803, -13138.14059861776, -15734.160229529716, -18330.179860441673, -16931.414123948416  …  7917.446636196438, 11214.74575150765, 14512.044866818862, -5271.7498250484105, -1974.4507097371986, 1322.8484055740134, 4620.147520885225, 7917.446636196438, -1974.4507097371986, 1322.8484055740134], [14849.799476742286, 15930.750176086378, 13924.57490146147, 11918.399626836563, 9912.224352211653, 7906.049077586746, 5899.873802961838, 3893.6985283369295, 1887.5232537120216, 2968.473953056113  …  6620.7549527512165, 6620.7549527512165, 6620.7549527512165, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 13082.032523775917, 13082.032523775917], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 15, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 67567.56940096187), … (4) … , (100000.0, 67567.56940096187)])]), [94831.18867433425, 73451.81728375112], [[[100000.0, 67567.56940096187], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [93611.10211788886, 80173.12477744419], [100000.0, 73544.51422997881], [100000.0, 67567.56940096187]]], 2.0, 8.539646907159662e7, 1.571295030917378e11, 7832.048497452011, 2.4418149897730017e18, [77.18179333098198, 137.42632491867568, 97.73618770158065, 91.60073242094602, 136.0549616278157], [5168.8113256657525, 6892.5235894391235, 4237.792527854595, 1583.0614662700682, -1071.6695953144595, -3726.400656898987, -2002.6883931256161, -998.6668209792022, -2587.3824831875054, -4176.098145395808  …  -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -4581.740944699632, -1920.8125055666267, 740.1159335663784, -1920.8125055666267], [-5884.247882789248, -6276.441752855163, -5672.414461926796, -5068.3871709984305, -4464.359880070064, -3860.3325891416976, -4252.526459207613, -5707.497483889769, -3405.220999899324, -1102.944515908879  …  -4116.480929822879, -4116.480929822879, -4116.480929822879, -1093.1404396396877, -1093.1404396396877, -1093.1404396396877, 1930.2000505435037, 1930.2000505435037, 1930.2000505435037, 4953.540540726695], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 16, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(53932.759746554315, 0.0), … (6) … , (53932.759746554315, 0.0)])]), [58425.22292764254, 20191.465410638164], [[[53932.759746554315, 0.0], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [71716.46244051738, 22675.03020194567], [68671.97979189342, 11357.646577717289], [60505.550527273575, 0.0], [53932.759746554315, 0.0]]], 2.0, 6.839042436264261e8, 1.258383808272624e12, 20685.200044051562, 1.565246038690451e20, [107.66171295242654, 120.03122187956625, 99.65053593600663, 144.92665928164806, 142.67320840749525, 159.33959078788573, 125.71707075497156], [-4492.463181088227, -3956.1290668852325, -4973.423632261886, -5990.71819763854, -7008.012763015193, -8025.307328391847, -9042.6018937685, -10059.896459145153, -11077.191024521808, -12094.48558989846  …  -8489.481337283707, -5153.989021085012, -1818.496704886315, 1516.9956113123812, 4852.487927511078, 8187.980243709774, -5153.989021085012, -1818.496704886315, 1516.9956113123812, -1818.496704886315], [-20191.465410638164, -21875.90766786816, -18680.932915017933, -15485.958162167708, -12290.983409317481, -9096.008656467255, -5901.033903617029, -2706.0591507668028, 488.91560208342327, 3683.890354933649  …  11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 14545.034352869088, 14545.034352869088, 14545.034352869088, 17841.907633923176], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 17, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(17908.287958045414, 100000.0), … (6) … , (17908.287958045414, 100000.0)])]), [33613.80575505652, 83906.94777996265], [[[17908.287958045414, 100000.0], [18547.875224604184, 80169.53795708867], [26710.280846516955, 68796.48411425078], [39281.31096026047, 60492.26994952193], [39877.20185576307, 60461.079867567416], [50658.43944984143, 89053.30431356005], [43884.94062412003, 100000.0], [17908.287958045414, 100000.0]]], 2.0, 9.177048547118509e8, 1.6885769326698057e12, 24268.06240472919, 2.6469913646102482e20, [88.15269295890626, 146.1804201768482, 159.11497680870042, 149.54814772225257, 107.66365980992538, 127.59207823601393, 121.74802428735329], [-15705.517797011107, -15648.532053579796, -15543.408897640567, -15438.285741701337, -15333.162585762107, -15228.039429822877, -15122.916273883648, -15065.930530452337, -14035.195499735168, -12001.550312909023  …  8784.729500515592, 12030.798121167067, -13937.750844044738, -10691.682223393262, -7445.613602741787, -4199.544982090311, -953.4763614388353, 2292.5922592126403, 5538.660879864116, 8784.729500515592], [16093.052220037353, 14326.204002151735, 11066.850880723716, 7807.4977592956975, 4548.144637867679, 1288.7915164396604, -1970.561604988358, -3737.409822873975, -5173.5801806144755, -8007.151223066773  …  11052.250155570997, 11052.250155570997, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 18, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 28356.11795020577), … (3) … , (100000.0, 28356.11795020577)])]), [94387.04219792974, 42856.57222855929], [[[100000.0, 28356.11795020577], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [100000.0, 53790.42066586246], [100000.0, 28356.11795020577]]], 2.0, 2.3895189855925107e8, 4.39671493349022e11, 15548.905735338496, 2.075893184387935e19, [49.74575563575111, 127.4091240568404, 104.1820578473644, 78.66306246004409], [5612.957802070261, 6962.090324090315, 4341.691870581147, 1721.2934170719775, -899.1050364371913, -3519.5034899463603, -6139.901943455529, -8760.300396964698, -7411.167874944644, -7323.422442346604  …  -2480.536998284219, 682.3269254098365, 3845.190849103892, -5643.400921978275, -2480.536998284219, 682.3269254098365, 3845.190849103892, -2480.536998284219, 682.3269254098365, 3845.190849103892], [-14500.454278353522, -15642.750841435021, -13424.086552721725, -11205.422264008428, -8986.757975295131, -6768.093686581835, -4549.429397868538, -2330.7651091552416, -3473.0616722367413, -1707.4737420961371  …  2909.290395844241, 2909.290395844241, 2909.290395844241, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 9166.081484336799, 9166.081484336799, 9166.081484336799], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 19, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 53790.42066586246), … (4) … , (100000.0, 53790.42066586246)])]), [93357.10571133424, 59553.21688748059], [[[100000.0, 53790.42066586246], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [100000.0, 67567.56940096187], [100000.0, 53790.42066586246]]], 2.0, 1.9638967050932837e8, 3.613569937371642e11, 10409.509630176246, 1.272724957257055e19, [101.3369375399559, 122.25235323322298, 104.25995041206583, 109.33255214573725, 102.81820666901802], [6642.8942886657605, 8376.168805988078, 5196.457351863297, 2016.7458977385152, -1162.9655563862664, -4342.677010511048, -7522.388464635829, -5789.113947313512, -7593.368392765311, -9397.62283821711  …  -4503.610080013208, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392], [-5762.796221618126, -5415.291620783709, -6052.792690920758, -6690.293761057807, -7327.794831194857, -7965.295901331906, -8602.796971468955, -8255.29237063454, -6359.929304333658, -4464.566238032778  …  2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 5005.14050374361, 5005.14050374361, 5005.14050374361, 7878.3069840965545, 7878.3069840965545, 7878.3069840965545], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 20, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)
 Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(40564.17237212511, 0.0), … (3) … , (40564.17237212511, 0.0)])]), [42958.001115460844, 13654.895707898613], [[[40564.17237212511, 0.0], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [53932.759746554315, 0.0], [40564.17237212511, 0.0]]], 2.0, 3.6688421815656054e8, 6.750669614080714e11, 17518.604505049727, 5.57569422307862e19, [108.1017472199212, 84.07101639433667, 95.48894933816871, 72.33828704757347], [-2393.8287433357327, -1844.5739782078895, -2905.0871299171754, -3965.600281626461, -5026.113433335747, -6086.626585045033, -7147.139736754319, -8207.652888463605, -9268.16604017289, -10328.679191882176  …  67.41052631492676, 3113.937576918985, -9072.170625497249, -6025.64357489319, -2979.1165242891316, 67.41052631492676, 3113.937576918985, -6025.64357489319, -2979.1165242891316, 67.41052631492676], [-13654.895707898613, -15335.169257944664, -12090.860220628758, -8846.551183312851, -5602.242145996945, -2357.933108681038, 886.3759286348686, 4130.684965950775, 7374.994003266682, 10619.303040582588  …  6896.868076149904, 6896.868076149904, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 13158.200353177286, 13158.200353177286, 13158.200353177286], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 21, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

Create Model

model = Model(grid, ocean, atmos, domain, floe_arr)

Model{Float64, RegRectilinearGrid{Float64}, Domain{Float64, PeriodicBoundary{North, Float64}, PeriodicBoundary{South, Float64}, PeriodicBoundary{East, Float64}, PeriodicBoundary{West, Float64}, StructArrays.StructVector{TopographyElement{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, rmax::Vector{Float64}}, Int64}}, StructArrays.StructVector{Floe{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, coords::Vector{Vector{Vector{Vector{Float64}}}}, height::Vector{Float64}, area::Vector{Float64}, mass::Vector{Float64}, rmax::Vector{Float64}, moment::Vector{Float64}, angles::Vector{Vector{Float64}}, x_subfloe_points::Vector{Vector{Float64}}, y_subfloe_points::Vector{Vector{Float64}}, α::Vector{Float64}, u::Vector{Float64}, v::Vector{Float64}, ξ::Vector{Float64}, status::Vector{Subzero.Status}, id::Vector{Int64}, ghost_id::Vector{Int64}, parent_ids::Vector{Vector{Int64}}, ghosts::Vector{Vector{Int64}}, fxOA::Vector{Float64}, fyOA::Vector{Float64}, trqOA::Vector{Float64}, hflx_factor::Vector{Float64}, overarea::Vector{Float64}, collision_force::Vector{Matrix{Float64}}, collision_trq::Vector{Float64}, interactions::Vector{Matrix{Float64}}, num_inters::Vector{Int64}, stress_accum::Vector{Matrix{Float64}}, stress_instant::Vector{Matrix{Float64}}, strain::Vector{Matrix{Float64}}, damage::Vector{Float64}, p_dxdt::Vector{Float64}, p_dydt::Vector{Float64}, p_dudt::Vector{Float64}, p_dvdt::Vector{Float64}, p_dξdt::Vector{Float64}, p_dαdt::Vector{Float64}}, Int64}}(RegRectilinearGrid{Float64}
  ⊢x extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m
  ∟y extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m, Ocean{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.90909 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C, Atmos{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.0 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C, Domain
  ⊢Northern boundary of type PeriodicBoundary{North, Float64}
  ⊢Southern boundary of type PeriodicBoundary{South, Float64}
  ⊢Eastern boundary of type PeriodicBoundary{East, Float64}
  ⊢Western boundary of type PeriodicBoundary{West, Float64}
  ∟0-element TopograpahyElement{Float64} list, Floe{Float64}[Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 15631.577920172793), … (4) … , (100000.0, 15631.577920172793)])]), [85711.61922005899, 27312.213227680793], [[[100000.0, 15631.577920172793], [71716.46244051738, 22675.03020194567], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [100000.0, 28356.11795020577], [100000.0, 15631.577920172793]]], 2.0, 4.69398005598528e8, 8.636923303012916e11, 18455.21786648942, 7.819040385550921e19, [76.01607055889471, 126.25405949075277, 90.67870268750109, 116.79692289860256, 130.2542443642489], [14288.380779941013, 16003.757384271776, 12832.328307457361, 9660.899230642946, 6489.470153828533, 3318.0410770141184, 146.6120001997042, -3024.81707661471, -6196.246153429124, -9367.675230243538  …  6226.312199831872, 9373.463013403258, -9509.441868025058, -6362.291054453672, -3215.1402408822855, -67.98942731089967, 3079.1613862604863, 6226.312199831872, -67.98942731089967, 3079.1613862604863], [-11680.635307508, -12107.815720582692, -11318.034409790465, -10528.25309899824, -9738.471788206014, -8948.69047741379, -8158.909166621563, -7369.127855829337, -6579.346545037111, -5789.565234244886  …  3564.730798936382, 3564.730798936382, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 10303.53037567539, 10303.53037567539], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 1, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 73544.51422997881), … (2) … , (100000.0, 73544.51422997881)])]), [97870.36737262971, 80465.95765510257], [[[100000.0, 73544.51422997881], [93611.10211788886, 80173.12477744419], [100000.0, 87680.23395788489], [100000.0, 73544.51422997881]]], 2.0, 4.5155834915867805e7, 8.308673624519676e10, 7522.042123017303, 8.810650693918952e17, [43.9450383721843, 95.65563805510841, 40.39932357270729], [2129.632627370287, 3356.4063526085815, 1145.79501946165, -1064.8163136852816, -3275.427646832213, -5486.038979979145, -4259.265254740851, -3113.5562055044256, -1064.8163136852818, 983.9235781338623, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 0.0, 0.0], [-6921.443425123754, -8194.245934181998, -5900.6920427865325, -3607.138151391068, -1313.5842599956027, 979.9696313998622, -292.8328776583803, 1053.4029020242754, 3460.721712561972, 5868.0405230996685, 7214.276302782324, 5446.509349815955, 1913.1140758248414, -1620.281198166272, -5153.676472157385, -1620.281198166272, 1913.1140758248414], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 2, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(76898.32736748773, 77881.43845066903), … (5) … , (76898.32736748773, 77881.43845066903)])]), [80594.00652552379, 67170.8397752156], [[[76898.32736748773, 77881.43845066903], [71678.17724422461, 74000.83776654628], [69731.37176368426, 60788.69094406669], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [76898.32736748773, 77881.43845066903]]], 2.0, 3.549290796181326e8, 6.53069506497364e11, 12598.756195436905, 3.7710184125653115e19, [131.74832702379143, 135.00889306377988, 103.44974080292391, 100.31882102703666, 113.24112293558707, 136.23309514688106], [-3695.6791580360587, -6305.754219667615, -8915.829281299171, -8658.131933678354, -9150.57196883475, -9643.012003991149, -10135.452039147545, -10627.892074303943, -11120.33210946034, -10862.634761839523  …  10466.726802718476, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663, 7206.461034119871, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663], [10710.598675453424, 8770.29833339205, 6829.997991330674, 8578.88108548797, 5236.898483329135, 1894.9158811703005, -1447.0667209885346, -4789.049323147369, -8131.031925306204, -6382.148831148908  …  2523.747090724229, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 8942.831722487055, 8942.831722487055, 8942.831722487055, 8942.831722487055], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 3, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (3) … , (0.0, 81608.99581128299)])]), [1992.8823393862158, 72479.44017635356], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [3770.057775250367, 67773.38988696749], [0.0, 67022.0962985524], [0.0, 81608.99581128299]]], 2.0, 4.472412766760963e7, 8.229239490840172e10, 9344.536698515838, 1.0370393406838369e18, [25.720422825686196, 108.82356664552104, 146.72621968139975, 78.72979084739306], [-1992.8823393862158, -1225.7064094673674, 174.08558838475847, 1573.8775862368843, 2973.66958408901, 3740.8455140078586, 1777.175435864151, -107.85345176103237, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, 0.0, 0.0, 0.0, 0.0], [9129.555634929435, 7536.934809825998, 4631.033047725278, 1725.131285624558, -1180.7704764761622, -2773.391301579599, -4706.0502893860685, -5081.697083593615, -5457.343877801162, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601, 7361.788681963066, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 4, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 2108.7909233259197), … (3) … , (100000.0, 2108.7909233259197)])]), [85349.20058304284, 12759.257210627215], [[[100000.0, 2108.7909233259197], [68671.97979189342, 11357.646577717289], [71716.46244051738, 22675.03020194567], [100000.0, 15631.577920172793], [100000.0, 2108.7909233259197]]], 2.0, 3.8259072868817246e8, 7.039669407862373e11, 18112.933381781473, 6.5333077688377205e19, [73.55201323481629, 91.39132522232647, 91.072732101752, 103.98392944110529], [14650.799416957161, 16346.224342551182, 13189.963428069854, 10033.702513588523, 6877.441599107194, 3721.1806846258646, 564.919770144535, -2591.3411443367945, -5747.602058818124, -8903.862973299454  …  12883.032463990792, -14909.45383818305, -11435.39305041132, -7961.332262639589, -4487.271474867859, -1013.2106870961288, 2460.8501006756014, 5934.910888447332, -11435.39305041132, -7961.332262639589], [-10650.466287301295, -11151.00032378084, -10219.18907583989, -9287.377827898938, -8355.566579957987, -7423.755332017036, -6491.944084076085, -5560.132836135135, -4628.321588194184, -3696.510340253233  …  1335.7238892772807, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 8148.006038352085, 8148.006038352085], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 5, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (4) … , (0.0, 81608.99581128299)])]), [8962.33309769216, 87293.65732020503], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [18547.875224604184, 80169.53795708867], [17908.287958045414, 100000.0], [0.0, 100000.0], [0.0, 81608.99581128299]]], 2.0, 4.543822952720574e8, 8.360634233005857e11, 17881.494978141724, 6.7991925253237375e19, [154.27957717431383, 76.48680563883177, 127.38631014576067, 91.84730704109374, 90.0], [-8962.33309769216, -8195.157167773312, -6795.365169921186, -5395.57317206906, -3995.7811742169347, -3228.6052442980863, -1859.3397652074657, 659.5643380497515, 3178.468441306969, 5697.372544564186  …  -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654, -7194.566144725792, -4192.097880991502, -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654], [-5684.661508922043, -7277.28233402548, -10183.1840961262, -13089.08585822692, -15994.98762032764, -17587.608445431077, -16469.52434851764, -14412.69412639568, -12355.86390427372, -10299.033682151758  …  7593.773574416933, 7593.773574416933, 7593.773574416933, 7593.773574416933, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 6, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(60505.550527273575, 0.0), … (3) … , (60505.550527273575, 0.0)])]), [78082.51384882341, 3876.119291709452], [[[60505.550527273575, 0.0], [68671.97979189342, 11357.646577717289], [100000.0, 2108.7909233259197], [100000.0, 0.0], [60505.550527273575, 0.0]]], 2.0, 2.5731412177668297e8, 4.734579840690967e11, 22257.594208507348, 4.574450329071635e19, [54.28292924502844, 109.26908398978784, 106.44798676518374, 90.0], [-17576.963321549832, -16544.970781905584, -14510.822720128468, -12476.674658351354, -10442.526596574238, -9410.534056929988, -7715.109131335967, -4610.979091678351, -1506.8490520207365, 1597.2809876368785  …  10342.742225448308, 13611.73454970228, -12540.204044329492, -9271.21172007552, -6002.219395821549, -2733.2270715675772, 535.7652526863941, 3804.7575769403657, -9271.21172007552, -6002.219395821549], [-3876.119291709452, -2440.8522788792225, 388.18523847305397, 3217.2227558253303, 6046.260273177607, 7481.527286007837, 6980.993249528292, 6064.572407146928, 5148.151564765563, 4231.730722384199  …  499.0185518517671, 499.0185518517671, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 5713.760333041468, 5713.760333041468], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 7, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(25482.063933135447, 48800.33213577528), … (8) … , (25482.063933135447, 48800.33213577528)])]), [39324.736398807945, 42187.873742157804], [[[25482.063933135447, 48800.33213577528], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [54356.76582028708, 48907.36542038052], [39877.20185576307, 60461.079867567416], [39281.31096026047, 60492.26994952193], [25482.063933135447, 48800.33213577528]]], 2.0, 7.194006581353397e8, 1.3236972109690251e12, 19610.54645883592, 1.5746305872438765e20, [124.64444697758991, 153.14239702539754, 118.89501290853333, 171.41958894582058, 144.47982878226506, 133.37283990700973, 132.9075986376685, 144.40873395402195, 136.72955286169338], [-13842.672465672498, -13669.254773411065, -13378.448592749266, -13087.642412087465, -12796.836231425666, -12506.030050763868, -12215.223870102067, -11924.417689440268, -11750.999997178835, -10509.063520727177  …  -2314.7503588359277, 938.6346924541394, 4192.019743744207, 7445.404795034274, 10698.78984632434, -5568.135410125995, -2314.7503588359277, 938.6346924541394, 4192.019743744207, 938.6346924541394], [6612.458393617475, 4853.2181303189345, 1903.1273110670156, -1046.9635081849035, -3997.0543274368224, -6947.145146688741, -9897.235965940661, -12847.32678519258, -14606.56704849112, -16556.965794839445  …  9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 13055.981688145439, 13055.981688145439, 13055.981688145439, 13055.981688145439, 16536.629254397754], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 8, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(75500.23553787598, 100000.0), … (3) … , (75500.23553787598, 100000.0)])]), [81812.79465207118, 93051.11247855646], [[[75500.23553787598, 100000.0], [73438.92036209526, 95968.05099128203], [77330.45737587528, 79599.16016136229], [94990.72859285733, 100000.0], [75500.23553787598, 100000.0]]], 2.0, 2.2352717453713274e8, 4.112900011483242e11, 14897.817985644646, 1.8595750879226687e19, [117.07817808302875, 139.54859765859973, 54.25478202082039, 49.118442237551136], [-6312.559114195203, -7343.216702085563, -8373.874289975924, -7965.001680463247, -7196.55373177458, -6428.105783085914, -5659.657834397247, -4891.20988570858, -4482.3372761959035, -3325.3382472102185  …  2402.0298254051086, 5404.742212876663, 8407.454600348217, -6606.107337009555, -3603.394949538, -600.6825620664458, 2402.0298254051086, 5404.742212876663, 8407.454600348217, 11410.166987819772], [6948.887521443539, 4932.913017084553, 2916.9385127255664, 1197.1061984605117, -2035.2003518868946, -5267.506902234301, -8499.813452581708, -11732.120002929114, -13451.95231719417, -12115.407051620952  …  1808.0593819361761, 1808.0593819361761, 1808.0593819361761, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 9, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 25731.58204587066), … (2) … , (0.0, 25731.58204587066)])]), [4111.848557080553, 18450.34292670981], [[[0.0, 25731.58204587066], [12335.54567124166, 21299.167510608822], [0.0, 8320.279223649946], [0.0, 25731.58204587066]]], 2.0, 1.0738896057966119e8, 1.975956874665766e11, 10932.771340534675, 4.3667696262100403e18, [70.23559643234292, 66.22020874088905, 43.54419482676801], [-4111.848557080553, -2448.218922594947, 554.5432114952022, 3557.305345585351, 6560.0674796755, 8223.697114161107, 9441.536313379724, 7331.361446282739, 5221.1865791857535, 3111.011712088769  …  589.2556509887896, -2344.0816041141843, 589.2556509887896, 3522.592906091764, -2344.0816041141843, 589.2556509887896, 3522.592906091764, 6455.930161194738, -2344.0816041141843, 589.2556509887896], [7281.239119160848, 6683.462874014483, 5604.508859024781, 4525.5548440350785, 3446.6008290453765, 2848.8245838990115, 4130.178432313476, 1909.950434629504, -310.27756305446826, -2530.5055607384406  …  -4893.354521021502, -1424.4122919495085, -1424.4122919495085, -1424.4122919495085, 2044.529937122485, 2044.529937122485, 2044.529937122485, 2044.529937122485, 5513.472166194479, 5513.472166194479], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 10, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)  …  Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 32207.178831597743), … (3) … , (0.0, 32207.178831597743)])]), [7157.5241000569895, 25954.12483010887], [[[0.0, 32207.178831597743], [17962.641911132076, 23027.3563608404], [12335.54567124166, 21299.167510608822], [0.0, 25731.58204587066], [0.0, 32207.178831597743]]], 2.0, 8.128930093923163e7, 1.4957231372818622e11, 11194.48724077687, 3.866678601257591e18, [62.93065978349994, 44.142063486863776, 143.16287316197923, 109.76440356765707], [-7157.5241000569895, -5583.404622328917, -2620.524031193731, 342.35655994145526, 3305.2371510766416, 6268.117742211827, 9230.998333347014, 10805.117811075086, 7991.569691129878, 5178.02157118467  …  6151.929257068849, 9037.350858108717, -5389.7571470906205, -2504.335546050753, 381.0860549891146, 3266.507656028982, -5389.7571470906205, -2504.335546050753, 381.0860549891146, -5389.7571470906205], [6253.054001488872, 5448.598987310266, 3934.41649883024, 2420.234010350214, 906.0515218701876, -608.1309666098387, -2122.313455089865, -2926.7684692684707, -3790.8628943842596, -4654.9573195000485  …  -2887.1903665336795, -2887.1903665336795, -429.69789484828533, -429.69789484828533, -429.69789484828533, -429.69789484828533, 2027.794576837109, 2027.794576837109, 2027.794576837109, 4485.287048522503], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 12, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 40825.89817681637), … (6) … , (0.0, 40825.89817681637)])]), [12697.167115723327, 31418.32758700577], [[[0.0, 40825.89817681637], [4420.1954736259295, 41993.70052733918], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [26292.921098828716, 22649.172216637835], [17962.641911132076, 23027.3563608404], [0.0, 32207.178831597743], [0.0, 40825.89817681637]]], 2.0, 2.7430885631593776e8, 5.047282956213255e11, 17126.016877750342, 3.914516269592779e19, [104.79922854795807, 133.2997300542375, 154.38840686641623, 107.27913607125417, 127.63412732875173, 155.53003091488236, 117.06934021650007], [-12697.167115723327, -10487.06937891036, -8276.971642097396, -6776.2045279990125, -3897.346428026668, -1018.4883280543233, 1860.3697719180213, 4739.227871890366, 7618.085971862711, 10496.944071835054  …  -7769.382791238454, -4609.36541971995, -1449.3480482014465, 1710.6693233170572, -10929.400162756958, -7769.382791238454, -4609.36541971995, -1449.3480482014465, -10929.400162756958, -7769.382791238454], [9407.570589810603, 9991.471765072005, 10575.372940333407, 9641.18983770266, 7849.185890786803, 6057.181943870945, 4265.177996955089, 2473.174050039232, 681.170103123375, -1110.8338437924817  …  2484.008225459597, 2484.008225459597, 2484.008225459597, 2484.008225459597, 5645.807106413317, 5645.807106413317, 5645.807106413317, 5645.807106413317, 8807.605987367038, 8807.605987367038], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 13, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(18708.048610863716, 0.0), … (4) … , (18708.048610863716, 0.0)])]), [29692.671074503003, 10025.935716776967], [[[18708.048610863716, 0.0], [26292.921098828716, 22649.172216637835], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [40564.17237212511, 0.0], [18708.048610863716, 0.0]]], 2.0, 3.4089674709385103e8, 6.272500146526859e11, 15997.228292747279, 4.233716017246853e19, [71.48506477943155, 158.2814367604344, 133.8258510202125, 104.5093946598428, 71.8982527800788], [-10984.622463639287, -10423.264807987005, -9346.238611876932, -8269.21241576686, -7192.186219656787, -6115.160023546715, -5038.133827436642, -3961.1076313265694, -3399.7499756742873, -2138.3740860482612  …  2996.8710595461866, 6050.302702100963, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -56.560583008589674, 2996.8710595461866], [-10025.935716776967, -8349.666769847336, -5133.5610493842405, -1917.4553289211447, 1298.6503915419507, 4514.756112005046, 7730.861832468142, 10946.967552931237, 12623.236499860868, 14114.10436520113  …  4486.231424657013, 4486.231424657013, 7672.331471773916, 7672.331471773916, 7672.331471773916, 10858.431518890819, 10858.431518890819, 10858.431518890819, 14044.531566007721, 14044.531566007721], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 14, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(69731.37176368426, 60788.69094406669), … (6) … , (69731.37176368426, 60788.69094406669)])]), [71288.1799442355, 45938.89146732441], [[[69731.37176368426, 60788.69094406669], [54356.76582028708, 48907.36542038052], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [69731.37176368426, 60788.69094406669]]], 2.0, 6.692501784369707e8, 1.2314203283240261e12, 17364.42744277742, 1.3688265153115973e20, [120.47160929676905, 123.37635799052929, 126.97662415698366, 124.39463803085083, 115.79395304455706, 133.56558891941265, 155.4212285608975], [-1556.8081805512338, -158.04244405797817, -2754.0620749699347, -5350.081705881891, -7946.101336793848, -10542.120967705803, -13138.14059861776, -15734.160229529716, -18330.179860441673, -16931.414123948416  …  7917.446636196438, 11214.74575150765, 14512.044866818862, -5271.7498250484105, -1974.4507097371986, 1322.8484055740134, 4620.147520885225, 7917.446636196438, -1974.4507097371986, 1322.8484055740134], [14849.799476742286, 15930.750176086378, 13924.57490146147, 11918.399626836563, 9912.224352211653, 7906.049077586746, 5899.873802961838, 3893.6985283369295, 1887.5232537120216, 2968.473953056113  …  6620.7549527512165, 6620.7549527512165, 6620.7549527512165, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 13082.032523775917, 13082.032523775917], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 15, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 67567.56940096187), … (4) … , (100000.0, 67567.56940096187)])]), [94831.18867433425, 73451.81728375112], [[[100000.0, 67567.56940096187], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [93611.10211788886, 80173.12477744419], [100000.0, 73544.51422997881], [100000.0, 67567.56940096187]]], 2.0, 8.539646907159662e7, 1.571295030917378e11, 7832.048497452011, 2.4418149897730017e18, [77.18179333098198, 137.42632491867568, 97.73618770158065, 91.60073242094602, 136.0549616278157], [5168.8113256657525, 6892.5235894391235, 4237.792527854595, 1583.0614662700682, -1071.6695953144595, -3726.400656898987, -2002.6883931256161, -998.6668209792022, -2587.3824831875054, -4176.098145395808  …  -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -4581.740944699632, -1920.8125055666267, 740.1159335663784, -1920.8125055666267], [-5884.247882789248, -6276.441752855163, -5672.414461926796, -5068.3871709984305, -4464.359880070064, -3860.3325891416976, -4252.526459207613, -5707.497483889769, -3405.220999899324, -1102.944515908879  …  -4116.480929822879, -4116.480929822879, -4116.480929822879, -1093.1404396396877, -1093.1404396396877, -1093.1404396396877, 1930.2000505435037, 1930.2000505435037, 1930.2000505435037, 4953.540540726695], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 16, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(53932.759746554315, 0.0), … (6) … , (53932.759746554315, 0.0)])]), [58425.22292764254, 20191.465410638164], [[[53932.759746554315, 0.0], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [71716.46244051738, 22675.03020194567], [68671.97979189342, 11357.646577717289], [60505.550527273575, 0.0], [53932.759746554315, 0.0]]], 2.0, 6.839042436264261e8, 1.258383808272624e12, 20685.200044051562, 1.565246038690451e20, [107.66171295242654, 120.03122187956625, 99.65053593600663, 144.92665928164806, 142.67320840749525, 159.33959078788573, 125.71707075497156], [-4492.463181088227, -3956.1290668852325, -4973.423632261886, -5990.71819763854, -7008.012763015193, -8025.307328391847, -9042.6018937685, -10059.896459145153, -11077.191024521808, -12094.48558989846  …  -8489.481337283707, -5153.989021085012, -1818.496704886315, 1516.9956113123812, 4852.487927511078, 8187.980243709774, -5153.989021085012, -1818.496704886315, 1516.9956113123812, -1818.496704886315], [-20191.465410638164, -21875.90766786816, -18680.932915017933, -15485.958162167708, -12290.983409317481, -9096.008656467255, -5901.033903617029, -2706.0591507668028, 488.91560208342327, 3683.890354933649  …  11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 14545.034352869088, 14545.034352869088, 14545.034352869088, 17841.907633923176], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 17, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(17908.287958045414, 100000.0), … (6) … , (17908.287958045414, 100000.0)])]), [33613.80575505652, 83906.94777996265], [[[17908.287958045414, 100000.0], [18547.875224604184, 80169.53795708867], [26710.280846516955, 68796.48411425078], [39281.31096026047, 60492.26994952193], [39877.20185576307, 60461.079867567416], [50658.43944984143, 89053.30431356005], [43884.94062412003, 100000.0], [17908.287958045414, 100000.0]]], 2.0, 9.177048547118509e8, 1.6885769326698057e12, 24268.06240472919, 2.6469913646102482e20, [88.15269295890626, 146.1804201768482, 159.11497680870042, 149.54814772225257, 107.66365980992538, 127.59207823601393, 121.74802428735329], [-15705.517797011107, -15648.532053579796, -15543.408897640567, -15438.285741701337, -15333.162585762107, -15228.039429822877, -15122.916273883648, -15065.930530452337, -14035.195499735168, -12001.550312909023  …  8784.729500515592, 12030.798121167067, -13937.750844044738, -10691.682223393262, -7445.613602741787, -4199.544982090311, -953.4763614388353, 2292.5922592126403, 5538.660879864116, 8784.729500515592], [16093.052220037353, 14326.204002151735, 11066.850880723716, 7807.4977592956975, 4548.144637867679, 1288.7915164396604, -1970.561604988358, -3737.409822873975, -5173.5801806144755, -8007.151223066773  …  11052.250155570997, 11052.250155570997, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 18, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 28356.11795020577), … (3) … , (100000.0, 28356.11795020577)])]), [94387.04219792974, 42856.57222855929], [[[100000.0, 28356.11795020577], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [100000.0, 53790.42066586246], [100000.0, 28356.11795020577]]], 2.0, 2.3895189855925107e8, 4.39671493349022e11, 15548.905735338496, 2.075893184387935e19, [49.74575563575111, 127.4091240568404, 104.1820578473644, 78.66306246004409], [5612.957802070261, 6962.090324090315, 4341.691870581147, 1721.2934170719775, -899.1050364371913, -3519.5034899463603, -6139.901943455529, -8760.300396964698, -7411.167874944644, -7323.422442346604  …  -2480.536998284219, 682.3269254098365, 3845.190849103892, -5643.400921978275, -2480.536998284219, 682.3269254098365, 3845.190849103892, -2480.536998284219, 682.3269254098365, 3845.190849103892], [-14500.454278353522, -15642.750841435021, -13424.086552721725, -11205.422264008428, -8986.757975295131, -6768.093686581835, -4549.429397868538, -2330.7651091552416, -3473.0616722367413, -1707.4737420961371  …  2909.290395844241, 2909.290395844241, 2909.290395844241, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 9166.081484336799, 9166.081484336799, 9166.081484336799], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 19, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 53790.42066586246), … (4) … , (100000.0, 53790.42066586246)])]), [93357.10571133424, 59553.21688748059], [[[100000.0, 53790.42066586246], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [100000.0, 67567.56940096187], [100000.0, 53790.42066586246]]], 2.0, 1.9638967050932837e8, 3.613569937371642e11, 10409.509630176246, 1.272724957257055e19, [101.3369375399559, 122.25235323322298, 104.25995041206583, 109.33255214573725, 102.81820666901802], [6642.8942886657605, 8376.168805988078, 5196.457351863297, 2016.7458977385152, -1162.9655563862664, -4342.677010511048, -7522.388464635829, -5789.113947313512, -7593.368392765311, -9397.62283821711  …  -4503.610080013208, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392], [-5762.796221618126, -5415.291620783709, -6052.792690920758, -6690.293761057807, -7327.794831194857, -7965.295901331906, -8602.796971468955, -8255.29237063454, -6359.929304333658, -4464.566238032778  …  2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 5005.14050374361, 5005.14050374361, 5005.14050374361, 7878.3069840965545, 7878.3069840965545, 7878.3069840965545], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 20, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(40564.17237212511, 0.0), … (3) … , (40564.17237212511, 0.0)])]), [42958.001115460844, 13654.895707898613], [[[40564.17237212511, 0.0], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [53932.759746554315, 0.0], [40564.17237212511, 0.0]]], 2.0, 3.6688421815656054e8, 6.750669614080714e11, 17518.604505049727, 5.57569422307862e19, [108.1017472199212, 84.07101639433667, 95.48894933816871, 72.33828704757347], [-2393.8287433357327, -1844.5739782078895, -2905.0871299171754, -3965.600281626461, -5026.113433335747, -6086.626585045033, -7147.139736754319, -8207.652888463605, -9268.16604017289, -10328.679191882176  …  67.41052631492676, 3113.937576918985, -9072.170625497249, -6025.64357489319, -2979.1165242891316, 67.41052631492676, 3113.937576918985, -6025.64357489319, -2979.1165242891316, 67.41052631492676], [-13654.895707898613, -15335.169257944664, -12090.860220628758, -8846.551183312851, -5602.242145996945, -2357.933108681038, 886.3759286348686, 4130.684965950775, 7374.994003266682, 10619.303040582588  …  6896.868076149904, 6896.868076149904, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 13158.200353177286, 13158.200353177286, 13158.200353177286], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 21, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)])

Create Outout Writers

initwriter = InitialStateOutputWriter(dir = dirs[1], overwrite = true)
checkpointer = CheckpointOutputWriter(
    250,
    dir = dirs[1],
    filename = "checkpoint.jld2",
    overwrite = true,
    jld2_kw = Dict{Symbol, Any}(),
)
floewriter = FloeOutputWriter(50, dir = dirs[1], overwrite = true)
writers = OutputWriters(initwriter, floewriter, checkpointer)

OutputWriters{StructArrays.StructVector{InitialStateOutputWriter, @NamedTuple{filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{FloeOutputWriter, @NamedTuple{Δtout::Vector{Int64}, outputs::Vector{Vector{Symbol}}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{GridOutputWriter, @NamedTuple{outputs::Vector{Vector{Symbol}}, Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}, xg::Vector{Vector{FT} where FT<:AbstractFloat}, yg::Vector{Vector{FT} where FT<:AbstractFloat}, data::Vector{Array{FT, 3} where FT<:AbstractFloat}, average::Vector{Bool}}, Int64}, StructArrays.StructVector{CheckpointOutputWriter, @NamedTuple{Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}}(InitialStateOutputWriter[InitialStateOutputWriter("restart_sim/run1/initial_state.jld2", true)], FloeOutputWriter[FloeOutputWriter(50, [:poly, :centroid, :coords, :height, :area, :mass, :rmax, :moment, :angles, :x_subfloe_points  …  :stress_accum, :stress_instant, :strain, :damage, :p_dxdt, :p_dydt, :p_dudt, :p_dvdt, :p_dξdt, :p_dαdt], "restart_sim/run1/floes.jld2", true)], GridOutputWriter[], CheckpointOutputWriter[CheckpointOutputWriter(250, "restart_sim/run1/checkpoint.jld2", true)])

Create Simulation and Constants

modulus = 1.5e3*(mean(sqrt.(floe_arr.area)) + minimum(sqrt.(floe_arr.area)))
consts = Constants(E = modulus, f = 0, turnθ = 0)
simulation = Simulation(
    model = model,
    consts = consts,
    Δt = Δt,
    nΔt = nΔt,
    verbose = true,
    writers = writers,
    rng = Xoshiro(1),
)

Simulation{Float64, Model{Float64, RegRectilinearGrid{Float64}, Domain{Float64, PeriodicBoundary{North, Float64}, PeriodicBoundary{South, Float64}, PeriodicBoundary{East, Float64}, PeriodicBoundary{West, Float64}, StructArrays.StructVector{TopographyElement{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, rmax::Vector{Float64}}, Int64}}, StructArrays.StructVector{Floe{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, coords::Vector{Vector{Vector{Vector{Float64}}}}, height::Vector{Float64}, area::Vector{Float64}, mass::Vector{Float64}, rmax::Vector{Float64}, moment::Vector{Float64}, angles::Vector{Vector{Float64}}, x_subfloe_points::Vector{Vector{Float64}}, y_subfloe_points::Vector{Vector{Float64}}, α::Vector{Float64}, u::Vector{Float64}, v::Vector{Float64}, ξ::Vector{Float64}, status::Vector{Subzero.Status}, id::Vector{Int64}, ghost_id::Vector{Int64}, parent_ids::Vector{Vector{Int64}}, ghosts::Vector{Vector{Int64}}, fxOA::Vector{Float64}, fyOA::Vector{Float64}, trqOA::Vector{Float64}, hflx_factor::Vector{Float64}, overarea::Vector{Float64}, collision_force::Vector{Matrix{Float64}}, collision_trq::Vector{Float64}, interactions::Vector{Matrix{Float64}}, num_inters::Vector{Int64}, stress_accum::Vector{Matrix{Float64}}, stress_instant::Vector{Matrix{Float64}}, strain::Vector{Matrix{Float64}}, damage::Vector{Float64}, p_dxdt::Vector{Float64}, p_dydt::Vector{Float64}, p_dudt::Vector{Float64}, p_dvdt::Vector{Float64}, p_dξdt::Vector{Float64}, p_dαdt::Vector{Float64}}, Int64}}, NoFracture, MonteCarloPointsGenerator{Float64}, DecayAreaScaledCalculator{Float64}, Random.Xoshiro, OutputWriters{StructArrays.StructVector{InitialStateOutputWriter, @NamedTuple{filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{FloeOutputWriter, @NamedTuple{Δtout::Vector{Int64}, outputs::Vector{Vector{Symbol}}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{GridOutputWriter, @NamedTuple{outputs::Vector{Vector{Symbol}}, Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}, xg::Vector{Vector{FT} where FT<:AbstractFloat}, yg::Vector{Vector{FT} where FT<:AbstractFloat}, data::Vector{Array{FT, 3} where FT<:AbstractFloat}, average::Vector{Bool}}, Int64}, StructArrays.StructVector{CheckpointOutputWriter, @NamedTuple{Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}}}(Model{Float64, RegRectilinearGrid{Float64}, Domain{Float64, PeriodicBoundary{North, Float64}, PeriodicBoundary{South, Float64}, PeriodicBoundary{East, Float64}, PeriodicBoundary{West, Float64}, StructArrays.StructVector{TopographyElement{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, rmax::Vector{Float64}}, Int64}}, StructArrays.StructVector{Floe{Float64}, @NamedTuple{poly::Vector{GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}}, centroid::Vector{Vector{Float64}}, coords::Vector{Vector{Vector{Vector{Float64}}}}, height::Vector{Float64}, area::Vector{Float64}, mass::Vector{Float64}, rmax::Vector{Float64}, moment::Vector{Float64}, angles::Vector{Vector{Float64}}, x_subfloe_points::Vector{Vector{Float64}}, y_subfloe_points::Vector{Vector{Float64}}, α::Vector{Float64}, u::Vector{Float64}, v::Vector{Float64}, ξ::Vector{Float64}, status::Vector{Subzero.Status}, id::Vector{Int64}, ghost_id::Vector{Int64}, parent_ids::Vector{Vector{Int64}}, ghosts::Vector{Vector{Int64}}, fxOA::Vector{Float64}, fyOA::Vector{Float64}, trqOA::Vector{Float64}, hflx_factor::Vector{Float64}, overarea::Vector{Float64}, collision_force::Vector{Matrix{Float64}}, collision_trq::Vector{Float64}, interactions::Vector{Matrix{Float64}}, num_inters::Vector{Int64}, stress_accum::Vector{Matrix{Float64}}, stress_instant::Vector{Matrix{Float64}}, strain::Vector{Matrix{Float64}}, damage::Vector{Float64}, p_dxdt::Vector{Float64}, p_dydt::Vector{Float64}, p_dudt::Vector{Float64}, p_dvdt::Vector{Float64}, p_dξdt::Vector{Float64}, p_dαdt::Vector{Float64}}, Int64}}(RegRectilinearGrid{Float64}
  ⊢x extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m
  ∟y extent (0.0 to 100000.0) with 10 grid cells of size 10000.0 m, Ocean{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.90909 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C, Atmos{Float64}
  ⊢Vector fields of dimension (11, 11)
  ⊢Tracer fields of dimension (11, 11)
  ⊢Average u-velocity of: 0.0 m/s
  ⊢Average v-velocity of: 0.0 m/s
  ∟Average temperature of: 0.0 C, Domain
  ⊢Northern boundary of type PeriodicBoundary{North, Float64}
  ⊢Southern boundary of type PeriodicBoundary{South, Float64}
  ⊢Eastern boundary of type PeriodicBoundary{East, Float64}
  ⊢Western boundary of type PeriodicBoundary{West, Float64}
  ∟0-element TopograpahyElement{Float64} list, Floe{Float64}[Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 15631.577920172793), … (4) … , (100000.0, 15631.577920172793)])]), [85711.61922005899, 27312.213227680793], [[[100000.0, 15631.577920172793], [71716.46244051738, 22675.03020194567], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [100000.0, 28356.11795020577], [100000.0, 15631.577920172793]]], 2.0, 4.69398005598528e8, 8.636923303012916e11, 18455.21786648942, 7.819040385550921e19, [76.01607055889471, 126.25405949075277, 90.67870268750109, 116.79692289860256, 130.2542443642489], [14288.380779941013, 16003.757384271776, 12832.328307457361, 9660.899230642946, 6489.470153828533, 3318.0410770141184, 146.6120001997042, -3024.81707661471, -6196.246153429124, -9367.675230243538  …  6226.312199831872, 9373.463013403258, -9509.441868025058, -6362.291054453672, -3215.1402408822855, -67.98942731089967, 3079.1613862604863, 6226.312199831872, -67.98942731089967, 3079.1613862604863], [-11680.635307508, -12107.815720582692, -11318.034409790465, -10528.25309899824, -9738.471788206014, -8948.69047741379, -8158.909166621563, -7369.127855829337, -6579.346545037111, -5789.565234244886  …  3564.730798936382, 3564.730798936382, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 6934.1305873058855, 10303.53037567539, 10303.53037567539], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 1, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 73544.51422997881), … (2) … , (100000.0, 73544.51422997881)])]), [97870.36737262971, 80465.95765510257], [[[100000.0, 73544.51422997881], [93611.10211788886, 80173.12477744419], [100000.0, 87680.23395788489], [100000.0, 73544.51422997881]]], 2.0, 4.5155834915867805e7, 8.308673624519676e10, 7522.042123017303, 8.810650693918952e17, [43.9450383721843, 95.65563805510841, 40.39932357270729], [2129.632627370287, 3356.4063526085815, 1145.79501946165, -1064.8163136852816, -3275.427646832213, -5486.038979979145, -4259.265254740851, -3113.5562055044256, -1064.8163136852818, 983.9235781338623, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 2129.632627370287, 0.0, 0.0], [-6921.443425123754, -8194.245934181998, -5900.6920427865325, -3607.138151391068, -1313.5842599956027, 979.9696313998622, -292.8328776583803, 1053.4029020242754, 3460.721712561972, 5868.0405230996685, 7214.276302782324, 5446.509349815955, 1913.1140758248414, -1620.281198166272, -5153.676472157385, -1620.281198166272, 1913.1140758248414], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 2, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(76898.32736748773, 77881.43845066903), … (5) … , (76898.32736748773, 77881.43845066903)])]), [80594.00652552379, 67170.8397752156], [[[76898.32736748773, 77881.43845066903], [71678.17724422461, 74000.83776654628], [69731.37176368426, 60788.69094406669], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [76898.32736748773, 77881.43845066903]]], 2.0, 3.549290796181326e8, 6.53069506497364e11, 12598.756195436905, 3.7710184125653115e19, [131.74832702379143, 135.00889306377988, 103.44974080292391, 100.31882102703666, 113.24112293558707, 136.23309514688106], [-3695.6791580360587, -6305.754219667615, -8915.829281299171, -8658.131933678354, -9150.57196883475, -9643.012003991149, -10135.452039147545, -10627.892074303943, -11120.33210946034, -10862.634761839523  …  10466.726802718476, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663, 7206.461034119871, -5834.602040274549, -2574.336271675944, 685.9294969226612, 3946.1952655212663], [10710.598675453424, 8770.29833339205, 6829.997991330674, 8578.88108548797, 5236.898483329135, 1894.9158811703005, -1447.0667209885346, -4789.049323147369, -8131.031925306204, -6382.148831148908  …  2523.747090724229, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 5733.289406605642, 8942.831722487055, 8942.831722487055, 8942.831722487055, 8942.831722487055], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 3, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (3) … , (0.0, 81608.99581128299)])]), [1992.8823393862158, 72479.44017635356], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [3770.057775250367, 67773.38988696749], [0.0, 67022.0962985524], [0.0, 81608.99581128299]]], 2.0, 4.472412766760963e7, 8.229239490840172e10, 9344.536698515838, 1.0370393406838369e18, [25.720422825686196, 108.82356664552104, 146.72621968139975, 78.72979084739306], [-1992.8823393862158, -1225.7064094673674, 174.08558838475847, 1573.8775862368843, 2973.66958408901, 3740.8455140078586, 1777.175435864151, -107.85345176103237, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, -1992.8823393862158, 0.0, 0.0, 0.0, 0.0], [9129.555634929435, 7536.934809825998, 4631.033047725278, 1725.131285624558, -1180.7704764761622, -2773.391301579599, -4706.0502893860685, -5081.697083593615, -5457.343877801162, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601, 7361.788681963066, -3689.5769248347933, -926.7355231353285, 1836.1058785641362, 4598.947280263601], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 4, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 2108.7909233259197), … (3) … , (100000.0, 2108.7909233259197)])]), [85349.20058304284, 12759.257210627215], [[[100000.0, 2108.7909233259197], [68671.97979189342, 11357.646577717289], [71716.46244051738, 22675.03020194567], [100000.0, 15631.577920172793], [100000.0, 2108.7909233259197]]], 2.0, 3.8259072868817246e8, 7.039669407862373e11, 18112.933381781473, 6.5333077688377205e19, [73.55201323481629, 91.39132522232647, 91.072732101752, 103.98392944110529], [14650.799416957161, 16346.224342551182, 13189.963428069854, 10033.702513588523, 6877.441599107194, 3721.1806846258646, 564.919770144535, -2591.3411443367945, -5747.602058818124, -8903.862973299454  …  12883.032463990792, -14909.45383818305, -11435.39305041132, -7961.332262639589, -4487.271474867859, -1013.2106870961288, 2460.8501006756014, 5934.910888447332, -11435.39305041132, -7961.332262639589], [-10650.466287301295, -11151.00032378084, -10219.18907583989, -9287.377827898938, -8355.566579957987, -7423.755332017036, -6491.944084076085, -5560.132836135135, -4628.321588194184, -3696.510340253233  …  1335.7238892772807, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 4741.864963814683, 8148.006038352085, 8148.006038352085], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 5, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 81608.99581128299), … (4) … , (0.0, 81608.99581128299)])]), [8962.33309769216, 87293.65732020503], [[[0.0, 81608.99581128299], [5733.727853394074, 69706.04887477396], [18547.875224604184, 80169.53795708867], [17908.287958045414, 100000.0], [0.0, 100000.0], [0.0, 81608.99581128299]]], 2.0, 4.543822952720574e8, 8.360634233005857e11, 17881.494978141724, 6.7991925253237375e19, [154.27957717431383, 76.48680563883177, 127.38631014576067, 91.84730704109374, 90.0], [-8962.33309769216, -8195.157167773312, -6795.365169921186, -5395.57317206906, -3995.7811742169347, -3228.6052442980863, -1859.3397652074657, 659.5643380497515, 3178.468441306969, 5697.372544564186  …  -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654, -7194.566144725792, -4192.097880991502, -1189.6296172572133, 1812.838646477076, 4815.306910211365, 7817.775173945654], [-5684.661508922043, -7277.28233402548, -10183.1840961262, -13089.08585822692, -15994.98762032764, -17587.608445431077, -16469.52434851764, -14412.69412639568, -12355.86390427372, -10299.033682151758  …  7593.773574416933, 7593.773574416933, 7593.773574416933, 7593.773574416933, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596, 10938.575726828596], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 6, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(60505.550527273575, 0.0), … (3) … , (60505.550527273575, 0.0)])]), [78082.51384882341, 3876.119291709452], [[[60505.550527273575, 0.0], [68671.97979189342, 11357.646577717289], [100000.0, 2108.7909233259197], [100000.0, 0.0], [60505.550527273575, 0.0]]], 2.0, 2.5731412177668297e8, 4.734579840690967e11, 22257.594208507348, 4.574450329071635e19, [54.28292924502844, 109.26908398978784, 106.44798676518374, 90.0], [-17576.963321549832, -16544.970781905584, -14510.822720128468, -12476.674658351354, -10442.526596574238, -9410.534056929988, -7715.109131335967, -4610.979091678351, -1506.8490520207365, 1597.2809876368785  …  10342.742225448308, 13611.73454970228, -12540.204044329492, -9271.21172007552, -6002.219395821549, -2733.2270715675772, 535.7652526863941, 3804.7575769403657, -9271.21172007552, -6002.219395821549], [-3876.119291709452, -2440.8522788792225, 388.18523847305397, 3217.2227558253303, 6046.260273177607, 7481.527286007837, 6980.993249528292, 6064.572407146928, 5148.151564765563, 4231.730722384199  …  499.0185518517671, 499.0185518517671, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 3106.3894424466175, 5713.760333041468, 5713.760333041468], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 7, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(25482.063933135447, 48800.33213577528), … (8) … , (25482.063933135447, 48800.33213577528)])]), [39324.736398807945, 42187.873742157804], [[[25482.063933135447, 48800.33213577528], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [54356.76582028708, 48907.36542038052], [39877.20185576307, 60461.079867567416], [39281.31096026047, 60492.26994952193], [25482.063933135447, 48800.33213577528]]], 2.0, 7.194006581353397e8, 1.3236972109690251e12, 19610.54645883592, 1.5746305872438765e20, [124.64444697758991, 153.14239702539754, 118.89501290853333, 171.41958894582058, 144.47982878226506, 133.37283990700973, 132.9075986376685, 144.40873395402195, 136.72955286169338], [-13842.672465672498, -13669.254773411065, -13378.448592749266, -13087.642412087465, -12796.836231425666, -12506.030050763868, -12215.223870102067, -11924.417689440268, -11750.999997178835, -10509.063520727177  …  -2314.7503588359277, 938.6346924541394, 4192.019743744207, 7445.404795034274, 10698.78984632434, -5568.135410125995, -2314.7503588359277, 938.6346924541394, 4192.019743744207, 938.6346924541394], [6612.458393617475, 4853.2181303189345, 1903.1273110670156, -1046.9635081849035, -3997.0543274368224, -6947.145146688741, -9897.235965940661, -12847.32678519258, -14606.56704849112, -16556.965794839445  …  9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 9575.334121893125, 13055.981688145439, 13055.981688145439, 13055.981688145439, 13055.981688145439, 16536.629254397754], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 8, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(75500.23553787598, 100000.0), … (3) … , (75500.23553787598, 100000.0)])]), [81812.79465207118, 93051.11247855646], [[[75500.23553787598, 100000.0], [73438.92036209526, 95968.05099128203], [77330.45737587528, 79599.16016136229], [94990.72859285733, 100000.0], [75500.23553787598, 100000.0]]], 2.0, 2.2352717453713274e8, 4.112900011483242e11, 14897.817985644646, 1.8595750879226687e19, [117.07817808302875, 139.54859765859973, 54.25478202082039, 49.118442237551136], [-6312.559114195203, -7343.216702085563, -8373.874289975924, -7965.001680463247, -7196.55373177458, -6428.105783085914, -5659.657834397247, -4891.20988570858, -4482.3372761959035, -3325.3382472102185  …  2402.0298254051086, 5404.742212876663, 8407.454600348217, -6606.107337009555, -3603.394949538, -600.6825620664458, 2402.0298254051086, 5404.742212876663, 8407.454600348217, 11410.166987819772], [6948.887521443539, 4932.913017084553, 2916.9385127255664, 1197.1061984605117, -2035.2003518868946, -5267.506902234301, -8499.813452581708, -11732.120002929114, -13451.95231719417, -12115.407051620952  …  1808.0593819361761, 1808.0593819361761, 1808.0593819361761, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717, 5181.12056847717], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 9, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 25731.58204587066), … (2) … , (0.0, 25731.58204587066)])]), [4111.848557080553, 18450.34292670981], [[[0.0, 25731.58204587066], [12335.54567124166, 21299.167510608822], [0.0, 8320.279223649946], [0.0, 25731.58204587066]]], 2.0, 1.0738896057966119e8, 1.975956874665766e11, 10932.771340534675, 4.3667696262100403e18, [70.23559643234292, 66.22020874088905, 43.54419482676801], [-4111.848557080553, -2448.218922594947, 554.5432114952022, 3557.305345585351, 6560.0674796755, 8223.697114161107, 9441.536313379724, 7331.361446282739, 5221.1865791857535, 3111.011712088769  …  589.2556509887896, -2344.0816041141843, 589.2556509887896, 3522.592906091764, -2344.0816041141843, 589.2556509887896, 3522.592906091764, 6455.930161194738, -2344.0816041141843, 589.2556509887896], [7281.239119160848, 6683.462874014483, 5604.508859024781, 4525.5548440350785, 3446.6008290453765, 2848.8245838990115, 4130.178432313476, 1909.950434629504, -310.27756305446826, -2530.5055607384406  …  -4893.354521021502, -1424.4122919495085, -1424.4122919495085, -1424.4122919495085, 2044.529937122485, 2044.529937122485, 2044.529937122485, 2044.529937122485, 5513.472166194479, 5513.472166194479], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 10, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)  …  Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 32207.178831597743), … (3) … , (0.0, 32207.178831597743)])]), [7157.5241000569895, 25954.12483010887], [[[0.0, 32207.178831597743], [17962.641911132076, 23027.3563608404], [12335.54567124166, 21299.167510608822], [0.0, 25731.58204587066], [0.0, 32207.178831597743]]], 2.0, 8.128930093923163e7, 1.4957231372818622e11, 11194.48724077687, 3.866678601257591e18, [62.93065978349994, 44.142063486863776, 143.16287316197923, 109.76440356765707], [-7157.5241000569895, -5583.404622328917, -2620.524031193731, 342.35655994145526, 3305.2371510766416, 6268.117742211827, 9230.998333347014, 10805.117811075086, 7991.569691129878, 5178.02157118467  …  6151.929257068849, 9037.350858108717, -5389.7571470906205, -2504.335546050753, 381.0860549891146, 3266.507656028982, -5389.7571470906205, -2504.335546050753, 381.0860549891146, -5389.7571470906205], [6253.054001488872, 5448.598987310266, 3934.41649883024, 2420.234010350214, 906.0515218701876, -608.1309666098387, -2122.313455089865, -2926.7684692684707, -3790.8628943842596, -4654.9573195000485  …  -2887.1903665336795, -2887.1903665336795, -429.69789484828533, -429.69789484828533, -429.69789484828533, -429.69789484828533, 2027.794576837109, 2027.794576837109, 2027.794576837109, 4485.287048522503], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 12, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(0.0, 40825.89817681637), … (6) … , (0.0, 40825.89817681637)])]), [12697.167115723327, 31418.32758700577], [[[0.0, 40825.89817681637], [4420.1954736259295, 41993.70052733918], [27573.73640162911, 27581.306693666684], [28815.672878080768, 25630.90794731836], [26292.921098828716, 22649.172216637835], [17962.641911132076, 23027.3563608404], [0.0, 32207.178831597743], [0.0, 40825.89817681637]]], 2.0, 2.7430885631593776e8, 5.047282956213255e11, 17126.016877750342, 3.914516269592779e19, [104.79922854795807, 133.2997300542375, 154.38840686641623, 107.27913607125417, 127.63412732875173, 155.53003091488236, 117.06934021650007], [-12697.167115723327, -10487.06937891036, -8276.971642097396, -6776.2045279990125, -3897.346428026668, -1018.4883280543233, 1860.3697719180213, 4739.227871890366, 7618.085971862711, 10496.944071835054  …  -7769.382791238454, -4609.36541971995, -1449.3480482014465, 1710.6693233170572, -10929.400162756958, -7769.382791238454, -4609.36541971995, -1449.3480482014465, -10929.400162756958, -7769.382791238454], [9407.570589810603, 9991.471765072005, 10575.372940333407, 9641.18983770266, 7849.185890786803, 6057.181943870945, 4265.177996955089, 2473.174050039232, 681.170103123375, -1110.8338437924817  …  2484.008225459597, 2484.008225459597, 2484.008225459597, 2484.008225459597, 5645.807106413317, 5645.807106413317, 5645.807106413317, 5645.807106413317, 8807.605987367038, 8807.605987367038], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 13, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(18708.048610863716, 0.0), … (4) … , (18708.048610863716, 0.0)])]), [29692.671074503003, 10025.935716776967], [[[18708.048610863716, 0.0], [26292.921098828716, 22649.172216637835], [28815.672878080768, 25630.90794731836], [32118.063536997226, 25838.234235751057], [40564.17237212511, 0.0], [18708.048610863716, 0.0]]], 2.0, 3.4089674709385103e8, 6.272500146526859e11, 15997.228292747279, 4.233716017246853e19, [71.48506477943155, 158.2814367604344, 133.8258510202125, 104.5093946598428, 71.8982527800788], [-10984.622463639287, -10423.264807987005, -9346.238611876932, -8269.21241576686, -7192.186219656787, -6115.160023546715, -5038.133827436642, -3961.1076313265694, -3399.7499756742873, -2138.3740860482612  …  2996.8710595461866, 6050.302702100963, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -3109.992225563366, -56.560583008589674, 2996.8710595461866, -56.560583008589674, 2996.8710595461866], [-10025.935716776967, -8349.666769847336, -5133.5610493842405, -1917.4553289211447, 1298.6503915419507, 4514.756112005046, 7730.861832468142, 10946.967552931237, 12623.236499860868, 14114.10436520113  …  4486.231424657013, 4486.231424657013, 7672.331471773916, 7672.331471773916, 7672.331471773916, 10858.431518890819, 10858.431518890819, 10858.431518890819, 14044.531566007721, 14044.531566007721], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 14, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(69731.37176368426, 60788.69094406669), … (6) … , (69731.37176368426, 60788.69094406669)])]), [71288.1799442355, 45938.89146732441], [[[69731.37176368426, 60788.69094406669], [54356.76582028708, 48907.36542038052], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [69731.37176368426, 60788.69094406669]]], 2.0, 6.692501784369707e8, 1.2314203283240261e12, 17364.42744277742, 1.3688265153115973e20, [120.47160929676905, 123.37635799052929, 126.97662415698366, 124.39463803085083, 115.79395304455706, 133.56558891941265, 155.4212285608975], [-1556.8081805512338, -158.04244405797817, -2754.0620749699347, -5350.081705881891, -7946.101336793848, -10542.120967705803, -13138.14059861776, -15734.160229529716, -18330.179860441673, -16931.414123948416  …  7917.446636196438, 11214.74575150765, 14512.044866818862, -5271.7498250484105, -1974.4507097371986, 1322.8484055740134, 4620.147520885225, 7917.446636196438, -1974.4507097371986, 1322.8484055740134], [14849.799476742286, 15930.750176086378, 13924.57490146147, 11918.399626836563, 9912.224352211653, 7906.049077586746, 5899.873802961838, 3893.6985283369295, 1887.5232537120216, 2968.473953056113  …  6620.7549527512165, 6620.7549527512165, 6620.7549527512165, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 9851.393738263567, 13082.032523775917, 13082.032523775917], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 15, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 67567.56940096187), … (4) … , (100000.0, 67567.56940096187)])]), [94831.18867433425, 73451.81728375112], [[[100000.0, 67567.56940096187], [92828.50028120863, 69199.29082454351], [88481.68077666825, 75498.45471114098], [93611.10211788886, 80173.12477744419], [100000.0, 73544.51422997881], [100000.0, 67567.56940096187]]], 2.0, 8.539646907159662e7, 1.571295030917378e11, 7832.048497452011, 2.4418149897730017e18, [77.18179333098198, 137.42632491867568, 97.73618770158065, 91.60073242094602, 136.0549616278157], [5168.8113256657525, 6892.5235894391235, 4237.792527854595, 1583.0614662700682, -1071.6695953144595, -3726.400656898987, -2002.6883931256161, -998.6668209792022, -2587.3824831875054, -4176.098145395808  …  -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -1920.8125055666267, 740.1159335663784, 3401.0443726993835, -4581.740944699632, -1920.8125055666267, 740.1159335663784, -1920.8125055666267], [-5884.247882789248, -6276.441752855163, -5672.414461926796, -5068.3871709984305, -4464.359880070064, -3860.3325891416976, -4252.526459207613, -5707.497483889769, -3405.220999899324, -1102.944515908879  …  -4116.480929822879, -4116.480929822879, -4116.480929822879, -1093.1404396396877, -1093.1404396396877, -1093.1404396396877, 1930.2000505435037, 1930.2000505435037, 1930.2000505435037, 4953.540540726695], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 16, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(53932.759746554315, 0.0), … (6) … , (53932.759746554315, 0.0)])]), [58425.22292764254, 20191.465410638164], [[[53932.759746554315, 0.0], [44832.48232119377, 28580.86301404227], [55044.67824938872, 39801.13999752771], [68140.10877192479, 31408.046754035153], [71716.46244051738, 22675.03020194567], [68671.97979189342, 11357.646577717289], [60505.550527273575, 0.0], [53932.759746554315, 0.0]]], 2.0, 6.839042436264261e8, 1.258383808272624e12, 20685.200044051562, 1.565246038690451e20, [107.66171295242654, 120.03122187956625, 99.65053593600663, 144.92665928164806, 142.67320840749525, 159.33959078788573, 125.71707075497156], [-4492.463181088227, -3956.1290668852325, -4973.423632261886, -5990.71819763854, -7008.012763015193, -8025.307328391847, -9042.6018937685, -10059.896459145153, -11077.191024521808, -12094.48558989846  …  -8489.481337283707, -5153.989021085012, -1818.496704886315, 1516.9956113123812, 4852.487927511078, 8187.980243709774, -5153.989021085012, -1818.496704886315, 1516.9956113123812, -1818.496704886315], [-20191.465410638164, -21875.90766786816, -18680.932915017933, -15485.958162167708, -12290.983409317481, -9096.008656467255, -5901.033903617029, -2706.0591507668028, 488.91560208342327, 3683.890354933649  …  11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 11248.161071814999, 14545.034352869088, 14545.034352869088, 14545.034352869088, 17841.907633923176], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 17, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(17908.287958045414, 100000.0), … (6) … , (17908.287958045414, 100000.0)])]), [33613.80575505652, 83906.94777996265], [[[17908.287958045414, 100000.0], [18547.875224604184, 80169.53795708867], [26710.280846516955, 68796.48411425078], [39281.31096026047, 60492.26994952193], [39877.20185576307, 60461.079867567416], [50658.43944984143, 89053.30431356005], [43884.94062412003, 100000.0], [17908.287958045414, 100000.0]]], 2.0, 9.177048547118509e8, 1.6885769326698057e12, 24268.06240472919, 2.6469913646102482e20, [88.15269295890626, 146.1804201768482, 159.11497680870042, 149.54814772225257, 107.66365980992538, 127.59207823601393, 121.74802428735329], [-15705.517797011107, -15648.532053579796, -15543.408897640567, -15438.285741701337, -15333.162585762107, -15228.039429822877, -15122.916273883648, -15065.930530452337, -14035.195499735168, -12001.550312909023  …  8784.729500515592, 12030.798121167067, -13937.750844044738, -10691.682223393262, -7445.613602741787, -4199.544982090311, -953.4763614388353, 2292.5922592126403, 5538.660879864116, 8784.729500515592], [16093.052220037353, 14326.204002151735, 11066.850880723716, 7807.4977592956975, 4548.144637867679, 1288.7915164396604, -1970.561604988358, -3737.409822873975, -5173.5801806144755, -8007.151223066773  …  11052.250155570997, 11052.250155570997, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984, 14325.285267070984], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 18, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 28356.11795020577), … (3) … , (100000.0, 28356.11795020577)])]), [94387.04219792974, 42856.57222855929], [[[100000.0, 28356.11795020577], [86975.8743229851, 39383.51055632255], [87567.99176402073, 51297.92451684605], [100000.0, 53790.42066586246], [100000.0, 28356.11795020577]]], 2.0, 2.3895189855925107e8, 4.39671493349022e11, 15548.905735338496, 2.075893184387935e19, [49.74575563575111, 127.4091240568404, 104.1820578473644, 78.66306246004409], [5612.957802070261, 6962.090324090315, 4341.691870581147, 1721.2934170719775, -899.1050364371913, -3519.5034899463603, -6139.901943455529, -8760.300396964698, -7411.167874944644, -7323.422442346604  …  -2480.536998284219, 682.3269254098365, 3845.190849103892, -5643.400921978275, -2480.536998284219, 682.3269254098365, 3845.190849103892, -2480.536998284219, 682.3269254098365, 3845.190849103892], [-14500.454278353522, -15642.750841435021, -13424.086552721725, -11205.422264008428, -8986.757975295131, -6768.093686581835, -4549.429397868538, -2330.7651091552416, -3473.0616722367413, -1707.4737420961371  …  2909.290395844241, 2909.290395844241, 2909.290395844241, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 6037.6859400905205, 9166.081484336799, 9166.081484336799, 9166.081484336799], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 19, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(100000.0, 53790.42066586246), … (4) … , (100000.0, 53790.42066586246)])]), [93357.10571133424, 59553.21688748059], [[[100000.0, 53790.42066586246], [87567.99176402073, 51297.92451684605], [83959.48287311713, 55088.65064944781], [92828.50028120863, 69199.29082454351], [100000.0, 67567.56940096187], [100000.0, 53790.42066586246]]], 2.0, 1.9638967050932837e8, 3.613569937371642e11, 10409.509630176246, 1.272724957257055e19, [101.3369375399559, 122.25235323322298, 104.25995041206583, 109.33255214573725, 102.81820666901802], [6642.8942886657605, 8376.168805988078, 5196.457351863297, 2016.7458977385152, -1162.9655563862664, -4342.677010511048, -7522.388464635829, -5789.113947313512, -7593.368392765311, -9397.62283821711  …  -4503.610080013208, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392, -1377.3642747756749, 1748.8815304618583, 4875.127335699392], [-5762.796221618126, -5415.291620783709, -6052.792690920758, -6690.293761057807, -7327.794831194857, -7965.295901331906, -8602.796971468955, -8255.29237063454, -6359.929304333658, -4464.566238032778  …  2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 2131.9740233906646, 5005.14050374361, 5005.14050374361, 5005.14050374361, 7878.3069840965545, 7878.3069840965545, 7878.3069840965545], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 20, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0), Floe{Float64}(GeoInterface.Wrappers.Polygon{false, false}([GeoInterface.Wrappers.LinearRing([(40564.17237212511, 0.0), … (3) … , (40564.17237212511, 0.0)])]), [42958.001115460844, 13654.895707898613], [[[40564.17237212511, 0.0], [32118.063536997226, 25838.234235751057], [44832.48232119377, 28580.86301404227], [53932.759746554315, 0.0], [40564.17237212511, 0.0]]], 2.0, 3.6688421815656054e8, 6.750669614080714e11, 17518.604505049727, 5.57569422307862e19, [108.1017472199212, 84.07101639433667, 95.48894933816871, 72.33828704757347], [-2393.8287433357327, -1844.5739782078895, -2905.0871299171754, -3965.600281626461, -5026.113433335747, -6086.626585045033, -7147.139736754319, -8207.652888463605, -9268.16604017289, -10328.679191882176  …  67.41052631492676, 3113.937576918985, -9072.170625497249, -6025.64357489319, -2979.1165242891316, 67.41052631492676, 3113.937576918985, -6025.64357489319, -2979.1165242891316, 67.41052631492676], [-13654.895707898613, -15335.169257944664, -12090.860220628758, -8846.551183312851, -5602.242145996945, -2357.933108681038, 886.3759286348686, 4130.684965950775, 7374.994003266682, 10619.303040582588  …  6896.868076149904, 6896.868076149904, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 10027.534214663596, 13158.200353177286, 13158.200353177286, 13158.200353177286], 0.0, 0.0, 0.0, 0.0, Subzero.Status(Subzero.active, Int64[]), 21, 0, Int64[], Int64[], 0.0, 0.0, 0.0, 0.0, 0.0, [0.0 0.0], 0.0, Matrix{Float64}(undef, 0, 7), 0, [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], [0.0 0.0; 0.0 0.0], 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)]), Constants{Float64}(1027.0, 1.2, 0.003, 0.001, 0.00125, 0.0, 0.0, 293000.0, 2.14, 0.3, 0.2, 3.595334669716184e7), Random.Xoshiro(0xfff0241072ddab67, 0xc53bc12f4c3f0b4e, 0x56d451780b2dd4ba, 0x50a4aa153d208dd8, 0x3649a58b3b63d5db), true, "sim", 20, 2500, FloeSettings{Float64, MonteCarloPointsGenerator{Float64}, DecayAreaScaledCalculator{Float64}}(920.0, 1.0e6, 0.1, 10.0, 0.05, 1.0e-5, MonteCarloPointsGenerator{Float64}(1000, 10, 0.1), DecayAreaScaledCalculator{Float64}(0.2, 0.0)), CouplingSettings(true, 10, 1, false), CollisionSettings{Float64}(true, 0.55, 0.75), FractureSettings{NoFracture}(false, NoFracture(), 0, false, 3), SimplificationSettings{Float64}(true, 30, 100.0, 20), RidgeRaftSettings{Float64}(false, 0, 0.95, 0.95, 0.01, 0.2, 5.0, 1.25, 0.25, 0.25, 1.0), WeldSettings{Float64}(false, Int64[], Int64[], Int64[], 1.0e6, 2.0e9, 150.0), OutputWriters{StructArrays.StructVector{InitialStateOutputWriter, @NamedTuple{filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{FloeOutputWriter, @NamedTuple{Δtout::Vector{Int64}, outputs::Vector{Vector{Symbol}}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}, StructArrays.StructVector{GridOutputWriter, @NamedTuple{outputs::Vector{Vector{Symbol}}, Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}, xg::Vector{Vector{FT} where FT<:AbstractFloat}, yg::Vector{Vector{FT} where FT<:AbstractFloat}, data::Vector{Array{FT, 3} where FT<:AbstractFloat}, average::Vector{Bool}}, Int64}, StructArrays.StructVector{CheckpointOutputWriter, @NamedTuple{Δtout::Vector{Int64}, filepath::Vector{String}, overwrite::Vector{Bool}}, Int64}}(InitialStateOutputWriter[InitialStateOutputWriter("restart_sim/run1/initial_state.jld2", true)], FloeOutputWriter[FloeOutputWriter(50, [:poly, :centroid, :coords, :height, :area, :mass, :rmax, :moment, :angles, :x_subfloe_points  …  :stress_accum, :stress_instant, :strain, :damage, :p_dxdt, :p_dydt, :p_dudt, :p_dvdt, :p_dξdt, :p_dαdt], "restart_sim/run1/floes.jld2", true)], GridOutputWriter[], CheckpointOutputWriter[CheckpointOutputWriter(250, "restart_sim/run1/checkpoint.jld2", true)]))

Run the first part of the simulation

run!(simulation)

Run second part of the simulation

new_initwriter = InitialStateOutputWriter(initwriter; dir = dirs[2])
new_floewriter = FloeOutputWriter(floewriter; dir = dirs[2])
writers = OutputWriters(new_initwriter, new_floewriter)  # didn't save checkpoint since not restarting
Subzero.restart!(
    dirs[1] * "/initial_state.jld2",
    dirs[1] * "/checkpoint.jld2",
    nΔt, writers; start_tstep = nΔt)

Plot all simulation parts

for i in 1:2
    plot_sim(
        joinpath(dirs[i], "floes.jld2"),
        joinpath(dirs[i], "initial_state.jld2"),
        Δt,
        joinpath(dirs[i], "restart_floes.mp4"),
    )
end

First Part of Simulation

Second Part of Simulation

This page was generated using Literate.jl.