Index: /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp	(revision 19160)
+++ /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp	(revision 19161)
@@ -2077,4 +2077,5 @@
 	/*Intermediaries*/
 	int      domaintype,dim;
+
 	Element* basalelement;
 
@@ -2141,7 +2142,6 @@
 	Input* cinput = basalelement->GetInput(FrictionCEnum);
 	Input* Asinput = basalelement->GetInput(FrictionAsEnum);
-	Input* Ninput = basalelement->GetInput(FrictionEffectivePressureEnum);
 	Input* nInput =basalelement->GetInput(MaterialsRheologyNEnum);
-	
+	Input* Ninput = basalelement->GetInput(FrictionEffectivePressureEnum);	
 	/* Start  looping on the number of gaussian points: */
 	Gauss* gauss=basalelement->NewGauss(4);
Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 19160)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 19161)
@@ -109,5 +109,6 @@
 	bool dakota_analysis,islevelset,isenthalpy;
 	int frictionlaw;
-
+	int FrictionCoupling;
+	
 	/*Now, is the model 3d? otherwise, do nothing: */
 	if(iomodel->domaintype==Domain2DhorizontalEnum)return;
@@ -175,8 +176,11 @@
 			break;
 		case 3:
+			iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
 			iomodel->FetchDataToInput(elements,FrictionCEnum);
 			iomodel->FetchDataToInput(elements,FrictionAsEnum);
 			iomodel->FetchDataToInput(elements,FrictionQEnum);
-			iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			if (FrictionCoupling==0){
+				iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			}
 			break;
 		case 4:
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 19160)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 19161)
@@ -673,4 +673,5 @@
 	int    materials_type,finiteelement,fe_FS;
 	int    approximation,frictionlaw;
+	int    FrictionCoupling;
 	int*   finiteelement_list=NULL;
 	bool   isSSA,isL1L2,isHO,isFS,iscoupling;
@@ -800,8 +801,11 @@
 			break;
 		case 3:
+			iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
 			iomodel->FetchDataToInput(elements,FrictionCEnum);
 			iomodel->FetchDataToInput(elements,FrictionAsEnum);
 			iomodel->FetchDataToInput(elements,FrictionQEnum);
-			iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			if (FrictionCoupling==0){
+				iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			}
 			break;
 		case 4:
Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 19160)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 19161)
@@ -62,4 +62,5 @@
 
 	int finiteelement = P1Enum;
+	
 	if(iomodel->domaintype==Domain3DEnum) iomodel->FetchData(2,MeshVertexonbaseEnum,MeshVertexonsurfaceEnum);
 	::CreateNodes(nodes,iomodel,ThermalAnalysisEnum,finiteelement);
@@ -72,5 +73,5 @@
 
 	int frictionlaw;
-
+	int FrictionCoupling;
 	/*Now, is the model 3d? otherwise, do nothing: */
 	if(iomodel->domaintype==Domain2DhorizontalEnum)return;
@@ -130,8 +131,11 @@
 			break;
 		case 3:
+			iomodel->Constant(&FrictionCoupling,FrictionCouplingEnum);
 			iomodel->FetchDataToInput(elements,FrictionCEnum);
 			iomodel->FetchDataToInput(elements,FrictionAsEnum);
 			iomodel->FetchDataToInput(elements,FrictionQEnum);
-			iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			if (FrictionCoupling==0){
+				iomodel->FetchDataToInput(elements,FrictionEffectivePressureEnum);
+			}
 			break;
 		case 4: