compare_cond

 COMPARE_COND -  Diagnostics plots to compare CTD and samples conductivity values.

 CTD Calibration toolbox.

 USAGE:
 compare_cond(cruiseid);

INPUT: cruiseid: Cruise identification name.
   
 OUTPUT: A series of plots to compare CTD and bottle samples conductivity values 
          
 DESCRIPTION: This function uses the summary conductivity calibration file created by  
 make_cond_sum and do a series of diagnostics plots to compare CTD and
 bottle samples conductivity values to evaluate the quality of the data
 that we are coleecting, as well, any potential problems with
 niskin/sample bottles or sensors.
 AUTHORS:Carlos Fonseca 
         UM/CIMAS
         Derrick Snowden
         NOAA/AOML/PhOD
         Tue May 11 11:23:37 EDT 2004
 CHANGELOG: 
   23-Jul-2004, Version 1.0
        * Initial version.

0001 function compare_cond(cruiseid)
0002 % COMPARE_COND -  Diagnostics plots to compare CTD and samples conductivity values.
0003 %
0004 % CTD Calibration toolbox.
0005 %
0006 % USAGE:
0007 % compare_cond(cruiseid);
0008 %
0009 %INPUT: cruiseid: Cruise identification name.
0010 %
0011 % OUTPUT: A series of plots to compare CTD and bottle samples conductivity values
0012 %
0013 % DESCRIPTION: This function uses the summary conductivity calibration file created by
0014 % make_cond_sum and do a series of diagnostics plots to compare CTD and
0015 % bottle samples conductivity values to evaluate the quality of the data
0016 % that we are coleecting, as well, any potential problems with
0017 % niskin/sample bottles or sensors.
0018 % AUTHORS:Carlos Fonseca
0019 %         UM/CIMAS
0020 %         Derrick Snowden
0021 %         NOAA/AOML/PhOD
0022 %         Tue May 11 11:23:37 EDT 2004
0023 % CHANGELOG:
0024 %   23-Jul-2004, Version 1.0
0025 %        * Initial version.
0026 
0027 % Retrieve the necessary preferences for this cruise
0028 group = group_name(cruiseid);
0029 
0030 % cruise_dir is the main calibration data directory
0031 if ispref(group,'cal_data_dir')
0032     cruise_dir=getpref(group,'cal_data_dir');
0033 else
0034     error(['The cal_data_dir preference was not set for cruise: ' cruiseid '.  Run register_cruise.m']);
0035 end
0036 
0037 % Load data from the cruise database
0038 db_file = fullfile(cruise_dir,[cruiseid '_db.mat']);
0039 if exist(db_file)==2
0040     load(db_file);
0041 else
0042     error(['The cruise data base file was not found on the path.  ' db_file]);
0043 end
0044 n=who('cond_cal');
0045  if isempty(n)==1
0046  error(['cond_cal not found. Run make_cond_sum.m ']);
0047  return;
0048  end
0049 
0050 figure;
0051 plot(cond_cal.btlsal,cond_cal.ctdsal1,'r+')
0052 hold on;
0053 xlabel('Bottle Salinity (psu)')
0054 ylabel('CTD Salinity (psu)')
0055 title('Primary Sensor')
0056 
0057 figure
0058 plot(cond_cal.btlsal,cond_cal.ctdsal1,'r+')
0059 text(cond_cal.btlsal,cond_cal.ctdsal1,int2str(cond_cal.stnnbr))
0060 hold on;
0061 xlabel('Bottle Salinity')
0062 ylabel('CTD Salinity')
0063 title('Primary Sensor - Station Number')
0064 
0065 figure
0066 plot(cond_cal.stnnbr,cond_cal.btlsal-cond_cal.ctdsal1,'r+')
0067 xlabel('Station number')
0068 ylabel('Bottle Salinity - CTD Salinity')
0069 title('Primary Sensor') 
0070  
0071 figure
0072 plot(cond_cal.niskinbtlpos,cond_cal.btlsal-cond_cal.ctdsal1,'r+')
0073 xlabel('Niskin number')
0074 ylabel('Bottle Salinity - CTD Salinity')
0075 title('Primary Sensor') 
0076  
0077 figure
0078 plot(cond_cal.sampnbr,cond_cal.btlsal-cond_cal.ctdsal1,'r+')
0079 xlabel('Sample Bottle number')
0080 ylabel('Bottle Salinity - CTD Salinity')
0081 title('Primary Sensor') 
0082   
0083 figure
0084 plot(cond_cal.ctdprs,cond_cal.btlsal-cond_cal.ctdsal1,'r+')
0085 xlabel('Pressure')
0086 ylabel('Bottle Salinity - CTD Salinity')
0087 title('Primary Sensor')
0088   
0089 pause
0090 
0091 figure;
0092 plot(cond_cal.btlsal,cond_cal.ctdsal2,'r+')
0093 hold on;
0094 xlabel('Bottle Salinity (psu)')
0095 ylabel('CTD Salinity (psu)')
0096 title('Secondary Sensor')
0097 
0098 figure
0099 plot(cond_cal.btlsal,cond_cal.ctdsal2,'r+')
0100 text(cond_cal.btlsal,cond_cal.ctdsal2,int2str(cond_cal.stnnbr))
0101 hold on;
0102 xlabel('Bottle Salinity')
0103 ylabel('CTD Salinity')
0104 title('Secondary Sensor - Station Number')
0105 
0106 figure
0107 plot(cond_cal.stnnbr,cond_cal.btlsal-cond_cal.ctdsal2,'r+')
0108 xlabel('Station number')
0109 ylabel('Bottle Salinity - CTD Salinity')
0110 title('Secondary Sensor') 
0111  
0112 figure
0113 plot(cond_cal.niskinbtlpos,cond_cal.btlsal-cond_cal.ctdsal2,'r+')
0114 xlabel('Niskin number')
0115 ylabel('Bottle Salinity - CTD Salinity')
0116 title('Secondary Sensor') 
0117  
0118 figure
0119 plot(cond_cal.sampnbr,cond_cal.btlsal-cond_cal.ctdsal2,'r+')
0120 xlabel('Sample Bottle number')
0121 ylabel('Bottle Salinity - CTD Salinity')
0122 title('Secondary Sensor') 
0123   
0124 figure
0125 plot(cond_cal.ctdprs,cond_cal.btlsal-cond_cal.ctdsal2,'r+')
0126 xlabel('Pressure')
0127 ylabel('Bottle Salinity - CTD Salinity')
0128 title('Secondary Sensor')
0129   
0130 figure
0131 plot(cond_cal.ctdprs,cond_cal.ctdsal1-cond_cal.ctdsal2,'r+')
0132   xlabel('Pressure, dbar')
0133   ylabel('Primary - Secondary Sensor Salinity')
0134   title('Primary  - Secondary Sensor')