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 expressions derived for use with NEMO.

Definition at line 28 of file MOM_EOS_NEMO.F90.

Private functions
subroutine	calculate_density_scalar_nemo (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 expressions derived for use with NEMO. More...

subroutine	calculate_density_array_nemo (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 expressions derived for use with NEMO. More...

Functions and subroutines

◆ calculate_density_array_nemo()

subroutine mom_eos_nemo::calculate_density_nemo::calculate_density_array_nemo	(	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 expressions derived for use with NEMO.

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 206 of file MOM_EOS_NEMO.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 :: zp, zt, zh, zs, zr0, zn, zn0, zn1, zn2, zn3, zs0
   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 potential temp to conservative temp
     zp = pressure(j)* pa2db         !Convert pressure from Pascal to decibar
  
     !The following algorithm was provided by Roquet in a private communication.
     !It is not necessarily the algorithm used in NEMO ocean!
     zp  = zp * r1_p0 !pressure
     zt  = zt * r1_t0 !temperature
     zs  = sqrt( abs( zs + rdeltas ) * r1_s0 )   ! square root salinity
  
     zn3 = eos013*zt   &
        &   + eos103*zs+eos003
  
     zn2 = (eos022*zt   &
        &   + eos112*zs+eos012)*zt   &
        &   + (eos202*zs+eos102)*zs+eos002
  
     zn1 = (((eos041*zt   &
        &   + eos131*zs+eos031)*zt   &
        &   + (eos221*zs+eos121)*zs+eos021)*zt   &
        &   + ((eos311*zs+eos211)*zs+eos111)*zs+eos011)*zt   &
        &   + (((eos401*zs+eos301)*zs+eos201)*zs+eos101)*zs+eos001
  
     zn0 = (((((eos060*zt   &
        &   + eos150*zs+eos050)*zt   &
        &   + (eos240*zs+eos140)*zs+eos040)*zt   &
        &   + ((eos330*zs+eos230)*zs+eos130)*zs+eos030)*zt   &
        &   + (((eos420*zs+eos320)*zs+eos220)*zs+eos120)*zs+eos020)*zt   &
        &   + ((((eos510*zs+eos410)*zs+eos310)*zs+eos210)*zs+eos110)*zs+eos010)*zt
  
     zs0 = (((((eos600*zs+eos500)*zs+eos400)*zs+eos300)*zs+eos200)*zs+eos100)*zs + eos000
  
     zr0 = (((((r05 * zp+r04) * zp+r03 ) * zp+r02 ) * zp+r01) * zp+r00) * zp
  
     if (present(rho_ref)) then
       zn  = ( ( zn3 * zp + zn2 ) * zp + zn1 ) * zp + (zn0 + (zs0 - rho_ref))
       rho(j) =  ( zn + zr0 ) ! density
     else
       zn  = ( ( zn3 * zp + zn2 ) * zp + zn1 ) * zp + (zn0 + zs0)
       rho(j) =  ( zn + zr0 ) ! density
     endif
  
  enddo

◆ calculate_density_scalar_nemo()

subroutine mom_eos_nemo::calculate_density_nemo::calculate_density_scalar_nemo	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	rho,
		real, intent(in), optional	rho_ref
	)