Changeset 20226

Timestamp:

02/21/16 23:27:58 (9 years ago)

Author:

dlcheng

Message:

CHG: test2002.py syntax fixes, ready for solve validation

File:

• issm/trunk-jpl/test/NightlyRun/test2002.py (modified) (3 diffs)

issm/trunk-jpl/test/NightlyRun/test2002.py ¶

@@ -9 +9 @@
 from solve import *
 from gmshplanet import *
+from maskpsl import *
+from gmtmask import *
+from paterson import *
 
 #mesh earth: 
@@ -26 +29 @@
 
 #elastic loading from love numbers: 
-love = dlmread('../Data/love_numbers_10k.txt')
+
+love = loadtxt('../Data/love_numbers_10k.txt')
 nlov=101
 print love
@@ -37 +41 @@
 
 icemask=ones((md.mesh.numberofvertices,1))
-pos=find(mask==0)
+pos=nonzero(mask==0)
 icemask[pos]=-1
-pos=find(sum(mask(md.mesh.elements),2)<3)
-icemask[md.mesh.elements[pos]]=-1
+pos=nonzero(sum(mask[md.mesh.elements.astype(int)-1],2)<3)
+icemask[md.mesh.elements[pos].astype(int)-1]=-1
 
 md.mask.ice_levelset=icemask