Index: ../trunk-jpl/src/m/classes/SMBautoregression.py
===================================================================
--- ../trunk-jpl/src/m/classes/SMBautoregression.py	(revision 26811)
+++ ../trunk-jpl/src/m/classes/SMBautoregression.py	(revision 26812)
@@ -132,7 +132,7 @@
                 indices = np.where(md.smb.basin_id==ii)[0]
                 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],:]-1])
                 temprefelevation[ii] = np.sum(areas[indices]*elemsh)/np.sum(areas[indices])
             if(np.any(templapserate_pos!=0) or np.any(templapserate_neg!=0)):
                 print('      smb.refelevation not specified: Reference elevations set to mean surface elevation of basins')