Index: ../trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (revision 17923) +++ ../trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (revision 17924) @@ -31,8 +31,6 @@ int stabilization,finiteelement,smb_model; bool dakota_analysis; - bool issmbgradients; - bool ispdd; bool isdelta18o; bool isgroundingline; bool islevelset; Index: ../trunk-jpl/src/c/analyses/SmoothedSurfaceSlopeYAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/SmoothedSurfaceSlopeYAnalysis.cpp (revision 17923) +++ ../trunk-jpl/src/c/analyses/SmoothedSurfaceSlopeYAnalysis.cpp (revision 17924) @@ -118,7 +118,6 @@ } /*Intermediaries */ - int input_enum; IssmDouble Jdet,thickness,slope[2]; IssmDouble taud_y,norms,normv,vx,vy; IssmDouble *xyz_list = NULL; Index: ../trunk-jpl/src/c/analyses/LsfReinitializationAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/LsfReinitializationAnalysis.cpp (revision 17923) +++ ../trunk-jpl/src/c/analyses/LsfReinitializationAnalysis.cpp (revision 17924) @@ -111,7 +111,7 @@ /*Intermediaries */ const int dim = 2; - int i,row,col,stabilization; + int row,col,stabilization; IssmDouble Jdet,D_scalar; IssmDouble dtau = 1.; IssmDouble mu = 1.; @@ -398,7 +398,7 @@ void LsfReinitializationAnalysis::SetDistanceOnIntersectedElements(FemModel* femmodel){/*{{{*/ /* Intermediaries */ - int i,k; + int i; IssmDouble dmaxp,dmaxm,val; Element* element; Index: ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 17923) +++ ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 17924) @@ -3197,7 +3197,7 @@ /*Intermediaries*/ bool mainlyfloating; - int i,dim,domaintype,epssize; + int dim,domaintype; int migration_style,point1; IssmDouble alpha2,Jdet,fraction1,fraction2; IssmDouble gllevelset,phi=1.; @@ -5259,7 +5259,7 @@ ElementVector* StressbalanceAnalysis::CreatePVectorCouplingHOFSFriction(Element* element){/*{{{*/ /*Intermediaries*/ - int i,j,approximation; + int i,approximation; int dim=3; IssmDouble Jdet,Jdet2d,FSreconditioning; IssmDouble bed_normal[3]; Index: ../trunk-jpl/src/c/modules/ModelProcessorx/CreateEdges.cpp =================================================================== --- ../trunk-jpl/src/c/modules/ModelProcessorx/CreateEdges.cpp (revision 17923) +++ ../trunk-jpl/src/c/modules/ModelProcessorx/CreateEdges.cpp (revision 17924) @@ -148,7 +148,7 @@ /*Intermediaries*/ bool isv1,isv2; int facenbv,v1,v2; - int id_edge,id_element,id_face; + int id_edge,id_element; int elementnbe; /*Mesh dependent variables*/ Index: ../trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp =================================================================== --- ../trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp (revision 17923) +++ ../trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp (revision 17924) @@ -18,10 +18,10 @@ void CreateParameters(Parameters* parameters,IoModel* iomodel,char* rootpath,FILE* toolkitsoptionsfid,const int solution_type){ int i,j,m,k; - int numoutputs,domaintype,materialtype,smb_model; + int numoutputs,materialtype,smb_model; char** requestedoutputs = NULL; IssmDouble time; - bool ispdd,isdelta18o; + bool isdelta18o; /*parameters for mass flux:*/ int mass_flux_num_profiles = 0; Index: ../trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp =================================================================== --- ../trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp (revision 17923) +++ ../trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp (revision 17924) @@ -32,7 +32,6 @@ void IoModelToConstraintsx(Constraints* constraints,IoModel* iomodel,IssmDouble* spcdata,int M,int N,int analysis_type,int finite_element,int dof){ /*intermediary: */ - bool isnan; int i,j,count,elementnbv,numfacevertices; IssmDouble value; IssmDouble *times = NULL; @@ -40,7 +39,7 @@ bool spcpresent = false; /*Higher-order finite elements*/ - int v1,v2,v3,v4; + int v1,v2; bool *my_edges = NULL; bool *my_faces = NULL; bool *boundaryedge = NULL; Index: ../trunk-jpl/src/c/toolkits/mumps/MumpsSolve.cpp =================================================================== --- ../trunk-jpl/src/c/toolkits/mumps/MumpsSolve.cpp (revision 17923) +++ ../trunk-jpl/src/c/toolkits/mumps/MumpsSolve.cpp (revision 17924) @@ -199,7 +199,7 @@ /*Variables*/ ISSM_MPI_Comm comm; int num_procs; - int i,j; + int i; int nnz,local_nnz; int *irn_loc = NULL; int *jcn_loc = NULL; @@ -216,9 +216,9 @@ comm = IssmComm::GetComm(); /*First, some checks*/ - if (Kff_M!=Kff_N)_error_("stiffness matrix Kff should be square"); - if (uf_M!=Kff_M || uf_M!=pf_M)_error_("solution vector should be the same size as stiffness matrix Kff and load vector pf"); - if (uf_m!=Kff_m || uf_m!=pf_m)_error_("solution vector should be locally the same size as stiffness matrix Kff and load vector pf"); + if(Kff_M!=Kff_N)_error_("stiffness matrix Kff should be square"); + if(uf_M!=Kff_M || uf_M!=pf_M)_error_("solution vector should be the same size as stiffness matrix Kff and load vector pf"); + if(uf_m!=Kff_m || uf_m!=pf_m)_error_("solution vector should be locally the same size as stiffness matrix Kff and load vector pf"); /*Initialize matrix */ /*figure out number of non-zero entries: */ Index: ../trunk-jpl/src/c/toolkits/issm/IssmMpiSparseMat.h =================================================================== --- ../trunk-jpl/src/c/toolkits/issm/IssmMpiSparseMat.h (revision 17923) +++ ../trunk-jpl/src/c/toolkits/issm/IssmMpiSparseMat.h (revision 17924) @@ -140,15 +140,12 @@ /*FUNCTION Echo{{{*/ void Echo(void){ - int my_rank; - int i,j,k; - /*Do a synchronized dump across all the rows: */ - my_rank=IssmComm::GetRank(); - for(i=0;im << "\n"); - for (j=0;jm;j++){ + for(int j=0;jm;j++){ _printf_("row " << j << ":"); this->matrix[j]->Echo(); _printf_("\n"); @@ -189,7 +186,6 @@ int this_row_numvalues; int* this_row_cols = NULL; int* this_row_mods = NULL; - doubletype* this_row_values = NULL; int* numvalues_perrow = NULL; int row; @@ -377,7 +373,7 @@ doubletype norm,local_norm; doubletype absolute; - int i,j; + int i; switch(mode){ case NORM_INF: @@ -420,9 +416,6 @@ /*FUNCTION MatMult{{{*/ void MatMult(IssmAbsVec* Xin,IssmAbsVec* AXin){ - int i,j; - doubletype *X_serial = NULL; - /*A check on the types: */ if(IssmVecTypeFromToolkitOptions()!=MpiEnum)_error_("MatMult operation only possible with 'mpi' vectors"); @@ -431,11 +424,11 @@ IssmMpiVec* AX=(IssmMpiVec*)AXin; /*Serialize input Xin: */ - X_serial=X->ToMPISerial(); + doubletype* X_serial=X->ToMPISerial(); /*Every cpu has a serial version of the input vector. Use it to do the Matrix-Vector multiply *locally and plug it into AXin: */ - for(i=0;im;i++){ + for(int i=0;im;i++){ AX->vector[i]=this->matrix[i]->Mult(X_serial); } Index: ../trunk-jpl/src/c/cores/transient_core.cpp =================================================================== --- ../trunk-jpl/src/c/cores/transient_core.cpp (revision 17923) +++ ../trunk-jpl/src/c/cores/transient_core.cpp (revision 17924) @@ -19,16 +19,15 @@ void transient_core(FemModel* femmodel){ /*parameters: */ - int i; IssmDouble starttime,finaltime,dt,yts; - bool isstressbalance,ismasstransport,isFS,isthermal,isgroundingline,isgia,islevelset,isdamageevolution,ishydrology; - bool save_results,dakota_analysis; - bool time_adapt=false; - int output_frequency; - int domaintype,groundingline_migration,smb_model; - int numoutputs = 0; - Analysis *analysis = NULL; - char** requested_outputs = NULL; + bool isstressbalance,ismasstransport,isFS,isthermal,isgroundingline,isgia,islevelset,isdamageevolution,ishydrology; + bool save_results,dakota_analysis; + bool time_adapt=false; + int output_frequency; + int domaintype,groundingline_migration; + int numoutputs = 0; + Analysis *analysis = NULL; + char** requested_outputs = NULL; /*intermediary: */ Index: ../trunk-jpl/src/c/cores/masstransport_core.cpp =================================================================== --- ../trunk-jpl/src/c/cores/masstransport_core.cpp (revision 17923) +++ ../trunk-jpl/src/c/cores/masstransport_core.cpp (revision 17924) @@ -12,10 +12,9 @@ void masstransport_core(FemModel* femmodel){ /*parameters: */ - int i; int numoutputs,domaintype; bool save_results; - bool issmbgradients,ispdd,isdelta18o,isFS,isfreesurface,dakota_analysis; + bool isFS,isfreesurface,dakota_analysis; int solution_type; char** requested_outputs = NULL; Index: ../trunk-jpl/src/c/cores/steadystate_core.cpp =================================================================== --- ../trunk-jpl/src/c/cores/steadystate_core.cpp (revision 17923) +++ ../trunk-jpl/src/c/cores/steadystate_core.cpp (revision 17924) @@ -20,7 +20,6 @@ void steadystate_core(FemModel* femmodel){ /*intermediary: */ - int i; int step; Vector* ug = NULL; Vector* ug_old = NULL; Index: ../trunk-jpl/src/c/cores/controlm1qn3_core.cpp =================================================================== --- ../trunk-jpl/src/c/cores/controlm1qn3_core.cpp (revision 17923) +++ ../trunk-jpl/src/c/cores/controlm1qn3_core.cpp (revision 17924) @@ -127,8 +127,7 @@ void simul(long* indic,long* n,double* X,double* pf,double* G,long izs[1],float rzs[1],void* dzs){ /*Recover Femmodel*/ - double f; - int intn,solution_type; + int solution_type; FemModel *femmodel = (FemModel*)dzs; /*Recover responses*/ Index: ../trunk-jpl/src/c/classes/Materials/Matpar.cpp =================================================================== --- ../trunk-jpl/src/c/classes/Materials/Matpar.cpp (revision 17923) +++ ../trunk-jpl/src/c/classes/Materials/Matpar.cpp (revision 17924) @@ -20,7 +20,6 @@ /*FUNCTION Matpar::Matpar(int matpar_mid,IoModel* iomodel){{{*/ Matpar::Matpar(int matpar_mid, IoModel* iomodel){ - bool ispdd; bool isefficientlayer; int hydrology_model,smb_model; iomodel->Constant(&hydrology_model,HydrologyModelEnum); Index: ../trunk-jpl/src/c/classes/kriging/Observations.cpp =================================================================== --- ../trunk-jpl/src/c/classes/kriging/Observations.cpp (revision 17923) +++ ../trunk-jpl/src/c/classes/kriging/Observations.cpp (revision 17924) @@ -352,7 +352,6 @@ /*Intermediaries*/ int i,j,n_obs; IssmPDouble prediction,error; - IssmPDouble numerator,denominator,ratio; IssmPDouble *x = NULL; IssmPDouble *y = NULL; IssmPDouble *obs = NULL; Index: ../trunk-jpl/src/c/classes/Loads/Riftfront.cpp =================================================================== --- ../trunk-jpl/src/c/classes/Loads/Riftfront.cpp (revision 17923) +++ ../trunk-jpl/src/c/classes/Loads/Riftfront.cpp (revision 17924) @@ -38,11 +38,8 @@ int riftfront_node_ids[2]; int riftfront_elem_ids[2]; int riftfront_matpar_id; - int riftfront_type; - int riftfront_fill; IssmDouble riftfront_friction; IssmDouble riftfront_fractionincrement; - bool riftfront_shelf; int penalty_lock; /*intermediary: */ Index: ../trunk-jpl/src/c/classes/Elements/Element.cpp =================================================================== --- ../trunk-jpl/src/c/classes/Elements/Element.cpp (revision 17923) +++ ../trunk-jpl/src/c/classes/Elements/Element.cpp (revision 17924) @@ -592,7 +592,7 @@ void Element::InputCreate(IssmDouble* vector,IoModel* iomodel,int M,int N,int vector_type,int vector_enum,int code){/*{{{*/ /*Intermediaries*/ - int i,t,row; + int i,t; IssmDouble time; /*Branch on type of vector: nodal or elementary: */ Index: ../trunk-jpl/src/c/classes/IoModel.cpp =================================================================== --- ../trunk-jpl/src/c/classes/IoModel.cpp (revision 17923) +++ ../trunk-jpl/src/c/classes/IoModel.cpp (revision 17924) @@ -285,7 +285,6 @@ /*FUNCTION IoModel::StartTrace{{{*/ void IoModel::StartTrace(void){ - int i; bool autodiff = false; bool keep=false; IssmDouble gcTriggerRatio; Index: ../trunk-jpl/src/c/classes/Massfluxatgate.h =================================================================== --- ../trunk-jpl/src/c/classes/Massfluxatgate.h (revision 17923) +++ ../trunk-jpl/src/c/classes/Massfluxatgate.h (revision 17924) @@ -30,9 +30,9 @@ /*Massfluxatgate constructors, destructors :*/ /*FUNCTION Massfluxatgate() {{{*/ Massfluxatgate(){ - this->name = 0; - this->numsegments = 0; - this->segments = 0; + this->name = 0; + this->numsegments = 0; + this->segments = 0; } /*}}}*/ /*FUNCTION Massfluxatgate(char* name, int numsegments, doubletype* segments) {{{*/