anchor.diff:      return
     end subroutine
   
  -  subroutine outputdyna_resfile(hecMESH, istp, dispre, dispim, velre, velim, accre, accim, iout)
  +  subroutine outputdyna_resfile(hecMESH, time, istp, dispre, dispim, velre, velim, accre, accim, iout)
       !---- args
       type(hecmwST_local_mesh), intent(in) :: hecMESH
  +    real(kind=kreal), intent(in)   :: time
       integer(kind=kint), intent(in) :: istp
       real(kind=kreal), intent(in)   :: dispre(:) !intend (numnodeDOF)
       real(kind=kreal), intent(in)   :: dispim(:) !intend (numnodeDOF)
  
positive.code: 
    end subroutine
  
    subroutine output_resfile(hecMESH, ifreq, disp, vel, acc, iout)
      !---- args
      type(hecmwST_local_mesh), intent(in) :: hecMESH
      integer(kind=kint), intent(in)       :: ifreq
      real(kind=kreal), intent(in)         :: disp(:) !intend (numnodeDOF)
      real(kind=kreal), intent(in)         :: vel(:) !intend (numnodeDOF)
  
negative.code:     comment='frequency_result'
  
      call hecmw_result_init(hecMESH, istp, header, comment)
  
      label = "time"
      absval(1) = time
      call hecmw_result_add(3, 1, label, absval)
  
      if(iout(1) == 1) then
        label='displacement_real'
  