Index: ../trunk-jpl/test/NightlyRun/test431.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test431.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test431.m	(revision 16039)
@@ -7,6 +7,7 @@
 md.timestepping.time_step=0.;
 md.thermal.isenthalpy=1;
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md=solve(md,SteadystateSolutionEnum());
 
 %Fields and tolerances to track changes
Index: ../trunk-jpl/test/NightlyRun/test121.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test121.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test121.m	(revision 16039)
@@ -5,6 +5,7 @@
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md.transient.isstressbalance=0;
 md.transient.ismasstransport=0;
 md.transient.isthermal=1;
Index: ../trunk-jpl/test/NightlyRun/test325.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test325.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test325.m	(revision 16039)
@@ -6,6 +6,7 @@
 md.timestepping.time_step=0.;
 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());
 
 %Fields and tolerances to track changes
Index: ../trunk-jpl/test/NightlyRun/test432.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test432.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test432.m	(revision 16039)
@@ -7,6 +7,7 @@
 md.timestepping.time_step=0.;
 md.thermal.isenthalpy=1;
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md=solve(md,SteadystateSolutionEnum());
 
 %Fields and tolerances to track changes
Index: ../trunk-jpl/test/NightlyRun/test122.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test122.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test122.m	(revision 16039)
@@ -4,6 +4,7 @@
 md=extrude(md,3,1.);
 md=setflowequation(md,'HO','all');
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md.thermal.isenthalpy=1;
 md.thermal.stabilization=2;
 md.cluster=generic('name',oshostname(),'np',3);
Index: ../trunk-jpl/test/NightlyRun/test326.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test326.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test326.m	(revision 16039)
@@ -5,6 +5,7 @@
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md.transient.isstressbalance=0;
 md.transient.ismasstransport=0;
 md.transient.isthermal=1;
Index: ../trunk-jpl/test/NightlyRun/test327.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test327.m	(revision 16038)
+++ ../trunk-jpl/test/NightlyRun/test327.m	(revision 16039)
@@ -5,6 +5,7 @@
 md=setflowequation(md,'HO','all');
 md.cluster=generic('name',oshostname(),'np',3);
 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
+md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
 md.initialization.temperature(:)=272.;
 md.thermal.spctemperature(find(md.mesh.vertexonsurface))=272.;
 md.thermal.isenthalpy=1;