autoReduction

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 function autoReduction(field,script,rn)

 Mass reduction of data: auto flagging, posting of plots on web, and
 saving of reduced data in data archive

 field - see getTracks.m
 script - data reduction script
 rn - mass reduction number 

 SJCM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function autoReduction(field,script,rn)
0002 
0003 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0004 % function autoReduction(field,script,rn)
0005 %
0006 % Mass reduction of data: auto flagging, posting of plots on web, and
0007 % saving of reduced data in data archive
0008 %
0009 % field - see getTracks.m
0010 % script - data reduction script
0011 % rn - mass reduction number
0012 %
0013 % SJCM
0014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0015 
0016 
0017 if(nargin==1)
0018     eval(sprintf('load autoreduce_%s.mat', field));
0019     start=resume;
0020 elseif(nargin==0)
0021   eval(sprintf('load autoreduce.mat', field));
0022   start=resume;
0023 else
0024   if (~exist('script'))
0025     disp('autoReduction:: Script does not exist!')
0026     return
0027   elseif (~exist('rn'))
0028     disp('autoReduction:: rn does not exist!')
0029     return
0030   else
0031     readscript(script);
0032   end
0033   list=getTracks(field);
0034   start=1;
0035 end 
0036 
0037 % cmb is different
0038 if strmatch('cmb',field)
0039   field='cmb';
0040 end
0041 
0042 if ~exist('oldFlagsRn')
0043   oldFlagsRn = 0;
0044 end
0045 
0046 
0047 for i=start:length(list)
0048    d = read_arc(list{i}{1}, list{i}{2});
0049 
0050    dcal=reduceData(d,script,1,1, field);
0051    webpipe(dcal,field,rn);
0052    archivedata(dcal,field,rn,list{i});
0053 
0054    resume=i+1;
0055   eval(sprintf('save autoreduce_%s.mat list resume script rn field', field));
0056 end
0057 
0058 
0059 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0060 function archivedata(dcal,field,rn,list)
0061 
0062 keyboard;
0063 [dir,fn]=archivepath(field,rn,list)
0064 eval(sprintf('%s=dcal;',fn));
0065 eval(sprintf('save %s/%s %s',dir,fn,fn));
0066 
0067 
0068 % write to fits file
0069 display('autoReduction::  Writing fits file not done yet');
0070 
0071 return;