Vacuum/positivity guard for the FVS split-flux precompute. The FVS split routines evaluate sqrt(gam*p/rho); a non-positive density or pressure would produce NaN that then propagates silently through the reconstructed face fluxes and the residual. Detect it here and fail clearly instead. Mirrors the 1D FVS guard in spatial_discretization.f90 (block-guarded check + parallel_fatal). Operands are guarded with nested ifs because Fortran .and. does not short-circuit (rho is a divisor in p).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | q(neq2d) | |||
| real(kind=wp), | intent(in) | :: | gam | |||
| character(len=*), | intent(in) | :: | ctx |