0001 function df = calcRADECJ2000(df)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 disp('Getting Equatorial co-ordinates in J2000');
0015
0016
0017 display('Calculating RA/DEC for J2000');
0018
0019
0020 long=df.antenna0.tracker.siteActual(2,1);
0021 lat =df.antenna0.tracker.siteActual(2,2);
0022
0023 az = df.antenna0.servo.apparent(:,1);
0024 el = df.antenna0.servo.apparent(:,2);
0025 jd=mjd2jd(df.antenna0.receiver.utc);
0026 [equa] = horiz_coo([pi/180*(az) pi/180*(el)],jd,[pi/180*(long) ...
0027 pi/180*(lat)],'e');
0028 df.antenna0.servo.equa=equa;
0029
0030
0031 yearoff = (mean(df.array.frame.utc) - date2mjd(2000, 01, ...
0032 01,0,0,0))/365.25;
0033 yearoff = round( (2000+yearoff)*10)/10;
0034 txt = sprintf('[equa2000] = coco(df.antenna0.servo.equa, ''j%4.1f'', ''j2000.0'', ''r'', ''r'');', yearoff);
0035 eval(txt);
0036 df.antenna0.servo.equa=equa2000;
0037
0038
0039
0040
0041 df.antenna0.servo.epoch = 'J2000'
0042
0043 clear az;
0044 clear el;
0045 clear equa;
0046 clear equa2000;
0047 end