[求助] abaqus调用fortran子程序读写txt时出现steptime重复，为什么？

C
C User subroutine VDLOAD
            subroutine vdload (
C Read only -
     *     nblock, ndim, stepTime, totalTime, amplitude,
     *     curCoords, velocity, dircos, jltyp, sname,
C Write only -
     *     value )
C
      include 'vaba_param.inc'
      parameter ( const = 1.d5 ) 
C
      dimension curCoords(nblock,ndim), velocity(nblock,ndim),
     *     dircos(nblock,ndim,ndim), value(nblock)
*
      character*80 sname
*
      integer ::i
c      character(len=20)::filename="duxie.txt"
      REAL::ta
      real ::ta1
      real::tz

      open(93,FILE="D:\abaqus2016\SIMULIA\duxie.txt",status="old")      
      do
      read(93,*,iostat = ios) ta,ta1
      if ( ios == -1 ) exit 
       end do
       backspace(93)
       backspace(93)

              ta=ta+2*stepTime
              ta1=ta1+5*stepTime
              tz=0.5*ta
              tx=0.5*ta1
      write(93,*) ta,ta1
      close(93) 
      open(89,FILE="D:\abaqus2016\SIMULIA\duxiet.txt",position="append")
      write(89,*) ta,ta1,stepTime, tz
      close(89) 
      
      if( jltyp .eq. 2 ) then
        do k = 1, nblock        
        value(k) =tz
        end do
      end if
        if( jltyp .eq. 1 ) then
        do k = 1, nblock        
        value(k) =tx
        end do
        end if
    
      return
      end

            subroutine vdload (
C Read only -
     *     nblock, ndim, stepTime, totalTime, amplitude,
     *     curCoords, velocity, dircos, jltyp, sname,
C Write only -
     *     value )
C
      include 'vaba_param.inc'
      parameter ( const = 1.d5 ) 
C
      dimension curCoords(nblock,ndim), velocity(nblock,ndim),
     *     dircos(nblock,ndim,ndim), value(nblock)
*
      character*80 sname
*
      REAL::ta
      open(93,FILE="D:\abaqus2016\SIMULIA\duxie.txt",status="old")      
       read(93,*) ta
       rewind(93)
       
       ta=ta+4*stepTime
       
      write(93,*) ta
      close(93) 
      open(89,FILE="D:\abaqus2016\SIMULIA\duxiet.txt",position="append")
      write(89,*) ta,stepTime
      close(89) 
      
      if( jltyp .eq. 2 ) then
        do k = 1, nblock        
        value(k) =-ta/4    ！体积是4， value(k)表示体力，体力X体积等于施加的反力（ value(k)是给到abaqus里的体力，abaqus能识别，cae里面后处理得到的反力等于 value(k)X4）
        end do
      end if    
      return
      end