Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 16747)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 16748)
@@ -237,15 +237,15 @@
 				break;
 			case PatersonEnum:
-				for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]);
+				for(i=0;i<numnodes;i++) B[i]=Paterson(temperature[i]);
 				element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum);
 				break;
 			case ArrheniusEnum:
 				element->GetVerticesCoordinates(&xyz_list);
-				for(i=0;i<numnodes;i++) B[i]=Arrhenius(values[i],surface[i]-xyz_list[i*3+2],element->GetMaterialParameter(MaterialsRheologyNEnum));
+				for(i=0;i<numnodes;i++) B[i]=Arrhenius(temperature[i],surface[i]-xyz_list[i*3+2],element->GetMaterialParameter(MaterialsRheologyNEnum));
 				element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum);
 				break;
 			case LliboutryDuvalEnum:
 				//for(i=0;i<numnodes;i++) B[i]=LliboutryDuval(values[i],pressure[i],material->GetN(),matpar->GetBeta(),matpar->GetReferenceTemperature(),matpar->GetHeatCapacity(),matpar->GetLatentHeat()); for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]); 
-				for(i=0;i<numnodes;i++) B[i]=Paterson(values[i]);
+				for(i=0;i<numnodes;i++) B[i]=Paterson(temperature[i]);
 				element->AddMaterialInput(MaterialsRheologyBEnum,&B[0],P1Enum); 
 				break;