logical function isnan(a)
        logical isnand
        double precision a
        isnan=isnand(a)
        return
        end