source: LMDZ5/branches/LF-private/tools/Max_diff_nc_with_lib/Jumble/Numerical/averge.f90 @ 2942

module AVERGE_m

  implicit none

contains

  pure real FUNCTION AVERGE(X,Y)

    real, intent(in):: x, y

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

    AVERGE = .5 * SQRT(X) * sqrt(Y) + 0.25 * (X + Y)
    ! (the square roots should be separated to avoid underflow)

  END FUNCTION AVERGE

end module AVERGE_m