FTN95

Question

我使用的是使用Silverfrost Plato的FORTRAN 95。在这里，我试图创建一个条件语句来比较双精度变量和双精度值。不知何故编译器不理解条件。

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

我无法自行解决此问题。任何帮助或相关文档的链接将不胜感激。