anchor.diff:      return
     end subroutine
     
  -  subroutine calcVelVector(freqData, bjRe, bjIm, velRe, velIm)
  +  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)
  
positive.code:     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)
      end do
    
    end subroutine 
  
negative.code:       write(ilog, *) freq, "[Hz] : ", val
        disp(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
        
        call calcVelVector(freqData, omega, bjRe, bjIm, dvaRe, dvaIm)
        vel(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
      
        call calcAccVector(freqData, omega, bjRe, bjIm, dvaRe, dvaIm)
        acc(:) = abs(cmplx(dvaRe(:), dvaIm(:)))
        
        if(IRESULT==1) then
  