Fix initatlization bug

examples/many_floes.jl

@@ -27,7 +27,7 @@ ocean = Ocean(grid, -0.2, 0.0, -1.0)
 atmos = Atmos(grid, 0.0, 0.0, -3.0)
 
 # Domain creation - boundaries and topography
-nboundary = (North, grid)
+nboundary = OpenBoundary(North, grid)
 sboundary = OpenBoundary(South, grid)
 eboundary = OpenBoundary(East, grid)
 wboundary = OpenBoundary(West, grid)
@@ -48,11 +48,11 @@ model = Model(grid, ocean, atmos, domain, floe_arr)
 modulus = 1.5e3*(mean(sqrt.(floe_arr.area)) + minimum(sqrt.(floe_arr.area)))
 consts = Constants(E = modulus)
 #consts = Constants(E = modulus, Cd_io = 0.0, Cd_ia = 0.0, Cd_ao = 0.0, f = 0.0, μ = 0.0) # collisions without friction 
-simulation = Simulation(model = model, consts = consts, Δt = Δt, nΔt = 4000, COLLISION = true, verbose = true)
+simulation = Simulation(; model, consts, Δt, nΔt = 4000, verbose = true)
 
 # Output setup
-initwriter = InitialStateOutputWriter(dir = "output/voronoi", filename = "initial_state.jld2", overwrite = true)
-floewriter = FloeOutputWriter(10, dir = "output/voronoi", filename = "f.jld2", overwrite = true)
+# initwriter = InitialStateOutputWriter(dir = "output/voronoi", filename = "initial_state.jld2", overwrite = true)
+# floewriter = FloeOutputWriter(10, dir = "output/voronoi", filename = "f.jld2", overwrite = true)
 
-# Run simulation
-run!(simulation, [initwriter, floewriter])
+# # Run simulation
+# run!(simulation, [initwriter, floewriter])

src/simulation_components/floe.jl

@@ -354,15 +354,15 @@ Output:
 """
 function initialize_floe_field(
  ::Type{FT},
- coords,
+ coords::V,
  domain,
  hmean,
  Δh,
  Δt;
  floe_settings = FloeSettings(min_floe_area = 0.0),
  rng = Xoshiro(),
  supress_warnings = false,
-) where {FT <: AbstractFloat}
+) where {FT <: AbstractFloat, V <: AbstractVector}
  floe_arr = StructArray{Floe{FT}}(undef, 0)
  floe_polys = [make_polygon(valid_polyvec!(c)) for c in coords]
  # Remove overlaps with topography