Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 22849)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 22850)
@@ -687,4 +687,16 @@
 		case 0:{
 					 element->GetInputValue(&thickness, gauss,ThicknessEnum);
+					 element->GetInputValue(&base, gauss,BaseEnum);
+					 element->GetInputValue(&sealevel, gauss,SealevelEnum);
+					 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum);
+					 IssmDouble rho_ice   = element->GetMaterialParameter(MaterialsRhoIceEnum);
+					 IssmDouble gravity   = element->GetMaterialParameter(ConstantsGEnum);
+					 p_ice   = gravity*rho_ice*thickness;
+					 p_water = rho_water*gravity*(sealevel-base);
+					 Neff = p_ice - p_water;
+				 }
+			  break;
+		case 1:{
+					 element->GetInputValue(&thickness, gauss,ThicknessEnum);
 					 IssmDouble rho_ice   = element->GetMaterialParameter(MaterialsRhoIceEnum);
 					 IssmDouble gravity   = element->GetMaterialParameter(ConstantsGEnum);
@@ -694,5 +706,5 @@
 				 }
 			  break;
-		case 1:{
+		case 2:{
 					 element->GetInputValue(&thickness, gauss,ThicknessEnum);
 					 element->GetInputValue(&base, gauss,BaseEnum);
@@ -706,16 +718,4 @@
 				 }
 			  break;
-		case 2:{
-					 element->GetInputValue(&thickness, gauss,ThicknessEnum);
-					 element->GetInputValue(&base, gauss,BaseEnum);
-					 element->GetInputValue(&sealevel, gauss,SealevelEnum);
-					 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoSeawaterEnum);
-					 IssmDouble rho_ice   = element->GetMaterialParameter(MaterialsRhoIceEnum);
-					 IssmDouble gravity   = element->GetMaterialParameter(ConstantsGEnum);
-					 p_ice   = gravity*rho_ice*thickness;
-					 p_water = rho_water*gravity*(sealevel-base);
-					 Neff = p_ice - p_water;
-				 }
-			  break;
 		case 3:
 			element->GetInputValue(&Neff,gauss,FrictionEffectivePressureEnum);