Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 20162)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 20163)
@@ -26,7 +26,7 @@
 love = dlmread('../Data/love_numbers_10k.txt')
 nlov=101
-md.slr.love_h = love(1:nlov,2)  # radial displacement (height) 
-md.slr.love_k = love(1:nlov,3)  # gravitational potential (phi) 
-
+print love
+md.slr.love_h = love[1][1:nlov]  # radial displacement (height) 
+md.slr.love_k = love[2][1:nlov]  # gravitational potential (phi) 
 #}}}
 #mask:  {{{
@@ -36,12 +36,12 @@
 icemask=ones(md.mesh.numberofvertices,1)
 pos=find(mask==0)
-icemask(pos)=-1
+icemask[pos]=-1
 pos=find(sum(mask(md.mesh.elements),2)<3)
-icemask(md.mesh.elements(pos))=-1
+icemask[md.mesh.elements[pos]]=-1
 
 md.mask.ice_levelset=icemask
 md.mask.ocean_levelset=zeros(md.mesh.numberofvertices,1)
 pos=numpy.nonzero(md.mask.ice_levelset==1)
-md.mask.ocean_levelset(pos)=1
+md.mask.ocean_levelset[pos]=1
 
 #make sure that the ice level set is all inclusive: