I\'ve been trying to find a standards-compliant way to check for Infinite and NaN values in Fortran 90/95 but it proved harder than I thought.
for testing NaN none of the things worked eg.if testing real s2p to see if it is NaN then
if(isnan(s2p))
did not work in gfortran nor did
if(s2p.ne.s2p).
The only thing that worked was
if(.not.s2p<1.and..not.s2p>1)
though to make real sure u may want to add
if(.not.s2p<1.and..not.s2p>1.and..not.s2p==1)