REALTYPE function eqstate1(S,T,p,g,rho_0)DESCRIPTION:
Calculates the in-situ buoyancy according to the selected method. S is salinity in psu, T is potential temperature in C (ITS-90), p is gauge pressure (absolute pressure - 10.1325 bar), g is the gravitational acceleration in ms and rho_0 the reference density in kgm. eqstate1 is the in-situ-density in kgm. For eq_state_method=1, the UNESCO equation of state is used, for eq_state_method=2, the Jackett et al. (2005) equation of state is used. Here, some care is needed, since the UNESCO equation used bar for pressure and the Jackett et al. (2005) uses dbar for pressure. For values of eq_state_method ranging from 1 to 4, one of the following methods will be used.
IMPLICIT NONEINPUT PARAMETERS:
REALTYPE,intent(in) :: S,T,p REALTYPE,optional,intent(in) :: g,rho_0REVISION HISTORY:
Original author(s): Hans Burchard & Karsten Bolding
Karsten Bolding 2012-12-28