anchor.diff:  
      ErrStat    = ErrID_None
      ErrMsg     = ""
  -   StifK(:,:) = 0.0D0
  -   RHS(:)     = 0.0D0
  +   StifK(:,:) = 0.0_BDKi
  +   RHS(:)     = 0.0_BDKi
   
      dof_elem = dof_node * node_elem
      rot_elem = (dof_node/2) * node_elem
  
positive.code: 
  END SUBROUTINE BD_InputGlobalLocal
  !-----------------------------------------------------------------------------------------------------------------------------------
  SUBROUTINE BD_CalcIC( u, p, x,  ErrStat, ErrMsg)
  !----------------------------------------------------------------------------
  ! This subroutine computes the initial states
  ! Rigid body assumption is used in initialization of the states.
  
negative.code: 
     DO i=1, node_total
         temp_id   = (i - 1) * dof_node
         rotf_temp = 0.0_BDKi
         roti_temp = 0.0_BDKi
         rot_temp  = 0.0_BDKi
         DO j=1,3
             uf(temp_id+j) = uf(temp_id+j) + uinc(temp_id+j)
             rotf_temp(j)  = uf(temp_id+3+j)
  