Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16784)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16785)
@@ -8611,6 +8611,4 @@
 	IssmDouble*    basis  = xNew<IssmDouble>(numnodes);
 	IssmDouble*    dbasis = xNew<IssmDouble>(numnodes*2);
-	IssmDouble*    HNx    = xNew<IssmDouble>(numnodes);
-	IssmDouble*    HNy    = xNew<IssmDouble>(numnodes);
 	IssmDouble*    H      = xNew<IssmDouble>(numnodes);
 	IssmDouble*    Nx     = xNew<IssmDouble>(numnodes);
@@ -8651,8 +8649,8 @@
 		GetNodalFunctionsDerivatives(dbasis,&xyz_list[0][0],gauss);
 
-		TriaRef::GetInputDerivativeValue(&dhnx[0],HNx,&xyz_list[0][0],gauss);
-		TriaRef::GetInputDerivativeValue(&dhny[0],HNy,&xyz_list[0][0],gauss);
-		TriaRef::GetInputValue(&hnx,HNx,gauss);
-		TriaRef::GetInputValue(&hny,HNy,gauss);
+		TriaRef::GetInputDerivativeValue(&dhnx[0],Nx,&xyz_list[0][0],gauss);
+		TriaRef::GetInputDerivativeValue(&dhny[0],Ny,&xyz_list[0][0],gauss);
+		TriaRef::GetInputValue(&hnx,Nx,gauss);
+		TriaRef::GetInputValue(&hny,Ny,gauss);
 
 		gamma=h/(2.*thickness+1.e-10);
@@ -8671,6 +8669,4 @@
 	xDelete<IssmDouble>(Nx);
 	xDelete<IssmDouble>(Ny);
-	xDelete<IssmDouble>(HNx);
-	xDelete<IssmDouble>(HNy);
 	delete gauss;
 	return pe;