[求助] Unclassifiable statement 错误

98 帖子	0 主题	0 精华

大师

F 币: 658 元
贡献: 293 点

发消息

楼主

发表于 2021-11-26 22:28:55 | 显示全部楼层

本帖最后由胡文刚于 2021-11-26 22:44 编辑

这样是你想要的吗？

[Fortran] 纯文本查看 复制代码

Program workfile
  Implicit None
  integer , parameter :: n = 1000
  real , parameter :: PI = acos(-1.0) , a = 1.42e-10
  real :: b1(2), b2(2)
  real , allocatable :: k(:,:,:)
  integer :: x , y
  allocate(k(2,n-1,n-1))
  b1 = [ 2*PI/a*sqrt(3.0) , 2*PI/a ]
  b2 = [ 2*PI/a*sqrt(3.0) ,-2*PI/a ]
  Forall(x=1:size(k,2) , y=1:size(k,3))
    k(:,x,y) = (b1*x+b2*y)/n
  End Forall
End Program workfile