Index: /issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyShreveAnalysis.cpp	(revision 17294)
@@ -117,4 +117,5 @@
 	element->GetVerticesCoordinates(&xyz_list);
 	element->FindParam(&dt,TimesteppingTimeStepEnum);
+	element->FindParam(&diffusivity,HydrologyshreveStabilizationEnum);
 	Input* vx_input=element->GetInput(HydrologyWaterVxEnum); _assert_(vx_input);
 	Input* vy_input=element->GetInput(HydrologyWaterVyEnum); _assert_(vy_input);
Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 17294)
@@ -814,5 +814,14 @@
 }/*}}}*/
 void Element::TransformSolutionCoord(IssmDouble* values,Node** nodes_list,int numnodes,int transformenum){/*{{{*/
-	this->TransformSolutionCoord(values,nodes_list,numnodes,transformenum);
+	/*NOT NEEDED*/
+	/*All nodes have the same Coordinate System*/
+	int* cs_array = xNew<int>(numnodes);
+	for(int i=0;i<numnodes;i++) cs_array[i]=transformenum;
+
+	/*Call core*/
+	this->TransformSolutionCoord(values,nodes_list,numnodes,cs_array);
+
+	/*Clean-up*/
+	xDelete<int>(cs_array);
 }/*}}}*/
 void Element::TransformSolutionCoord(IssmDouble* solution,Node** nodes_list,int numnodes,int* cs_array){/*{{{*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17294)
@@ -200,5 +200,5 @@
 	IssmDouble  pressure,viscosity,Jdet2d;
 	IssmDouble  bed_normal[3];
-	IssmDouble  basalforce[3];
+	IssmDouble  basalforce[3] = {0.};
 	IssmDouble  epsilon[6]; /* epsilon=[exx,eyy,ezz,exy,exz,eyz];*/
 	IssmDouble  stresstensor[6]={0.0};
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17294)
@@ -117,5 +117,5 @@
 	this->parameters->FindParam(&waitonlock,SettingsWaitonlockEnum);
 	this->parameters->FindParam(&lockfilename,LockFileNameEnum);
-	if (waitonlock>0){
+	if(waitonlock){
 		_printf0_("write lock file:\n");
 		WriteLockFile(lockfilename);
Index: /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp	(revision 17294)
@@ -589,5 +589,5 @@
 
 	IssmDouble  penetration;
-	int         activate;
+	bool        activate;
 	int         unstable;
 	IssmDouble  vx1;
@@ -634,6 +634,6 @@
 
 	/*activation: */
-	if(penetration<0)activate=1;
-	else  activate=0;
+	if(penetration<0)activate=true;
+	else  activate=false;
 
 	/*Here, we try to avoid zigzaging. When a penalty activates and deactivates for more than penalty_lock times, 
Index: /issm/trunk-jpl/src/c/classes/Misfit.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Misfit.h	(revision 17293)
+++ /issm/trunk-jpl/src/c/classes/Misfit.h	(revision 17294)
@@ -113,8 +113,8 @@
 
 			 int i;
-			 IssmDouble misfit_t=0;
-			 IssmDouble all_misfit_t=0;
+			 IssmDouble misfit_t=0.;
+			 IssmDouble all_misfit_t=0.;
 			 IssmDouble dt;
-			 IssmDouble area_t;
+			 IssmDouble area_t=0.;
 			 IssmDouble all_area_t;
 			 IssmDouble time,starttime,finaltime;
Index: /issm/trunk-jpl/src/c/shared/FSanalyticals/FSanalytical3.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/FSanalyticals/FSanalytical3.cpp	(revision 17293)
+++ /issm/trunk-jpl/src/c/shared/FSanalyticals/FSanalytical3.cpp	(revision 17294)
@@ -4,5 +4,5 @@
  */
 
-#include <cmath.h>
+#include <math.h>
 #include "../Numerics/types.h"
 #include "../Numerics/constants.h"