Details of Rotang

Rotang can be called either through the IDL command line or through the afla_ang GUI. The full cummand line syntax is:

rotang, Date, RA, Dec, Time, accS, ast=ast, ha=ha, lst=lst, jd=jd, az=az, za=za, sep=sep, rotang=rotang, alfabeams=alfabeams

date, RA, Dec and Time are covered in the basic rotang documentation; the fifth input parameter is accS, which is the accuracy to which the requested time is required. This defaults to half a second: if this requirement is not met then the calculation will be iterated until it is. It is not normally necessary to modify this parameter.

ast, ha and lst are used to specify which of AST, LST or HA is wanted, as described in the basic documentation. The precedence sequence for this is that if both LST takes precedence over HA but AST takes precedence over LST. If none are specified, the program defaults to HA. If these parameters are specified with a variable name, e.g. ha=haHr then the value of the variable on return will be that of the HA, LST or AST of the requested observation (in hours). It is thus possible to specify, for instance:

IDL> haHr = 1
IDL> lstHr = 0
IDL> astHr = 0
IDL> rotang, ... , ha=haHr, lst=lstHr, ast=astHr

In which case, the variables haHr, lstHr and astHr will, upon return, contain the requested HA and the corresponding LST and AST in hours.

jd, az, za, sep and rotang are return-parameters, returning the Julian Date, Azimuth and Zenith Angle of the observation and the mean declination separation of the ALFA beams and optimum rotation angle of the array respectively. alfabeams is a structure which contains the elements {angle, sep, dev, bnum, azoff, zaoff, az, za, ra, dec, raoff, decoff}. alfabeams.angle and alfabeams.sep should be identical to rotang and sep, and alfabeams.dev contains the error on the mean declination separation of the ALFA beams. The other elements are 7-element arrays containing the beam number (0-6) and the corresponding azimuth offset, zenith angle offset, azimuth, zenith angle, right ascension, declination, right ascension offset and declination offset. All offsets are with respect to the central beam (beam 0).

The algorithm works quite simply. Firstly the Julian Date of the requested observations is calculated as the next time after midnight AST on the requested date that meets the specified time criterion. This is used to calculate the Az and ZA of the source from its RA and Dec. From these, the parallactic angle of the source is calculated and used to make an inital estimate of the source position. This should ensure that the rotation of the source is correct to within 4 degrees, thus the general orientation of the beams will be correct. The minimising algorithm then starts from this initial guess and finds the slope of the error on the mean vs angle, which falls almost linearly to close to zero at its minimum. By moving down this slope in one degree steps until the slop reverses, a position accurate to one degree is found. The algorithm then double-checks that the beams are in the correct places, not rotated to a minimum different from the standard at meridian crossing, and that ALFA is within its rotation limits. If ALFA is outside its rotation limits, then hte angle is increased or decreased by 180 degrees to bring it inside. A second minimization then takes place using the same algorithm, this time with 0.1 degree steps. Finally, a check is made to see whether the source is at a Zenith Angle of less than two or greater than twenty degrees, and if it is a warning is printed.

Rotang consists of the following files in addition to the standard aoidl installation:

Original page written by Robert Minchin, 22nd November 2004. Last modified 22nd November 2004 by Robert Minchin.