Detailed Description

Compute the freezing point potential temperature [degC] from salinity [PSU] and pressure [Pa] using the expression from Millero (1978) (and in appendix A of Gill 1982), but with the of the pressure dependence changed from 7.53e-8 to 7.75e-8 to make this an expression for potential temperature (not in situ temperature), using a value that is correct at the freezing point at 35 PSU and 5e6 Pa (500 dbar).

Definition at line 27 of file MOM_TFreeze.F90.

Private functions
subroutine	calculate_tfreeze_millero_scalar (S, pres, T_Fr)
	This subroutine computes the freezing point potential temperature [degC] from salinity [ppt], and pressure [Pa] using the expression from Millero (1978) (and in appendix A of Gill 1982), but with the of the pressure dependence changed from 7.53e-8 to 7.75e-8 to make this an expression for potential temperature (not in situ temperature), using a value that is correct at the freezing point at 35 PSU and 5e6 Pa (500 dbar). More...

subroutine	calculate_tfreeze_millero_array (S, pres, T_Fr, start, npts)
	This subroutine computes the freezing point potential temperature [degC] from salinity [ppt], and pressure [Pa] using the expression from Millero (1978) (and in appendix A of Gill 1982), but with the of the pressure dependence changed from 7.53e-8 to 7.75e-8 to make this an expression for potential temperature (not in situ temperature), using a value that is correct at the freezing point at 35 PSU and 5e6 Pa (500 dbar). More...

Functions and subroutines

◆ calculate_tfreeze_millero_array()

subroutine mom_tfreeze::calculate_tfreeze_millero::calculate_tfreeze_millero_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 potential temperature [degC] from salinity [ppt], and pressure [Pa] using the expression from Millero (1978) (and in appendix A of Gill 1982), but with the of the pressure dependence changed from 7.53e-8 to 7.75e-8 to make this an expression for potential temperature (not in situ temperature), using a value that is correct at the freezing point at 35 PSU and 5e6 Pa (500 dbar).

Parameters

[in]	s	Salinity [PSU].
[in]	pres	Pressure [Pa].
[out]	t_fr	Freezing point potential temperature [degC].
[in]	start	The starting point in the arrays.
[in]	npts	The number of values to calculate.

Definition at line 106 of file MOM_TFreeze.F90.

   real,  dimension(:), intent(in)  :: S     !< Salinity [PSU].
   real,  dimension(:), intent(in)  :: pres  !< Pressure [Pa].
   real,  dimension(:), intent(out) :: T_Fr  !< Freezing point potential 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 :: cS1 = -0.0575, cs3_2 = 1.710523e-3, cs2 = -2.154996e-4
   real, parameter :: dTFr_dp = -7.75e-8
   integer :: j
  
   do j=start,start+npts-1
     t_fr(j) = s(j)*(cs1 + (cs3_2 * sqrt(max(s(j),0.0)) + cs2 * s(j))) + &
               dtfr_dp*pres(j)
   enddo
  

◆ calculate_tfreeze_millero_scalar()

subroutine mom_tfreeze::calculate_tfreeze_millero::calculate_tfreeze_millero_scalar	(	real, intent(in)	S,
		real, intent(in)	pres,
		real, intent(out)	T_Fr
	)