Index: /issm/trunk-jpl/test/NightlyRun/test2002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24889)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.py	(revision 24890)
@@ -40,7 +40,5 @@
 
 #mask:
-md.mask = maskpsl()  # use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
 mask = gmtmask(md.mesh.lat, md.mesh.long)
-
 icemask = np.ones((md.mesh.numberofvertices))
 pos = np.where(mask == 0)[0]
@@ -49,7 +47,5 @@
 icemask[md.mesh.elements[pos, :].astype(int) - 1] = -1
 md.mask.ice_levelset = icemask
-md.mask.ocean_levelset = np.zeros((md.mesh.numberofvertices))
-pos = np.where(md.mask.ice_levelset == 1)
-md.mask.ocean_levelset[pos] = 1
+md.mask.ocean_levelset = -icemask
 
 #make sure that the ice level set is all inclusive:
Index: /issm/trunk-jpl/test/NightlyRun/test2003.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24889)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.py	(revision 24890)
@@ -42,7 +42,5 @@
 
 #mask:  {{{
-md.mask = maskpsl()  # use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
 mask = gmtmask(md.mesh.lat, md.mesh.long)
-
 icemask = np.ones((md.mesh.numberofvertices, ))
 pos = np.where(mask == 0)
@@ -52,7 +50,5 @@
 icemask[md.mesh.elements[pos, :] - 1] = -1
 md.mask.ice_levelset = icemask
-md.mask.ocean_levelset = np.zeros((md.mesh.numberofvertices, ))
-pos = np.where(md.mask.ice_levelset == 1)
-md.mask.ocean_levelset[pos] = 1
+md.mask.ocean_levelset = -icemask
 
 #make sure that the ice level set is all inclusive:
Index: /issm/trunk-jpl/test/NightlyRun/test2010.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24889)
+++ /issm/trunk-jpl/test/NightlyRun/test2010.py	(revision 24890)
@@ -43,7 +43,5 @@
 #}}}
 #mask:  {{{
-md.mask = maskpsl()  # use maskpsl class (instead of mask) to store the ocean function as a ocean_levelset
 mask = gmtmask(md.mesh.lat, md.mesh.long)
-
 icemask = np.ones((md.mesh.numberofvertices, ))
 pos = np.where(mask == 0)
@@ -52,7 +50,5 @@
 icemask[md.mesh.elements[pos, :] - 1] = -1
 md.mask.ice_levelset = icemask
-md.mask.ocean_levelset = np.zeros((md.mesh.numberofvertices, ))
-pos = np.where(md.mask.ice_levelset == 1)
-md.mask.ocean_levelset[pos] = 1
+md.mask.ocean_levelset = -icemask
 
 #make sure that the ice level set is all inclusive: