Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 27473)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 27474)
@@ -19,8 +19,8 @@
 Friction::Friction(){/*{{{*/
 	this->element=NULL;
-	this->law=0;
-	this->linearize=false;
-	this->apply_dim = 1.;
-	this->domaintype=-1;
+	this->law        = 0;
+	this->linearize  = 0;
+	this->apply_dim  = 1.;
+	this->domaintype = -1;
 	this->vx_input=NULL;
 	this->vy_input=NULL;
@@ -33,5 +33,7 @@
 	/* Determine the dimension according to the domain type automatically. 
 	 * There are exceptions, e.g. HO, which needs the user to specify the dimension used in Friciton.*/
+
 	this->element=element_in;
+	this->linearize  = 0;
 
 	/* Load necessary parameters */
@@ -65,5 +67,5 @@
 		element_in->FindParam(&this->linearize,FrictionLinearizeEnum);
 		if(this->linearize){
-			this->linearize = false; /*Change to make sure we do the calculation once*/
+			this->linearize = 0; /*Change to make sure we do the calculation once*/
 			int numvertices = this->element->GetNumberOfVertices();
 			this->alpha2_list            = xNew<IssmDouble>(numvertices);
@@ -77,5 +79,5 @@
 				this->GetAlphaComplement(&this->alpha2_complement_list[iv], gauss);
 			}
-			this->linearize = true ; /*Change back, we are now all set!*/
+			this->linearize = 1; /*Change back, we are now all set!*/
 		}
 	}
@@ -92,4 +94,5 @@
 Friction::~Friction(){/*{{{*/
 	if(this->linearize){
+		printf("linearize = %i\n",this->linearize);
 		xDelete<IssmDouble>(this->alpha2_list);
 		xDelete<IssmDouble>(this->alpha2_complement_list);