Program Listing for File UnescoFreezing.hpp

↰ Return to documentation for file (src/include/UnescoFreezing.hpp)

#ifndef SRC_INCLUDE_UNESCOFREEZING_HPP
#define SRC_INCLUDE_UNESCOFREEZING_HPP

#include <cmath>

#include "IFreezingPoint.hpp"

namespace Nextsim {

class UnescoFreezing : public IFreezingPoint {
    inline double operator()(double sss) const override
    {
        // Fofonoff and Millard, Unesco technical papers in marine science 44, (1983)
        const double a0 = -0.0575;
        const double a1 = +1.710523e-3;
        const double a2 = -2.154996e-4;
        const double b = -7.53e-4;
        const double p0 = 0; // Zero hydrostatic pressure

        return sss * (a0 + a1 * std::sqrt(sss) + a2 * sss) + b * p0;
    }
};
}

#endif /* SRC_INCLUDE_UNESCOFREEZING_HPP */