Index: ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp	(revision 18497)
+++ ../trunk-jpl/src/c/analyses/SeaiceAnalysis.cpp	(revision 18498)
@@ -130,6 +130,7 @@
 	/*Clean up and return*/
 	xDelete<IssmDouble>(xyz_list);
 	xDelete<IssmDouble>(B);
+	xDelete<IssmDouble>(M);
 	delete gauss;
 	return Ke;
 }
@@ -206,7 +207,7 @@
 		windvy_input->GetInputValue(&windvy,gauss);
 		oceanvx_input->GetInputValue(&oceanvx,gauss);
 		oceanvy_input->GetInputValue(&oceanvy,gauss);
-		oceanssh_input->GetInputDerivativeValue(&dssh[2],xyz_list,gauss);
+		oceanssh_input->GetInputDerivativeValue(&dssh[0],xyz_list,gauss);
 
 		/*Previous speed (finite differencing)*/
 		for(int i=0;i<numnodes;i++){
@@ -214,7 +215,7 @@
 			pe->values[i*2+1]+=Jdet*gauss->weight*rho_ice*thickness*vy*basis[i];
 		}
 
-		/*Atmespheric forcing*/
+		/*Atmospheric forcing*/
 		vnorm    = sqrt(windvx*windvx + windvy*windvy);
 		constant_part = concentration*air_coef*rho_air*(air_lin_drag_coef+air_quad_drag_coef*vnorm);
 		for(int i=0;i<numnodes;i++){
