Detailed Description

For a given thermodynamic state, return the second derivatives of density with various combinations of conservative temperature, absolute salinity, and pressure, using the TEOS10 expressions.

Definition at line 47 of file MOM_EOS_TEOS10.F90.

Private functions
subroutine	calculate_density_second_derivs_scalar_teos10 (T, S, pressure, drho_dS_dS, drho_dS_dT, drho_dT_dT, drho_dS_dP, drho_dT_dP)
	Calculate the 5 second derivatives of the equation of state for scalar inputs. More...

subroutine	calculate_density_second_derivs_array_teos10 (T, S, pressure, drho_dS_dS, drho_dS_dT, drho_dT_dT, drho_dS_dP, drho_dT_dP, start, npts)
	Calculate the 5 second derivatives of the equation of state for scalar inputs. More...

Functions and subroutines

◆ calculate_density_second_derivs_array_teos10()

subroutine mom_eos_teos10::calculate_density_second_derivs_teos10::calculate_density_second_derivs_array_teos10	(	real, dimension(:), intent(in)	T,
		real, dimension(:), intent(in)	S,
		real, dimension(:), intent(in)	pressure,
		real, dimension(:), intent(out)	drho_dS_dS,
		real, dimension(:), intent(out)	drho_dS_dT,
		real, dimension(:), intent(out)	drho_dT_dT,
		real, dimension(:), intent(out)	drho_dS_dP,
		real, dimension(:), intent(out)	drho_dT_dP,
		integer, intent(in)	start,
		integer, intent(in)	npts
	)

private

Calculate the 5 second derivatives of the equation of state for scalar inputs.

Parameters

[in]	t	Conservative temperature [degC]
[in]	s	Absolute Salinity [g kg-1]
[in]	pressure	pressure [Pa].
[out]	drho_ds_ds	Partial derivative of beta with respect to S
[out]	drho_ds_dt	Partial derivative of beta with resepct to T
[out]	drho_dt_dt	Partial derivative of alpha with respect to T
[out]	drho_ds_dp	Partial derivative of beta with respect to pressure
[out]	drho_dt_dp	Partial derivative of alpha with respect to pressure
[in]	start	The starting point in the arrays.
[in]	npts	The number of values to calculate.

Definition at line 277 of file MOM_EOS_TEOS10.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)    :: drho_dS_dS !< Partial derivative of beta with respect to S
   real, dimension(:), intent(out)    :: drho_dS_dT !< Partial derivative of beta with resepct to T
   real, dimension(:), intent(out)    :: drho_dT_dT !< Partial derivative of alpha with respect to T
   real, dimension(:), intent(out)    :: drho_dS_dP !< Partial derivative of beta with respect to pressure
   real, dimension(:), intent(out)    :: drho_dT_dP !< Partial derivative of alpha with respect to pressure
   integer, intent(in)  :: start    !< The starting point in the arrays.
   integer, intent(in)  :: npts     !< The number of values to calculate.
  
   ! Local variables
   real :: zs, zt, zp
   integer :: j
  
   do j=start,start+npts-1
     !Conversions
     zs = s(j) !gsw_sr_from_sp(S)       !Convert practical salinity to absolute salinity
     zt = t(j) !gsw_ct_from_pt(S,T)  !Convert potantial temp to conservative temp
     zp = pressure(j)* pa2db         !Convert pressure from Pascal to decibar
     if (s(j) < -1.0e-10) then ; !Can we assume safely that this is a missing value?
       drho_ds_ds(j) = 0.0 ; drho_ds_dt(j) = 0.0 ; drho_dt_dt(j) = 0.0
       drho_ds_dp(j) = 0.0 ; drho_dt_dp(j) = 0.0
     else
       call gsw_rho_second_derivatives(zs, zt, zp, rho_sa_sa=drho_ds_ds(j), rho_sa_ct=drho_ds_dt(j), &
                                       rho_ct_ct=drho_dt_dt(j), rho_sa_p=drho_ds_dp(j), rho_ct_p=drho_dt_dp(j))
     endif
   enddo
  

◆ calculate_density_second_derivs_scalar_teos10()

subroutine mom_eos_teos10::calculate_density_second_derivs_teos10::calculate_density_second_derivs_scalar_teos10	(	real, intent(in)	T,
		real, intent(in)	S,
		real, intent(in)	pressure,
		real, intent(out)	drho_dS_dS,
		real, intent(out)	drho_dS_dT,
		real, intent(out)	drho_dT_dT,
		real, intent(out)	drho_dS_dP,
		real, intent(out)	drho_dT_dP
	)