Index: /issm/trunk-jpl/src/c/classes/Vertex.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Vertex.cpp	(revision 17006)
+++ /issm/trunk-jpl/src/c/classes/Vertex.cpp	(revision 17007)
@@ -132,17 +132,24 @@
 
 	/*sigma remains constant. z=bed+sigma*thickness*/
-	if(this->meshtype==Mesh2DverticalEnum){
-		oldy = this->y;
-		newy = bed[this->pid]+sigma*(surface[this->pid] - bed[this->pid]);
-		vely = (newy-oldy)/dt;
-		this->y = newy;
-		vy->SetValue(this->pid,vely,INS_VAL);
-	}
-	else{
-		oldz = this->z;
-		newz = bed[this->pid]+sigma*(surface[this->pid] - bed[this->pid]);
-		velz = (newz-oldz)/dt;
-		this->z = newz;
-		vz->SetValue(this->pid,velz,INS_VAL);
+	switch(this->meshtype){
+		case Mesh2DhorizontalEnum:
+			/*Nothing*/
+			return;
+		case Mesh2DverticalEnum:
+			oldy = this->y;
+			newy = bed[this->pid]+sigma*(surface[this->pid] - bed[this->pid]);
+			vely = (newy-oldy)/dt;
+			this->y = newy;
+			vy->SetValue(this->pid,vely,INS_VAL);
+			return;
+		case Mesh3DEnum:
+			oldz = this->z;
+			newz = bed[this->pid]+sigma*(surface[this->pid] - bed[this->pid]);
+			velz = (newz-oldz)/dt;
+			this->z = newz;
+			vz->SetValue(this->pid,velz,INS_VAL);
+			return;
+		default:
+			_error_("not implemented");
 	}
 }