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:     velIm(:) = 0.0D0
      
      do imode=1, freqData%numMode
        velRe(:) = velRe(:) - freqData%eigOmega(imode) * bjIm(imode) * freqData%eigVector(:,imode)
        velIm(:) = velIm(:) + freqData%eigOmega(imode) * bjRe(imode) * freqData%eigVector(:,imode)
      end do
    end subroutine
    
    subroutine calcAccVector(freqData, bjRe, bjIm, accRe, accIm)
    !---- args
      type(fstr_freqanalysis_data), intent(in) :: freqData
      real(kind=kreal),             intent(in) :: bjRe(:)       !intend (numMode)
      real(kind=kreal),             intent(in) :: bjIm(:)       !intend (numMode)
      real(kind=kreal),          intent(inout) :: accRe(:)      !intend (numNodeDOF)
  
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
  