REALTYPE function rho_feistel(s,th,p,UNPress)DESCRIPTION:
Computes the in-situ density in (244) according to the Jackett et al. (2005) equation of state for sea water, which is based on the Gibbs potential developed by Feistel (2003). The pressure dependence can be switched on (UNPress=.true.) or off (UNPress=.false.). s is salinity in psu, th is potential temperature in C (ITS-90), p is gauge pressure (absolute pressure - 10.1325 dbar) and rho_feistel is the in-situ density in kgm. The check value is rho_feistel(20,20,1000) = 1017.728868019642 .
IMPLICIT NONEINPUT PARAMETERS:
REALTYPE, intent(in) :: s,th,p LOGICAL, intent(in) :: UNPressREVISION HISTORY:
Original author(s): Hans Burchard & Karsten Bolding