Limitations Include:

• When edgesOnVertex contains zeros (i.e. the vertex is along a culled boundary) an BoundsError is raised in the signIndexField! method, which tries to index arrays at 0. For example:

  ERROR: BoundsError: attempt to access 2×555 Matrix{Int32} at index [1, 0]
  Stacktrace:
   [1] getindex
     @ ./essentials.jl:14 [inlined]
   [2] getindex
     @ ./abstractarray.jl:1294 [inlined]
   [3] signIndexField!(dualMesh::MOKA.MPASMesh.DualCells{Int64, Vector{Float64}, Matrix{Int32}}, edges::MOKA.MPASMesh.Edges{Int64, Vector{Float64}, Vector{Int32}, Matrix{Float64}, Matrix{Int32}})
     @ MOKA.MPASMesh /global/u2/a/anolan/Moka.jl/src/infra/MPASMesh/HorzMesh.jl:323
   [4] ReadHorzMesh(meshPath::String; backend::KernelAbstractions.CPU)
     @ MOKA.MPASMesh /global/u2/a/anolan/Moka.jl/src/infra/MPASMesh/HorzMesh.jl:345
   [5] ReadHorzMesh(meshPath::String)
     @ MOKA.MPASMesh /global/u2/a/anolan/Moka.jl/src/infra/MPASMesh/HorzMesh.jl:334
   [6] top-level scope
     @ REPL[2]:1
  

• Horizontal mesh structures do not yet have *Mask arrays (i.e. edgeMask) arrays, meaning horizontal boundary conditions (i.e. no slip) will not be correctly applied.