Index: /issm/trunk-jpl/test/NightlyRun/test330.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test330.m	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test330.m	(revision 22900)
@@ -16,4 +16,5 @@
 md.hydrology.isefficientlayer=0;
 md.hydrology.sedimentlimit_flag=0;
+md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
 md.hydrology.rel_tol=1.0e-6;
 md.hydrology.penalty_lock=0;
Index: /issm/trunk-jpl/test/NightlyRun/test330.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test330.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test330.py	(revision 22900)
@@ -28,4 +28,5 @@
 md.hydrology.isefficientlayer=0
 md.hydrology.sedimentlimit_flag=0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.hydrology.rel_tol=1.0e-6
 md.hydrology.penalty_lock=0
Index: /issm/trunk-jpl/test/NightlyRun/test3300.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3300.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test3300.py	(revision 22900)
@@ -26,4 +26,5 @@
 md.hydrology.sedimentlimit_flag=1
 md.hydrology.sedimentlimit=400.0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.hydrology.sediment_thickness=20.0
 md.initialization.sediment_head=np.zeros((md.mesh.numberofvertices))
@@ -58,14 +59,14 @@
 md=solve(md,'Transient')
 
-sedvol=np.zeros(4000)
-eplvol=np.zeros(4000)
-totvol=np.zeros(4001)
-time=np.arange(0.002,8.001,0.002)
-store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)
-sedstore=20.0*store
-for i in range(0,4000):
-	sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore
-	eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThicknessHydrostep)
-	totvol[i+1]=totvol[i]+md.basalforcings.groundedice_melting_rate[0,i]*0.002
+# sedvol=np.zeros(4000)
+# eplvol=np.zeros(4000)
+# totvol=np.zeros(4001)
+# time=np.arange(0.002,8.001,0.002)
+# store=md.constants.g*md.hydrology.sediment_porosity*md.materials.rho_freshwater*((md.hydrology.sediment_compressibility/md.hydrology.sediment_porosity)+md.hydrology.water_compressibility)
+# sedstore=20.0*store
+# for i in range(0,4000):
+# 	sedvol[i]=np.mean(md.results.TransientSolution[i].SedimentHead)*sedstore
+# 	eplvol[i]=np.mean(md.results.TransientSolution[i].EplHead)*store*np.mean(md.results.TransientSolution[i].HydrologydcEplThicknessHydrostep)
+# 	totvol[i+1]=totvol[i]+md.basalforcings.groundedice_melting_rate[0,i]*0.002
 
 field_names     =['SedimentWaterHead5','EplWaterHead5','SedimentWaterHead40','EplWaterHead40']
Index: /issm/trunk-jpl/test/NightlyRun/test332.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test332.m	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test332.m	(revision 22900)
@@ -12,4 +12,5 @@
 md.hydrology.sedimentlimit_flag=1;
 md.hydrology.sedimentlimit=8000.0;
+md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/NightlyRun/test332.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test332.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test332.py	(revision 22900)
@@ -26,4 +26,5 @@
 
 md.hydrology.isefficientlayer=0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.hydrology.sedimentlimit_flag=1
 md.hydrology.sedimentlimit=8000.0
Index: /issm/trunk-jpl/test/NightlyRun/test333.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test333.m	(revision 22900)
@@ -18,4 +18,5 @@
 md.basalforcings.floatingice_melting_rate = 0.0*ones(md.mesh.numberofvertices,1);
 md.hydrology.sediment_transmitivity=3*ones(md.mesh.numberofvertices,1);
+md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
 
 md.initialization.epl_head=0.0*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/NightlyRun/test333.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test333.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test333.py	(revision 22900)
@@ -26,4 +26,5 @@
 md.hydrology.sedimentlimit_flag=1
 md.hydrology.sedimentlimit=800.0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.initialization.sediment_head=np.zeros((md.mesh.numberofvertices))
 md.hydrology.spcsediment_head=np.nan*np.ones((md.mesh.numberofvertices))
Index: /issm/trunk-jpl/test/NightlyRun/test334.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test334.m	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test334.m	(revision 22900)
@@ -12,5 +12,7 @@
 md.hydrology.sedimentlimit_flag=1;
 md.hydrology.sedimentlimit=8000.0;
-md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
+md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices, ...
+                                         1);
+md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
 pos=find(md.mesh.y==0);
Index: /issm/trunk-jpl/test/NightlyRun/test334.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test334.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test334.py	(revision 22900)
@@ -26,4 +26,5 @@
 md.hydrology.sedimentlimit_flag=1
 md.hydrology.sedimentlimit=8000.0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.initialization.sediment_head=np.zeros((md.mesh.numberofvertices))
 md.hydrology.spcsediment_head=np.nan*np.ones((md.mesh.numberofvertices))
Index: /issm/trunk-jpl/test/NightlyRun/test335.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test335.m	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test335.m	(revision 22900)
@@ -14,4 +14,5 @@
 md.hydrology.sedimentlimit=800.0;
 md.initialization.sediment_head=0.0*ones(md.mesh.numberofvertices,1);
+md.hydrology.mask_thawed_node=ones(md.mesh.numberofvertices,1);
 md.hydrology.spcsediment_head=NaN*ones(md.mesh.numberofvertices,1);
 md.basalforcings.groundedice_melting_rate = 2.0*ones(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/test/NightlyRun/test335.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test335.py	(revision 22899)
+++ /issm/trunk-jpl/test/NightlyRun/test335.py	(revision 22900)
@@ -27,4 +27,5 @@
 md.hydrology.sedimentlimit_flag=1
 md.hydrology.sedimentlimit=800.0
+md.hydrology.mask_thawed_node=np.ones((md.mesh.numberofvertices))
 md.initialization.sediment_head=np.zeros((md.mesh.numberofvertices))
 md.hydrology.spcsediment_head=np.nan*np.ones((md.mesh.numberofvertices))