Adding Mrho to setRHS_PBC, instead of M.
Adding initJet() for class Helmholtz. Code ran fine, but the proposed changes didn't take the expected effect of refinement.