Changeset 26813

Timestamp:

01/25/22 10:10:20 (3 years ago)

Author:

vverjans

Message:

BUG: fixing Python syntax wrt revision 26810

File:

• issm/trunk-jpl/src/m/classes/SMBautoregression.py (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/SMBautoregression.py

@@ -119 +119 @@
         templapserate_pos = np.copy(md.smb.lapserate_pos)
         templapserate_neg = np.copy(md.smb.lapserate_neg)
-        temprefelevation  = np.copy(md.smb.lapserate_neg)
+        temprefelevation  = np.copy(md.smb.refelevation)
         if(np.any(np.isnan(md.smb.lapserate_pos))):
             templapserate_pos = np.zeros((md.smb.num_basins))
@@ -130 +130 @@
             areas = GetAreas(md.mesh.elements, md.mesh.x, md.mesh.y)
             for ii in range(int(md.smb.num_basins)):
-                indices = np.where(md.smb.basin_id==ii)[0]
+                indices = np.where(md.smb.basin_id==ii+1)[0]
                 elemsh  = np.zeros((len(indices)))
                 for jj in range(len(indices)):
                     elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj],:]-1])
+                    if(ii==0 and jj==0):
+                        print(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)):