Index: /issm/trunk-jpl/src/m/classes/SMBautoregression.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBautoregression.py	(revision 26810)
+++ /issm/trunk-jpl/src/m/classes/SMBautoregression.py	(revision 26811)
@@ -133,5 +133,5 @@
                 elemsh  = np.zeros((len(indices)))
                 for jj in range(len(indices)):
-                    elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj,:]]])
+                    elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj],:]])
                 temprefelevation[ii] = np.sum(areas[indices]*elemsh)/np.sum(areas[indices])
             if(np.any(templapserate_pos!=0) or np.any(templapserate_neg!=0)):