Index: ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18569)
+++ ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18570)
@@ -247,7 +247,6 @@
 	/*LATH parameters*/
 	iomodel->Constant(&fe_FS,FlowequationFeFSEnum);
 	if(fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){
-		iomodel->FetchDataToInput(elements,PressureEnum,0.);
 		InputUpdateFromConstantx(elements,0.,SigmaNNEnum);
 	}
 
@@ -3028,6 +3027,11 @@
 					BtBUzawa,1);
 	}
 
+	/*The pressure augmentation should not be transformed*/
+	MatrixMultiply(BtBUzawa,pnumdof,numdof,1,
+				BtBUzawa,pnumdof,numdof,0,
+				&Ke->values[0],1);
+
 	if(element->IsOnBase() && 0){ 
 		element->FindParam(&rl,AugmentedLagrangianRlambdaEnum);
 		element->GetVerticesCoordinatesBase(&xyz_list_base);
@@ -3069,11 +3073,6 @@
 	/*Transform Coordinate System*/
 	element->TransformStiffnessMatrixCoord(Ke,cs_list);
 
-	/*The pressure augmentation should not be transformed*/
-	MatrixMultiply(BtBUzawa,pnumdof,numdof,1,
-				BtBUzawa,pnumdof,numdof,0,
-				&Ke->values[0],1);
-
 	/*Clean up and return*/
 	delete gauss;
 	xDelete<IssmDouble>(xyz_list);
@@ -3444,20 +3443,13 @@
 		delete pe2;
 		return pe3;
 	}
-	else if(fe_FS==LATaylorHoodEnum){
-		ElementVector* pe2=CreatePVectorFSViscousLATH(element);
+	else if(fe_FS==LATaylorHoodEnum || fe_FS==LACrouzeixRaviartEnum){
+		ElementVector* pe2=CreatePVectorFSViscousLA(element);
 		ElementVector* pe3 = new ElementVector(pe,pe2);
 		delete pe;
 		delete pe2;
 		return pe3;
 	}
-	else if(fe_FS==LACrouzeixRaviartEnum){
-		ElementVector* pe2=CreatePVectorFSViscousLACR(element);
-		ElementVector* pe3 = new ElementVector(pe,pe2);
-		delete pe;
-		delete pe2;
-		return pe3;
-	}
 	return pe;
 }/*}}}*/
 ElementVector* StressbalanceAnalysis::CreatePVectorFSFriction(Element* element){/*{{{*/
@@ -3774,7 +3766,7 @@
 ElementVector* StressbalanceAnalysis::CreatePVectorFSViscousLA(Element* element){/*{{{*/
 
 	int         i,dim;
-	IssmDouble  Jdet,r,pressure;
+	IssmDouble  Jdet,pressure;
 	IssmDouble  bed_normal[3];
 	IssmDouble *xyz_list      = NULL;
 	IssmDouble *xyz_list_base = NULL;
@@ -3796,7 +3788,6 @@
 	IssmDouble*    dbasis  = xNew<IssmDouble>(3*numnodes);
 
 	/*Retrieve all inputs and parameters*/
-	element->FindParam(&r,AugmentedLagrangianREnum);
 	element->GetVerticesCoordinates(&xyz_list);
 
 	/*Get pressure and sigmann*/
@@ -3845,7 +3836,7 @@
 	}
 
 	/*Transform coordinate system*/
-	//element->TransformLoadVectorCoord(pe,cs_list); Do not transform augmentation
+	element->TransformLoadVectorCoord(pe,cs_list);
 
 	/*Clean up and return*/
 	delete gauss;
