raw_autosal

 Creating a new data set for the autosal values

0001 % Creating a new data set for the autosal values
0002 function [output]=raw_autosal(datfile,rawfile,outfile)
0003 addpath C:\WWARD2\CTD\Calibration\SeaBird_tools;
0004 addpath C:\WWARD2\CTD\Calibration\Seawater;
0005 load_dat=['load ',datfile,'.dat']
0006 eval(load_dat)
0007 aux_dat=eval(datfile)
0008 load_raw=['load ',rawfile,'.raw']
0009 eval(load_raw)
0010 aux_raw=eval(rawfile)
0011 bottle=aux_dat(:,2);
0012 year=input('Enter the year: ');
0013 month=input('Enter the month: ');
0014 day=input('Enter the day: ');
0015 hour=aux_raw(:,11);
0016 minute=aux_raw(:,12);
0017 second=aux_raw(:,13);
0018 k=1;
0019 limite=1;
0020 while limite <= length(bottle)
0021     bottle_idx=find(bottle==bottle(limite,1));
0022     limite=max(bottle_idx)+1
0023     data=aux_raw(min(bottle_idx):max(bottle_idx),1:10);
0024     d=size(data);
0025     data=reshape(data,d(1)*d(2),1)
0026     mean_rt=mean(data);
0027     std_rt=std(data);
0028     mean_rt_new=0;
0029     std_rt_new=0;
0030     last_std=std_rt;
0031     while std_rt_new < last_std
0032      if std_rt < 0.000024
0033      aux_idx=find(data >= mean_rt-1.8*std_rt & data <= mean_rt+1.8*std_rt);
0034      else
0035      aux_idx=find(data >= mean_rt-1.45*std_rt & data <= mean_rt+1.45*std_rt);
0036      end
0037      new_data=data(aux_idx);
0038      mean_rt_new=mean(new_data);
0039      std_rt_new=std(new_data);
0040      data=new_data;
0041      mean_rt=mean_rt_new;
0042      last_std=std_rt;
0043      std_rt=std_rt_new;
0044      end
0045      data
0046      length(data)
0047      pause
0048      output(k,1)=aux_dat(max(bottle_idx),1);
0049      output(k,2)=bottle(max(bottle_idx));
0050      output(k,3)=mean_rt;
0051      output(k,4)=std_rt;
0052      output(k,5)=length(data);
0053      sw_sals(max(data)/2,24)
0054      sw_sals(min(data)/2,24)
0055      output(k,6)=sw_sals(max(data)/2,24)-sw_sals(min(data)/2,24);
0056      output(k,7)=year;
0057      output(k,8)=month;
0058      output(k,9)=day;
0059      output(k,10)=hour(max(bottle_idx));
0060      output(k,11)=minute(max(bottle_idx));
0061      output(k,12)=second(max(bottle_idx));
0062      sw_sals(max(data)/2,24)
0063      sw_sals(min(data)/2,24)
0064      output(k,6)=sw_sals(max(data)/2,24)-sw_sals(min(data)/2,24);
0065      k=k+1;
0066      fid=fopen(outfile,'w')
0067      fprintf(fid,'%%Station\t Bottle\t Mean RT\t Std Dev RT\t Data Used\t QC  year month day hour min sec\n')
0068      fprintf(fid,'%d\t %d\t %1.5f\t %1.6f\t %d\t %1.4f %d\t %d\t %d\t %d\t %d\t %d\n',output')
0069      fclose(fid)
0070  end
0071  return