
function c = svel(T,S)

%  Function SVEL (WILSON) from RDSUBS.FOR used. Units = m/sec 

   P = 1.0332272; 	% Pressure = surface value ie:depth of 0.0 dBar 
   SD = S - 35.0;      	% SD is (Sal - 35.0). Assume sal=35. 

   V00 = (((7.9851e-6*T - 2.6045e-4).*T - 4.4532e-2).*T + 4.5721).*T + 1449.14;
   V01 = (7.7711e-7*T - 1.1244e-2).*T + 1.39799;
   V0 = (1.69202e-3*SD + V01).*SD + V00;
   V10 = ((4.5283e-8*T + 7.4812e-6).*T - 1.8607e-4).*T + 0.160272;
   V11 = (1.579e-9*T + 3.158e-8).*T + 7.7016e-5;
   V1 = V11*SD + V10;
   V20 = (1.8563e-9*T - 2.5294e-7).*T + 1.0268e-5;
   V2 = -1.2943e-7*SD + V20;
   V3 = -1.9646e-10*T + 3.5216e-9;
   c = (((-3.3603e-12*P + V3).*P + V2).*P + V1).*P + V0;
