Detailed Description

Compute the in situ density of sea water (in [kg m-3]), or its anomaly with respect to a reference density, from salinity (in psu), potential temperature (in deg C), and pressure [Pa], using the expressions from Wright, 1997, J. Atmos. Ocean. Tech., 14, 735-740.

Definition at line 32 of file MOM_EOS_Wright.F90.

Private functions
subroutine	calculate_density_scalar_wright (T, S, pressure, rho, rho_ref)
	This subroutine computes the in situ density of sea water (rho in [kg m-3]) from salinity (S [PSU]), potential temperature (T [degC]), and pressure [Pa]. It uses the expression from Wright, 1997, J. Atmos. Ocean. Tech., 14, 735-740. More...

subroutine	calculate_density_array_wright (T, S, pressure, rho, start, npts, rho_ref)
	This subroutine computes the in situ density of sea water (rho in [kg m-3]) from salinity (S [PSU]), potential temperature (T [degC]), and pressure [Pa]. It uses the expression from Wright, 1997, J. Atmos. Ocean. Tech., 14, 735-740. More...

Functions and subroutines

◆ calculate_density_array_wright()

subroutine mom_eos_wright::calculate_density_wright::calculate_density_array_wright	(	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 salinity (S [PSU]), potential temperature (T [degC]), and pressure [Pa]. It uses the expression from Wright, 1997, J. Atmos. Ocean. Tech., 14, 735-740.

Parameters

[in]	t	potential temperature relative to the surface [degC].
[in]	s	salinity [PSU].
[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 111 of file MOM_EOS_Wright.F90.

   real, dimension(:), intent(in)  :: T        !< potential temperature relative to the surface [degC].
   real, dimension(:), intent(in)  :: S        !< salinity [PSU].
   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].
  
   ! Original coded by R. Hallberg, 7/00, anomaly coded in 3/18.
   ! Local variables
   real :: al0, p0, lambda
   real :: al_TS, p_TSp, lam_TS, pa_000
   integer :: j
  
   if (present(rho_ref)) pa_000 = (b0*(1.0 - a0*rho_ref) - rho_ref*c0)
   if (present(rho_ref)) then ; do j=start,start+npts-1
     al_ts = a1*t(j) +a2*s(j)
     al0 = a0 + al_ts
     p_tsp = pressure(j) + (b4*s(j) + t(j) * (b1 + (t(j)*(b2 + b3*t(j)) + b5*s(j))))
     lam_ts = c4*s(j) + t(j) * (c1 + (t(j)*(c2 + c3*t(j)) + c5*s(j)))
  
     ! The following two expressions are mathematically equivalent.
     ! rho(j) = (b0 + p0_TSp) / ((c0 + lam_TS) + al0*(b0 + p0_TSp)) - rho_ref
     rho(j) = (pa_000 + (p_tsp - rho_ref*(p_tsp*al0 + (b0*al_ts + lam_ts)))) / &
              ( (c0 + lam_ts) + al0*(b0 + p_tsp) )
   enddo ; else ; do j=start,start+npts-1
     al0 = (a0 + a1*t(j)) +a2*s(j)
     p0 = (b0 + b4*s(j)) + t(j) * (b1 + t(j)*(b2 + b3*t(j)) + b5*s(j))
     lambda = (c0 +c4*s(j)) + t(j) * (c1 + t(j)*(c2 + c3*t(j)) + c5*s(j))
     rho(j) = (pressure(j) + p0) / (lambda + al0*(pressure(j) + p0))
   enddo ; endif
  

◆ calculate_density_scalar_wright()

subroutine mom_eos_wright::calculate_density_wright::calculate_density_scalar_wright	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	rho,
		real, intent(in), optional	rho_ref
	)