Takuya Miyashita

NaNの有無を確認 †

isnan †

isnanを使ってNaNがあるかどうかを確認できる．
NaNが発生した時点で計算を中止したい場合などに．

do I = 1,NY
 do J = 1,NX
  if(isnan(VEL(I,J))) then
    write(*,*)'NaN発生'
    stop
  end if
 end do
end do

上記の方法だとループの途中で終わるので，NaN が複数ある場合にも最初にある I,J しかわからない．
配列すべてを確認したい場合は logical 変数を使用し，ループを抜け出した後にもう一度配列に NaN があったかどうか確かめる．

hoge = .false.
do I = 1,NY
 do J = 1,NX
  if(isnan(VEL(I,J))) then
    hoge = .true.
  end if
 end do
end do
if(hoge)then
 write(*,*) 'VEL(I,J)をファイルに出力して確認したりできる'
 stop
end if

その他 †

a が NaN の場合, a==a が常に偽になるので, isnan (Intel の ifort だけ対応？)が使えなければ.