Index: /issm/trunk-jpl/test/NightlyRun/test243.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test243.py	(revision 22868)
+++ /issm/trunk-jpl/test/NightlyRun/test243.py	(revision 22869)
@@ -1,6 +1,4 @@
 #Test Name: SquareShelfSMBGemb
 import numpy as np
-import scipy as sci
-import sci.io as spio
 from model import *
 from socket import gethostname
@@ -24,5 +22,5 @@
 
 #load hourly surface forcing date from 1979 to 2009:
-inputs = spio.loadmat('../Data/gemb_input.mat',squeeze_me=True)
+inputs = np.load('../Data/gemb_input.npy').item()
 
 #setup the inputs:
Index: /issm/trunk-jpl/test/NightlyRun/test701.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test701.py	(revision 22868)
+++ /issm/trunk-jpl/test/NightlyRun/test701.py	(revision 22869)
@@ -1,6 +1,4 @@
 #Test Name: FlowbandFSshelf
 import numpy as np
-import scipy as sci
-from sci.interpolate import interp1d
 from model import *
 from solve import *
@@ -16,6 +14,6 @@
 
 #Geometry	    #interp1d returns a function to be called on md.mesh.x
-md.geometry.surface = interp1d(x,b+h)(md.mesh.x)
-md.geometry.base = interp1d(x,b)(md.mesh.x)
+md.geometry.surface = np.interp(md.mesh.x,x,b+h)
+md.geometry.base = np.interp(md.mesh.x,x,b)
 md.geometry.thickness = md.geometry.surface-md.geometry.base
 
Index: /issm/trunk-jpl/test/NightlyRun/test703.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test703.py	(revision 22868)
+++ /issm/trunk-jpl/test/NightlyRun/test703.py	(revision 22869)
@@ -1,6 +1,4 @@
 #Test Name: FlowbandFSsheetshelfTrans
 import numpy as np
-import scipy as sci
-from sci.interpolate import interp1d
 from model import  *
 from setflowequation import  *
@@ -35,6 +33,6 @@
 
 #geometry
-md.geometry.surface = interp1d(x,b+h)(md.mesh.x)
-md.geometry.base = interp1d(x,b)(md.mesh.x)
+md.geometry.surface = np.interp(md.mesh.x,x,b+h)
+md.geometry.base = np.interp(md.mesh.x,x,b)
 md.geometry.thickness = md.geometry.surface - md.geometry.base