Detailed Description

Compute the in situ density of sea water ([kg m-3]), or its anomaly with respect to a reference density, from absolute salinity (g/kg), conservative temperature (in deg C), and pressure [Pa], using the TEOS10 expressions.

Definition at line 28 of file MOM_EOS_TEOS10.F90.

Private functions
subroutine	calculate_density_scalar_teos10 (T, S, pressure, rho, rho_ref)
	This subroutine computes the in situ density of sea water (rho in [kg m-3]) from absolute salinity (S [g kg-1]), conservative temperature (T [degC]), and pressure [Pa]. It uses the expression from the TEOS10 website. More...

subroutine	calculate_density_array_teos10 (T, S, pressure, rho, start, npts, rho_ref)
	This subroutine computes the in situ density of sea water (rho in [kg m-3]) from absolute salinity (S [g kg-1]), conservative temperature (T [degC]), and pressure [Pa]. It uses the expression from the TEOS10 website. More...

Functions and subroutines

◆ calculate_density_array_teos10()

subroutine mom_eos_teos10::calculate_density_teos10::calculate_density_array_teos10	(	real, dimension(:), intent(in)	T,
		real, dimension(:), intent(in)	S,
		real, dimension(:), intent(in)	pressure,
		real, dimension(:), intent(out)	rho,
		integer, intent(in)	start,
		integer, intent(in)	npts,
		real, intent(in), optional	rho_ref
	)

private

This subroutine computes the in situ density of sea water (rho in [kg m-3]) from absolute salinity (S [g kg-1]), conservative temperature (T [degC]), and pressure [Pa]. It uses the expression from the TEOS10 website.

Parameters

[in]	t	Conservative temperature [degC].
[in]	s	Absolute salinity [g kg-1]
[in]	pressure	pressure [Pa].
[out]	rho	in situ density [kg m-3].
[in]	start	the starting point in the arrays.
[in]	npts	the number of values to calculate.
[in]	rho_ref	A reference density [kg m-3].

Definition at line 84 of file MOM_EOS_TEOS10.F90.

   real, dimension(:), intent(in)  :: T        !< Conservative temperature [degC].
   real, dimension(:), intent(in)  :: S        !< Absolute salinity [g kg-1]
   real, dimension(:), intent(in)  :: pressure !< pressure [Pa].
   real, dimension(:), intent(out) :: rho      !< in situ density [kg m-3].
   integer,            intent(in)  :: start    !< the starting point in the arrays.
   integer,            intent(in)  :: npts     !< the number of values to calculate.
   real,     optional, intent(in)  :: rho_ref  !< A reference density [kg m-3].
  
   ! Local variables
   real :: zs, zt, zp
   integer :: j
  
   do j=start,start+npts-1
     !Conversions
     zs = s(j) !gsw_sr_from_sp(S(j))       !Convert practical salinity to absolute salinity
     zt = t(j) !gsw_ct_from_pt(S(j),T(j))  !Convert potantial temp to conservative temp
     zp = pressure(j)* pa2db         !Convert pressure from Pascal to decibar
  
     if (s(j) < -1.0e-10) then !Can we assume safely that this is a missing value?
       rho(j) = 1000.0
     else
       rho(j) = gsw_rho(zs,zt,zp)
     endif
     if (present(rho_ref)) rho(j) = rho(j) - rho_ref
   enddo

◆ calculate_density_scalar_teos10()

subroutine mom_eos_teos10::calculate_density_teos10::calculate_density_scalar_teos10	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	rho,
		real, intent(in), optional	rho_ref
	)