Index: /issm/trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 25180)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 25181)
@@ -3,5 +3,5 @@
 %mesh earth:
 md=model;
-md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',700.); %500 km resolution mesh
+md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',700.); %700 km resolution mesh
 
 %parameterize solidearth solution:
Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 25180)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 25181)
@@ -25,6 +25,6 @@
 
 #antarctica
-late = np.sum(md.mesh.lat[md.mesh.elements - 1], axis=1) / 3
-longe = np.sum(md.mesh.long[md.mesh.elements - 1], axis=1) / 3
+late = md.mesh.lat[md.mesh.elements - 1].sum(axis=1) / 3
+longe = md.mesh.long[md.mesh.elements - 1].sum(axis=1) / 3
 pos = np.where(late < -80)[0]
 md.solidearth.surfaceload.icethicknesschange[pos] = -100
@@ -42,12 +42,8 @@
 pos = np.where(mask == 0)[0]
 icemask[pos] = -1
-pos = np.where(np.sum(mask[md.mesh.elements - 1], axis=1) < 3)
+pos = np.where(mask[md.mesh.elements - 1].sum(axis=1) < 3)[0]
 icemask[md.mesh.elements[pos, :] - 1] = -1
 md.mask.ice_levelset = icemask
 md.mask.ocean_levelset = -icemask
-
-#make sure that the ice level set is all inclusive:
-md.mask.land_levelset = np.zeros((md.mesh.numberofvertices, 1))
-md.mask.ocean_levelset = -np.ones((md.mesh.numberofvertices, 1))
 
 #make sure that the elements that have loads are fully grounded
@@ -57,5 +53,5 @@
 #make sure wherever there is an ice load, that the mask is set to ice:
 #pos = np.nonzero(md.solidearth.surfaceload.icethicknesschange)[0] # Do we need to do this twice?
-md.mask.ice_levelset[md.mesh.elements[pos, :] - 1] = 1
+md.mask.ice_levelset[md.mesh.elements[pos, :] - 1] = -1
 # }}}
 
Index: /issm/trunk-jpl/test/NightlyRun/test2010.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 25180)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 25181)
@@ -41,6 +41,6 @@
 mask = gmtmask(md.mesh.lat, md.mesh.long)
 icemask = np.ones((md.mesh.numberofvertices, 1))
-pos = np.where(mask == 0)
-icemask[pos[0]] = -1
+pos = np.where(mask == 0)[0]
+icemask[pos] = -1
 pos = np.where(mask[md.mesh.elements - 1].sum(axis=1) < 3)[0]
 icemask[md.mesh.elements[pos, :] - 1] = -1