59 const PetscInt PRES = 1;  variable
824   for (c = 0; c < Nc; ++c) f1[c * dim + c] -= u[uOff[PRES]];  in f1_conduct_v()
1225     PetscCall(PetscDSSetResidual(ds, PRES, f0_q, NULL));  in SetupProblem()
1229     PetscCall(PetscDSSetJacobian(ds, VEL, PRES, NULL, NULL, g2_vp, NULL));  in SetupProblem()
1230     PetscCall(PetscDSSetJacobian(ds, PRES, VEL, NULL, g1_qu, NULL, NULL));  in SetupProblem()
1240       exactFuncs[PRES]   = quadratic_p;  in SetupProblem()
1243       exactFuncs_t[PRES] = NULL;  in SetupProblem()
1253       exactFuncs[PRES]   = cubic_p;  in SetupProblem()
1256       exactFuncs_t[PRES] = NULL;  in SetupProblem()
1266       exactFuncs[PRES]   = cubic_trig_p;  in SetupProblem()
1269       exactFuncs_t[PRES] = NULL;  in SetupProblem()
1278       exactFuncs[PRES]   = taylor_green_p;  in SetupProblem()
1281       exactFuncs_t[PRES] = taylor_green_p_t;  in SetupProblem()
1292     PetscCall(PetscDSSetResidual(ds, PRES, f0_conduct_q, NULL));  in SetupProblem()
1296     PetscCall(PetscDSSetJacobian(ds, VEL, PRES, NULL, NULL, g2_conduct_vp, NULL));  in SetupProblem()
1298     PetscCall(PetscDSSetJacobian(ds, PRES, VEL, g0_conduct_qu, g1_conduct_qu, NULL, NULL));  in SetupProblem()
1299     PetscCall(PetscDSSetJacobian(ds, PRES, TEMP, g0_conduct_qT, g1_conduct_qT, NULL, NULL));  in SetupProblem()
1306 …PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, PRES, 0, 1, f0_conduct_quadratic_q, 0, NULL));  in SetupProblem()
1310       exactFuncs[PRES]   = quadratic_p;  in SetupProblem()
1313       exactFuncs_t[PRES] = NULL;  in SetupProblem()
1321       PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, PRES, 0, 1, f0_conduct_pipe_q, 0, NULL));  in SetupProblem()
1325       exactFuncs[PRES]   = pipe_p;  in SetupProblem()
1328       exactFuncs_t[PRES] = pipe_p_t;  in SetupProblem()
1352 …PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, PRES, 0, 1, f0_conduct_pipe_wiggly_q, 0, NULL…  in SetupProblem()
1356       exactFuncs[PRES]   = pipe_wiggly_p;  in SetupProblem()
1359       exactFuncs_t[PRES] = pipe_wiggly_p_t;  in SetupProblem()
1412   PetscCall(PetscDSSetExactSolution(ds, PRES, exactFuncs[PRES], ctx));  in SetupProblem()
1415   PetscCall(PetscDSSetExactSolutionTimeDerivative(ds, PRES, exactFuncs_t[PRES], ctx));  in SetupProblem()
1486   PetscCall(DMSetField(dm, PRES, NULL, (PetscObject)fe[PRES]));  in SetupDiscretization()
1496   PetscCall(PetscFEDestroy(&fe[PRES]));  in SetupDiscretization()
1503     PetscCall(DMGetField(dm, PRES, NULL, &pressure));  in SetupDiscretization()
1517 …PRES, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_WRONG, "Nullspace must be for pressure field…  in CreatePressureNullSpace()