%FRICTIONWEERTMAN class definition
%
%   Usage:
%      friction=frictionhydro();

classdef frictionhydro
	properties (SetAccess=public) 
		q                  = NaN;
		C                  = NaN;
		As                 = NaN;
		effective_pressure = NaN; 
	end
	methods
		function obj = frictionhydro(varargin) % {{{
			switch nargin
				case 0
					obj=setdefaultparameters(obj);
				otherwise
					error('constructor not supported');
			end
		end % }}}
		function obj = setdefaultparameters(obj) % {{{

		end % }}}
		function md = checkconsistency(obj,md,solution,analyses) % {{{

			%Early return
			if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
			md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
			md = checkfield(md,'fieldname','friction.C','NaN',1,'size',[md.mesh.numberofelements 1]);
			md = checkfield(md,'fieldname','friction.As','NaN',1,'size',[md.mesh.numberofelements 1]);
		end % }}}
		function disp(obj) % {{{
			disp(sprintf('Effective Pressure based friction law described in Gagliardini 2007'));
			fielddisplay(obj,'q','friction law exponent q>=1');
			fielddisplay(obj,'C','friction law max value [SI]');
			fielddisplay(obj,'As','Sliding Parameter without cavitation [m Pa^-n s^-1]');
			fielddisplay(obj,'effective_pressure','Effective Pressure [Pa]');
		end % }}}
		function marshall(obj,md,fid) % {{{
			WriteData(fid,'enum',FrictionLawEnum,'data',3,'format','Integer');
			WriteData(fid,'class','friction','object',obj,'fieldname','q','format','DoubleMat','mattype',2);
			WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',2);
			WriteData(fid,'class','friction','object',obj,'fieldname','As','format','DoubleMat','mattype',2);
			WriteData(fid,'class','friction','object',obj,'fieldname','effective_pressure','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1);
		end % }}}
	end
end