Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 22046)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 22047)
@@ -289,4 +289,5 @@
 
 	/*diverse: */
+	int         CoupledFlag;
 	IssmDouble  r,s;
 	IssmDouble  drag_p, drag_q;
@@ -308,4 +309,5 @@
 	IssmDouble rho_ice          = element->GetMaterialParameter(MaterialsRhoIceEnum);
 	IssmDouble gravity          = element->GetMaterialParameter(ConstantsGEnum);
+	element->parameters->FindParam(&CoupledFlag,FrictionCouplingEnum);
 
 	//compute r and q coefficients: */
@@ -314,5 +316,18 @@
 
 	//From base and thickness, compute effective pressure when drag is viscous:
-	Neff=gravity*(rho_ice*thickness+rho_water*(base-sealevel));
+	switch(CoupledFlag){
+		case 0:
+			Neff=gravity*(rho_ice*thickness+rho_water*(base-sealevel));
+			break;	
+		case 1:
+			element->GetInputValue(&Neff,gauss,FrictionEffectivePressureEnum);
+			break;
+		case 2:
+			element->GetInputValue(&Neff,gauss,EffectivePressureEnum);
+			break;
+		default:
+			_error_("not supported");
+	}
+
 	if(Neff<0)Neff=0;