Detailed Description

Compute the in situ specific volume of sea water (in [m3 kg-1]), or an anomaly with respect to a reference specific volume, from salinity [PSU], potential temperature [degC], and pressure [Pa], using the UNESCO (1981) equation of state.

Definition at line 29 of file MOM_EOS_UNESCO.F90.

Private functions
subroutine	calculate_spec_vol_scalar_unesco (T, S, pressure, specvol, spv_ref)
	This subroutine computes the in situ specific volume of sea water (specvol in [m3 kg-1]) from salinity (S [PSU]), potential temperature (T [degC]) and pressure [Pa], using the UNESCO (1981) equation of state. If spv_ref is present, specvol is an anomaly from spv_ref. More...

subroutine	calculate_spec_vol_array_unesco (T, S, pressure, specvol, start, npts, spv_ref)
	This subroutine computes the in situ specific volume of sea water (specvol in [m3 kg-1]) from salinity (S [PSU]), potential temperature (T [degC]) and pressure [Pa], using the UNESCO (1981) equation of state. If spv_ref is present, specvol is an anomaly from spv_ref. More...

Functions and subroutines

◆ calculate_spec_vol_array_unesco()

subroutine mom_eos_unesco::calculate_spec_vol_unesco::calculate_spec_vol_array_unesco	(	real, dimension(:), intent(in)	T,
		real, dimension(:), intent(in)	S,
		real, dimension(:), intent(in)	pressure,
		real, dimension(:), intent(out)	specvol,
		integer, intent(in)	start,
		integer, intent(in)	npts,
		real, intent(in), optional	spv_ref
	)

private

This subroutine computes the in situ specific volume of sea water (specvol in [m3 kg-1]) from salinity (S [PSU]), potential temperature (T [degC]) and pressure [Pa], using the UNESCO (1981) equation of state. If spv_ref is present, specvol is an anomaly from spv_ref.

Parameters

[in]	t	potential temperature relative to the surface [degC].
[in]	s	salinity [PSU].
[in]	pressure	pressure [Pa].
[out]	specvol	in situ specific volume [m3 kg-1].
[in]	start	the starting point in the arrays.
[in]	npts	the number of values to calculate.
[in]	spv_ref	A reference specific volume [m3 kg-1].

Definition at line 159 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) :: specvol  !< in situ specific volume [m3 kg-1].
   integer,            intent(in)  :: start    !< the starting point in the arrays.
   integer,            intent(in)  :: npts     !< the number of values to calculate.
   real,     optional, intent(in)  :: spv_ref  !< A reference specific volume [m3 kg-1].
  
   ! 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 :: 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?
       specvol(j) = 0.001
       if (present(spv_ref)) specvol(j) = 0.001 - spv_ref
       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) ).
  
     rho0 = r00 + 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
  
 !  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(spv_ref)) then
       specvol(j) = (ks*(1.0 - (rho0*spv_ref)) - p1) / (rho0*ks)
     else
       specvol(j) = (ks - p1) / (rho0*ks)
     endif
   enddo

◆ calculate_spec_vol_scalar_unesco()

subroutine mom_eos_unesco::calculate_spec_vol_unesco::calculate_spec_vol_scalar_unesco	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	specvol,
		real, intent(in), optional	spv_ref
	)