anchor.diff:          allocate( rbTorque(numrbs, 3) )
   c
           allocate( rbDisp(numrbs, 3) )
  +        allocate( rbTotalDisp(numrbs, 3) )
           allocate( rbForceOld(numrbs, 3)  )
           allocate( rbTorqueOld(numrbs, 3)  )
           allocate( rbVelOld(numrbs, 3)  )
  
positive.code:         rbTorque = 0.0
  c
          rbDisp = 0.0
          rbForceOld = 0.0
          rbTorqueOld = 0.0
          rbVelOld = 0.0
  
negative.code:           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(:,:)
  