Linux for ansys中的分段错误

我正在使用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,MaxEleNoc 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)' )  10000end program

输出为：

Element_   1Element_  10Element_ 100Element_1000Element_****

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

Element_0001Element_0010Element_0100Element_1000Element_****

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

最后,如果您希望数字直接在斜杠后开始而没有前导零,则可以使用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)') filenameend program

结果是

Element_1                       Element_10                      Element_100  

来源：https://www.icode9.com/content-3-505551.html