[文件读写] abaqus子程序txt文件读写，并行计算报错

SUBROUTINE FILM(H,SINK,TEMP,KSTEP,KINC,TIME,NOEL,NPT,
     1 COORDS,JLTYP,FIELD,NFIELD,SNAME,NODE,AREA)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION H(2),TIME(2),COORDS(3), FIELD(NFIELD)
      CHARACTER*80 SNAME
      real a(2,241)
      integer i,j
      
      data iread /1/
      save iread
      save a
      
      if(iread.eq.1) then
      !读入温度数据
          iread = 2
          open(12,file='D:\abaqusworks\abaqus\temperature.txt')
          do i=1,241
              read(12,*) a(1,i),a(2,i)
              print*, a(1,i),a(2,i)
          end do
      end if
      close(12)
      ! 对应时间步长
      do j = 1,241
          if (TIME(2).EQ.a(1,j)) then
              SINK = a(2,j)
              exit
          end if
      end do
         
      v = 1.0                               ! wind speed (m/s)工程估算风速
      hc=(2.6*(sqrt(sqrt(ABS(SINK-TEMP))))+4*v)*3600    !β，J/(m2*h*℃）经验公式
      H(1)=hc                             !β，J/(m2*h*℃）
      H(2)=0                              !换热系数对于时间的变化率
      
      FIELD(1)=NFIELD
      FIELD(2)=H(1)
      FIELD(3)=SINK
      FIELD(4)=TEMP
      


      RETURN
      END