MOM_intrinsic_functions.F90

!> A module with intrinsic functions that are used by MOM but are not supported
!!  by some compilers.
module mom_intrinsic_functions
 
! This file is part of MOM6. See LICENSE.md for the license.
 
implicit none ; private
 
public :: invcosh
 
contains
 
!> Evaluate the inverse cosh, either using a math library or an
!! equivalent expression
function invcosh(x)
  real, intent(in) :: x !< The argument of the inverse of cosh.  NaNs will
                        !! occur if x<1, but there is no error checking
  real :: invcosh
 
#ifdef __INTEL_COMPILER
  invcosh = acosh(x)
#else
  invcosh = log(x+sqrt(x*x-1))
#endif
 
end function invcosh
 
end module mom_intrinsic_functions