| 1 | 
greg | 
1.1 | 
{ SCCSid "$SunId$ LBL" } | 
| 2 | 
  | 
  | 
 | 
| 3 | 
  | 
  | 
{ | 
| 4 | 
  | 
  | 
        Source distribution coordinates (degrees). | 
| 5 | 
  | 
  | 
 | 
| 6 | 
  | 
  | 
        Theta is measured from the negative z-axis, | 
| 7 | 
  | 
  | 
        phi is measured from the positive x-axis | 
| 8 | 
  | 
  | 
        according to the right-hand rule. | 
| 9 | 
  | 
  | 
 | 
| 10 | 
  | 
  | 
        srcB_vert and srcB_horiz are angles | 
| 11 | 
  | 
  | 
        used in type B photometry. | 
| 12 | 
  | 
  | 
 | 
| 13 | 
  | 
  | 
        A1      - optional multipier | 
| 14 | 
  | 
  | 
} | 
| 15 | 
  | 
  | 
 | 
| 16 | 
  | 
  | 
flatcorr(v) = corr(v) / Rdot;           { correction for flat sources } | 
| 17 | 
  | 
  | 
corr(v) = if(AC-.5, A1*v, v);           { multiplier correction } | 
| 18 | 
  | 
  | 
 | 
| 19 | 
  | 
  | 
src_theta = Acos(Dz) / DEGREE;                  { 0-180 } | 
| 20 | 
  | 
  | 
src_phi = norm_deg( atan2(Dy, Dx) / DEGREE );   { 0-360 } | 
| 21 | 
  | 
  | 
                                { bilateral symmetry } | 
| 22 | 
  | 
  | 
src_phi2 = tri( src_phi, 180 );         { 0-180 } | 
| 23 | 
  | 
  | 
                                { quadrilateral symmetry } | 
| 24 | 
  | 
  | 
src_phi4 = tri( src_phi, 90 );          { 0-90 } | 
| 25 | 
  | 
  | 
 | 
| 26 | 
  | 
  | 
norm_deg(d) : if( d, d, d+360 ); | 
| 27 | 
  | 
  | 
 | 
| 28 | 
  | 
  | 
                                { Type B photometry coordinates } | 
| 29 | 
  | 
  | 
srcB_vert = atan( -Dx/Dz ) / DEGREE; | 
| 30 | 
  | 
  | 
srcB_horiz = Asin(Dy) / DEGREE; | 
| 31 | 
  | 
  | 
                                { w/ bilateral symmetry } | 
| 32 | 
  | 
  | 
srcB_horiz2 = abs( srcB_horiz ); |