Index: /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/analyses/hydrology_core.cpp	(revision 15037)
@@ -63,5 +63,5 @@
 
 			/*transfer water column thickness to old water column thickness: */
-			
+
 			InputDuplicatex(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,WatercolumnEnum,WaterColumnOldEnum);
 
@@ -83,5 +83,4 @@
 
 			if(VerboseSolution()) _pprintLine_("   computing water transfer");
-			
 
 			if(VerboseSolution()) _pprintLine_("   computing water head");
Index: /issm/trunk-jpl/src/c/classes/Constraints/Constraints.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Constraints/Constraints.h	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Constraints/Constraints.h	(revision 15037)
@@ -28,5 +28,4 @@
 		/*}}}*/
 
-
 		/*numerics*/
 		int   NumberOfConstraints(void);
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15037)
@@ -9356,5 +9356,4 @@
 }
 /*}}}*/
-	
 
 #endif
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15037)
@@ -1914,5 +1914,5 @@
 		/*Get dof list: */
 		GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
-		
+
 		/*Use the dof list to index into the vector: */
 		for(int i=0;i<numdof;i++){
@@ -1922,5 +1922,5 @@
 		/*Add input to the element: */
 		this->inputs->AddInput(new TriaP1Input(name,values));
-		
+
 		/*Free ressources:*/
 		xDelete<int>(doflist);
@@ -6173,5 +6173,5 @@
 	/* Start  looping on the number of gaussian points: */
 	gauss=new GaussTria(2);
-	
+
 	for(int ig=gauss->begin();ig<gauss->end();ig++){
 
@@ -6245,5 +6245,5 @@
 		if(reCast<bool,IssmDouble>(dt)) scalar = scalar*dt;
 		for(int i=0;i<numdof;i++) pe->values[i]+=scalar*basis[i];
-		
+
 		/*Transient term*/
 		if(reCast<bool,IssmDouble>(dt)){
@@ -6336,5 +6336,4 @@
 	GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
 
-
 	/*Use the dof list to index into the solution vector: */
 	for(int i=0;i<numdof;i++){
@@ -6352,5 +6351,5 @@
 		this->parameters->FindParam(&penalty_factor,HydrologydcPenaltyFactorEnum);
 		kappa=kmax*pow(10.,penalty_factor);
-		
+
 		for(int i=0;i<NUMVERTICES;i++){
 			this->GetHydrologyDCInefficientHmax(&h_max,nodes[i]);
@@ -6474,5 +6473,5 @@
 /*FUNCTION Tria::GetHydrologyTransfer{{{*/
 void  Tria::GetHydrologyTransfer(Vector<IssmDouble>* transfer){
-	
+
 	const int  numdof         = NDOF1 *NUMVERTICES;
 	int        *doflist       = NULL;
@@ -6485,10 +6484,9 @@
 	IssmDouble epl_head[numdof],sed_head[numdof];
 
-
 	GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
-	
+
 	/*Get the flag to know if the efficient layer is present*/
 	this->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
-	
+
 	if(!isefficientlayer){
 		transfer->SetValues(numdof,doflist,&wh_trans[0],ADD_VAL);
@@ -6497,5 +6495,5 @@
 	/*Also get the flag to the transfer method (TO DO)*/
 	this->parameters->FindParam(&transfermethod,HydrologydcTransferFlagEnum);
-	
+
 	/*Switch between the different transfer methods cases(TO DO)*/
 	switch(transfermethod){
@@ -6505,11 +6503,11 @@
 	case 1:
 		this->parameters->FindParam(&leakage,HydrologydcLeakageFactorEnum);
-		
+
 		sed_trans = matpar->GetSedimentTransmitivity();
 		sed_thick = matpar->GetSedimentThickness();
-		
+
 		GetInputListOnVertices(&sed_head[0],SedimentHeadEnum);
 		GetInputListOnVertices(&epl_head[0],EplHeadEnum);
-		
+
 		for(int i=0;i<numdof;i++){
 			if(sed_head[i]>epl_head[i]){
@@ -6526,5 +6524,5 @@
 		_error_("no case higher than 1 for the Transfer method");
 	}
-	
+
 	/*Assign output pointer*/
 	transfer->SetValues(numdof,doflist,&wh_trans[0],ADD_VAL);
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 15037)
@@ -1458,5 +1458,5 @@
 #ifdef  _HAVE_DAKOTA_
 void FemModel::DakotaResponsesx(double* d_responses,char** responses_descriptors,int numresponsedescriptors,int d_numresponses){/*{{{*/
-	
+
 	int        i,j;
 	int        my_rank;
@@ -1562,8 +1562,8 @@
 
   int      i;
-	
+
 	/*intermediary: */
 	Element *element     = NULL;
-	
+
 	/*Go through elements, and add contribution from each element to the deflection vector wg:*/
 	for (i=0;i<elements->Size();i++){
@@ -1578,8 +1578,8 @@
 
 	Vector<IssmDouble>* basisg=NULL; 
-	
+
 	/*Vector allocation*/
 	basisg=new Vector<IssmDouble>(nodes->NumberOfNodes());
-	
+
 	for (int i=0;i<elements->Size();i++){
 		Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
@@ -1588,8 +1588,8 @@
 	/*Assemble*/
 	basisg->Assemble();
-	
+
 	/*Update Inputs*/
 	InputUpdateFromVectorx(elements,nodes,vertices,loads,materials,parameters,basisg,BasisIntegralEnum,NodesEnum);
-	
+
 }
 /*}}}*/
@@ -1598,8 +1598,8 @@
 
 	Vector<IssmDouble>* transferg=NULL; 
-	
+
 	/*Vector allocation*/
 	transferg=new Vector<IssmDouble>(nodes->NumberOfNodes());
-	
+
 	for (int i=0;i<elements->Size();i++){
 		Element* element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
@@ -1608,8 +1608,8 @@
 	/*Assemble*/
 	transferg->Assemble();
-	
+
 	/*Update Inputs*/
 	InputUpdateFromVectorx(elements,nodes,vertices,loads,materials,parameters,transferg,WaterTransferEnum,NodesEnum);
-	
-}
-/*}}}*/
+
+}
+/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Inputs/Input.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/Input.h	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Inputs/Input.h	(revision 15037)
@@ -17,5 +17,4 @@
 #include "../../shared/shared.h"
 /*}}}*/
-
 
 class Input: public Object{
Index: /issm/trunk-jpl/src/c/classes/Inputs/Inputs.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/Inputs.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Inputs/Inputs.cpp	(revision 15037)
@@ -428,5 +428,5 @@
 /*Methods relating to inputs: */
 void IsInputConverged(IssmDouble* peps, Input** new_inputs,Input** old_inputs,int num_inputs,int criterion_enum){ /*{{{*/
- 
+
 	/*figure out convergence at the input level. 
 	  We are given a list of inputs, new and old, and a criterion, and using the
Index: /issm/trunk-jpl/src/c/classes/Node.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15037)
@@ -926,5 +926,4 @@
 /*}}}*/
 
-
 /*Methods inherent to Node: */
 int* GetLocalDofList(Node** nodes,int numnodes,int setenum,int approximation){ /*{{{*/
@@ -987,5 +986,4 @@
 int* GetGlobalDofList(Node** nodes,int numnodes,int setenum,int approximation){/*{{{*/
 
-
 	int  i,numdof,count;
 	int* ndof_list=NULL;
@@ -1025,5 +1023,4 @@
 int GetNumberOfDofs(Node** nodes,int numnodes,int setenum,int approximation){/*{{{*/
 
-
 	/*output: */
 	int numberofdofs=0;
@@ -1038,5 +1035,4 @@
 #ifdef _HAVE_DIAGNOSTIC_
 void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum){/*{{{*/
-
 
 	int* cs_array=NULL;
@@ -1054,5 +1050,4 @@
 /*}}}*/
 void TransformInvStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array){/*{{{*/
-
 
 	int     i,j;
@@ -1105,5 +1100,4 @@
 void TransformLoadVectorCoord(ElementVector* pe,Node** nodes,int numnodes,int* cs_array){/*{{{*/
 
-
 	int     i;
 	int     numdofs   = 0;
@@ -1139,5 +1133,4 @@
 void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int cs_enum){/*{{{*/
 
-
 	int* cs_array=NULL;
 
@@ -1154,5 +1147,4 @@
 /*}}}*/
 void TransformSolutionCoord(IssmDouble* solution,Node** nodes,int numnodes,int* cs_array){/*{{{*/
-
 
 	int     i;
@@ -1189,5 +1181,4 @@
 void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int cs_enum){/*{{{*/
 
-
 	int* cs_array=NULL;
 
@@ -1204,5 +1195,4 @@
 /*}}}*/
 void TransformStiffnessMatrixCoord(ElementMatrix* Ke,Node** nodes,int numnodes,int* cs_array){/*{{{*/
-
 
 	int     i,j;
@@ -1239,5 +1229,4 @@
 /*}}}*/
 void CoordinateSystemTransform(IssmDouble** ptransform,Node** nodes,int numnodes,int* cs_array){/*{{{*/
-
 
 	int     i,counter;
Index: /issm/trunk-jpl/src/c/classes/Nodes.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Nodes.h	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Nodes.h	(revision 15037)
@@ -9,5 +9,4 @@
 class Nodes;
 class Materials;
-
 
 /*!\brief Declaration of Nodes class.
Index: /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Params/Parameters.cpp	(revision 15037)
@@ -571,5 +571,4 @@
 	 * and parameters
 	 */ 
-
 
 	char* options=NULL;
Index: /issm/trunk-jpl/src/c/classes/Params/Parameters.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Params/Parameters.h	(revision 15036)
+++ /issm/trunk-jpl/src/c/classes/Params/Parameters.h	(revision 15037)
@@ -62,4 +62,3 @@
 void        ToolkitsOptionsFromAnalysis(Parameters* parameters,int analysis_type);
 
-
 #endif //ifndef _PARAMETERS_H_
Index: /issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/GiaDeflectionCorex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/GiaDeflectionCorex.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/GiaDeflectionCorex.cpp	(revision 15037)
@@ -16,5 +16,5 @@
 	double pset[7];
 };
-	
+
 struct blocko{
 	double rhoi;
@@ -154,5 +154,5 @@
 	/*dmi: */
 	blockt_dmi=xNew<IssmDouble>(Ntimm);
-	
+
 	/*irate: */
 	blockn_.irate=irate; 
@@ -187,3 +187,2 @@
 
 }
-      
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.h
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.h	(revision 15036)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.h	(revision 15037)
@@ -9,5 +9,4 @@
 
 #include "../../classes/classes.h"
-
 
 void ModelProcessorx(Elements** pelements, Nodes** pnodes, Vertices** pvertices, Materials** pmaterials, Constraints** pconstraints, Loads** ploads, Parameters** pparameters, FILE* iomodel_handle,char* rootpath,const int solution_type,const int nummodels,const int* analysis_type_listh);
Index: /issm/trunk-jpl/src/c/shared/Numerics/extrema.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Numerics/extrema.cpp	(revision 15036)
+++ /issm/trunk-jpl/src/c/shared/Numerics/extrema.cpp	(revision 15037)
@@ -38,4 +38,2 @@
 }
 #endif
-
-