lookAtScan

d = read_arc('03-nov-2009:01:09:01' , '03-nov-2009:01:34:55');

0001 %d = read_arc('03-nov-2009:01:09:01' , '03-nov-2009:01:34:55');
0002 
0003 %d = read_arc('09-nov-2009:21:01:54', '09-nov-2009:21:45:00');
0004 %d = read_arc('09-nov-2009:22:34:37', '09-nov-2009:23:09:00');
0005 
0006 %d = read_arc('10-nov-2009:01:49:02', '10-nov-2009:02:25:00');
0007 d = pipe_read('02-Mar-2010:09:18:24',      '02-Mar-2010:12:13:25');
0008 
0009 time = (d.antenna0.servo.utc - d.antenna0.servo.utc(1))*24*60;
0010 
0011 azpos = d.antenna0.servo.fast_az_pos;
0012 temp = d.antenna0.thermal.ccTemperatureLoad;
0013 temp(length(temp)-4:length(temp)) = [];
0014 time(1:5) = [];
0015 azpos(1:5) = [];
0016 
0017 
0018 % for temps, want to show from 30 to 1030
0019 
0020   subplot(2,2,1)
0021   plot(time(30:1030), azpos(30:1030));
0022 xlabel('Time (minutes)');
0023 ylabel('Degrees');
0024 title('Azimuth during scan');
0025 
0026 subplot(2,2,3)
0027   plot(time(30:1030), temp(30:1030));
0028 xlabel('Time (minutes)');
0029 ylabel('Temperature [K]');
0030 title('Cold Load during scan');
0031 
0032 subplot(2,2,2)
0033   plot(time(480:565), azpos(480:565));
0034 xlabel('Time (minutes)');
0035 ylabel('Degrees');
0036 title('Azimuth at turnaround');
0037 
0038 subplot(2,2,4)
0039   plot(time(480:565), temp(480:565));
0040 xlabel('Time (minutes)');
0041 ylabel('Temperature [K]');
0042 title('Cold Load at turnaround');
0043 
0044 
0045 
0046 dt = time;
0047 dt(1) = [];
0048 daz = diff(azpos);
0049 
0050 azErr = d.antenna0.servo.fast_az_err;
0051 % 60 to 510 is the first scan
0052 
0053 
0054 % 4 panel plot.  [az     zoom(dAz)   ]
0055 %         [dAz    zoom(azErr)]
0056 
0057 figure(1)
0058   subplot(2,2,1)
0059   plot(time, azpos);
0060 xlabel('Time (Minutes)');
0061 ylabel('Azimuth Position (Degrees)');
0062 title('Antenna Azimuth in 30 minute scan');
0063 
0064 
0065 subplot(2,2,3)
0066   plot(dt, daz*5);
0067 xlabel('Time (Minutes)');
0068 ylabel('Velocity (degrees/second)');
0069 title('Velocity in Azimuth');
0070 
0071 subplot(2,2,2);
0072 plot(60*time(55:505), daz(55:505)*5);
0073 xlabel('Time (Seconds)');
0074 ylabel('Velocity (degrees/second)');
0075 title('Velocity in Azimuth Over One Scan');
0076 
0077 
0078 subplot(2,2,4);
0079 plot(60*time(75:505), azErr(75:505)*60*60);
0080 xlabel('Time (Seconds)');
0081 ylabel('Azimuth Error (Arcseconds)');
0082 title('Azimuth Error over 1 Scan');
0083 
0084 print -djpeg -r250 az_scan.jpg
0085 print -dpng -r250 az_scan.png
0086