Index: /issm/trunk/src/m/model/marshall.m
===================================================================
--- /issm/trunk/src/m/model/marshall.m	(revision 7645)
+++ /issm/trunk/src/m/model/marshall.m	(revision 7646)
@@ -82,4 +82,5 @@
 WriteData(fid,md.spctemperature,'Mat','spctemperature');
 WriteData(fid,md.spcthickness,'Mat','spcthickness');
+WriteData(fid,md.spcwatercolumn,'Mat','spcwatercolumn');
 
 WriteData(fid,md.pressureload,'Mat','pressureload');
@@ -95,4 +96,5 @@
 end
 WriteData(fid,md.dhdt,'Mat','dhdt');
+WriteData(fid,md.watercolumn,'Mat','watercolumn');
 
 %Get materials
@@ -151,4 +153,6 @@
 %Thermal parameters
 WriteData(fid,md.beta,'Scalar','beta');
+WriteData(fid,md.hydro_gamma,'Scalar','gamma'); 
+WriteData(fid,md.hydro_kn,'Scalar','kn');
 WriteData(fid,md.meltingpoint,'Scalar','meltingpoint');
 WriteData(fid,md.latentheat,'Scalar','latentheat');
Index: /issm/trunk/src/m/model/process_solve_options.m
===================================================================
--- /issm/trunk/src/m/model/process_solve_options.m	(revision 7645)
+++ /issm/trunk/src/m/model/process_solve_options.m	(revision 7646)
@@ -22,5 +22,5 @@
 if ~ismember(solution_type,[DiagnosticSolutionEnum,PrognosticSolutionEnum,ThermalSolutionEnum,...
 		SteadystateSolutionEnum,ParametersSolutionEnum,Transient2DSolutionEnum,Transient3DSolutionEnum...
-		BalancedthicknessSolutionEnum,BalancedvelocitiesSolutionEnum,BedSlopeSolutionEnum,SurfaceSlopeSolutionEnum,GroundingLineMigration2DSolutionEnum]),
+		BalancedthicknessSolutionEnum,BalancedvelocitiesSolutionEnum,BedSlopeSolutionEnum,SurfaceSlopeSolutionEnum,GroundingLineMigration2DSolutionEnum,HydrologySolutionEnum]),
 	error(['process_solve_options error message: solution_type ' EnumToString(solution_type) ' not supported yet!']);
 end
Index: /issm/trunk/src/m/model/tres.m
===================================================================
--- /issm/trunk/src/m/model/tres.m	(revision 7645)
+++ /issm/trunk/src/m/model/tres.m	(revision 7646)
@@ -90,4 +90,7 @@
 	md.temperature=PatchToVec(md.results.ThermalSolution.Temperature);
 	md.melting_rate=PatchToVec(md.results.ThermalSolution.MeltingRate);
+elseif strcmpi(string,'hydrology'),
+	md.watercolumn=PatchToVec(md.results.HydrologySolution.WaterColumn);
+
 else 
 	error(['tres error message: analysis ' string ' not supported yet!']);