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 [PSU], potential temperature [degC], and pressure [Pa], using the UNESCO (1981) equation of state.

Definition at line 22 of file MOM_EOS_UNESCO.F90.

Private functions
subroutine	calculate_density_scalar_unesco (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], using the UNESCO (1981) equation of state. More...

subroutine	calculate_density_array_unesco (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], using the UNESCO (1981) equation of state. More...

Functions and subroutines

◆ calculate_density_array_unesco()

subroutine mom_eos_unesco::calculate_density_unesco::calculate_density_array_unesco	(	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], using the UNESCO (1981) equation of state.

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 83 of file MOM_EOS_UNESCO.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].
  
   ! Local variables
   real :: t_local, t2, t3, t4, t5  ! Temperature to the 1st - 5th power [degC^n].
   real :: s_local, s32, s2         ! Salinity to the 1st, 3/2, & 2nd power [PSU^n].
   real :: p1, p2      ! Pressure (in bars) to the 1st and 2nd power [bar] and [bar2].
   real :: rho0        ! Density at 1 bar pressure [kg m-3].
   real :: sig0        ! The anomaly of rho0 from R00 [kg m-3].
   real :: ks          ! The secant bulk modulus [bar].
   integer :: j
  
   do j=start,start+npts-1
     if (s(j) < -1.0e-10) then !Can we assume safely that this is a missing value?
       rho(j) = 1000.0
       cycle
     endif
  
     p1 = pressure(j)*1.0e-5; p2 = p1*p1
     t_local = t(j); t2 = t_local*t_local; t3 = t_local*t2; t4 = t2*t2; t5 = t3*t2
     s_local = s(j); s2 = s_local*s_local; s32 = s_local*sqrt(s_local)
  
 !  Compute rho(s,theta,p=0) - (same as rho(s,t_insitu,p=0) ).
  
     sig0 = r10*t_local + r20*t2 + r30*t3 + r40*t4 + r50*t5 + &
            s_local*(r01 + r11*t_local + r21*t2 + r31*t3 + r41*t4) + &
            s32*(r032 + r132*t_local + r232*t2) + r02*s2
     rho0 = r00 + sig0
  
 !  Compute rho(s,theta,p), first calculating the secant bulk modulus.
  
     ks = s00 + s10*t_local + s20*t2 + s30*t3 + s40*t4 + s_local*(s01 + s11*t_local + s21*t2 + s31*t3) + &
          s32*(s032 + s132*t_local + s232*t2) + &
          p1*(sp00 + sp10*t_local + sp20*t2 + sp30*t3 + &
              s_local*(sp01 + sp11*t_local + sp21*t2) + sp032*s32) + &
          p2*(sp000 + sp010*t_local + sp020*t2 + s_local*(sp001 + sp011*t_local + sp021*t2))
  
     if (present(rho_ref)) then
       rho(j) = ((r00 - rho_ref)*ks + (sig0*ks + p1*rho_ref)) / (ks - p1)
     else
       rho(j) = rho0*ks / (ks - p1)
     endif
   enddo

◆ calculate_density_scalar_unesco()

subroutine mom_eos_unesco::calculate_density_unesco::calculate_density_scalar_unesco	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	rho,
		real, intent(in), optional	rho_ref
	)