#ifndef __JACOUSTICS__JUNESCO__
#define __JACOUSTICS__JUNESCO__

/**
 * \file
 *
 * Sound velocity according UNESCO equation.
 * \author mdejong
 */
namespace JACOUSTICS {}
namespace JPP { using namespace JACOUSTICS; }

namespace JACOUSTICS {

  /**
   * Get sound velocity according UNESCO equation.
   *
   * \param  D            Depth [m]
   * \param  S            Salinity [ppk]
   * \param  T            Temperature [C]
   * \return              sound velocity [m/s]
   */
  inline double getVelocity(const double D, const double S, const double T)
  {
    return (1448.96
	    + 4.591 * T
	    - 5.304e-2  * T*T
	    + 2.374e-4  * T*T*T
	    + 1.340 * (S - 35)
	    + 1.630e-2  * D
	    + 1.675e-7  * D*D
	    - 1.025e-2  * T * (S - 35)
	    - 7.139e-13 * T*D*D*D);
  }
}

#endif