anchor.diff:            real*8, allocatable  :: rbTorque(:,:)
   c
             real*8, allocatable  :: rbDisp(:,:)
  +          real*8, allocatable  :: rbTotalDisp(:,:)
             real*8, allocatable  :: rbForceOld(:,:)
             real*8, allocatable  :: rbTorqueOld(:,:)
             real*8, allocatable  :: rbVelOld(:,:)
  
positive.code:         deallocate( rbTorque )
  c
          deallocate( rbDisp )
          deallocate( rbForceOld )
          deallocate( rbTorqueOld )
          deallocate( rbVelOld )
  
negative.code: c
  c----------------------------------------------------------------------
  c
          subroutine synchronize_rbForce ( )
  c
          use rigidBodyFlag
          use rigidBodyForce
          use core_rigid_body
  c
          include "common.h"
  c
          real*8, dimension(numrbs) :: ax, ay, az, px, py, pz, ag, sc
  c
  c.... XXX need to add rotation motion here
          do j = 1,numrbs
            if (rbsMM(j) .ne. 1)
       &      call error('rigidBodyBCElas','not support mode',rbsMM(j))
          enddo
  c
          ax = 1.0
          ay = 0.0
          az = 0.0
          px = 0.0
          py = 0.0
          pz = 0.0
          ag = 0.0
          sc = 0.0
  c
          call core_update_rbms(rbTotalDisp(:,1),
       &                        rbTotalDisp(:,2),
       &                        rbTotalDisp(:,3),
       &                        ax, ay, az,
       &                        px, py, pz,
       &                        ag, sc,
       &                        rbsTags(:), numrbs)
  c
          return
          end
  c
  c
  c----------------------------------------------------------------------
  c
  
  c
  c----------------------------------------------------------------------
  c
  
  