anchor.diff:  c
           allocate( rbForce(numrbs, 3)  )
           allocate( rbTorque(numrbs, 3) )
  +c
           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
  +c
           rbDisp = 0.0
  -c
           rbForceOld = 0.0
           rbTorqueOld = 0.0
           rbVelOld = 0.0
           rbAccOld = 0.0
  +c
  +        return
  +        end
  +c
  +c----------------------------------------------------------------------
  +c
  +        subroutine init_rbForce ( )
  +c
  +        use rigidBodyForce
  +c
  +        rbForce  = 0.0
  +        rbTorque = 0.0
   c
           return
           end
  
positive.code: c
          deallocate( rbForce )
          deallocate( rbTorque )
          deallocate( rbDisp )
  c
          deallocate( rbForceOld )
          deallocate( rbTorqueOld )
          deallocate( rbVelOld )
  
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(:,:)
  