JUNESCO.hh

Go to the documentation of this file.

#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