Just discovered that WCS_ROTATE is still not working properly for 180 degree rotations; specifically, it fails at longitude 180 degrees:
HIDL> wcs_rotate, 180d0+[1,0,-1], 7d0+[0,0,0], phi, theta, [180d0,0d0], theta0=0d0
HIDL> print, phi
1.0000000 180.00000 -1.0000000
HIDL> print, theta
7.0000000 7.0000000 7.0000000
Just discovered that WCS_ROTATE is still not working properly for 180 degree rotations; specifically, it fails at longitude 180 degrees:
HIDL> wcs_rotate, 180d0+[1,0,-1], 7d0+[0,0,0], phi, theta, [180d0,0d0], theta0=0d0
HIDL> print, phi
1.0000000 180.00000 -1.0000000
HIDL> print, theta
7.0000000 7.0000000 7.0000000