Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15229)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15230)
@@ -5827,4 +5827,17 @@
 
 #ifdef _HAVE_HYDROLOGY_
+/*FUNCTION Tria::AllActive{{{*/
+bool Tria::AllActive(void){
+
+	/*Intermediaries*/
+	const int  numnodes = NUMVERTICES;
+
+	for(int i=0;i<numnodes;i++){
+		if(!this->nodes[i]->IsActive()) return false;
+	}
+
+	return true;
+}
+/*}}}*/
 /*FUNCTION Tria::CreateHydrologyWaterVelocityInput {{{*/
 void Tria::CreateHydrologyWaterVelocityInput(void){
@@ -6061,4 +6074,7 @@
 	GaussTria   *gauss = NULL;
 
+	/*Check that all nodes are active, else return empty matrix*/
+	if(!this->AllActive()) return NULL;
+
 	/*Initialize Element matrix*/
 	ElementMatrix* Ke=new ElementMatrix(nodes,NUMVERTICES,this->parameters,NoneApproximationEnum);
@@ -6229,4 +6245,7 @@
 	IssmDouble basis[numdof];
 	GaussTria* gauss=NULL;
+
+	/*Check that all nodes are active, else return empty matrix*/
+	if(!this->AllActive()) return NULL;
 
 	/*Initialize Element vector*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 15229)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 15230)
@@ -259,4 +259,5 @@
 		void    GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode);
 		void    GetHydrologyTransfer(Vector<IssmDouble>* transfer);
+		bool    AllActive(void);
 		#endif
 		#ifdef _HAVE_BALANCED_
Index: /issm/trunk-jpl/src/c/classes/Node.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15229)
+++ /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15230)
@@ -503,4 +503,11 @@
 
 	indexing.Activate();
+
+}
+/*}}}*/
+/*FUNCTION Node::IsActive{{{*/
+bool Node::IsActive(void){
+
+	return indexing.active;
 
 }
Index: /issm/trunk-jpl/src/c/classes/Node.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Node.h	(revision 15229)
+++ /issm/trunk-jpl/src/c/classes/Node.h	(revision 15230)
@@ -87,4 +87,5 @@
 		void   GetLocalDofList(int* poutdoflist,int approximation_enum,int setenum);
 		void   FreezeDof(int dof);
+		bool   IsActive(void);
 		void   Activate(void);
 		void   Deactivate(void);