Бесконечность в Фортране

Какой самый безопасный способ установить для переменной значение + Infinity в Фортране? В настоящий момент я использую:

program test
  implicit none
  print *,infinity()
contains
  real function infinity()
    implicit none
    real :: x
    x = huge(1.)
    infinity = x + x
  end function infinity
end program test

, но мне интересно, есть ли способ лучше?