Index: /issm/trunk-jpl/test/NightlyRun/runme.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 21613)
+++ /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 21614)
@@ -53,5 +53,5 @@
 	#Process options
 	#GET benchmark {{{
-	if not benchmark in ['all','nightly','ismip','eismint','thermal','mesh','validation','tranforcing','adolc','slr']:
+	if not benchmark in ['all','nightly','ismip','eismint','thermal','mesh','validation','tranforcing','adolc','slr','referential']:
 		print("runme warning: benchmark '{}' not supported, defaulting to test 'nightly'.".format(benchmark))
 		benchmark='nightly'
Index: /issm/trunk-jpl/test/NightlyRun/test1601.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1601.py	(revision 21613)
+++ /issm/trunk-jpl/test/NightlyRun/test1601.py	(revision 21614)
@@ -18,5 +18,5 @@
 md.initialization.vx[:]=0.
 md.initialization.vy[:]=0.
-md.initialization.vel=np.zeros_like(md.initialization.vx)
+md.initialization.vel=np.zeros((md.mesh.numberofvertices))
 
 md.cluster=generic('name',gethostname(),'np',2)
@@ -30,6 +30,7 @@
 md.mesh.y=np.sin(theta)*x+np.cos(theta)*y
 
-md.stressbalance.referential[:,0:3]=np.tile([np.cos(theta),np.sin(theta),0],(md.mesh.numberofvertices))
-md.stressbalance.referential[:,3: ]=np.tile([0,0,1],(md.mesh.numberofvertices))
+rotation_array=np.array([np.cos(theta),np.sin(theta),0])
+md.stressbalance.referential[:,0:3]=(np.tile(rotation_array,(md.mesh.numberofvertices,1)))
+md.stressbalance.referential[:,3: ]=np.tile([0,0,1],(md.mesh.numberofvertices,1))
 md=solve(md,'Stressbalance')
 vel1=md.results.StressbalanceSolution.Vel
@@ -39,8 +40,8 @@
 
 #Now, put CS back to normal except on the side where the spc are applied
-pos=np.nonzero(np.logical_or.reduce(x==0.,x==1000000.))[0]
+pos=np.where(np.logical_or(x==0.,x==1000000.))[0]
 md.stressbalance.referential[:]=np.nan
-md.stressbalance.referential[pos,0:3]=np.tile([np.cos(theta),np.sin(theta),0],(len(pos)))
-md.stressbalance.referential[pos,3: ]=np.tile([0,0,1],(len(pos)))
+md.stressbalance.referential[pos,0:3]=np.tile([np.cos(theta),np.sin(theta),0],(len(pos),1))
+md.stressbalance.referential[pos,3: ]=np.tile([0,0,1],(len(pos),1))
 md=solve(md,'Stressbalance')
 vel2=md.results.StressbalanceSolution.Vel
@@ -52,6 +53,3 @@
 field_names     =['vel1','vel2']
 field_tolerances=[1e-11,1e-11]
-field_values=[\
-	vel1, \
-	vel2, \
-	]
+field_values=[vel1,vel2]
Index: /issm/trunk-jpl/test/NightlyRun/test1602.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1602.py	(revision 21613)
+++ /issm/trunk-jpl/test/NightlyRun/test1602.py	(revision 21614)
@@ -31,6 +31,6 @@
 md.mesh.y=np.sin(theta)*x+np.cos(theta)*y
 
-md.stressbalance.referential[:,0:3]=np.tile([np.cos(theta),np.sin(theta),0],(md.mesh.numberofvertices))
-md.stressbalance.referential[:,3: ]=np.tile([0,0,1],(md.mesh.numberofvertices))
+md.stressbalance.referential[:,0:3]=np.tile([np.cos(theta),np.sin(theta),0],(md.mesh.numberofvertices,1))
+md.stressbalance.referential[:,3: ]=np.tile([0,0,1],(md.mesh.numberofvertices,1))
 md=solve(md,'Stressbalance')
 vel1=md.results.StressbalanceSolution.Vel
@@ -42,5 +42,3 @@
 field_names     =['vel1']
 field_tolerances=[1e-9]
-field_values=[\
-	vel1, \
-	]
+field_values=[vel1]
Index: /issm/trunk-jpl/test/NightlyRun/test201.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test201.py	(revision 21613)
+++ /issm/trunk-jpl/test/NightlyRun/test201.py	(revision 21614)
@@ -20,8 +20,6 @@
 field_names     =['Vx','Vy','Vel','Pressure']
 field_tolerances=[1e-13,1e-13,1e-13,1e-13]
-field_values=[\
-	md.results.StressbalanceSolution.Vx,\
-	md.results.StressbalanceSolution.Vy,\
-	md.results.StressbalanceSolution.Vel,\
-	md.results.StressbalanceSolution.Pressure,\
-	]
+field_values=[md.results.StressbalanceSolution.Vx,
+							md.results.StressbalanceSolution.Vy,
+							md.results.StressbalanceSolution.Vel,
+							md.results.StressbalanceSolution.Pressure]