Fortran如何处理NaN异常？

在Fortran中，处理NaN（非数字）异常通常涉及检查和标记。由于Fortran本身并不直接提供类似于其他高级语言（如Python或Java）的内置函数来判断一个值是否为NaN，你需要借助于一些库或者自定义函数来进行判断。这里我们假设你在使用支持IEEE 754浮点数规范的编译器，如gfortran或Intel Fortran。

一种常见的做法是使用库函数`ieee754.isnan()`，虽然这在Fortran 90/95中不是标准的一部分，但在一些现代工具链中可用。例如：

use ieee_arithmetic

real :: x = some_value_where_nans_might_occur
logical :: nan_flag

nan_flag = ieee754.isnan(x)
if (nan_flag) then
   write(*,*) 'Value is NaN'
else
   do_something_with_x
end if

另外，如果没有现成的库，你也可以创建自己的函数，利用浮点数比较的特性来检查：

pure function is_nan(value) result(is_nan)
    real :: value
    logical :: is_nan
    
    is_nan = any([value .ne. value, value /= value])
end function

这样，当`value`是NaN时，上述函数将返回`.true.`。