write_out_fits_and_detrendfits

0001 addpath([home,'/',installeddir,'/Angelas_Raster_Code/']);
0002 addpath([home,'/',installeddir,'/Raster_analysis_newcal2//']);
0003 
0004 % 19/3/2014 - look at data in n_raster_elliptical_newcal directories
0005 datadir = 'n_raster_elliptical_newcal';
0006 datadir_fits = 'n_raster_elliptical_newcal_fits';
0007 %Load in the pipelined data
0008 %TauA
0009 %load('./Raster_analysis_newcal2/17-Jul-2013:13:53:34_taua_I_pipelined.mat')
0010 load(['./',datadir,'/17-Jul-2013:13:53:34_taua_I_pipelined.mat'])
0011 %M42
0012 load(['./',datadir,'/16-Jul-2013:13:48:20_m42_I_pipelined.mat'])
0013 % Run reduceData to make the .fits file
0014 
0015 dred = reduceData(d,'redScript_newcal_raster_2014_fits.red',0,datadir_fits)
0016 clear
0017 
0018 %load in detrended data
0019 load(['./',datadir,'/17-Jul-2013:13:53:34_taua_I_dgood.mat'])
0020 %load('./Raster_analysis_newcal2/17-Jul-2013:13:53:34_taua_I_dgood.mat')
0021 load(['./',datadir,'/16-Jul-2013:13:48:20_m42_I_dgood.mat'])
0022 % gets loaded as dgood - lets set to d and remove data(:,9)
0023 % then write out fits
0024 
0025 d = dgood;
0026 d.antenna0.receiver.data(:,9)=[];
0027 
0028 %  For old dtaa, need to add day night flag
0029 % Set up a daynight flag
0030 % sun within  sunlim degrees
0031 [dist az el]      = calcSourceDist(d, 'sun');
0032 d.antenna0.servo.solarDist = dist;
0033 az = wrap360(az);
0034 raz=round(60*az);
0035 raz(raz >= 21600) = 0;
0036 d.flags.dayNight = ones(size(az),'uint8');
0037 % make sure that we don't overwrite the last .fits file
0038 system(['mv /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/17-Jul-2013:13:53:34/fits/17-Jul-2013:13:53:34.fits   /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/17-Jul-2013:13:53:34/fits/17-Jul-2013:13:53:34_nodetrend.fits']);
0039   
0040 % make new fits
0041 dred = reduceData(d,'redScript_newcal_raster_2014_fits.red',0,datadir_fits)
0042 system(['mv /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/17-Jul-2013:13:53:34/fits/17-Jul-2013:13:53:34.fits   /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/17-Jul-2013:13:53:34/fits/17-Jul-2013:13:53:34_detrend.fits']);
0043 
0044 %%
0045 %M42 - load data
0046 load('./Raster_analysis_newcal2/16-Jul-2013:13:48:20_m42_I_pipelined.mat')
0047 
0048 % write out fits of data (no detrending)
0049 dred = reduceData(d,'redScript_newcal_raster_2014_fits.red')
0050 
0051 clear
0052 
0053 %reload detrended data - has 9 columns in data: I1,Q1,U1,Q2,U2,Q3,U3,I2,V
0054 load('./Raster_analysis_newcal2/16-Jul-2013:13:48:20_m42_I_dgood.mat')
0055 
0056 % gets loaded as dgood - lets set to d and remove data(:,9)
0057 % then write out fits
0058 
0059 d = dgood;
0060 d.antenna0.receiver.data(:,9)=[];
0061 
0062 % make sure that we don't overwrite the last .fits file
0063 system(['mv /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/16-Jul-2013:13:49:14/fits/16-Jul-2013:13:49:14.fits   /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/16-Jul-2013:13:49:14/fits/16-Jul-2013:13:49:14_nodetrend.fits']);
0064   
0065 % make new fits
0066 dred = reduceData(d,'redScript_newcal_raster_2014_fits.red')
0067 system(['mv /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/16-Jul-2013:13:49:14/fits/16-Jul-2013:13:49:14.fits   /usersVol1/act/CBASS/cbass_analysis/reduc/raster_tests/16-Jul-2013:13:49:14/fits/16-Jul-2013:13:49:14_detrend.fits']);