BUG FIX: p%mesh should be p_theta%mesh.  This was only a problem when using the viscous fs while testing the divergence using cv.  This is tested but apparently not in parallel when you can't get away with using the wrong global node numbering.