gsw_fdelta.f90

Go to the documentation of this file.

!==========================================================================
elemental function gsw_fdelta (p, long, lat)
!==========================================================================
!
! Calculates fdelta. 
!
! p      : sea pressure                                    [dbar]
! long   : longitude                                       [deg E]     
! lat    : latitude                                        [deg N]
!
! gsw_fdelta : Absolute Salinty Anomaly                    [unitless]
!--------------------------------------------------------------------------

use gsw_mod_toolbox, only : gsw_saar

use gsw_mod_error_functions, only : gsw_error_code, gsw_error_limit

use gsw_mod_kinds

implicit none

real (r8), intent(in) :: p, long, lat

real (r8) :: gsw_fdelta

real (r8) ::  saar

character (*), parameter :: func_name = "gsw_fdelta"

saar = gsw_saar(p,long,lat)

if (saar.gt.gsw_error_limit) then
   gsw_fdelta = gsw_error_code(1,func_name)
else
   gsw_fdelta = ((1.0_r8 + 0.35_r8)*saar)/(1.0_r8 - 0.35_r8*saar)
end if

return
end function

!--------------------------------------------------------------------------