Added SQ(Field3D) function

Just calls f*f. Useful shorthand sometimes, and matches
existing SQ(BoutReal) and SQ(Field2D) functions