Index: /issm/trunk-jpl/src/c/analyses/ad_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ad_core.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/analyses/ad_core.cpp	(revision 16144)
@@ -34,5 +34,4 @@
 	double     *xp  = NULL;
 	int my_rank=IssmComm::GetRank();
-
 
 	/*AD mode on?: */
Index: /issm/trunk-jpl/src/c/analyses/balancevelocity_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/balancevelocity_core.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/analyses/balancevelocity_core.cpp	(revision 16144)
@@ -14,5 +14,4 @@
 	/*parameters: */
 	bool save_results;
-
 
 	/*recover parameters: */
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16144)
@@ -5417,5 +5417,5 @@
 /*FUNCTION Penta::DrainWaterfraction{{{*/
 void Penta::DrainWaterfraction(void){
-    
+
     /*Intermediaries*/
 	bool isenthalpy;
@@ -5754,5 +5754,4 @@
 /*FUNCTION Penta::CreatePVectorAdjointHO{{{*/
 ElementVector* Penta::CreatePVectorAdjointHO(void){
-
 
 	/*Nothing to be done if not on surface*/
@@ -7592,5 +7591,5 @@
 	ElementMatrix* Ke=new ElementMatrix(Ke1,Ke2);
 	delete Ke1; delete Ke2;
-	
+
 	/*Compute HO Matrix with P1 element type\n");*/
 	this->element_type=P1Enum;
@@ -8282,5 +8281,4 @@
 		D_scalar=gauss->weight*Jdet;
 
-
 		/*Add stabilization*/
 		GetNodalFunctionsP1Derivatives(&dbasis[0][0],&xyz_list[0][0],gauss);
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16144)
@@ -1766,5 +1766,4 @@
 /*FUNCTION Tria::InputUpdateFromVector(IssmDouble* vector, int name, int type);{{{*/
 void  Tria::InputUpdateFromVector(IssmDouble* vector, int name, int type){
-
 
 	/*Check that name is an element input*/
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 16144)
@@ -189,5 +189,5 @@
 	/*Open input file on cpu 0: */
 	if(my_rank==0) IOMODEL = pfopen0(inputfilename ,"rb");
-	
+
 	/*Open toolkits file: */
 	toolkitsoptionsfid=pfopen(toolkitsfilename,"r");
Index: /issm/trunk-jpl/src/c/classes/IndependentObject.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/IndependentObject.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/classes/IndependentObject.cpp	(revision 16144)
@@ -106,10 +106,9 @@
 		if(my_rank==0){
 			if(fread(&pscalar,sizeof(IssmPDouble),1,fid)!=1)_error_("could not read scalar ");
-			
+
 			/*Now, before we even broadcast this to other nodes, declare the scalar  as an independent variable!: */
 			scalar<<=pscalar;
 		}
 
-		
 		ISSM_MPI_Bcast(&scalar,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm()); 
 
Index: /issm/trunk-jpl/src/c/modules/AllocateSystemMatricesx/AllocateSystemMatricesx.h
===================================================================
--- /issm/trunk-jpl/src/c/modules/AllocateSystemMatricesx/AllocateSystemMatricesx.h	(revision 16143)
+++ /issm/trunk-jpl/src/c/modules/AllocateSystemMatricesx/AllocateSystemMatricesx.h	(revision 16144)
@@ -11,4 +11,3 @@
 void MatrixNonzeros(int** pd_nnz,int** po_nnz,FemModel* femmodel,int set1enum,int set2enum);
 
-
 #endif
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp	(revision 16144)
@@ -245,5 +245,5 @@
 	CreateParametersDakota(&parameters,iomodel,rootpath,solution_type,analysis_type);
 	#endif
-	
+
 	/*Now, deal with toolkits options, which need to be put into the parameters dataset: */
 	ParseToolkitsOptionsx(parameters,toolkitsoptionsfid);
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 16144)
@@ -38,5 +38,4 @@
 	iomodel->FetchData(1,FlowequationElementEquationEnum);
 	finiteelement_list=xNewZeroInit<int>(iomodel->numberofelements);
-
 
 	/*Do we have coupling*/
Index: /issm/trunk-jpl/src/c/modules/ResetConstraintsx/ResetConstraintsx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ResetConstraintsx/ResetConstraintsx.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/modules/ResetConstraintsx/ResetConstraintsx.cpp	(revision 16144)
@@ -15,5 +15,4 @@
 
 void ResetConstraintsx(FemModel* femmodel){
-
 
 	/*Display message*/
Index: /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp	(revision 16144)
@@ -27,5 +27,5 @@
     else 
         Dret=D2;
-    
+
     /*check if dt*Dret>waterfraction. If so, drain whole waterfraction*/
     if(dt==0.){
Index: /issm/trunk-jpl/src/c/shared/Elements/LliboutryDuval.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/LliboutryDuval.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/shared/Elements/LliboutryDuval.cpp	(revision 16144)
@@ -6,5 +6,4 @@
 #include "../Numerics/types.h"
 #include "../Exceptions/exceptions.h"
-
 
 /* get ice stiffness B from enthalpy, pressure and flow law exponent*/
@@ -35,5 +34,5 @@
    *  
    *  Convert A to B :  B = A^(-1/n) */
-  
+
   /*Some physical constants (Aschwanden 2012)*/
   /*TODO: get those constants from model*/
@@ -46,5 +45,5 @@
   /*Intermediaries*/
   IssmDouble A,B,Tstar,Tpmp,H_sp,waterfraction;
-	
+
   _assert_(pressure>0);
   _assert_(enthalpy>0);
@@ -81,5 +80,5 @@
 	 * get enthalpy from temperature and water fraction,
 	 * and use LliboutryDuval(IssmDouble enthalpy, IssmDouble pressure,IssmDouble n) */
-  
+
   IssmDouble rho_ice=910; // kg/m^3
   IssmDouble g=9.81; //kg*m/s^2
@@ -103,8 +102,7 @@
     enthalpy=heatcapacity*(273.15 - Tref) + waterfraction*latentheat;
   }
-	
+
   B=LliboutryDuval(enthalpy, pressure, n);
 
   return B;
 }
-
Index: /issm/trunk-jpl/src/c/shared/MemOps/MemOps.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/MemOps/MemOps.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/shared/MemOps/MemOps.cpp	(revision 16144)
@@ -5,5 +5,4 @@
  *      Author: utke
  */
-
 
 #ifdef HAVE_CONFIG_H
Index: /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h	(revision 16143)
+++ /issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h	(revision 16144)
@@ -513,5 +513,4 @@
 			pf=(IssmMpiVec<IssmDouble>*)pfin;
 
-						
 			switch(IssmSolverTypeFromToolkitOptions()){
 				case MumpsEnum: {
Index: /issm/trunk-jpl/src/c/toolkits/issm/IssmToolkitUtils.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/issm/IssmToolkitUtils.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/toolkits/issm/IssmToolkitUtils.cpp	(revision 16144)
@@ -65,5 +65,5 @@
 	 *we try and stick with the Petsc vector types: */
 	vec_type=ToolkitOptions::GetToolkitOptionValue("vec_type");
-	
+
 	if (strcmp(vec_type,"mpi")==0){
 		vec_type_enum=MpiEnum;
Index: /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.cpp	(revision 16143)
+++ /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.cpp	(revision 16144)
@@ -110,5 +110,5 @@
 }/*}}}*/
 int ISSM_MPI_Allreduce(void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, ISSM_MPI_Comm comm){/*{{{*/
- 
+
   int rc=0;
 #ifdef _HAVE_MPI_
@@ -178,5 +178,5 @@
 }/*}}}*/
 int ISSM_MPI_Comm_free(ISSM_MPI_Comm *comm){ /*{{{*/
- 
+
   int rc=0;
 #ifdef _HAVE_MPI_
@@ -203,5 +203,5 @@
 }/*}}}*/
 int ISSM_MPI_Comm_size( ISSM_MPI_Comm comm, int *size){ /*{{{*/
- 
+
   int rc=0;
 #ifdef _HAVE_MPI_
@@ -263,5 +263,5 @@
 }/*}}}*/
 int ISSM_MPI_Gatherv(void *sendbuf, int sendcnt, ISSM_MPI_Datatype sendtype, void *recvbuf, int *recvcnts, int *displs, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm){/*{{{*/
- 
+
   int rc=0;
   assert(sendtype==recvtype); // we handle only identical representations
Index: /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h	(revision 16143)
+++ /issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h	(revision 16144)
@@ -51,5 +51,5 @@
 	#define ISSM_MPI_ANY_TAG       MPI_ANY_TAG
 	#define ISSM_MPI_ANY_SOURCE    MPI_ANY_SOURCE
-	
+
     /*other include files: */
 	#include "./commops/commops.h"
Index: /issm/trunk-jpl/src/m/classes/snowpack.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/snowpack.m	(revision 16143)
+++ /issm/trunk-jpl/src/m/classes/snowpack.m	(revision 16144)
@@ -124,5 +124,5 @@
 		end % }}}
 		function obj = setdefaultparameters(obj) % {{{
-	
+
 		%snowpack:  %{{{
 		obj.snowpack_meas_tss = 1;
@@ -279,5 +279,4 @@
 			%filters {{{
 			filter_values={'MIN_MAX','RATE_FILTER1','RATE_FILTER2','UNHEATED_RAIN_GAUGE_FILTER','WMO_UNDERCATCH_FILTER','WMO_UNDERCATCH_FILTER-SIMPLIFIED','UNVENTILLATED_TEMPERATURE_SENSOR','ADD_AN_OFFSET'};
-
 
 			md=checkfield(md,'snowpack.filters_ta_filter1','values',{filter_values});
Index: /issm/trunk-jpl/src/m/classes/thermal.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.m	(revision 16143)
+++ /issm/trunk-jpl/src/m/classes/thermal.m	(revision 16144)
@@ -40,5 +40,5 @@
 			%Should we use cold ice (default) or enthalpy formulation
 			obj.isenthalpy=0;
-            
+
             		%will basal boundary conditions be set dynamically
             		obj.isdynamicbasalspc=0;
@@ -69,5 +69,5 @@
 			fielddisplay(obj,'isenthalpy','use an enthalpy formulation to include temperate ice (default is 0)');
             		fielddisplay(obj,'isdynamicbasalspc',['enable dynamic setting of basal forcing. required for enthalpy formulation (default is 0)']);
-            
+
 		end % }}}
 		function marshall(obj,md,fid) % {{{
Index: /issm/trunk-jpl/src/m/exp/Shp2Exp.m
===================================================================
--- /issm/trunk-jpl/src/m/exp/Shp2Exp.m	(revision 16143)
+++ /issm/trunk-jpl/src/m/exp/Shp2Exp.m	(revision 16144)
@@ -9,5 +9,4 @@
 %
 %   See also EXPMASTER, EXPDOC
-
 
 	if ~exist(shapefilename,'file'),
@@ -22,5 +21,5 @@
 			ids=find(isnan(x));
 			x(ids)=[]; y(ids)=[];
-			
+
 			exp(end+1).x=x;
 			exp(end).y=y;
@@ -33,3 +32,2 @@
 
 	expwrite(exp,expfilename);
-
Index: /issm/trunk-jpl/src/wrappers/python/include/wrapper_macros.h
===================================================================
--- /issm/trunk-jpl/src/wrappers/python/include/wrapper_macros.h	(revision 16143)
+++ /issm/trunk-jpl/src/wrappers/python/include/wrapper_macros.h	(revision 16144)
@@ -43,5 +43,4 @@
 /* WRAPPER 3.2 {{{*/
 #define WRAPPER(modulename,...)  \
-\
 static PyObject* modulename(PyObject* self,PyObject* args);\
 static PyMethodDef modulename##_funcs[] = {\
@@ -49,5 +48,4 @@
 	{NULL,NULL,0,NULL}\
 };\
-\
 static struct PyModuleDef modulename##module= {\
 	PyModuleDef_HEAD_INIT,\
@@ -58,11 +56,8 @@
 	modulename##_funcs\
 };\
-\
 PyMODINIT_FUNC PyInit_##modulename(void){\
-\
 	import_array();\
 	return PyModule_Create(&modulename##module);\
 }\
-\
 static PyObject* modulename(PyObject* self,PyObject* args)
 /*}}}*/
@@ -70,5 +65,4 @@
 /* WRAPPER 2.7 {{{*/
 #define WRAPPER(modulename,...)  \
-\
 static PyObject* modulename(PyObject* self,PyObject* args);\
 static PyMethodDef modulename##_funcs[] = {\
@@ -76,11 +70,8 @@
 	{NULL,NULL,0,NULL}\
 };\
-\
 PyMODINIT_FUNC init##modulename(void){\
-\
 	import_array();\
 	(void) Py_InitModule(#modulename, modulename##_funcs);\
 }\
-\
 static PyObject* modulename(PyObject* self,PyObject* args)
 /*}}}*/