Index: /issm/trunk-jpl/src/c/analyses/ResetBoundaryConditions.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ResetBoundaryConditions.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/analyses/ResetBoundaryConditions.cpp	(revision 16158)
@@ -10,6 +10,5 @@
 
 	/*variables: */
-	Vector<IssmDouble>*    yg    = NULL;
-	Nodes *nodes = NULL;
+	Vector<IssmDouble>* yg = NULL;
 
 	if(VerboseSolution()) _printf0_("   updating boundary conditions...\n");
Index: /issm/trunk-jpl/src/c/analyses/controltao_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16158)
@@ -29,5 +29,4 @@
 	AppCtx              user;
 	TaoSolver           tao = 0;
-	IssmDouble         *dummy        = NULL;
 	int                *control_list = NULL;
 	Vector<IssmDouble> *X            = NULL;
@@ -163,5 +162,4 @@
 	AppCtx   *user      = (AppCtx *)userCtx;
 	FemModel *femmodel  = user->femmodel;
-	Element  *element   = NULL;
 	int      *responses = NULL;
 
Index: /issm/trunk-jpl/src/c/analyses/dakota_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/dakota_core.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/analyses/dakota_core.cpp	(revision 16158)
@@ -161,7 +161,5 @@
 	char              *dakota_output_file = NULL;
 	char              *dakota_error_file  = NULL;
-	int                status             = 0;
 	Dakota::ModelLIter ml_iter;
-	Parameters        *parameters         = NULL;
 
 	/*Recover dakota_input_file, dakota_output_file and dakota_error_file, in the parameters dataset in parallel */
Index: /issm/trunk-jpl/src/c/analyses/gia_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/gia_core.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/analyses/gia_core.cpp	(revision 16158)
@@ -11,9 +11,8 @@
 void gia_core(FemModel* femmodel){
 
-	int i;
-	Vector<IssmDouble>*  wg  = NULL;
-	Vector<IssmDouble>*  dwdtg  = NULL;
-	IssmDouble*          x   = NULL;
-	IssmDouble*          y   = NULL;
+	Vector<IssmDouble> *wg    = NULL;
+	Vector<IssmDouble> *dwdtg = NULL;
+	IssmDouble          *x    = NULL;
+	IssmDouble          *y    = NULL;
 
 	/*parameters: */
Index: /issm/trunk-jpl/src/c/analyses/stressbalance_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/stressbalance_core.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/analyses/stressbalance_core.cpp	(revision 16158)
@@ -78,5 +78,5 @@
 	}
 
-	if (dim==3 & (isSIA || isSSA || isL1L2 || isHO)){
+	if (dim==3 && (isSIA || isSSA || isL1L2 || isHO)){
 		if(VerboseSolution()) _printf0_("   computing vertical velocities\n");
 		femmodel->SetCurrentConfiguration(StressbalanceVerticalAnalysisEnum);
Index: /issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/bamg/BamgQuadtree.cpp	(revision 16158)
@@ -535,5 +535,4 @@
 		register int l=0; // level
 		register BamgQuadtreeBox * b;
-		Icoor1 h=MaxISize;
 		Icoor1 hb =  MaxISize;
 		Icoor1 i0=0,j0=0;
@@ -549,5 +548,4 @@
 		ii[0]=i0;
 		jj[0]=j0;
-		h=hb;
 		do {    
 			b= pb[l];
@@ -560,6 +558,5 @@
 					  {
 						R2 XY(X,b->v[k]->r);
-						double dd;
-						if( (dd= LengthInterpole(Mx(XY), b->v[k]->m(XY)))  < seuil ){
+						if(LengthInterpole(Mx(XY), b->v[k]->m(XY)) < seuil){
 							return b->v[k]; 
 						}
Index: /issm/trunk-jpl/src/c/bamg/Direction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/Direction.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/bamg/Direction.cpp	(revision 16158)
@@ -27,5 +27,5 @@
 		int r =1; 
 		if (dir!= MaxICoor) {
-			Icoor2 x(dir/2),y1(MaxICoor/2-Abs(x)),y(dir%2?-y1:y1);
+			Icoor2 x(dir/2),y1(MaxICoor/2-Abs(x)),y((dir%2)?-y1:y1);
 			r = (x*i + y*j) >=0;
 		}
Index: /issm/trunk-jpl/src/c/bamg/Geometry.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/Geometry.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/bamg/Geometry.cpp	(revision 16158)
@@ -62,6 +62,4 @@
 
 		int verbose;
-		nbv=0;
-		nbe=0;
 		nbcurves=0;
 
@@ -286,5 +284,4 @@
 		int nbreqv=0;
 		int nbtan=0;
-		int nbcracked=0;
 		int i,count;
 
@@ -378,5 +375,4 @@
 		if (nbtan){
 			bamggeom->TangentAtEdges=xNew<double>(4*nbtan);
-			count=0;
 			for (i=0;i<nbe;i++){
 				if (edges[i].TgA() && edges[i][0].Corner()){
@@ -392,5 +388,4 @@
 					bamggeom->TangentAtEdges[4*i+3]=edges[i].tg[1].y;
 				}
-				count=count+1;
 			}
 		}
Index: /issm/trunk-jpl/src/c/bamg/ListofIntersectionTriangles.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/ListofIntersectionTriangles.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/bamg/ListofIntersectionTriangles.cpp	(revision 16158)
@@ -48,5 +48,4 @@
 
 		SegInterpolation* SegI=lSegsI;
-		SegI=lSegsI;
 		lSegsI[NbSeg].last=Size;
 		int EndSeg=Size;
@@ -225,8 +224,7 @@
 		Triangle *tbegin, *t;
 
-		long int verbose=2;
 		Icoor2 deta[3], deti,detj;
 		double ba[3];
-		int nbt =0,ifirst=-1,ilast;
+		int ifirst=-1,ilast;
 		int i0,i1,i2;
 		int ocut,i,j,k=-1;
@@ -288,5 +286,4 @@
 				if (t == tbegin) { // 
 					double ba,bb;
-					long int verbose=2;
 					AdjacentTriangle edge=CloseBoundaryEdge(a,t,ba,bb);
 					BamgVertex & v0 = *edge.EdgeVertex(0), & v1 = *edge.EdgeVertex(1);
Index: /issm/trunk-jpl/src/c/bamg/Mesh.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 16158)
@@ -269,5 +269,4 @@
 	void Mesh::ReadMesh(int* index,double* x,double* y,int nods,int nels){
 
-		double Hmin = HUGE_VAL;// the infinie value 
 		long i1,i2,i3;
 		long i;
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16158)
@@ -43,15 +43,11 @@
 Penta::Penta(int penta_id, int penta_sid, int index, IoModel* iomodel,int nummodels)
 	:PentaRef(nummodels)
-	,ElementHook(nummodels,index+1,6,iomodel) //index+1: material id, iomodel->numberofelements+1: matpar id
-                                                                      { //i is the element index
-
-	int i;
+	,ElementHook(nummodels,index+1,6,iomodel){
+
 	int penta_elements_ids[2];
 
 	/*Checks in debugging mode*/
-	/*{{{*/
 	_assert_(iomodel->Data(MeshUpperelementsEnum));
 	_assert_(iomodel->Data(MeshLowerelementsEnum));
-	/*}}}*/
 
 	/*id: */
@@ -4544,5 +4540,5 @@
 
 	/*Intermediaries*/
-	int        i,j,found=0;
+	int        i,found=0;
 	int        stabilization;
 	IssmDouble Jdet,phi,dt;
@@ -6526,5 +6522,4 @@
 	int*         pdoflist=NULL;
 	IssmDouble   FSreconditioning;
-	GaussPenta  *gauss;
 
 	/*Fetch number of nodes and dof for this finite element*/
@@ -8586,5 +8581,4 @@
 
 	/*Intermediaries */
-	int         i,j;
 	int         approximation;
 	IssmDouble  viscosity,Jdet;
@@ -8632,5 +8626,5 @@
 		material->GetViscosity3dFS(&viscosity,&epsilon[0]);
 
-		for(i=0;i<NUMVERTICES;i++){
+		for(int i=0;i<NUMVERTICES;i++){
 			pe->values[i*NDOF3+0]+=-Jdet*gauss->weight*viscosity*dw[0]*dbasis[2][i];
 			pe->values[i*NDOF3+1]+=-Jdet*gauss->weight*viscosity*dw[1]*dbasis[2][i];
Index: /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 16158)
@@ -846,7 +846,6 @@
 	 * where h is the interpolation function for node i.
 	 */
-	int i;
-	int num_dof=4;
-
+
+	int        num_dof=4;
 	IssmDouble L1L2l3[NUMNODESP1_2d];
 	IssmDouble dbasis[3][NUMNODESP1];
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16157)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16158)
@@ -24,7 +24,4 @@
 /*FUNCTION Tria::Tria(){{{*/
 Tria::Tria(){
-
-	int i;
-
 	this->nodes    = NULL;
 	this->vertices = NULL;
Index: /issm/trunk-jpl/src/c/toolkits/issm/Bucket.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/issm/Bucket.h	(revision 16157)
+++ /issm/trunk-jpl/src/c/toolkits/issm/Bucket.h	(revision 16158)
@@ -131,9 +131,7 @@
 		void SpawnBucketsPerCpu(DataSet* bucketsofcpu_i,int rank_i,int* rowranks){ /*{{{*/
 
-			int i,j;
-
 			/*go through our idxm index of rows this bucket owns, and spawn buckets  
 			 *if these rows belong to cpu rank_i. Use rowranks to determine this.*/
-			for(i=0;i<m;i++){
+			for(int i=0;i<m;i++){
 				if (rowranks[idxm[i]]==rank_i){
 					/*This row belongs to cpu rank_i, so spawn a bucket with this row, and add it to the bucketsofcpu_i dataset: */
@@ -147,6 +145,5 @@
 			}
 
-		};
-		/*}}}*/
+		}; /*}}}*/
 		int BucketType(void){ /*{{{*/
 
Index: /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h	(revision 16157)
+++ /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiVec.h	(revision 16158)
@@ -149,5 +149,4 @@
 
 			int        *row_indices_forcpu = NULL;
-			int        *col_indices_forcpu = NULL;
 			int        *modes_forcpu       = NULL;
 			doubletype *values_forcpu      = NULL;