aeattcor -- correct Euler angles in attitude file for Suzaku
aeattcor ae123456789.att ae123456789.attcor \ ae123456789.orb ae123456789.hk
This tool calculates the effects of the "thermal wobbling" caused by thermal distortions of the satellite bodies, and correct the Euler angles in attitude file for Suzaku.
The attitude correction function is described in 'Empirical Correction of Thermal Wobbling', JX-ISAS-SUZAKU-MEMO-2007-04, as
XOFFSET = -38.5*(b/60) - 9.25 + 0.4*T_86 [pixel] YOFFSET = 10.0 + ( 7.0/1000)*T_DY [pixel] for T_DY < 0 [s] = 10.0 + ( 5.0/2000)*T_DY [pixel] for 0 < T_DY < 2000 [s] = 45.0 - (15.0/1000)*T_DY [pixel] for 2000 < T_DY [s] b: Ecliptic latitude [deg] T_86: Temperature difference, HK_XIS_RAD8_T1_CAL - HK_XIS_RAD6_T1_CAL T_DY: Time after night-day transision [s], defined as T_DY = T_DY_NT, for day-time T_DY = -TN_DY_NT, for night-time
Following columns are also added.
EULER_OLD: contains original Z-Y-Z euler angles [deg] XOFFSET: calculated XOFFSET by the above formula [pixel] YOFFSET: calculated YOFFSET by the above formula [pixel] T_86: HK_XIS_RAD8_T1_CAL - HK_XIS_RAD6_T1_CAL [K] T_DY: T_DY_NT for day-time, -TN_DY_NT for night-time [s] in EHK FLAG_ATTCOR: attitude correction was conducted (=1) or not (=0)
If the EULER_OLD column already exists in the input attitude file, the aeattcor task reads the Euler angles from this column. Therefore, applying the aeattcor to the corrected products gives the same output with the input attitude file, except for history.
% aeattcor ae400006010.att.gz ae400006010.attcor \ ae400006010.orb.gz ae400006010.hk aeattcor version 2007-04-08 Written by Y.ISHISAKI (TMU) Built on ANL HEADAS converter 1.73 for ANL version 1.73 aeattcor: *** show parameter *** INFILE 'ae400006010.att.gz' OUTFILE 'ae400006010.attcor' ORBIT 'ae400006010.orb.gz' HKFILE 'ae400006010.hk' HK_TIME_MARGIN 60.0 (s) CLOBBER YES aste_orbit: reading 'ae400006010.orb.gz[EXTNAME=PAR_ORBIT,EXTVER=0]' ... nrow=28805, nkp=25921, tstart=183340800.0, tstop=184896000.0 aste_gethk-2.5: found 'HK_XIS_RAD6_T1_CAL' at hdu=3, col=134 \ in 'ae400006010.hk' aste_gethk-2.5: t=184176126.248 < TSTART=184209153.963 \ for 'HK_XIS_RAD6_T1_CAL' in 'ae400006010.hk' aeattcor: INFO: TSTART=184209154.0 for ae400006010.hk aste_gethk-2.5: found 'HK_XIS_RAD6_T1_CAL' at hdu=3, col=134 \ in 'ae400006010.hk' aste_gethk-2.5: found 'HK_XIS_RAD8_T1_CAL' at hdu=3, col=138 \ in 'ae400006010.hk' aste_gethk-2.5: t=184289964.983 > TSTOP=184289962.826 \ for 'HK_XIS_RAD6_T1_CAL' aeattcor: INFO: TSTOP=184289962.8 for ae400006010.hk aste_gethk-2.5: found 'HK_XIS_RAD6_T1_CAL' at hdu=3, col=134 \ in 'ae400006010.hk' NUM_CORR = 150819 / number of corrected Euler angles AVG_LAMB = 350.4655 / average ecliptic longitude (deg) AVG_BETA = +52.6553 / average ecliptic latitude (deg) AVG_XOFF = -42.4934 / average DETX offset (pixel) AVG_YOFF = 6.1975 / average DETY offset (pixel) SGM_XOFF = 5.1810 / 1 sigma DETX offset standard deviation (pixel) SGM_YOFF = 18.5209 / 1 sigma DETY offset standard deviation (pixel) aeattcor: INFO: created 'ae400006010.attcor'
none
This program was developed primarily in the ASTRO-E ANL environment by Yoshitaka Ishisaki (TMU).
Questions and comments should be addressed to astroehelp@olegacy.gsfc.nasa.gov.