176   const SWNode *uL = (const SWNode *)xL, *uR = (const SWNode *)xR;  in PhysicsRiemann_SW_Rusanov()  local
180   const SWNode *uR = &uRreal.swnode;  in PhysicsRiemann_SW_Rusanov()  local
194   if (uL->h < 0 || uR->h < 0) {  in PhysicsRiemann_SW_Rusanov()
211   ierr = SWFlux(phys, nn, uR, &fR.swnode);  in PhysicsRiemann_SW_Rusanov()
218   cR    = PetscSqrtReal(sw->gravity * uR->h); /* gravity wave speed */  in PhysicsRiemann_SW_Rusanov()
219 …(PetscAbsReal(Dot2Real(uL->uh, nn) / uL->h) + cL, PetscAbsReal(Dot2Real(uR->uh, nn) / uR->h) + cR);  in PhysicsRiemann_SW_Rusanov()
276   const SWNode *uL = (const SWNode *)xL, *uR = (const SWNode *)xR;  in PhysicsRiemann_SW_HLL()  local
280   const SWNode *uR = &uRreal.swnode;  in PhysicsRiemann_SW_HLL()  local
293   if (uL->h <= 0 || uR->h <= 0) {  in PhysicsRiemann_SW_HLL()
308   ierr = SWFlux(phys, nn, uR, &fR.swnode);  in PhysicsRiemann_SW_HLL()
316   aR = PetscSqrtReal(sw->gravity * uR->h);  in PhysicsRiemann_SW_HLL()
320   u_R = Dot2Real(uR->uh, nn) / uR->h;  in PhysicsRiemann_SW_HLL()
801     const EulerNode *uL = (const EulerNode *)xL, *uR = (const EulerNode *)xR;  in PhysicsRiemann_Euler_Godunov()  local
809     ierr = EulerFlux(phys, nn, uR, &fR.eulernode);  in PhysicsRiemann_Euler_Godunov()
821     ierr = SpeedOfSound_PG(gamma, uR, &cR);  in PhysicsRiemann_Euler_Godunov()
828     velR  = DotDIMReal(uR->ru, nn) / uR->r;  in PhysicsRiemann_Euler_Godunov()