Index: /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/Elements/Tria.cpp	(revision 22233)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/c/classes/Elements/Tria.cpp	(revision 22234)
@@ -4030,4 +4030,5 @@
 	/*early return if we are not on an ice cap:*/
 	if(!(this->inputs->Max(MaskIceLevelsetEnum)<0)){
+		this->inputs->AddInput(new TriaInput(SealevelEustaticMaskEnum,0,P0Enum));
 		*peustatic=0; //do not forget to assign this pointer, otherwise, global eustatic will be garbage!
 		return;
@@ -4036,4 +4037,5 @@
 	/*early return if we are fully floating: */
 	if (this->inputs->Max(MaskGroundediceLevelsetEnum)<=0){
+		this->inputs->AddInput(new TriaInput(SealevelEustaticMaskEnum,0,P0Enum));
 		*peustatic=0; //do not forget to assign this pointer, otherwise, global eustatic will be garbage!
 		return;
@@ -4044,4 +4046,7 @@
 		notfullygrounded=true; //used later on.
 	}
+		
+	/*Inform mask: */
+	this->inputs->AddInput(new TriaInput(SealevelEustaticMaskEnum,1,P0Enum));
 
 	/*recover material parameters: */