Index: /issm/trunk/src/c/objects/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Tria.cpp	(revision 3089)
+++ /issm/trunk/src/c/objects/Tria.cpp	(revision 3090)
@@ -4604,16 +4604,22 @@
 
 	/*Update internal data if inputs holds new values: */
-	if(inputs->Recover("thickness",&new_h[0],1,dofs,3,(void**)nodes)){
-		/*density, needed later: */
-		double di=(this->matpar->GetRhoIce()/this->matpar->GetRhoWater());
-		/*Go through grids: */
-		for (i=0;i<3;i++){
-			if(this->nodes[i]->IsOnShelf()){
-				this->b[i]=this->b[i]-di*(new_h[i]-h[i]); //hydrostatic equilibrium;
-			}
-			this->s[i]=this->b[i]+new_h[i];
-			this->h[i]=new_h[i];
-		}
-	}
+	if (id==1) printf("WARNING if QMU: no hydrostatic equilibrium is applied here (conflict with prognostic, which do not have matpar)\n");
+	//For now
+	inputs->Recover("thickness",&h[0],1,dofs,3,(void**)nodes);
+	//Later
+	/*
+		if(inputs->Recover("thickness",&new_h[0],1,dofs,3,(void**)nodes)){
+	//density, needed later:
+	double di=(this->matpar->GetRhoIce()/this->matpar->GetRhoWater());
+	//Go through grids:
+	for (i=0;i<3;i++){
+	if(this->nodes[i]->IsOnShelf()){
+	this->b[i]=this->b[i]-di*(new_h[i]-h[i]); //hydrostatic equilibrium;
+	}
+	this->s[i]=this->b[i]+new_h[i];
+	this->h[i]=new_h[i];
+	}
+	}
+	*/
 	inputs->Recover("surface",&s[0],1,dofs,3,(void**)nodes);
 	inputs->Recover("bed",&b[0],1,dofs,3,(void**)nodes);