Index: /issm/trunk/src/c/objects/Loads/Numericalflux.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Numericalflux.cpp	(revision 5455)
+++ /issm/trunk/src/c/objects/Loads/Numericalflux.cpp	(revision 5456)
@@ -723,5 +723,7 @@
 		vyaverage_input=tria->inputs->GetInput(VyEnum);
 	}
-	thickness_input=tria->inputs->GetInput(ThicknessEnum);
+	/*Here, as it is a forcing, we have H=Hobs by default (for control methods)*/
+	thickness_input=tria->inputs->GetInput(ThicknessObsEnum);
+	if (!thickness_input) thickness_input=tria->inputs->GetInput(ThicknessEnum);
 
 	/* Get node coordinates and normal vector: */