logical function sisnan( sin )
real, intent(in) :: sin

sisnan = (sin.ne.sin)
return
end
