Question

我使用标记f进行错误检查。当我想再次检查时，Fortran（或者gfortran）不会让我重新分配它的价值。

integer, dimension(:,:), allocatable :: A
integer :: f, n        

write (*, *) "Give an integer n > 0. n = "

   read (*, IOSTAT=f) n

   do while(f /= 0)
      print *, "Error with input. Please try again."
      read (*, IOSTAT=f) n
   end do

   write (*, "(a, i5)") "You have entered n = ", n

   allocate(A(n), STAT=f)
   if (f /= 0) 
      print *, "Error: not enough memory for A."
   end if

注意：我认为复制粘贴可能会弄乱我的间距。

f已被声明为整数（而不是参数整数）：integer :: f。

我是Fortran的初学者，所以我很可能犯了一些不可思议的错误！

Answer 1

此错误消息令人困惑，但问题是

   if (f /= 0) 
      print *, "Error: not enough memory for A."
   end if

应该是

   if (f /= 0) then
      print *, "Error: not enough memory for A."
   end if

Answer 2

   implicit none
   integer :: f, n
   integer, dimension(:,:), allocatable :: A


   write (*, *) "Give an integer n > 0. n = "
   read (*, *, IOSTAT=f) n

   do while (f /= 0)
      print *, "Error with input. Please try again."
      read (*, IOSTAT=f) n
   end do

   write (*, "(a, i5)") "You have entered n = ", n

   allocate(A(n,n), STAT=f)
   if (f /= 0) then
      print *, "Error: not enough memory for A."
      !exit program. How do I do this?
   end if

这似乎有效。

（1）正如Vladimir F指出的那样，Fortran想要if (<condition>) then <stuff> endif。

（2）正如我在上述评论中所提到的，我应该写allocate(A(n,n), STAT=f)。

感谢您的帮助！这个答案只是为了完整性 - 弗拉基米尔实际上回答了这个问题。

＆＃34;无法分配命名常量＆＃34; （重新分配变量）

2 个答案: