anchor.diff:     INTEGER(IntKi)                   :: j ! Index counter
      INTEGER(IntKi)                   :: temp_id ! Index counter
      INTEGER(IntKi)                   :: allo_stat ! Allows for an error code return
  +   INTEGER(IntKi)          :: ErrStat2                     ! Temporary Error status
  +   CHARACTER(ErrMsgLen)    :: ErrMsg2                      ! Temporary Error message
  +   CHARACTER(*), PARAMETER :: RoutineName = 'BD_GenerateDynamicElementForce'
   
      dof_elem = dof_node * node_elem
      rot_elem = (dof_node/2) * node_elem
  
positive.code: !   for uuN0, uuN, and vvn (initial nodal configuration, nodal displacements, and velocity of mass) these values
  !   are then passed back to GenerateDynamicElement.
  !**********************************************************************************************************************************
  SUBROUTINE BD_ElemNodalDisp(uu,node_elem,dof_node,nelem,Nu)
  
     REAL(ReKi),INTENT(IN):: uu(:) ! Initial position vector, nodal displacements, and velocity of mass
     INTEGER(IntKi),INTENT(IN):: node_elem ! Node per element
     INTEGER(IntKi),INTENT(IN):: dof_node ! Degrees of freedom per node
     INTEGER(IntKi),INTENT(IN):: nelem ! Number of elements, being looped from 1, elem_total
     REAL(ReKi),INTENT(INOUT):: Nu(:) ! Initial nodal configuration, nodal displacements, and velocity of mass at the nodes of each element
  
     INTEGER(IntKi):: i ! Index counter
     INTEGER(IntKi):: j ! Index counter
     INTEGER(IntKi):: temp_id1 ! Counter
     INTEGER(IntKi):: temp_id2 ! Counter
  
     DO i=1,node_elem
         DO j=1,dof_node
  
negative.code:    INTEGER(IntKi)                  :: j
     INTEGER(IntKi)                  :: temp_id
     INTEGER(IntKi)                  :: allo_stat
     INTEGER(IntKi)          :: ErrStat2                     ! Temporary Error status
     CHARACTER(ErrMsgLen)    :: ErrMsg2                      ! Temporary Error message
     CHARACTER(*), PARAMETER :: RoutineName = 'BD_GenerateStaticElement'
  
     dof_elem = dof_node * node_elem
     rot_elem = (dof_node/2) * node_elem
  