Ошибка при сравнении между переменной двойной точности и значением двойной точности в FTN95

Я использую FORTRAN 95 с Silverfrost Plato. Здесь я пытаюсь создать условный оператор для сравнения переменной двойной точности со значением двойной точности. Почему-то компилятор не понимает условия.

program db
implicit none
integer :: j,n
double precision :: t,t_max,dt
n = 20
dt = 0.000020d0
t_max = 0.00632d0
t = 0.0d0
while (t .LE. t_max) do
    if ( t == 0.00158d0) then
        do j = 0,n
            print  *, j
        end do
    end if
print *, t 
t = t + dt    
end do
end program db

Я не могу решить эту проблему самостоятельно. Буду очень признателен за любую помощь или ссылки на соответствующую документацию.