Index: ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp	(revision 18519)
+++ ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp	(revision 18520)
@@ -440,7 +440,7 @@
 	IssmDouble f = exp(ridging_exponent*(1.-concentration));
 	IssmDouble E = young_modulus*(1.-damage)*f;
 
-	/*Let's build C = hE*D */
+	/*Let's build C = E*D */
 	for(int i=0;i<3;i++){
 		for(int j=0;j<3;j++){
 			C[i*3+j] = E*D[i][j];
@@ -471,8 +471,10 @@
 		B[2*numnodes*0+2*i+1] = 0.;
 		B[2*numnodes*1+2*i+0] = 0.;
 		B[2*numnodes*1+2*i+1] = dbasis[1*numnodes+i];
-		B[2*numnodes*2+2*i+0] = .5*dbasis[1*numnodes+i];
-		B[2*numnodes*2+2*i+1] = .5*dbasis[0*numnodes+i];
+		//B[2*numnodes*2+2*i+0] = .5*dbasis[1*numnodes+i];
+		//B[2*numnodes*2+2*i+1] = .5*dbasis[0*numnodes+i];
+		B[2*numnodes*2+2*i+0] = dbasis[1*numnodes+i];
+		B[2*numnodes*2+2*i+1] = dbasis[0*numnodes+i];
 	}
 
 	/*Clean-up*/
