anchor.diff:  c
           deallocate( rbForce )
           deallocate( rbTorque )
  +c
           deallocate( rbDisp )
  -c
           deallocate( rbForceOld )
           deallocate( rbTorqueOld )
           deallocate( rbVelOld )
  
positive.code: c
          allocate( rbForce(numrbs, 3)  )
          allocate( rbTorque(numrbs, 3) )
          allocate( rbDisp(numrbs, 3) )
  c
          allocate( rbForceOld(numrbs, 3)  )
          allocate( rbTorqueOld(numrbs, 3)  )
          allocate( rbVelOld(numrbs, 3)  )
          allocate( rbAccOld(numrbs, 3)  )
  c
          return
          end
  c
  c----------------------------------------------------------------------
  c
          subroutine init_rbForce ( )
  c
          use rigidBodyForce
  c
          rbForce  = 0.0
          rbTorque = 0.0
          rbDisp = 0.0
  c
          rbForceOld = 0.0
          rbTorqueOld = 0.0
          rbVelOld = 0.0
          rbAccOld = 0.0
  c
          return
          end
  
negative.code:           integer, allocatable :: rbIndex(:)
            real*8, allocatable  :: rbForce(:,:)
            real*8, allocatable  :: rbTorque(:,:)
  c
            real*8, allocatable  :: rbDisp(:,:)
            real*8, allocatable  :: rbForceOld(:,:)
            real*8, allocatable  :: rbTorqueOld(:,:)
            real*8, allocatable  :: rbVelOld(:,:)
  