Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21828)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21829)
@@ -2419,11 +2419,12 @@
 		element->GetVerticesCoordinates(&xyz_list);
 		for(int i=0;i<numvertices;i++){
-			x = *(xyz_list+3*i+0);
-			y = *(xyz_list+3*i+1);
+			x = xyz_list[3*i+0];
+			y = xyz_list[3*i+1];
 			
 			bx=-150.-728.8*std::pow(x/300000.,2)+343.91*std::pow(x/300000.,4)-50.57*std::pow(x/300000.,6);
 			by=500./(1.+std::exp((-2./4000.)*(y-80000./2.-24000.)))+500./(1.+std::exp((2./4000.)*(y-80000./2.+24000.)));
 			
-			r[i]=std::max(bx+by,-720.);
+			r[i] = -720.;
+			if(bx+by>-720.) r[i] = bx+by;
 		}	
 
@@ -2437,6 +2438,4 @@
    /*Delete*/
    xDelete<IssmDouble>(r);
-	
-	return;
 }
 /*}}}*/