Skip to content

HorzMesh structure does not support non-periodic meshes #2

@andrewdnolan

Description

@andrewdnolan

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.

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions