anchor.diff:           if (allocated(indx       )) deallocate(indx       )
   
         end subroutine Cleanup
  -   END SUBROUTINE BD_DynamicSolutionGA2
  +END SUBROUTINE BD_DynamicSolutionGA2
   
      SUBROUTINE BD_GenerateDynamicElementGA2(uuN0,uuNf,vvNf,aaNf,                 &
                                              Stif0,Mass0,gravity,u,damp_flag,beta,&
  
positive.code: 
     END SUBROUTINE lubksb
  
     SUBROUTINE BD_UpdateDynamicGA2(ainc,uf,vf,af,xf,coef,node_total,dof_node)
  
     REAL(ReKi), INTENT(IN):: ainc(:)
     REAL(DbKi),INTENT(IN)::coef(:)
     INTEGER(IntKi), INTENT(IN):: node_total, dof_node
     REAL(ReKi), INTENT(INOUT):: uf(:),vf(:),af(:),xf(:)
  
     REAL(ReKi):: rotf_temp(3), roti_temp(3), rot_temp(3)
     INTEGER(IntKi):: i, j, temp_id
     INTEGER(IntKi)          :: ErrStat2                     ! Temporary Error status
     CHARACTER(ErrMsgLen)    :: ErrMsg2                      ! Temporary Error message
     CHARACTER(*), PARAMETER :: RoutineName = 'BD_UpdateDynamicGA2'
  
     DO i=2, node_total
         temp_id = (i - 1) * dof_node
  
negative.code:        ENDDO
     ENDDO
  
  END SUBROUTINE BD_AssembleStiffK
  
  SUBROUTINE BD_AssembleRHS(nelem,dof_elem,node_elem,dof_node,ElemRHS,GlobalRHS,ErrStat,ErrMsg)
     !-------------------------------------------------------------------------------
     ! This subroutine assembles global force vector.
     !-------------------------------------------------------------------------------
  
     REAL(ReKi),    INTENT(IN   ):: ElemRHS(:)    ! Total element force (Fc, Fd, Fb)
     INTEGER(IntKi),INTENT(IN   ):: nelem         ! Number of elements
     INTEGER(IntKi),INTENT(IN   ):: dof_elem      ! Degrees of freedom per element
     INTEGER(IntKi),INTENT(IN   ):: node_elem     ! Nodes per element
     INTEGER(IntKi),INTENT(IN   ):: dof_node      ! Degrees of freedom per node
     REAL(ReKi),    INTENT(INOUT):: GlobalRHS(:)  ! Global force vector
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
     INTEGER(IntKi)              :: i
     INTEGER(IntKi)              :: temp_id
  
     ErrStat = ErrID_None
     ErrMsg  = ""
  
     DO i=1,dof_elem
         temp_id = (nelem-1)*(node_elem-1)*dof_node+i
         GlobalRHS(temp_id) = GlobalRHS(temp_id)+ElemRHS(i)
     ENDDO
  
  END SUBROUTINE BD_AssembleRHS
  
  SUBROUTINE ludcmp(a,n,indx,d)
  !***************************************************************************************
  