longTermOp

 records long term opacity data alongside weather info

0001 % records long term opacity data alongside weather info
0002 clear
0003 
0004 home = '/scratch/falcon_4/mirfan/cbass';
0005 logfile = ([home, '/cbass_analysis/log/obs_log.html']);
0006 start_scan= 1; %which scan to start from in case it stops half way through
0007 max_time= 1;%maximum length of each data section in hours
0008 
0009 t3=extract_obslog(logfile,'cbass_survey_0.2','10-Feb-2012:00:00:41','20-Feb-2012:23:59:51');
0010 
0011  t=t3;
0012 
0013 for Z=1:length(t)
0014     tstring=strtrim(t{Z,1});
0015     tstring2=strtrim(t{Z,2});
0016     tt(Z,1) = datenum(tstring,'dd-mmm-yyyy:HH:MM:ss');
0017     tt(Z,2) = datenum(tstring2,'dd-mmm-yyyy:HH:MM:ss');
0018     
0019 end
0020 
0021 aa=[];
0022 for Z=1:length(t)
0023    
0024     
0025     startt=tt(Z,1);
0026     endt=tt(Z,2);
0027     scan_time=endt-startt;
0028     while (scan_time > max_time/24)
0029       aa=[aa;[startt startt+max_time/24]];
0030       startt=startt+(max_time+0.00001)/24;
0031       scan_time=endt-startt;
0032     end
0033     aa=[aa;[startt endt]];
0034           
0035 end
0036 %make new time data section to avoid huge data sections
0037 for Z=1:length(aa)
0038     tsplit{Z,1}=datestr(aa(Z,1),'dd-mmm-yyyy:HH:MM:ss');
0039     tsplit{Z,2}=datestr(aa(Z,2),'dd-mmm-yyyy:HH:MM:ss');
0040 end
0041 
0042 numscans=length(tsplit);
0043 n=0;
0044 
0045 
0046 Z=start_scan;
0047 
0048 for Z=start_scan:numscans
0049     close all;
0050     try
0051         time1{Z}=tsplit{Z,1};
0052         time2{Z}=tsplit{Z,2};
0053         disp(['Time Start: ',time1{Z},' Time End: ',time2{Z}]);
0054         disp(['Total Reduction Run from ',tsplit{1},' to ',tsplit{length(tsplit),2}]);
0055         d = pipe_read(time1{Z},time2{Z});
0056         disp(['Total Reduction Run from ',tsplit{1},' to ',tsplit{length(tsplit),2}]);
0057     catch
0058         disp('Failed to Read in Data')
0059         clear flags
0060         continue
0061     end 
0062     
0063   [years,months,days,hours,mins,secs] = datevec(datenum(tsplit{Z}, 'dd-mmm-yyyy:HH:MM:ss'));
0064 
0065 
0066     try      
0067         Scans(Z) = date2mjd(years,months,days,hours,mins,secs);
0068         d = reduceData_simple(d);
0069         Temp(Z) = mean(d.array.weather.airTemperature);
0070         relH(Z) = mean(d.array.weather.relativeHumidity); 
0071         
0072         skydip(Z) = NaN;
0073         SkyDipMethod2v2
0074         skydip(Z) = OpacityOne;
0075         
0076         pipe = calculateTau(d);
0077         pipeTheo(Z) = last(pipe(:,4));
0078 
0079         myTheo(Z) = CalcTheoOp(d); 
0080         noiseD(Z) = NoiseDiodeOp(d);
0081        
0082     catch
0083         
0084          
0085       
0086      disp('No sky dips')
0087      clear flags
0088         continue
0089     
0090     end 
0091 
0092   
0093     clear d
0094     clear flags
0095    
0096 
0097 end
0098 
0099 All_Time = [time1;time2]';
0100 [rows,cols] = size(All_Time);
0101 
0102 fid3 = fopen('/scratch/falcon_4/mirfan/cbass/cbass_analysis/MelsTests/Feb2012OpV2Scheds.txt', 'a');
0103 
0104 for g = 1:rows
0105 fprintf(fid3, '%s %s \n', All_Time{g,1:2});
0106 end
0107 
0108 
0109 fid = fopen('/scratch/falcon_4/mirfan/cbass/cbass_analysis/MelsTests/Feb2012OpV2.txt','a');
0110 for i = 1 : Z
0111 fprintf(fid, '%10.4f %7.4f %7.4f %6.4f %6.4f %6.4f %6.4f \n',Scans(i), Temp(i), relH(i), pipeTheo(i), myTheo(i), skydip(i), noiseD(i));
0112 end
0113 fclose(fid);