anchor.diff:        if(IRESULT==1) then
           write(*,   *) freq, "[Hz] : ", im, ".res"
           write(ilog,*) freq, "[Hz] : ", im, ".res"
  -        call output_resfile(hecMESH, im, disp, vel, acc, freqiout)
  +        call output_resfile(hecMESH, freq, im, disp, vel, acc, freqiout)
         end if
         if(IVISUAL==1 .and. vistype==1) then
           write(*,   *) freq, "[Hz] : ", im, ".vis"
  
positive.code:     return
    end subroutine
  
    subroutine outputdyna_resfile(hecMESH, istp, dispre, dispim, velre, velim, accre, accim, iout)
      !---- args
      type(hecmwST_local_mesh), intent(in) :: hecMESH
      integer(kind=kint), intent(in) :: istp
      real(kind=kreal), intent(in)   :: dispre(:) !intend (numnodeDOF)
      real(kind=kreal), intent(in)   :: dispim(:) !intend (numnodeDOF)
  
negative.code: 
    end subroutine
  
    subroutine output_resfile(hecMESH, freq, ifreq, disp, vel, acc, iout)
      !---- args
      type(hecmwST_local_mesh), intent(in) :: hecMESH
      real(kind=kreal), intent(in)         :: freq
      integer(kind=kint), intent(in)       :: ifreq
      real(kind=kreal), intent(in)         :: disp(:) !intend (numnodeDOF)
      real(kind=kreal), intent(in)         :: vel(:) !intend (numnodeDOF)
  