diff --git a/tsg_util/calibration.m b/tsg_util/calibration.m
index 032fa26414555decb04ea880dde95182ba3e24b5..994512ec530b41e14651f2f5756295e3e9b704eb 100644
--- a/tsg_util/calibration.m
+++ b/tsg_util/calibration.m
@@ -7,12 +7,6 @@ function calibration( hMainFig )
 %
 % Compute salinity from calibrated conductivity and jacket temperature 
 %
-% Bug :
-% 1 - tsg.SSPS_CAL = sw_salt( cond/sw_c3515(),t90TOt68(temp),zeros(size(cond)));
-%     replaced by
-%     tsg.SSPS_CAL = sw_salt( cond/sw_c3515(),temp,zeros(size(cond)));
-%
-%     CSIRO toolbox modified sw_salt modified to use temprature ITS90
 %
 % $Id$
 
@@ -88,11 +82,7 @@ if ~isempty( indC ) || ~isempty( indT )
   end
   
   if ~isempty( temp ) && ~isempty( cond )
-
-% Modification du 13/02/2009
-% tsg.SSPS_CAL = sw_salt( ...
-%       cond/sw_c3515(), t90TOt68(temp), zeros(size(cond)));
-    tsg.SSPS_CAL = sw_salt( cond/sw_c3515(), temp, zeros(size(cond)));
+    tsg.SSPS_CAL = gsw_SP_from_R( cond/gsw_C3515, temp, zeros(size(cond)));
   end
 
 else
diff --git a/tsg_util/tsg_accuracy.m b/tsg_util/tsg_accuracy.m
index f8f835531743395c8bd0a94eca7fde09a8ed8429..b7167dab320597d5227af6316c567525989f8f8a 100644
--- a/tsg_util/tsg_accuracy.m
+++ b/tsg_util/tsg_accuracy.m
@@ -66,25 +66,27 @@ if strcmp( PARA, 'SSPS')
 
     if (~isempty( T) & ~isempty( C))
 
+        c68 = 1.00024;
+        c0 =  0.6766097;
+        c1 =  2.00564e-2;
+        c2 =  1.104259e-4;
+        c3 = -6.9698e-7;
+        c4 =  1.0031e-9;
+        
         % salinity error is computed from temperature/conductivity errors
         % by error propagation in the salinity equation
         % see Emery and Thomson, Data analysis methods in phys. oceano., p.273
 
-        R=C/sw_c3515();
-        rt = sw_salrt(T);
-        Rt=R./rt;
+        R = C/gsw_C3515;
+        T68 = T * c68;
+        rt = c0 + (c1 + (c2 + (c3 + c4.*T68).*T68).*T68).*T68;
+        Rt = R./rt;
 
         %Rt is a function of C,T
         %error on Rt : err(Rt)^2=(d(Rt)/dC*errC)^2+(d(Rt)/dT*errT)^2
 
-        dRtdC=1./(rt*sw_c3515());
+        dRtdC=1./(rt*gsw_C3515);
 
-        c68 = 1.00024;
-        c0 =  0.6766097;
-        c1 =  2.00564e-2;
-        c2 =  1.104259e-4;
-        c3 = -6.9698e-7;
-        c4 =  1.0031e-9;
         dRtdT=c68*R.*(c1+(2*c2+(3*c3+4*c4*T).*T*c68).*T*c68)./(rt.^2);
 
         errRt=sqrt((dRtdC*errC).^2+(dRtdT*errT).^2);