Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 16840)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 16841)
@@ -170,12 +170,12 @@
 
 		element->JacobianDeterminant(&Jdet,xyz_list,gauss);
-		D_scalar=gauss->weight*Jdet*kappa;
+		D_scalar=gauss->weight*Jdet;
 		if(dt!=0.) D_scalar=D_scalar*dt;
 
 		/*Conduction: */
 		GetBConduct(B,element,xyz_list,gauss); 
-		D[0][0]=D_scalar;
-		D[1][1]=D_scalar;
-		D[2][2]=D_scalar;
+		D[0][0]=D_scalar*kappa;
+		D[1][1]=D_scalar*kappa;
+		D[2][2]=D_scalar*kappa;
 		TripleMultiply(B,3,numnodes,1,
 					&D[0][0],3,3,0,