smurf: Allow calcqu to create maps of 2nd harmonic

New parameter HARMONIC allows the harmonic of the half-wave plate
rotation to be specified. This should always be 4 if *real* Q and U
values are required, but can be set to something else (e.g. 2) to
investigate instrumental effects.