anchor.diff:      character(len=HECMW_HEADER_LEN) :: header
       character(len=HECMW_MSG_LEN)    :: comment
       character(len=HECMW_NAME_LEN)   :: label, nameid
  +    real(kind=kreal)                :: freqval(1)
       !---- body
   
       nameid='fstrRES'
       header='*fstrresult'
       comment='frequency_result'
       call hecmw_result_init(hecMESH, ifreq, header, comment)
  +
  +    label = "frequency"
  +    freqval(1) = freq
  +    call hecmw_result_add(3, 1, label, freqval)
  +
       if(iout(1) == 1) then
         label='displacement'
         call hecmw_result_add(1, 3, label, disp) !mode=node, ndof=3
  
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:       if(IRESULT==1) then
          write(*,   *) "time=", time, " : ", im, ".res"
          write(ilog,*) "time=", time, " : ", im, ".res"
          call outputdyna_resfile(hecMESH, time, im, dvaRe, dvaIm, velRe, velIm, accRe, accIm, freqiout)
        end if
        if(IVISUAL==1 .and. vistype==2) then
          write(*,   *) "time=", time, " : ", im, ".vis"
  