Fortran Coder

标题: 新人代码求大佬指正Unexpected data declaration statement at [打印本页]

作者: Ice-Tea 时间: 2018-9-19 20:21
标题: 新人代码求大佬指正Unexpected data declaration statement at
刚刚学写fortran，读取一个三维周期图形的坐标后，以每一个坐标为原点，做周期小球，若小球包裹了4个点（不含小球中心点），则输出小球及四个点的坐标代码如下，不知为什么一直报错，用的是linux的f90

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

program up

implicit none



!read the file......................................

real :: x(2250),y(2250),z(2250)

integer :: i

open(1,file='XYZ1.txt')

loopi: do i = 1,2250

read(1,*) x(i),y(i),z(i)

end do loopi

close(1)



!........draw a ball and get 4 points..............................................

integer :: a(4)

integer :: j

integer :: n

integer :: k

integer :: l

integer :: m(2250,4)

loopj: do j = 1,2250

k = 0

loopn: do n = 1,2250

if (j == n) then

cycle

end if

if ((x(n)-x(j))**2+(y(n)-y(j))**2+(z(n)-z(j))**2 < 1.8**2) then

k = k + 1

a(k) = n

if (k == 4) then

loopl: do l = 1,4

m(j,l) = a(l)

print *,m(j,l)

end do loopl

end if

end if

end do loopn

end do loopj

end program up

报错如下：

integer :: a(4)
            1
Error: Unexpected data declaration statement at (1)
5.f90:15.12:

integer :: j
         1
Error: Unexpected data declaration statement at (1)
5.f90:16.12:

integer :: n
         1
Error: Unexpected data declaration statement at (1)
5.f90:17.12:

integer :: k
         1
Error: Unexpected data declaration statement at (1)
5.f90:18.12:

integer :: l
         1
Error: Unexpected data declaration statement at (1)
5.f90:19.20:

integer :: m(2250,4)
                  1
Error: Unexpected data declaration statement at (1)
5.f90:28.8:

a(k) = n
      1
Error: Unexpected STATEMENT FUNCTION statement at (1)
5.f90:31.13:

m(j,l) = a(l)
         1
Error: Unexpected STATEMENT FUNCTION statement at (1)
5.f90:20.4:

do j = 1,2250
1
Error: Symbol 'j' at (1) has no IMPLICIT type
5.f90:21.1:

k = 0
1
Error: Symbol 'k' at (1) has no IMPLICIT type
5.f90:30.4:

do l = 1,4
1
Error: Symbol 'l' at (1) has no IMPLICIT type
5.f90:22.4:

do n = 1,2250
1
Error: Symbol 'n' at (1) has no IMPLICIT type

XYZ1.txt

64.5 KB, 下载次数: 0

作者: fcode 时间: 2018-9-19 20:31
把所有定义语句放到前面即可。

作者: fcode 时间: 2018-9-19 20:33
忠告
1. 排版很重要，一定要做。
2. 对于代码中大量出现的同一个常数，可以定义为 Parameter，减少今后的修改量。比如你要把 2250 改成 2230，会很麻烦。

如下

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

program up

  implicit none

  Integer , parameter :: N = 2250 , S = 4

  real :: x(N),y(N),z(N)

  integer :: i , j , k , l

  integer :: a(S)

  integer :: m(N,S)

  !read the file......................................

  open(1,file='XYZ1.txt')

  do i = 1,N

    read(1,*) x(i),y(i),z(i)

  end do

  close(1)

  !........draw a ball and get 4 points..............................................

  do j = 1,N

    k = 0

    do i = 1,N

      if (j == i) cycle

      if ((x(i)-x(j))**2+(y(i)-y(j))**2+(z(i)-z(j))**2 < 1.8**2) then

        k = k + 1

        a(k) = i

        if (k == 4) then

          do l = 1,S

            m(j,l) = a(l)

            print *,m(j,l)

          end do

        end if

      end if

    end do

  end do

end program up

作者: Ice-Tea 时间: 2018-9-19 20:42

fcode 发表于 2018-9-19 20:31
把所有定义语句放到前面即可。

谢谢大佬是我学术不精刚开始学以为没用到的定义可以往后放一个小问题困扰了一天

，一直在网上搜索问题盲目学，终于找到组织了决定长期混论坛了

作者: Ice-Tea 时间: 2018-9-19 20:51

fcode 发表于 2018-9-19 20:31
把所有定义语句放到前面即可。

谢谢，我会改正的. 同时这里我有一个问题想请教下师兄说这里建议我用动态数组来赋值，增加程序的适用性，但导师说动态数组会增加程序负荷，这里静态数组就行，这段程序只是我要计算东西的一部分，还没完全写完，想知道如果改成动态数组会怎样

作者: fcode 时间: 2018-9-19 21:48
程序负荷这个问题可以忽略不计。

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)