Fortran：从读取的文件分配内存

Question

我尝试从文件中读取大量数据。 文件的结构是一致的，但是重复的次数却不一致。 输入文件的第一行告诉文件重复的次数。 我想分配许多与输入文件重复长度相同的向量。

有没有一种优雅的方法可以解决此问题，而无需分配每个变量？ 我当前的子例程不起作用（我是Fortran新手）。

subroutine geometry(N, Aw, Vp, Vtot,Ds, Spos, DIR)
implicit none

  integer, intent(out) :: N
  real, dimension(:), allocatable, intent(out):: Aw, Vp, Vtot, Ym, Zm, Zt, Ds
  real, dimension(:), allocatable, intent(out) :: Spos
  integer, dimension(:), allocatable, intent(out) :: Ns, DIR
  real, dimension(:), allocatable :: Hp
  real, parameter :: pi = 4.0*atan(1.0) 
  real :: Vs, Ds, Ns
  integer :: i

  open(unit = 11, file = "semisub.dat")
  read(11,*) N

  allocate(Aw(N))
  allocate(Vp(N))
  allocate(Vtot(N))
  allocate(Ym(N))
  allocate(Zm(N))
  allocate(Zt(N))
  allocate(Ds(N))
  allocate(Spos(N,4))
  allocate(Ns(N))
  allocate(DIR(N))
  allocate(Hp(N))

  do i = 1,N
     read(11,*) Bp, Hp, Lp, Ym(i), Zm(i), DIR(i)
     read(11,*) Ns(i), Ds(i)
     read(11,*) Spos(i,1:Ns(i))
  end do

Aw(:) = 2 * Ds(:) * Ns(:)   
Vp(:) = 2 * Lp(:) * Bp(:) * Hp(:)   
Vs(:) = 2 * Ns(:) * pi * (Ds/2)**2 * (-Zt(:))


 Vtot(:) = Vp(:) + Vs(:)


end subroutine geometry

错误日志：

geometry.f03:7.54:

  integer, dimension(:),allocatable, intent(out) :: Ns, DIR
                                                      1
Error: Symbol at (1) is not a DUMMY variable
geometry.f03:5.65:

  real, dimension(:), allocatable, intent(out):: Aw, Vp, Vtot, Ym, Zm, Zt, Ds
                                                                 1
Error: Symbol at (1) is not a DUMMY variable
geometry.f03:5.69:

  real, dimension(:), allocatable, intent(out):: Aw, Vp, Vtot, Ym, Zm, Zt, Ds
                                                                     1
Error: Symbol at (1) is not a DUMMY variable
geometry.f03:5.73:

  real, dimension(:), allocatable, intent(out):: Aw, Vp, Vtot, Ym, Zm, Zt, Ds
                                                                         1
Error: Symbol at (1) is not a DUMMY variable
geometry.f03:23.15:

  allocate(Spos(N,4))
               1
Error: Rank mismatch in array reference at (1) (2/1)
geometry.f03:32.20:

     read(11,*) Spos(i,1:Ns(i))
                    1
Error: Rank mismatch in array reference at (1) (2/1)

Answer 1

阅读之前，您必须分配要填充的每个数组。 您的方法还可以，但有些错误不会影响逻辑。

唯一的其他选择是通过分配使用重新分配：

do...
...
  read(11,*) Ns_tmp, Ds_tmp
  Ns = [Ns, Ns_tmp]
  Ds = [Ds, Ds_tmp]
...
end do...

但是如果您知道文件开头的部分数量，我就不会这样做。 它比较慢，而且过于复杂。

---

您的错误是参数列表(N, Aw, Vp, Vtot,Ds, Spos, DIR)

与您的伪参数声明不一致：

  integer, intent(out) :: N
  real, dimension(:), allocatable, intent(out):: Aw, Vp, Vtot, Ym, Zm, Zt, Ds
  real, dimension(:), allocatable, intent(out) :: Spos
  integer, dimension(:), allocatable, intent(out) :: Ns, DIR

只是解决它，它应该工作。 只需将所有参数添加到参数列表中，不要对非参数使用intent 。

另一个错误是您将Spos声明为按维，而将其用作2维。