% program to read compile_noref.txt file from Pete Strutton
% compile_noref.txt is a high resolution compilation of all sensors he has
  % access to.
% INPUTS:
  % date, exceldate, year_day, lon, lat, flr_2m, chl_2m, flr_10m, chl_10m, flr_15m, chl_15m,
  % b100_2m, b125_2m, b150_2m, b_back_2m, b100_10m, b125_10m, b150_10m, b_back_10m, 
  % b100_15m, b125_15m, b150_15m, b_back_15m, par_1m, flr_1m, par_5m, flr_5m, 
  % par_10m, c_10m, par_20m, c_20m, par_35m, c_35m
  %
  % Data notes:
  % flr_2m, flr_10m, flr_15m data, Wetlabs dfls
  % chl data in mg/m3
  % backscatter,vsf in 1/m.sr - scaling factor and blank offset applied
  % b_back calculated
  % par data in umol/m2.s, licor sensor
  % flr_1m, flr_5m data, chelsea fluorometer
  % c_10m, c_20m, c_35m in 1/m, CStar
%
close all
clear all
[date,exceldate,yearday,lon,lat,flr_2m,chl_2m,flr_10m,chl_10m,flr_15m,chl_15m,b100_2m,b125_2m,b150_2m,b_back_2m,b100_10m,b125_10m,b150_10m,b_back_10m,b100_15m,b125_15m,b150_15m,b_back_15m,par_1m,flr_1m,par_5m,flr_5m,par_10m,c_10m,par_20m,c_20m,par_35m,c_35m] = textread( ...
    'C:\Documents and Settings\Lisa Taylor\My Documents\Data\COASOSU\data\compile_noref.txt', ...
    '%f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f', ...
    'headerlines',2, 'delimiter','\t');
allDataCompileNoref=[date exceldate yearday lon lat flr_2m chl_2m flr_10m chl_10m flr_15m chl_15m b100_2m b125_2m b150_2m b_back_2m b100_10m b125_10m b150_10m b_back_10m b100_15m b125_15m b150_15m b_back_15m par_1m flr_1m par_5m flr_5m par_10m c_10m par_20m c_20m par_35m c_35m];
save compile-noref -v6
%
% figure
% subplot(3,3,1)
% plot(date,b100_2m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 100 2m')
% title({'COASOSU'})
% subplot(3,3,2)
% plot(date,b125_2m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 125 2m')
% title({'COASOSU'})
% subplot(3,3,3)
% plot(date,b150_2m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 150 2m')
% title({'COASOSU'})
% subplot(3,3,4)
% plot(date,b100_10m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 100 10m')
% title({'COASOSU'})
% subplot(3,3,5)
% plot(date,b125_10m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 125 10m')
% title({'COASOSU'})
% subplot(3,3,6)
% plot(date,b150_10m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 150 10m')
% title({'COASOSU'})
% subplot(3,3,7)
% plot(date,b100_15m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 100 15m')
% title({'COASOSU'})
% subplot(3,3,8)
% plot(date,b125_15m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 125 15m')
% title({'COASOSU'})
% subplot(3,3,9)
% plot(date,b150_15m,'.')
% axis([min(date),max(date),0,0.003])
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b 150 15m')
% title({'COASOSU'})
% figure
% plot(date,b_back_2m)
% datetick('x','mmmyy')
% xlabel('Date')
% ylabel('b back 2m')
% title({'COASOSU Data','b back 2m vs. Date'})
%
% figure % to determine what time zone was used:
% plot(date,chl_2m,'.')
% datetick('x','mmmyy')
