radpt

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
  function radpt(start, stop, filename)


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function radpt(start, stop, filename, doabs)
0002 
0003 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0004 %
0005 %  function radpt(start, stop, filename)
0006 %
0007 %
0008 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0009 
0010 % wrapper around radio pointing schedule that doesn't read in unecessary
0011 % registers.
0012 
0013   regs={'array.frame.received'...
0014     'array.frame.utc double',...
0015     'array.frame.features',...
0016     'array.weather.utc double',...
0017     'array.weather.pressure double',...
0018     'array.weather.airTemperature double',...
0019     'array.weather.windSpeed[0] double',...
0020     'array.weather.windDirection[0] double',...
0021     'array.weather.status double',...
0022     'array.weather.relativeHumidity double',...
0023     'antenna0.thermal.ccHeaterCurrent double',...
0024     'antenna0.thermal.lsTemperatureSensors double',...
0025     'antenna0.thermal.ccTemperatureLoad double',...
0026     'antenna0.thermal.ccColdPlate double', ...
0027     'antenna0.thermal.utc double',...
0028     'antenna0.thermal.dlpTemperatureSensors double',...
0029     'antenna0.tracker.lst double',...
0030     'antenna0.tracker.lacking double',...
0031     'antenna0.tracker.equat_geoc double',...
0032     'antenna0.tracker.horiz_topo double',...
0033     'antenna0.tracker.horiz_mount double',...
0034     'antenna0.tracker.flexure double',...
0035     'antenna0.tracker.horiz_off double',...
0036     'antenna0.tracker.tilts double',...
0037     'antenna0.tracker.fixedCollimation double',...
0038     'antenna0.tracker.encoder_off double',...
0039     'antenna0.tracker.sky_xy_off double',...
0040     'antenna0.tracker.source string',...
0041     'antenna0.tracker.scan_off double',...
0042     'antenna0.tracker.refraction double',...
0043     'antenna0.tracker.ut1utc double',...
0044     'antenna0.tracker.eqneqx double',...
0045     'antenna0.tracker.utc double',...
0046     'antenna0.tracker.time_diff double',...
0047     'antenna0.tracker.offSource double',...
0048     'antenna0.tracker.siteActual double',...
0049     'antenna0.tracker.siteFiducial double',...
0050     'antenna0.servo.utc double',...
0051     'antenna0.servo.fast_az_pos double',...
0052     'antenna0.servo.fast_el_pos double',...
0053     'antenna0.servo.fast_az_err double',...
0054     'antenna0.servo.fast_el_err double',...
0055     'antenna0.servo.slow_az_pos double',...
0056     'antenna0.servo.actual_current_az1[0] double',...
0057     'antenna0.servo.actual_current_az2[0] double',...
0058     'antenna0.servo.actual_current_el1[0] double',...
0059     'antenna0.frame.utc double',...
0060     'antenna0.frame.received double',...
0061     'antenna0.receiver.flags double',...
0062     'antenna0.receiver.utc double',...
0063     'antenna0.receiver.data double',...
0064     'antenna0.receiver.diagnostics double',...
0065     'antenna0.receiver.drainCurrent double',...
0066     'antenna0.receiver.drainVoltage double',...
0067     'antenna0.receiver.gateVoltage double',...
0068 };
0069 
0070 if(nargin<4)
0071   doabs = 0;
0072 end
0073 
0074 
0075 d = read_arc(start, stop, regs);
0076 d = pipe_read([],[],d);
0077 
0078 % remove the load
0079 plotparams.save = 0;
0080 plotparams.plot = 1;
0081 plotparams.interactive = 1;
0082 parm = [];
0083 field = '';
0084 
0085 % cold load needs alpha wrapper
0086 % make sure the alpha database is uptodate
0087 isup2date = checkAlphaDatabase((d.array.frame.utc(1)));
0088 if(~isup2date)
0089   display('Your alpha database is not up to date');
0090   display('You will update the database');
0091   display('Don''t fudge it up');
0092   updateAlphaDatabase(last(d.array.frame.utc));
0093 end
0094 d = alphaWrapper(d, plotparams, parm, field);
0095 
0096 % remove the cold load
0097 d = loadWrapper(d, plotparams, parm, field);
0098 
0099 % flag for RFI
0100 % fucking RFI flagging is awful and flagging EVERYTHING in pointing
0101 % scheds
0102 %parm.rfi.short = [inf inf inf inf];
0103 %parm.rfi.long =  [ inf inf inf inf inf];
0104 %parm.rfi.smooth = [inf];
0105 %parm.rfi.rejection = [inf];
0106 
0107 
0108 %d = rfiWrapper(d, plotparams, parm, field);
0109 
0110 switch(doabs)
0111   case 1
0112     d.antenna0.receiver.data(:,[1 8]) = d.antenna0.receiver.data(:,[1 8])+200;
0113   
0114   case 2
0115     d.antenna0.receiver.data(:,[1 8]) = 300 - d.antenna0.receiver.data(:,[1 ...
0116       8]);
0117 end
0118 
0119 
0120 
0121 
0122 radioPointing(d, 'scan', filename);
0123 
0124 return;