|
@inline function upwinded_divergence_flux_Uᶠᶜᶜ(i, j, k, grid, scheme::VectorInvariantCrossVerticalUpwinding, u, v) |
|
@inbounds û = u[i, j, k] |
|
δ_stencil = scheme.upwinding.divergence_stencil |
|
|
|
δᴿ = _biased_interpolate_xᶠᵃᵃ(i, j, k, grid, scheme, scheme.divergence_scheme, bias(û), flux_div_xyᶜᶜᶜ, δ_stencil, u, v) |
|
∂t_σ = _symmetric_interpolate_xᶠᵃᵃ(i, j, k, grid, scheme, cross_scheme, Az_Δr_∂t_σ) |
|
|
|
return û * (δᴿ + ∂t_σ) # For static grids, ∂t_σ == 0 |
|
end |
|
|
|
@inline function upwinded_divergence_flux_Vᶜᶠᶜ(i, j, k, grid, scheme::VectorInvariantCrossVerticalUpwinding, u, v) |
|
@inbounds v̂ = v[i, j, k] |
|
δ_stencil = scheme.upwinding.divergence_stencil |
|
|
|
δᴿ = _biased_interpolate_yᵃᶠᵃ(i, j, k, grid, scheme, scheme.divergence_scheme, bias(v̂), flux_div_xyᶜᶜᶜ, δ_stencil, u, v) |
|
∂t_σ = _symmetric_interpolate_yᵃᶠᵃ(i, j, k, grid, scheme, cross_scheme, Az_Δr_∂t_σ) |
|
|
|
return v̂ * (δᴿ + ∂t_σ) # For static grids, ∂t_σ == 0 |
|
end |
Oceananigans.jl/src/Advection/vector_invariant_cross_upwinding.jl
Lines 36 to 54 in 478f45b
cross_schemeis undefined in both methods. Probably should be a field ofschemeor something like that? It was introduced by @simone-silvestri in #3956.Detected with JETLS