0001
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