% datadir = 'C:\Documents and Settings\sethb\My Documents\BOG\CO2\Regional\Filtered Data';
datadir = '/ssdsdata/ocr3000cimt/TEMP';
files = dir([datadir '/*.mat']);

seas_col_num = 1;
year_col_num = 2;
box_col_num = 3;
pco2_col_num = 4;
tmp_col_num = 7;
chla_col_num = 8;
total_pco2 = zeros(2,11,4);
total_pco2_points = zeros(2,11,4);

for f = 1:length(files) 
   filename = files(f).name; 
   indata = open ([datadir '/' filename]);
   
   season = indata(1).outdata(1,seas_col_num);
   year = indata(1).outdata(1,year_col_num);
   
   for b = 1:4
       m = find(indata(1).outdata(:,box_col_num) == b);
       pco2_sum = sum(indata(1).outdata(m,pco2_col_num));
       pco2_points = length(m);

       n = find(~isnan(indata(1).outdata(m,tmp_col_num)));
       tmp_sum = sum(indata(1).outdata(m(n),tmp_col_num));
       tmp_points = length(n);

       o = find(~isnan(indata(1).outdata(m,chla_col_num)));
       chla_sum = sum(indata(1).outdata(m(o),chla_col_num));
       chla_points = length(n);

       total_pco2(season,year,b) = total_pco2(season,year,b) + pco2_sum;
       total_pco2_points(season,year,b) = total_pco2_points(season,year,b) + pco2_points;
       
       total_tmp(season,year,b) = total_pco2(season,year,b) + pco2_sum;
       total_pco2_points(season,year,b) = total_pco2_points(season,year,b) + pco2_points;
       
       total_pco2(season,year,b) = total_pco2(season,year,b) + pco2_sum;
       total_pco2_points(season,year,b) = total_pco2_points(season,year,b) + pco2_points;
   end
   
end