mom_eos_wright::calculate_density_wright Interface Reference

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...

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.

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(inout)  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].
[in,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 110 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(inout) :: 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  )

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]	rho_ref	A reference density [kg m-3].

Definition at line 80 of file MOM_EOS_Wright.F90.

  real,           intent(in)  :: T        !< Potential temperature relative to the surface [degC].
  real,           intent(in)  :: S        !< Salinity [PSU].
  real,           intent(in)  :: pressure !< pressure [Pa].
  real,           intent(out) :: rho      !< In situ density [kg m-3].
  real, optional, intent(in)  :: rho_ref  !< A reference density [kg m-3].
 
! *====================================================================*
! *  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.                *
! *  Coded by R. Hallberg, 7/00                                        *
! *====================================================================*
 
  real, dimension(1) :: T0, S0, pressure0, rho0
 
  t0(1) = t
  s0(1) = s
  pressure0(1) = pressure
 
  call calculate_density_array_wright(t0, s0, pressure0, rho0, 1, 1, rho_ref)
  rho = rho0(1)
 

---

The documentation for this interface was generated from the following file:

• /home/cermak/src/MOM6.devrob/src/equation_of_state/MOM_EOS_Wright.F90