Skip to content

z-boundary backward FD stencil uses fields(2) instead of fields(3) #1194

Description

@sbryngelson

Bug

In src/common/m_finite_differences.fpp, the backward finite-difference stencil for the z-direction upper boundary references fields(2) (the y-direction field) instead of fields(3) (the z-direction field) for the z - 2 term.

File: src/common/m_finite_differences.fpp

Impact: Corrupts z-direction derivative values at the upper z-boundary cells. The stencil silently mixes y-direction data into z-direction derivatives. This affects any 3D simulation using the backward FD stencil at the z-boundary (e.g., viscous stress, gradient computations).

Severity: Critical — produces silently wrong results in 3D simulations near the z upper boundary.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working or doesn't seem right

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions