Changeset 26816

Timestamp:

01/25/22 12:28:42 (3 years ago)

Author:

vverjans

Message:

BUG: working on Python syntax wrt revision 26810

File:

• issm/trunk-jpl/src/m/classes/SMBautoregression.py (modified) (1 diff)

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

@@ -133 +133 @@
                 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])
+                    #elemsh[jj] = np.mean(md.geometry.surface[md.mesh.elements[indices[jj],:]-1])
+                    sumval = 0
+                    cntval = 0
+                    elemindices = md.mesh.elements[indices[jj],:]-1
+                    for kk in range(len(elemindices)):
+                        sumval += md.geometry.surface[elemindices[kk]]
+                        cntval += 1
+                    elemsh[jj] = sumval/cntval
                 temprefelevation[ii] = np.sum(areas[indices]*elemsh)/np.sum(areas[indices])
+                print(indices)
+                print(areas[indices])
+                print(elemsh)
             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')