Detailed Description

Compute the freezing point conservative temperature [degC] from absolute salinity [g/kg] and pressure [Pa] using the TEOS10 package.

Definition at line 33 of file MOM_TFreeze.F90.

Private functions
subroutine	calculate_tfreeze_teos10_scalar (S, pres, T_Fr)
	This subroutine computes the freezing point conservative temperature [degC] from absolute salinity [g/kg], and pressure [Pa] using the TEOS10 package. More...

subroutine	calculate_tfreeze_teos10_array (S, pres, T_Fr, start, npts)
	This subroutine computes the freezing point conservative temperature [degC] from absolute salinity [g/kg], and pressure [Pa] using the TEOS10 package. More...

Functions and subroutines

◆ calculate_tfreeze_teos10_array()

subroutine mom_tfreeze::calculate_tfreeze_teos10::calculate_tfreeze_teos10_array	(	real, dimension(:), intent(in)	S,
		real, dimension(:), intent(in)	pres,
		real, dimension(:), intent(out)	T_Fr,
		integer, intent(in)	start,
		integer, intent(in)	npts
	)

private

This subroutine computes the freezing point conservative temperature [degC] from absolute salinity [g/kg], and pressure [Pa] using the TEOS10 package.

Parameters

[in]	s	absolute salinity [g/kg].
[in]	pres	pressure [Pa].
[out]	t_fr	Freezing point conservative temperature [degC].
[in]	start	the starting point in the arrays.
[in]	npts	the number of values to calculate.

Definition at line 148 of file MOM_TFreeze.F90.

   real, dimension(:), intent(in)  :: S     !< absolute salinity [g/kg].
   real, dimension(:), intent(in)  :: pres  !< pressure [Pa].
   real, dimension(:), intent(out) :: T_Fr  !< Freezing point conservative temperature [degC].
   integer,            intent(in)  :: start !< the starting point in the arrays.
   integer,            intent(in)  :: npts  !< the number of values to calculate.
  
   ! Local variables
   real, parameter :: Pa2db  = 1.e-4  ! The conversion factor from Pa to dbar.
   real :: zs,zp
   integer :: j
   ! Assume sea-water contains no dissolved air.
   real, parameter :: saturation_fraction = 0.0
  
   do j=start,start+npts-1
     !Conversions
     zs = s(j)
     zp = pres(j)* pa2db         !Convert pressure from Pascal to decibar
  
     if (s(j) < -1.0e-10) cycle !Can we assume safely that this is a missing value?
     t_fr(j) = gsw_ct_freezing_exact(zs,zp,saturation_fraction)
  enddo
  

◆ calculate_tfreeze_teos10_scalar()

subroutine mom_tfreeze::calculate_tfreeze_teos10::calculate_tfreeze_teos10_scalar	(	real, intent(in)	S,
		real, intent(in)	pres,
		real, intent(out)	T_Fr
	)