如果Fortran中存在NAN，如何引发错误？

Question

如果某些计算结果为NAN，我希望我的Fortran程序警告我。我已经天真地尝试过：

if(a==NAN) write(*,*) 'error'

和

if(a=='NAN') write(*,*) 'error'

但两者显然都失败了：第一个编译器询问NAN是什么类型，而第二个客栈我试图将数字与字符进行比较，显然编译器抱怨。

我该怎么办？