0001 function d = read_arcSouth_nojunk(start, finish, regs, arcdir, calfile)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014 if(~exist('regs'))
0015 regs=[];
0016 end
0017 if(~exist('arcdir'))
0018 arcdir=[];
0019 end
0020 if(~exist('calfile'))
0021 calfile=[];
0022 end
0023
0024 if(isempty(regs))
0025 regs={'array.frame.received'...
0026 'array.frame.utc double',...
0027 'array.frame.features',...
0028 'antenna0.roach1.LL double',...
0029 'antenna0.roach1.LLfreq double',...
0030 'antenna0.roach1.RR double',...
0031 'antenna0.roach1.RRfreq double',...
0032 'antenna0.roach1.U double',...
0033 'antenna0.roach1.Ufreq double',...
0034 'antenna0.roach1.Utime double',...
0035 'antenna0.roach1.Q double',...
0036 'antenna0.roach1.Qfreq double',...
0037 'antenna0.roach1.Qtime double',...
0038 'antenna0.roach1.load1 double',...
0039 'antenna0.roach1.load1freq double',...
0040 'antenna0.roach1.load2 double',...
0041 'antenna0.roach1.load2freq double',...
0042 'antenna0.roach1.utc double',...
0043 'antenna0.roach1.intCount double',...
0044 'antenna0.roach1.ntpSeconds double',...
0045 'antenna0.roach1.ntpUSeconds double',...
0046 'antenna0.roach2.LL double',...
0047 'antenna0.roach2.LLfreq double',...
0048 'antenna0.roach2.RR double',...
0049 'antenna0.roach2.RRfreq double',...
0050 'antenna0.roach2.U double',...
0051 'antenna0.roach2.Ufreq double',...
0052 'antenna0.roach2.Utime double',...
0053 'antenna0.roach2.Q double',...
0054 'antenna0.roach2.Qfreq double',...
0055 'antenna0.roach2.Qtime double',...
0056 'antenna0.roach2.load1 double',...
0057 'antenna0.roach2.load1freq double',...
0058 'antenna0.roach2.load2 double',...
0059 'antenna0.roach2.load2freq double',...
0060 'antenna0.roach2.utc double',...
0061 'antenna0.roach2.intCount double',...
0062 'antenna0.roach1.switchstatus double',...
0063 'antenna0.roach2.switchstatus double',...
0064 'antenna0.roach2.ntpSeconds double',...
0065 'antenna0.roach2.ntpUSeconds double',...
0066 'antenna0.thermal.ccTemperatureLoad double',...
0067 'antenna0.thermal.ccHeaterCurrent double',...
0068 'antenna0.thermal.lsTemperatureSensors double',...
0069 'antenna0.thermal.utc double',...
0070 'antenna0.roach1.version double',...
0071 'antenna0.servo.fast_az_pos double',...
0072 'antenna0.servo.fast_el_pos double',...
0073 'antenna0.servo.fast_az_err double',...
0074 'antenna0.servo.fast_el_err double',...
0075 'antenna0.servo.utc double',...
0076 'antenna0.servo.ntpSecond double',...
0077 'antenna0.servo.ntpUSecond double',...
0078 'antenna0.tracker.lst double',...
0079 'antenna0.tracker.lacking double',...
0080 'antenna0.tracker.equat_geoc double',...
0081 'antenna0.tracker.horiz_topo double',...
0082 'antenna0.tracker.horiz_mount double',...
0083 'antenna0.tracker.flexure double',...
0084 'antenna0.tracker.horiz_off double',...
0085 'antenna0.tracker.tilts double',...
0086 'antenna0.tracker.fixedCollimation double',...
0087 'antenna0.tracker.encoder_off double',...
0088 'antenna0.tracker.sky_xy_off double',...
0089 'antenna0.tracker.source string',...
0090 'antenna0.tracker.scan_off double',...
0091 'antenna0.tracker.refraction double',...
0092 'antenna0.tracker.ut1utc double',...
0093 'antenna0.tracker.eqneqx double',...
0094 'antenna0.tracker.utc double',...
0095 'antenna0.tracker.time_diff double',...
0096 'antenna0.tracker.offSource double',...
0097 'antenna0.tracker.siteActual double',...
0098 'antenna0.tracker.siteFiducial double',...
0099 'array.weather.airTemperature double',...
0100 'array.weather.pressure double',...
0101 'array.weather.relativeHumidity double',...
0102
0103 };
0104 end
0105
0106
0107 if(any(size(unique(regs))~=size(regs)))
0108 error('regs should be unique');
0109 end
0110
0111
0112
0113 [defaultDataDir defaultCalFile defaultReaderPath] = whichHostSa();
0114
0115 if(isempty(arcdir))
0116
0117 arcdir = defaultDataDir;
0118 end
0119 if(isempty(calfile))
0120 calfile = defaultCalFile;
0121 end
0122
0123 eval(sprintf('addpath %s', defaultReaderPath));
0124
0125
0126
0127 mjdstartorig = tstr2mjd(start);
0128 mjdstart = mjdstartorig - 5/60/60/24;
0129 mjdstoporig = tstr2mjd(finish);
0130 mjdstop = mjdstoporig + 5/60/60/24;
0131
0132
0133 d = cbassMatReadArcOpt(regs, mjd2string(mjdstart), mjd2string(mjdstop), arcdir, calfile);
0134
0135
0136
0137 d = massageDataSa(d);
0138
0139
0140 d = reshapeRegisters(d);
0141
0142
0143
0144
0145
0146
0147
0148
0149
0150
0151 return;
0152
0153
0154
0155
0156
0157
0158
0159
0160
0161
0162
0163
0164
0165
0166
0167
0168
0169
0170
0171
0172
0173
0174
0175
0176
0177
0178
0179
0180
0181
0182
0183
0184
0185
0186
0187
0188
0189
0190
0191
0192
0193
0194
0195
0196
0197
0198
0199
0200
0201
0202
0203
0204
0205
0206
0207
0208
0209
0210
0211
0212
0213
0214
0215
0216
0217
0218
0219
0220
0221
0222
0223
0224
0225
0226
0227
0228
0229
0230
0231
0232
0233
0234
0235
0236
0237
0238
0239
0240
0241
0242
0243
0244
0245 host = strfind(w, 'aslx5');
0246 if(~isempty(host))
0247 hostNum = 7;
0248 defaultDataDir = '/home/cbassuser/cbass_data/arc';
0249 defaultCalfile = '/home/cbassuser/cbass/gcpCbass/control/conf/cbass/cal';
0250 defaultReader = '/home/cbassuser/cbass/gcpCbass/matlab/common/';
0251 end
0252
0253
0254
0255
0256
0257
0258
0259
0260
0261
0262
0263
0264
0265
0266 function d = cutDesiredData(d, mjdstartorig, mjdstoporig)
0267
0268 if(~isfield(d, 'array'))
0269 return;
0270 end
0271
0272 if(~isfield(d, 'antenna0'))
0273 return;
0274 end
0275
0276 if(issubfield(d.antenna0, 'receiver', 'utc'))
0277 indFast = d.antenna0.receiver.utc>=mjdstartorig & ...
0278 d.antenna0.receiver.utc<mjdstoporig;
0279 noRx = 0;
0280 else
0281 indFast = ones(length(d.array.frame.utc)*100,1);
0282 indFast = logical(indFast);
0283 noRx = 1;
0284 end
0285
0286 if(isfield(d.antenna0, 'servo'))
0287 indMed = d.antenna0.servo.utc>=mjdstartorig & ...
0288 d.antenna0.servo.utc<mjdstoporig;
0289 noServo = 0;
0290 else
0291 indMed = ones(length(d.array.frame.utc)*5,1);
0292 indMed = logical(indMed);
0293 noServo = 1;
0294 end
0295
0296 indSlow = d.array.frame.utc>=mjdstartorig & d.array.frame.utc<mjdstoporig;
0297
0298
0299 indF = reshape(indFast, [5 20 length(indSlow)]);
0300 indF = permute(indF, [3 1 2]);
0301 indFs = mean(mean(indF,3),2)==1;
0302
0303 indM = reshape(indMed, [5 length(indSlow)]);
0304 indM = permute(indM, [2 1]);
0305 indMs = mean(indM,2)==1;
0306
0307 ind = indMs & indFs & indSlow;
0308
0309
0310
0311
0312 d = framecut(d, ind, 'regular');
0313
0314 return;
0315