quickLookAzScan

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