Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.naic.edu/~tghosh/a2186/may09/abber.pro Дата изменения: Wed May 6 22:36:25 2009 Дата индексирования: Mon Apr 11 01:36:19 2016 Кодировка: Поисковые слова: http astrokuban.info astrokuban

;
pro abber,c1,c10, abr_ra,abr_dec
;
readcol,'pos_diff.data',c1,c2,c3,c4,c5,c6,c7,c8,c9,c10
;
dr=3.141593/180.0
print,dr
;
epsi=23.45*dr
alfa=300.0*dr
delta=33.0*dr
kappa=20.5
xlam=c10*dr
;
; calculating using eqn. 5-30,31 , page 501, astrphysical formulae
;
abr_ra=(1./cos(delta))*(-1.0*kappa*sin(xlam)*sin(alfa)-kappa*cos(xlam)*cos(epsi) $
*cos(alfa))
;
abr_dec=-1.0*kappa*sin(xlam)*cos(alfa)*sin(delta) + kappa*cos(xlam)*(cos(epsi) $
*sin(alfa)*sin(delta) - sin(epsi)*cos(delta))

plot,c1-54674,abr_ra,psym=2,xtitle='MJD - 54674',ytitle='Delta_RA_aberration'
set_plot,'ps'
device,file='aberr_RA.ps',ysize=12,set_font='Times-roman'
plot,c1-54674,abr_ra,psym=2,xtitle='MJD - 54674',ytitle='Delta_RA_aberration'
device,/close
set_plot,'x'
;
print,'enter a number to continue'
read,an
plot,c1-54674,abr_dec,psym=4,xtitle='MJD - 54674',ytitle='Delta_Dec._aberration'
set_plot,'ps'
device,file='aberr_Dec.ps',ysize=12,set_font='Times-roman'
plot,c1-54674,abr_dec,psym=4,xtitle='MJD - 54674',ytitle='Delta_Dec._aberration'
device,/close
set_plot,'x'
end