我正在使用Fortran(Intel Composer(2011.1.107))编写用于ANSYS14.5.7的usersubroutine
我已经编辑了代码,将一些数据写入外部顺序文件,并将它们读入数组以进行计算,如下所示：

c Writes the array elements into the relevent elementID file and Integration point
   Do ElementNo=1,MaxEleNo
c writes the file name for the relevent element
   write (filename1, '( "Element_", I4)' )  ElementNo 
c opens the relevent file for the element data  
   OPEN(unit=ElementNo,status='unknown',ACCESS='APPEND'
 & ,file=filename1)    
   Write(ElementNo,fmt='(*(D))')(sthistory(ElementNo,:))
   close (ElementNo)
  end do
  sthistory=ZERO
  else
  endif

运行此命令时,我从HPC系统收到此错误：

Lyra: Ansys (v14.5.7) loaded.
Lyra: Intel Composer (2011.1.107) module loaded.
/pkg/suse11/ansys/v145/ansys/bin/ansys145: line 817: 50102 Segmentation fault           /pkg/suse11/ansys/v145/ansys/bin/linx64/ansys.e145 -np 4

但是,当我在子例程中使用带有save属性的静态数组运行时,代码没有问题.但这不足以将所有数据保存在我的计算中.
有人可以帮我找到可能的问题吗？

解决方法:

低单元号通常保留给“特殊单元”,例如STDOUT,STDERR,STDIN.不要使用这些(除非您知道自己在做什么),否则可能会发生意外情况.我有点意识到单位编号有一些上限,但是目前找不到参考.

因此,解决问题的最简单方法是在单元上添加偏移量(这又会导致大型数组出现问题),或者如果编译器支持,则使用newunit =.
但是,由于您在循环主体的末尾关闭了文件,因此不使用像1234这样的固定数字吗？

但是您的代码还有更多问题：该行

write (filename1, '( "Element_", I4)' )  ElementNo

将导致问题(对于大多数编译器).

考虑以下简单程序：

program test
  write (*, '( "Element_", I4)' )  1
  write (*, '( "Element_", I4)' )  10
  write (*, '( "Element_", I4)' )  100
  write (*, '( "Element_", I4)' )  1000
  write (*, '( "Element_", I4)' )  10000
end program

输出为：

Element_   1
Element_  10
Element_ 100
Element_1000
Element_****

导致文件名包含空格.这可能会导致错误！您可以做的是通过使用'((“ Element_”,I4.4)’将格式说明符更改为固定长度,这将为您提供：

Element_0001
Element_0010
Element_0100
Element_1000
Element_****

您会看到四位数仍然太小而无法容纳较大的元素,但不再有空格.

最后,如果您希望数字直接在斜杠后开始而没有前导零,则可以使用adjustl()和trim()的组合：

program test
  character(len=32)  :: filename

  write (filename, '(I4)')  1
  filename = "Element_" // adjustl(trim(filename))
  write(*,'(a)') filename
  write (filename, '(I4)')  10
  filename = "Element_" // adjustl(trim(filename))
  write(*,'(a)') filename
  write (filename, '(I4)')  100
  filename = "Element_" // adjustl(trim(filename))
  write(*,'(a)') filename
end program

结果是

Element_1                       
Element_10                      
Element_100