Index: /issm/trunk-jpl/src/m/classes/hydrologyglads.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologyglads.m	(revision 23935)
+++ /issm/trunk-jpl/src/m/classes/hydrologyglads.m	(revision 23935)
@@ -0,0 +1,92 @@
+%HYDROLOGYGLADS class definition
+%
+%   Usage:
+%      hydrologyglads=hydrologyglads();
+
+classdef hydrologyglads
+	properties (SetAccess=public) 
+		%Sheet
+		pressure_melt_coefficient = 0.; %c_t
+
+		%Channels
+		ischannels          = 0;
+
+		%Other
+		requested_outputs = {};
+	end
+	methods
+		function self = hydrologyglads(varargin) % {{{
+			switch nargin
+				case 0
+					self=setdefaultparameters(self);
+				case 1
+					self=structtoobj(self,varargin{1});
+				otherwise
+					error('constructor not supported');
+			end
+		end % }}}
+		function list = defaultoutputs(self,md) % {{{
+			list = {};
+		end % }}}    
+
+		function self = setdefaultparameters(self) % {{{
+
+			%Sheet parameters
+			self.pressure_melt_coefficient = 7.5e-8; %K/Pa (See table 1 in Erder et al. 2013)
+
+			%Channel parameters
+			self.ischannels=false;
+
+			%Other
+			self.requested_outputs={'default'};
+		end % }}}
+		function md = checkconsistency(self,md,solution,analyses) % {{{
+
+			%Early return
+			if ~ismember('HydrologyGladsAnalysis',analyses)
+				return;
+			end
+
+			%Sheet
+			md = checkfield(md,'fieldname','hydrology.pressure_melt_coefficient','numel',[1],'>=',0);	
+
+			%Channels
+			md = checkfield(md,'fieldname','hydrology.ischannels','numel',[1],'values',[0 1]);
+
+			%Other
+			md = checkfield(md,'fieldname','hydrology.requested_outputs','stringrow',1);
+		end % }}}
+		function disp(self) % {{{
+			disp(sprintf('   GLaDS (hydrologyglads) solution parameters:'));
+			disp(sprintf('      SHEET'));
+			fielddisplay(self,'pressure_melt_coefficient','Pressure melt coefficient (c_t) [K Pa^-1]');
+			disp(sprintf('      CHANNELS'));
+			fielddisplay(self,'ischannels','Do we allow for channels? 1: yes, 0: no');
+			disp(sprintf('      OTHER'));
+			fielddisplay(self,'requested_outputs','additional outputs requested');
+		end % }}}
+		function marshall(self,prefix,md,fid) % {{{
+
+			yts=md.constants.yts;
+
+			%Marshall model code first
+			WriteData(fid,prefix,'name','md.hydrology.model','data',5,'format','Integer');
+
+			%Sheet
+			WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','pressure_melt_coefficient','format','Double');
+
+			%Channels
+			WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','ischannels','format','Boolean');
+
+			%Others
+			outputs = self.requested_outputs;
+			pos  = find(ismember(outputs,'default'));
+			if ~isempty(pos),
+				outputs(pos) = [];  %remove 'default' from outputs
+				outputs      = [outputs defaultoutputs(self,md)]; %add defaults
+			end
+			WriteData(fid,prefix,'data',outputs,'name','md.hydrology.requested_outputs','format','StringArray');
+		end % }}}
+	end
+end
+
Index: /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.js
===================================================================
--- /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.js	(revision 23934)
+++ /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.js	(revision 23935)
@@ -79,5 +79,5 @@
 		
 	}else if(solutiontype ==='TransientSolution'){
-		analyses=['StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis','MasstransportAnalysis','HydrologyShaktiAnalysis'];
+		analyses=['StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis','MasstransportAnalysis','HydrologyShaktiAnalysis','HydrologyGladsAnalysis'];
 		
 	}else if(solutiontype ==='SealevelriseSolution'){
Index: /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.m
===================================================================
--- /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.m	(revision 23934)
+++ /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.m	(revision 23935)
@@ -71,5 +71,5 @@
 		analyses={'EsaAnalysis'};
 	elseif strcmp(solutiontype,'TransientSolution')
-		analyses={'StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis','MasstransportAnalysis','HydrologyShaktiAnalysis'};
+		analyses={'StressbalanceAnalysis','StressbalanceVerticalAnalysis','StressbalanceSIAAnalysis','L2ProjectionBaseAnalysis','ThermalAnalysis','MeltingAnalysis','EnthalpyAnalysis','MasstransportAnalysis','HydrologyShaktiAnalysis','HydrologyGladsAnalysis'};
 	elseif strcmp(solutiontype,'SealevelriseSolution')
 		analyses={'SealevelriseAnalysis'};
Index: /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.py
===================================================================
--- /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.py	(revision 23934)
+++ /issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.py	(revision 23935)
@@ -28,5 +28,5 @@
         analyses = ['LoveAnalysis']
     elif solutiontype == 'TransientSolution':
-        analyses = ['StressbalanceAnalysis', 'StressbalanceVerticalAnalysis', 'StressbalanceSIAAnalysis', 'L2ProjectionBaseAnalysis', 'ThermalAnalysis', 'MeltingAnalysis', 'EnthalpyAnalysis', 'MasstransportAnalysis']
+        analyses = ['StressbalanceAnalysis', 'StressbalanceVerticalAnalysis', 'StressbalanceSIAAnalysis', 'L2ProjectionBaseAnalysis', 'ThermalAnalysis', 'MeltingAnalysis', 'EnthalpyAnalysis', 'MasstransportAnalysis','HydrologyShaktiAnalysis','HydrologyGladsAnalysis']
     elif solutiontype == 'HydrologySolution':
         analyses = ['L2ProjectionBaseAnalysis', 'HydrologyShreveAnalysis', 'HydrologyDCInefficientAnalysis', 'HydrologyDCEfficientAnalysis']