Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 20225)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 20226)
@@ -9,4 +9,7 @@
 from solve import *
 from gmshplanet import *
+from maskpsl import *
+from gmtmask import *
+from paterson import *
 
 #mesh earth: 
@@ -26,5 +29,6 @@
 
 #elastic loading from love numbers: 
-love = dlmread('../Data/love_numbers_10k.txt')
+
+love = loadtxt('../Data/love_numbers_10k.txt')
 nlov=101
 print love
@@ -37,8 +41,8 @@
 
 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