Home > ctdcal > final_results.m

final_results

PURPOSE ^

SYNOPSIS ^

This is a script file.

DESCRIPTION ^





CROSS-REFERENCE INFORMATION ^

This function calls:

This function is called by:





SOURCE CODE ^


0001 sub=length(boxy0(ii_old));
0002 diary on
0003 if(pri_sec==0)
0004   disp([' Primary Sensor']);
0005  else
0006  disp([' Secondary Sensor']);
0007 end
0008 disp(['Number of points used   ',int2str(sub)]);
0009 disp(['Total number of points  ',int2str(tot)]);
0010 disp(['Percent of points used  ',num2str(100*sub/tot,4)]);
0011 disp(['fit standard deviation  ',num2str(s_old,4)]);
0012 disp(['fit mean residual  ',num2str(m_old,4)]);
0013 disp(['fit mean absolute residual  ',num2str(mean(abs(d_old)),4)]);
0014 disp(['fit mean absolute residual (%)  ',num2str(m_abs_old,4)]);
0015 disp(['Soc = ', num2str(Soc_old)]);
0016 disp(['Voffset = ', num2str(Voffset_old)]);
0017 disp(['Boc = ', num2str(Boc_old)]);
0018 disp(['tcor = ', num2str(tcor_old)]);
0019 disp(['pcor = ', num2str(pcor_old)]);
0020 disp(['cal = ', num2str(cal_old)]);
0021 disp(['gama = ', num2str(gama_old)]);
0022 disp(['p1 = ', num2str(p1_old)]);
0023 disp(['p2 = ', num2str(p2_old)]);
0024 disp(['p3 = ', num2str(p3_old)]);
0025 disp(['p4 = ', num2str(p4_old)]);
0026 disp(['p5 = ', num2str(p5_old)]);
0027 disp(['p6 = ', num2str(p6_old)]);
0028 diary off
0029 figure;
0030 subplot(2,1,1);
0031 deep_idx=find(prs(ii_old)>=H_COLOR);
0032 shallow_idx=find(prs(ii_old)<H_COLOR);
0033 plot(stn(ii_old(deep_idx)),d_old(deep_idx),'bo');
0034 hold on
0035 plot(stn(ii_old(shallow_idx)),d_old(shallow_idx),'ro');
0036 grid on
0037 title(['Dissolved Oxygen Calibration - Final Result : ',num2str(iter-1),' Data used: ',num2str(100*sub/tot,4),'%']);
0038 ylabel('Bottle - CTD Oxygen (ml/l)');
0039 xlabel('Station Number');
0040 %gtext({'mean residual= ',num2str(m_old,4),'std.dev= ',num2str(s_old,4)})
0041 %set(gca,'YLim', [-0.4 0.4]);
0042 %set(gca,'Ytick', [-0.4 -0.2 0 0.2 0.4]);
0043 hold on
0044 subplot(2,1,2);
0045 plot(prs(ii_old(deep_idx)),d_old(deep_idx),'bo');
0046 hold on
0047 plot(prs(ii_old(shallow_idx)),d_old(shallow_idx),'ro');
0048 grid on
0049 title([' mean residual = ',num2str(m_old,4),'  std.dev = ',num2str(s_old,4)]);
0050 ylabel('Bottle - CTD Oxygen (ml/l)');
0051 xlabel('Pressure');
0052 hold on
0053 
0054 figure;
0055 subplot(2,1,1);
0056 plot(nis(ii_old(deep_idx)),d_old(deep_idx),'bo');
0057 hold on
0058 plot(nis(ii_old(shallow_idx)),d_old(shallow_idx),'ro');
0059 grid on
0060 title(['Dissolved Oxygen Calibration - Final Result: ',num2str(iter-1)]);
0061 ylabel('Bottle - Primary CTD Oxygen (ml/l)');
0062 xlabel('Niskin number');
0063 %set(gca,'YLim', [-0.4 0.4]);
0064 %set(gca,'Ytick', [-0.4 -0.2 0 0.2 0.4]);
0065 hold on
0066 subplot(2,1,2);
0067 plot(o2_bot_num(ii_old(deep_idx)),d_old(deep_idx),'bo');
0068 hold on
0069 plot(o2_bot_num(ii_old(shallow_idx)),d_old(shallow_idx),'ro');
0070 grid on
0071 %title('Bottle-Secondary CTD Oxygen vs. pressure - all stations');
0072 ylabel('Bottle - CTD Oxygen (ml/l)');
0073 xlabel('Sample Bottle Number');
0074 hold on
0075 
0076 figure;
0077 subplot(2,1,1);
0078 plot(stn(ii_old(deep_idx)),ox_units(boxy0(ii_old(deep_idx)),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk')-ox_units(new_oxy_old(deep_idx),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk'),'bo');
0079 hold on
0080 plot(stn(ii_old(shallow_idx)),ox_units(boxy0(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk')-ox_units(new_oxy_old(shallow_idx),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk'),'ro');
0081 grid on
0082 title(['Dissolved Oxygen Calibration - Final Result']);
0083 ylabel('Bottle -  CTD Oxygen (\mu mol/kg)');
0084 xlabel('Station number');
0085 diff_o2=ox_units(boxy0(ii_old),csal(ii_old),sw_ptmp(ctem(ii_old),csal(ii_old),cprs(ii_old),0),'mll2mmk')-ox_units(new_oxy_old,csal(ii_old),sw_ptmp(ctem(ii_old),csal(ii_old),cprs(ii_old),0),'mll2mmk');
0086 mumol=nan_mean(diff_o2);
0087 sumol=nan_std(diff_o2);
0088 %gtext({'mean residual= ',num2str(mumol,4),'std.dev= ',num2str(sumol,4)})
0089 %set(gca,'YLim', [-0.4 0.4]);
0090 %set(gca,'Ytick', [-0.4 -0.2 0 0.2 0.4]);
0091 hold on
0092 subplot(2,1,2);
0093 plot(prs(ii_old(deep_idx)),ox_units(boxy0(ii_old(deep_idx)),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk')-ox_units(new_oxy_old(deep_idx),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk'),'bo');
0094 hold on
0095 plot(prs(ii_old(shallow_idx)),ox_units(boxy0(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk')-ox_units(new_oxy_old(shallow_idx),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk'),'ro');
0096 grid on
0097 title([' mean residual = ',num2str(mumol,4),'   std.dev = ',num2str(sumol,4)]);
0098 ylabel('Bottle - Primary CTD Oxygen (\mu mol/kg)');
0099 xlabel('Pressure');
0100 hold on
0101 
0102 aux_contour=[stn(ii_old) prs(ii_old) d_old round(d_old*100)];
0103 std_dev=std(d_old);
0104 pos_idx=find(aux_contour(:,4)>=0);
0105 neg_idx=find(aux_contour(:,4)<0);
0106 figure
0107 plot(aux_contour(pos_idx,1),aux_contour(pos_idx,2),'r+','MarkerSize',1)
0108 ax1=gca;
0109 set(ax1,'ydir','reverse')
0110 hold on
0111 plot(aux_contour(neg_idx,1),aux_contour(neg_idx,2),'k+','MarkerSize',1)
0112 text(aux_contour(pos_idx,1),aux_contour(pos_idx,2),num2str(aux_contour(pos_idx,4)),'Color','r')
0113 text(aux_contour(neg_idx,1),aux_contour(neg_idx,2),num2str(aux_contour(neg_idx,4)),'Color','k')
0114 title('Final Result : Bottle - CTD Oxygen (.10^{-2} ml/l)')
0115 xlabel('Station')
0116 ylabel('Pressure (dbar)')
0117 figure;
0118 subplot(2,1,1);
0119 plot(stn(ii_old(deep_idx)),ox_units(boxy0(ii_old(deep_idx)),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk')-ox_units(new_oxy_old(deep_idx),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk'),'bo');
0120 hold on
0121 plot(stn(ii_old(shallow_idx)),ox_units(boxy0(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk')-ox_units(new_oxy_old(shallow_idx),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk'),'ro');
0122 grid on
0123 title(['Dissolved Oxygen Calibration - Final Result']);
0124 ylabel('Bottle -  CTD Oxygen (\mu mol/kg)');
0125 xlabel('Station number');
0126 diff_o2=ox_units(boxy0(ii_old),csal(ii_old),sw_ptmp(ctem(ii_old),csal(ii_old),cprs(ii_old),0),'mll2mmk')-ox_units(new_oxy_old,csal(ii_old),sw_ptmp(ctem(ii_old),csal(ii_old),cprs(ii_old),0),'mll2mmk');
0127 mumol=nan_mean(diff_o2);
0128 sumol=nan_std(diff_o2);
0129 %gtext({'mean residual= ',num2str(mumol,4),'std.dev= ',num2str(sumol,4)})
0130 %set(gca,'YLim', [-0.4 0.4]);
0131 %set(gca,'Ytick', [-0.4 -0.2 0 0.2 0.4]);
0132 hold on
0133 subplot(2,1,2);
0134 plot(ox_units(boxy0(ii_old(deep_idx)),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk')-ox_units(new_oxy_old(deep_idx),csal(ii_old(deep_idx)),sw_ptmp(ctem(ii_old(deep_idx)),csal(ii_old(deep_idx)),cprs(ii_old(deep_idx)),0),'mll2mmk'),prs(ii_old(deep_idx)),'bo');
0135 hold on
0136 plot(ox_units(boxy0(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk')-ox_units(new_oxy_old(shallow_idx),csal(ii_old(shallow_idx)),sw_ptmp(ctem(ii_old(shallow_idx)),csal(ii_old(shallow_idx)),cprs(ii_old(shallow_idx)),0),'mll2mmk'),prs(ii_old(shallow_idx)),'ro');
0137 grid on
0138 title([' mean residual = ',num2str(mumol,4),'   std.dev = ',num2str(sumol,4)]);
0139 ylabel('Bottle - Primary CTD Oxygen (\mu mol/kg)');
0140 xlabel('Pressure');
0141 hold on
0142 
0143 figure;
0144 plot(boxy0(ii_old(deep_idx)),new_oxy_old(deep_idx),'b+');
0145 hold on
0146 plot(boxy0(ii_old(shallow_idx)),new_oxy_old(shallow_idx),'r+');
0147 legend('deep','shallow',0)
0148 grid on
0149 title('Bottle vs. CTD Oxygen - Final Result');
0150 %hold on
0151 ylabel('Bottle Oxygen (ml/l)');
0152 xlabel('CTD Oxygen (ml/l)');
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162


Generated on Fri 08-Oct-2004 11:57:17 by m2html © 2003