anchor.diff:  
      DO nelem=1,elem_total
          Nuu0(:) = uuN0(:,nelem)
  -       CALL BD_ElemNodalDisp(uuN,node_elem,dof_node,nelem,Nuuu)
  +       CALL BD_ElemNodalDisp(uuN,node_elem,dof_node,nelem,Nuuu,ErrStat2,ErrMsg2)
          temp_id = (nelem-1)*ngp
          DO j=1,ngp
              EStif0_GL(1:6,1:6,j) = Stif0(1:6,1:6,temp_id+j)
  
positive.code:    elk(:,:) = 0.0D0
  
     DO nelem=1,elem_total
         CALL BD_ElemNodalDisp(uuNf,node_elem,dof_node,nelem,Nuuu)
         CALL BD_NodalRelRot(uuN0(:,nelem),node_elem,dof_node,Nrr0)
         CALL BD_NodalRelRot(Nuuu,node_elem,dof_node,Nrrr)
  
  
negative.code:    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_GenerateStaticElementForce'
  
     dof_elem = dof_node * node_elem
     rot_elem = (dof_node/2) * node_elem
  