anchor.diff:      accIm(:) = 0.0D0
       
       do imode=1, freqData%numMode
  -      accRe(:) = accRe(:) - freqData%eigOmega(imode)**2 * bjRe(imode) * freqData%eigVector(:,imode)
  -      accIm(:) = accIm(:) - freqData%eigOmega(imode)**2 * bjIm(imode) * freqData%eigVector(:,imode)
  +      accRe(:) = accRe(:) - omega**2 * bjRe(imode) * freqData%eigVector(:,imode)
  +      accIm(:) = accIm(:) - omega**2 * bjIm(imode) * freqData%eigVector(:,imode)
       end do
     
     end subroutine 
  
positive.code:       write(ilog, *) freq, "[Hz] : ", val
        disp(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
        
        call calcVelVector(freqData, bjRe, bjIm, dvaRe, dvaIm)
        vel(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
      
        call calcAccVector(freqData, bjRe, bjIm, dvaRe, dvaIm)
        acc(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
        
        if(IRESULT==1) then
  
negative.code:     return
    end subroutine
    
    subroutine calcVelVector(freqData, omega, bjRe, bjIm, velRe, velIm)
    !---- args
      type(fstr_freqanalysis_data), intent(in) :: freqData
      real(kind=kreal),              intent(in) :: omega
      real(kind=kreal),              intent(in) :: bjRe(:)      !intend (numMode)
      real(kind=kreal),              intent(in) :: bjIm(:)      !intend (numMode)
      real(kind=kreal),           intent(inout) :: velRe(:)    !intend (numNodeDOF)
  