Очень долго разбираясь в проблеме вещественных чисел, и найдя экспертов в этом вопросе, я решил проблему.
Дело в том, что у вещественных чисел такая проблема, что их нельзя сравнивать на равенство, это в 99% случаев даёт абсолютно неправильный ответ.
Вместо:
IF ((DX+DY).EQ.90.0D0)
Нужно писать:
IF (ABS((DX+DY)-90.0D0) .LT. 1D-4)
А теперь по русски:
Нельзя сравнить два числа с плавающей точкой на равенство(==), нужно найти модуль разности текущего числа и числа, с которым мы сравниваем. Этот модуль разности называется абсолютной погрешностью измерения. Так вот, когда мы находим погрешность, мы должны её сравнить с малым числом. Это число нужно выбрать исходя из того, сколько знаков после запятой мы хотим сохранить.
На моём конкретном примере 89.99 складывается с 0.0099999999999 и получается 90.0000000001.
Естественно это число не равно 90.0
Но если найти абсолютную погрешность:
90.0000000001 - 90.0 = 0.0000000001
Её можно сравнить (как я это делаю в примере) с числом 0.0001.
IF (0.0000000001 < 0.0001) THEN
Наше число равно числу 90 с точность 4 знаков после запятой.
Но естественно 4 знака это мало. Исходя из официальных данных тип REAL(он же float, он же число с плавающей точкой одинарной точности), на который отводится 4 байта, имеет 7 значащих цифр. А тип DOUBLE PERCISION(он же double, REAL*8, число с плавающей точкой двойной точности), на который отводится 8 байт, имеет 15 значащих цифр.