Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 17452)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 17453)
@@ -108,5 +108,4 @@
 	case 334, name='SquareSheetConstrainedExtrudedHydrologyDC';
 	case 335, name='SquareSheetExtrudedHydrologyDCTwoLayers';
-	case 336, name='SquareSheetLevelsetSSAMasstransp2d';
 	case 401, name='SquareSheetShelfStressSSA2d';
 	case 402, name='SquareSheetShelfStressSSA3d';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 17452)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 17453)
@@ -107,5 +107,4 @@
 	332  : 'SquareSheetConstrainedHydrologyDC',
 	333  : 'SquareSheetHydrologyDCTwoLayers',
-	336  : 'SquareSheetLevelsetSSAMasstransp2d',
 	401  : 'SquareSheetShelfStressSSA2d',
 	402  : 'SquareSheetShelfStressSSA3d',
Index: sm/trunk-jpl/test/NightlyRun/test336.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test336.m	(revision 17452)
+++ 	(revision )
@@ -1,50 +1,0 @@
-md=triangle(model(),'../Exp/Square.exp',100000.);
-md=setmask(md,'','');
-md=parameterize(md,'../Par/SquareSheetConstrained.par');
-md=setflowequation(md,'SSA','all');
-md.cluster=generic('name',oshostname(),'np',3);
-
-% parameters
-md.transient.isstressbalance=1;
-md.transient.ismasstransport=1;
-md.transient.isthermal=0;
-md.transient.isgroundingline=0;
-md.transient.isgia=0;
-md.transient.islevelset=1;
-
-% init levelset function
-ymed = (max(md.mesh.y)+min(md.mesh.y))/2;
-md.mask.ice_levelset=md.mesh.y - ymed;
-
-md=solve(md,TransientSolutionEnum())
-
-%Fields and tolerances to track changes
-field_names     ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1',...
-                  'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2',...
-                  'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3'};
-field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
-		1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,...
-		1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
-field_values={...
-	md.results.TransientSolution(1).Vx,...
-	md.results.TransientSolution(1).Vy,...
-	md.results.TransientSolution(1).Vel,...
-	md.results.TransientSolution(1).Pressure,...
-	md.results.TransientSolution(1).Thickness,...
-	md.results.TransientSolution(1).Surface,...
-	md.results.TransientSolution(1).MaskIceLevelset,...
-	md.results.TransientSolution(2).Vx,...
-	md.results.TransientSolution(2).Vy,...
-	md.results.TransientSolution(2).Vel,...
-	md.results.TransientSolution(2).Pressure,...
-	md.results.TransientSolution(2).Thickness,...
-	md.results.TransientSolution(2).Surface,...
-	md.results.TransientSolution(2).MaskIceLevelset,...
-	md.results.TransientSolution(3).Vx,...
-	md.results.TransientSolution(3).Vy,...
-	md.results.TransientSolution(3).Vel,...
-	md.results.TransientSolution(3).Pressure,...
-	md.results.TransientSolution(3).Thickness,...
-	md.results.TransientSolution(3).Surface,...
-	md.results.TransientSolution(3).MaskIceLevelset...
-	};
Index: sm/trunk-jpl/test/NightlyRun/test336.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test336.py	(revision 17452)
+++ 	(revision )
@@ -1,60 +1,0 @@
-from model import *
-from triangle import *
-from setmask import *
-from parameterize import *
-from setflowequation import *
-from EnumDefinitions import *
-from solve import *
-from MatlabFuncs import *
-
-md=triangle(model(),'../Exp/Square.exp',100000.)
-md=setmask(md,'','')
-md=parameterize(md,'../Par/SquareSheetConstrained.py') 
-md=setflowequation(md,'SSA','all')
-md.cluster=generic('name',oshostname(),'np',3)
-
-# parameters
-md.transient.isstressbalance=True
-md.transient.ismasstransport=True
-md.transient.isthermal=False
-md.transient.isgroundingline=False
-md.transient.isgia=False
-md.transient.islevelset=True
-
-# init levelset function
-xmin, xmax=min(md.mesh.x), max(md.mesh.x)
-ymin, ymax=min(md.mesh.y), max(md.mesh.y)
-xmed, ymed=(xmax+xmin)/2.,(ymax+ymin)/2.  
-
-md.mask.ice_levelset=md.mesh.y.reshape(-1,1) - ymed
-
-md=solve(md,TransientSolutionEnum())
-
-#Fields and tolerances to track changes
-field_names     =['Vx','Vy','Vel','Pressure','Thickness','Surface','MaskIceLevelset']
-field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
-		1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
-		1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
-field_values=[\
-	md.results.TransientSolution[0].Vx,\
-	md.results.TransientSolution[0].Vy,\
-	md.results.TransientSolution[0].Vel,\
-	md.results.TransientSolution[0].Pressure,\
-	md.results.TransientSolution[0].Thickness,\
-	md.results.TransientSolution[0].Surface,\
-	md.results.TransientSolution[0].MaskIceLevelset,\
-	md.results.TransientSolution[1].Vx,\
-	md.results.TransientSolution[1].Vy,\
-	md.results.TransientSolution[1].Vel,\
-	md.results.TransientSolution[1].Pressure,\
-	md.results.TransientSolution[1].Thickness,\
-	md.results.TransientSolution[1].Surface,\
-	md.results.TransientSolution[1].MaskIceLevelset,\
-	md.results.TransientSolution[2].Vx,\
-	md.results.TransientSolution[2].Vy,\
-	md.results.TransientSolution[2].Vel,\
-	md.results.TransientSolution[2].Pressure,\
-	md.results.TransientSolution[2].Thickness,\
-	md.results.TransientSolution[2].Surface,\
-	md.results.TransientSolution[2].MaskIceLevelset,\
-	]