diff --git a/tsg_util/tsg_plot_SalTsgSample.m b/tsg_util/tsg_plot_SalTsgSample.m
new file mode 100644
index 0000000000000000000000000000000000000000..dd7ab959f0ebda5f8a5a409eaa866d305dc3dc91
--- /dev/null
+++ b/tsg_util/tsg_plot_SalTsgSample.m
@@ -0,0 +1,74 @@
+function tsg_plot_SalTsgSample( hGUI, hAxes )
+%
+% Function to plot the Salinity from TSG and smaple meausrements
+%
+% Input
+% -----
+% hGUI ............... Handel to the user interface
+% hPlotAxes .......... Handels to the graphic axes
+%
+% Output
+% ------
+
+% Retrieve named application data
+% -------------------------------
+tsg    = getappdata( hGUI, 'tsg_data');
+sample = getappdata( hGUI, 'sample' );
+qc     = getappdata( hGUI, 'qcColor');
+
+
+% Plot the Salinity with the right color code
+% -------------------------------------------
+axes( hAxes(1));
+
+ind = find(tsg.PSAL_QC <= qc.Code.NO_CONTROL );
+plot(hAxes(1), tsg.TIME(ind), tsg.PSAL(ind), ...
+               strcat('.',qc.Color.NO_CONTROL));
+
+hold on;
+
+ind = find(tsg.PSAL_QC == qc.Code.GOOD );
+plot(hAxes(1), tsg.TIME(ind), tsg.PSAL(ind), ...
+               strcat('.',qc.Color.GOOD));
+           
+ind = find(tsg.PSAL_QC == qc.Code.PROBABLY_GOOD );
+plot(hAxes(1), tsg.TIME(ind), tsg.PSAL(ind), ...
+               strcat('.',qc.Color.PROBABLY_GOOD));
+           
+ind = find(tsg.PSAL_QC == qc.Code.PROBABLY_BAD );
+plot(hAxes(1), tsg.TIME(ind), tsg.PSAL(ind), ...
+               strcat('.',qc.Color.PROBABLY_BAD));
+           
+ind = find(tsg.PSAL_QC == qc.Code.BAD );
+plot(hAxes(1), tsg.TIME(ind), tsg.PSAL(ind), ...
+               strcat('.',qc.Color.BAD));
+
+% Plot sample measurements
+% ------------------------
+plot(hAxes(2), sample.TIME, sample.PSAL, '.');
+
+%plot(hAxes(3), tsg.TIME, tsg.TEMP_TSG, '.');
+
+% Formatted x-TIME axes
+% ---------------------
+datetick(hAxes(1), 'x');
+datetick(hAxes(2), 'x');
+%datetick(hAxes(3), 'x');
+
+% Write some 'Y' label
+% ------------------
+set(get(hAxes(1), 'Ylabel'), 'String', 'Salinity');
+set(get(hAxes(2), 'Ylabel'), 'String', 'Salinity (samples');
+%set(get(hAxes(3), 'Ylabel'), 'String', 'Ship Velocity');
+
+% Make the axes visible
+% ---------------------
+set(hAxes(1), 'Visible', 'on' );
+set(hAxes(2), 'Visible', 'on' );
+%set(hAxes(3), 'Visible', 'on' );
+
+% The 3 axes will behave identically when zoomed and panned
+% ---------------------------------------------------------
+%linkaxes([hAxes(1),hAxes(2),hAxes(3)], 'x');
+
+end