Index: /issm/trunk-jpl/test/NightlyRun/test120.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test120.m	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test120.m	(revision 16041)
@@ -7,4 +7,5 @@
 md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md=solve(md,EnthalpySolutionEnum());
 
Index: /issm/trunk-jpl/test/NightlyRun/test120.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test120.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test120.py	(revision 16041)
@@ -17,4 +17,5 @@
 md.cluster=generic('name',oshostname(),'np',3)
 md.initialization.waterfraction=numpy.zeros(md.mesh.numberofvertices)
+md.initialization.watercolumn=numpy.zeros(md.mesh.numberofvertices)
 md=solve(md,EnthalpySolutionEnum())
 
Index: /issm/trunk-jpl/test/NightlyRun/test121.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test121.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test121.py	(revision 16041)
@@ -15,4 +15,5 @@
 md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md.transient.isstressbalance=False
 md.transient.ismasstransport=False
Index: /issm/trunk-jpl/test/NightlyRun/test122.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test122.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test122.py	(revision 16041)
@@ -14,4 +14,5 @@
 md=setflowequation(md,'HO','all')
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md.thermal.isenthalpy=1
 md.thermal.stabilization=2
Index: /issm/trunk-jpl/test/NightlyRun/test325.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test325.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test325.py	(revision 16041)
@@ -17,4 +17,5 @@
 md.cluster=generic('name',oshostname(),'np',3)
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md=solve(md,EnthalpySolutionEnum())
 
Index: /issm/trunk-jpl/test/NightlyRun/test326.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test326.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test326.py	(revision 16041)
@@ -16,4 +16,5 @@
 md.cluster=generic('name',oshostname(),'np',3)
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md.transient.isstressbalance=False
 md.transient.ismasstransport=False
Index: /issm/trunk-jpl/test/NightlyRun/test327.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test327.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test327.py	(revision 16041)
@@ -16,4 +16,5 @@
 md.cluster=generic('name',oshostname(),'np',3)
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md.initialization.temperature[:]=272.
 md.thermal.spctemperature[numpy.nonzero(md.mesh.vertexonsurface)[0]]=272.
Index: /issm/trunk-jpl/test/NightlyRun/test431.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test431.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test431.py	(revision 16041)
@@ -18,4 +18,5 @@
 md.thermal.isenthalpy=1
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md=solve(md,SteadystateSolutionEnum())
 
Index: /issm/trunk-jpl/test/NightlyRun/test432.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test432.py	(revision 16040)
+++ /issm/trunk-jpl/test/NightlyRun/test432.py	(revision 16041)
@@ -18,4 +18,5 @@
 md.thermal.isenthalpy=1
 md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
+md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
 md=solve(md,SteadystateSolutionEnum())