Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18195)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18196)
@@ -108,7 +108,7 @@
 	parameters->AddObject(iomodel->CopyConstantObject(FrictionLawEnum));
 
-	/*XTH parameters*/
+	/*XTH LATH parameters*/
 	iomodel->Constant(&fe_FS,FlowequationFeFSEnum);
-	if(fe_FS==XTaylorHoodEnum){
+	if(fe_FS==XTaylorHoodEnum || fe_FS==LATaylorHoodEnum){
 		parameters->AddObject(iomodel->CopyConstantObject(AugmentedLagrangianREnum));
 		parameters->AddObject(iomodel->CopyConstantObject(AugmentedLagrangianThetaEnum));
@@ -125,5 +125,5 @@
 
 	/*Intermediaries*/
-	int    materials_type,finiteelement;
+	int    materials_type,finiteelement,fe_FS;
 	int    approximation,frictionlaw;
 	int*   finiteelement_list=NULL;
@@ -236,4 +236,9 @@
 	if(islevelset){
 		iomodel->FetchDataToInput(elements,IceMaskNodeActivationEnum);
+	}
+	/*LATH parameters*/
+	iomodel->Constant(&fe_FS,FlowequationFeFSEnum);
+	if(fe_FS==LATaylorHoodEnum){
+		iomodel->FetchDataToInput(elements,PressureEnum);
 	}
 
@@ -4568,5 +4573,5 @@
 	element->InputChangeName(VxEnum,VxPicardEnum);
 	element->InputChangeName(VyEnum,VyPicardEnum);
-	element->InputChangeName(PressureEnum,PressurePicardEnum);
+	if(pnumdof>0) element->InputChangeName(PressureEnum,PressurePicardEnum);
 	if(dim==3) element->InputChangeName(VzEnum,VzPicardEnum);
 
@@ -4575,5 +4580,5 @@
 	element->AddInput(VyEnum,      vy,      element->VelocityInterpolation());
 	element->AddInput(VelEnum,     vel,     element->VelocityInterpolation());
-	element->AddInput(PressureEnum,pressure,element->PressureInterpolation());
+	if(pnumdof>0) element->AddInput(PressureEnum,pressure,element->PressureInterpolation());
 	if(dim==3) element->AddInput(VzEnum,vz, element->VelocityInterpolation());