Jan2014_M42_and_TauA_1

0001 Intensities = {'I'};%,'I1','I2'}
0002 cbass='N'
0003 save_path = 'Raster_analysis_newcal2'
0004 
0005 % 16-Jul M42
0006 
0007 reduced_data_file='16-Jul-2013:13:48:20_m42_I_pipelined.mat'
0008 start_date = '16-Jul-2013:13:48:20' % M42
0009 stop_date = '16-Jul-2013:21:08:10'
0010 
0011 run_lots_of_rasters2(cbass,Intensities,save_path,[],start_date,stop_date,reduced_data_file)
0012 
0013 % 17-Jul TauA
0014 reduced_data_file='17-Jul-2013:13:53:34_taua_I_pipelined.mat'
0015 start_date = '17-Jul-2013:13:53:34' % TauA
0016 stop_date = '17-Jul-2013:21:57:50'
0017 
0018 
0019 run_lots_of_rasters2(cbass,Intensities,save_path,[],start_date,stop_date,reduced_data_file)
0020 
0021 % Lokk at resulting calibration matrices
0022  
0023 m42_calmat.I  = load('./Raster_analysis_newcal/16-Jul-2013:13:48:20_m42_I_caldata.mat')
0024 m42_calmat.I1 =load('./Raster_analysis_newcal/16-Jul-2013:13:48:20_m42_I1_caldata.mat')
0025 m42_calmat.I2 =load('./Raster_analysis_newcal/16-Jul-2013:13:48:20_m42_I2_caldata.mat')
0026 
0027 tauA_calmat.I =load('./Raster_analysis_newcal/17-Jul-2013:13:53:34_taua_I_caldata.mat')
0028 tauA_calmat.I1 =load('./Raster_analysis_newcal/17-Jul-2013:13:53:34_taua_I1_caldata.mat')
0029 tauA_calmat.I2 =load('./Raster_analysis_newcal/17-Jul-2013:13:53:34_taua_I2_caldata.mat')
0030 %%
0031 %Run as far as saving .red files for m42 and ytaua since Nov2012-Oct 2013
0032 clear;
0033 close all;
0034 Intensities = {'I','I1','I2'}
0035 cbass='N'
0036 save_path = 'Raster_analysis_newcal2'
0037 input_filename = 'taua_Nov2012_Oc2013.txt'
0038 run_lots_of_rasters2(cbass,Intensities,save_path,input_filename)
0039 
0040 Intensities = {'I','I1','I2'}
0041 cbass='N'
0042 save_path = 'Raster_analysis_newcal2'
0043 input_filename = 'm42_Nov2012_Oc2013.txt'
0044 run_lots_of_rasters2(cbass,Intensities,save_path,input_filename
0045 
0046 %%
0047 % Make contour plots for each raster in RA and DEC
0048 Intensities = {'I'};%,'I1','I2'}
0049 IntensityType='I'
0050 cbass='N'
0051 save_path = 'Raster_analysis_newcal2'
0052 start_date='17-Jul-2013:13:53:34'
0053 source_name='TauA'
0054 load('17-Jul-2013:13:53:34_taua_I_dgood.mat')
0055 nchan=9;
0056 plot_labels = {'I1','Q1','U1','Q2','U2','Q3','U3','I2','V'};
0057 fit_IQUV2_radec(nchan, plot_labels, cbass,IntensityType,save_path,start_date,source_name,dgood)
0058