source: issm/oecreview/Archive/26740-27031/ISSM-26811-26812.diff

Last change on this file was 27032, checked in by Mathieu Morlighem, 3 years ago

CHG: added 26740-27031
File size: 932 bytes
RevLine 
[27032]1Index: ../trunk-jpl/src/m/classes/SMBautoregression.py
2===================================================================
3--- ../trunk-jpl/src/m/classes/SMBautoregression.py (revision 26811)
4+++ ../trunk-jpl/src/m/classes/SMBautoregression.py (revision 26812)
5@@ -132,7 +132,7 @@
6 indices = np.where(md.smb.basin_id==ii)[0]
7 elemsh = np.zeros((len(indices)))
8 for jj in range(len(indices)):
9- elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj],:]])
10+ elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj],:]-1])
11 temprefelevation[ii] = np.sum(areas[indices]*elemsh)/np.sum(areas[indices])
12 if(np.any(templapserate_pos!=0) or np.any(templapserate_neg!=0)):
13 print(' smb.refelevation not specified: Reference elevations set to mean surface elevation of basins')
Note: See TracBrowser for help on using the repository browser.