Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 28252)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.cpp	(revision 28253)
@@ -101,8 +101,8 @@
 /*}}}*/
 Friction::Friction(Element* element_in,int dim) : Friction(element_in) {/*{{{*/
-	this->apply_dim = reCast<IssmDouble>(dim);
+	this->apply_dim = reCast<IssmPDouble>(dim);
 }
 /*}}}*/
-Friction::Friction(Element* element_in,IssmDouble dim) : Friction(element_in) {/*{{{*/
+Friction::Friction(Element* element_in,IssmPDouble dim) : Friction(element_in) {/*{{{*/
 	this->apply_dim = dim;
 }
Index: /issm/trunk-jpl/src/c/classes/Loads/Friction.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Loads/Friction.h	(revision 28252)
+++ /issm/trunk-jpl/src/c/classes/Loads/Friction.h	(revision 28253)
@@ -21,5 +21,5 @@
 		int         domaintype;
 		int         linearize;
-		IssmDouble  apply_dim;
+		IssmPDouble apply_dim;
 		Input      *vx_input;
 		Input      *vy_input;
@@ -32,5 +32,5 @@
 		Friction(Element* element_in);
 		Friction(Element* element_in, int dim);
-		Friction(Element* element_in, IssmDouble dim);
+		Friction(Element* element_in, IssmPDouble dim);
 		~Friction();
 
Index: /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 28252)
+++ /issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/Gembx.cpp	(revision 28253)
@@ -228,5 +228,5 @@
 	IssmDouble z1=0.0;
 	int lice=0;
-	for(int l=0;(l<m & d[l]<dPHC-Dtol);l++){
+	for(int l=0;(l<m && d[l]<dPHC-Dtol);l++){
 		z1=z1+dz[l]*d[l]; //mm
 		lice=l+1;
@@ -665,5 +665,5 @@
 				lice=l+1;
 			}
-			if (depthsnow<=0.1+Dtol & lice<m & d[lice]>=dPHC-Dtol){
+			if (depthsnow<=0.1+Dtol && lice<m && d[lice]>=dPHC-Dtol){
 				aice = ai_max + (as_min - ai_max)*(d[lice]-dIce)/(dPHC-dIce);
 				a[0]= aice + max(a[0]-aice,0.0)*(depthsnow/0.1);