Index: /issm/trunk-jpl/src/m/classes/hydrologyshakti.m =================================================================== --- /issm/trunk-jpl/src/m/classes/hydrologyshakti.m (revision 23025) +++ /issm/trunk-jpl/src/m/classes/hydrologyshakti.m (revision 23025) @@ -0,0 +1,105 @@ +%HYDROLOGYSHAKTI class definition +% +% Usage: +% hydrologyshakti=hydrologyshakti(); + +classdef hydrologyshakti + properties (SetAccess=public) + head = NaN; + gap_height = NaN; + bump_spacing = NaN; + bump_height = NaN; + englacial_input = NaN; + moulin_input = NaN; + reynolds = NaN; + spchead = NaN; + neumannflux = NaN; + relaxation = 0; + storage = 0; + requested_outputs = {}; + end + methods + function self = extrude(self,md) % {{{ + end % }}} + function self = hydrologyshakti(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 = {'HydrologyHead','HydrologyGapHeight','EffectivePressure','HydrologyBasalFlux','DegreeOfChannelization'}; + end % }}} + + function self = setdefaultparameters(self) % {{{ + % Set under-relaxation parameter to be 1 (no under-relaxation of nonlinear iteration) + self.relaxation=1; + self.storage=0; + self.requested_outputs={'default'}; + end % }}} + function md = checkconsistency(self,md,solution,analyses) % {{{ + + %Early return + if ~ismember('HydrologyShaktiAnalysis',analyses) + return; + end + + md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'Inf',1,'timeseries',1); + md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'Inf',1,'timeseries',1); + md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1,'Inf',1); + md = checkfield(md,'fieldname','hydrology.spchead','size',[md.mesh.numberofvertices 1]); + md = checkfield(md,'fieldname','hydrology.relaxation','>=',0); + md = checkfield(md,'fieldname','hydrology.storage','>=',0); + md = checkfield(md,'fieldname','hydrology.requested_outputs','stringrow',1); + end % }}} + function disp(self) % {{{ + disp(sprintf(' hydrologyshakti solution parameters:')); + fielddisplay(self,'head','subglacial hydrology water head (m)'); + fielddisplay(self,'gap_height','height of gap separating ice to bed (m)'); + fielddisplay(self,'bump_spacing','characteristic bedrock bump spacing (m)'); + fielddisplay(self,'bump_height','characteristic bedrock bump height (m)'); + fielddisplay(self,'englacial_input','liquid water input from englacial to subglacial system (m/yr)'); + fielddisplay(self,'moulin_input','liquid water input from moulins (at the vertices) to subglacial system (m^3/s)'); + fielddisplay(self,'reynolds','Reynolds'' number'); + fielddisplay(self,'neumannflux','water flux applied along the model boundary (m^2/s)'); + fielddisplay(self,'spchead','water head constraints (NaN means no constraint) (m)'); + fielddisplay(self,'relaxation','under-relaxation coefficient for nonlinear iteration'); + fielddisplay(self,'storage','englacial storage coefficient (void ratio)'); + fielddisplay(self,'requested_outputs','additional outputs requested'); + end % }}} + function marshall(self,prefix,md,fid) % {{{ + + yts=md.constants.yts; + + WriteData(fid,prefix,'name','md.hydrology.model','data',3,'format','Integer'); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','head','format','DoubleMat','mattype',1); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','gap_height','format','DoubleMat','mattype',2); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_spacing','format','DoubleMat','mattype',2); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_height','format','DoubleMat','mattype',2); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','englacial_input','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','reynolds','format','DoubleMat','mattype',2); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','neumannflux','format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','spchead','format','DoubleMat','mattype',1); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','relaxation','format','Double'); + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','storage','format','Double'); + 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/classes/hydrologyshakti.py =================================================================== --- /issm/trunk-jpl/src/m/classes/hydrologyshakti.py (revision 23025) +++ /issm/trunk-jpl/src/m/classes/hydrologyshakti.py (revision 23025) @@ -0,0 +1,106 @@ +from fielddisplay import fielddisplay +from checkfield import checkfield +from WriteData import WriteData + +class hydrologyshakti(object): + """ + HYDROLOGYSHAKTI class definition + + Usage: + hydrologyshakti=hydrologyshakti() + """ + + def __init__(self): # {{{ + self.head = float('NaN') + self.gap_height = float('NaN') + self.bump_spacing = float('NaN') + self.bump_height = float('NaN') + self.englacial_input = float('NaN') + self.moulin_input = float('NaN') + self.reynolds = float('NaN') + self.spchead = float('NaN') + self.neumannflux = float('NaN') + self.relaxation = 0 + self.storage = 0 + self.requested_outputs = [] + + #set defaults + self.setdefaultparameters() + + #}}} + def __repr__(self): # {{{ + + string=' hydrologyshakti solution parameters:' + string="%s\n%s"%(string,fielddisplay(self,'head','subglacial hydrology water head (m)')) + string="%s\n%s"%(string,fielddisplay(self,'gap_height','height of gap separating ice to bed (m)')) + string="%s\n%s"%(string,fielddisplay(self,'bump_spacing','characteristic bedrock bump spacing (m)')) + string="%s\n%s"%(string,fielddisplay(self,'bump_height','characteristic bedrock bump height (m)')) + string="%s\n%s"%(string,fielddisplay(self,'englacial_input','liquid water input from englacial to subglacial system (m/yr)')) + string="%s\n%s"%(string,fielddisplay(self,'moulin_input','liquid water input from moulins (at the vertices) to subglacial system (m^3/s)')) + string="%s\n%s"%(string,fielddisplay(self,'reynolds','Reynolds'' number')) + string="%s\n%s"%(string,fielddisplay(self,'neumannflux','water flux applied along the model boundary (m^2/s)')) + string="%s\n%s"%(string,fielddisplay(self,'spchead','water head constraints (NaN means no constraint) (m)')) + string="%s\n%s"%(string,fielddisplay(self,'relaxation','under-relaxation coefficient for nonlinear iteration')) + string="%s\n%s"%(string,fielddisplay(self,'storage','englacial storage coefficient (void ratio)')) + string="%s\n%s"%(string,fielddisplay(self,'requested_outputs','additional outputs requested')) + return string + #}}} + def extrude(self,md): # {{{ + return self + #}}} + def setdefaultparameters(self): # {{{ + # Set under-relaxation parameter to be 1 (no under-relaxation of nonlinear iteration) + self.relaxation=1 + self.storage=0 + self.requested_outputs=['default'] + return self + #}}} + def defaultoutputs(self,md): # {{{ + list = ['HydrologyHead','HydrologyGapHeight','EffectivePressure','HydrologyBasalFlux','DegreeOfChannelization'] + return list + #}}} + + def checkconsistency(self,md,solution,analyses): # {{{ + + #Early return + if 'HydrologyShaktiAnalysis' not in analyses: + return md + + md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'Inf',1,'timeseries',1) + md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'Inf',1,'timeseries',1) + md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1,'Inf',1) + md = checkfield(md,'fieldname','hydrology.spchead','size',[md.mesh.numberofvertices]) + md = checkfield(md,'fieldname','hydrology.relaxation','>=',0) + md = checkfield(md,'fieldname','hydrology.storage','>=',0) + md = checkfield(md,'fieldname','hydrology.requested_outputs','stringrow',1) + + return md + # }}} + def marshall(self,prefix,md,fid): # {{{ + yts=md.constants.yts + + WriteData(fid,prefix,'name','md.hydrology.model','data',3,'format','Integer') + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','head','format','DoubleMat','mattype',1) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','gap_height','format','DoubleMat','mattype',2) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_spacing','format','DoubleMat','mattype',2) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_height','format','DoubleMat','mattype',2) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','englacial_input','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','reynolds','format','DoubleMat','mattype',2) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','neumannflux','format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','spchead','format','DoubleMat','mattype',1) + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','relaxation','format','Double') + WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','storage','format','Double') + #process requested outputs + outputs = self.requested_outputs + indices = [i for i, x in enumerate(outputs) if x == 'default'] + if len(indices) > 0: + outputscopy=outputs[0:max(0,indices[0]-1)]+self.defaultoutputs(md)+outputs[indices[0]+1:] + outputs =outputscopy + WriteData(fid,prefix,'data',outputs,'name','md.hydrology.requested_outputs','format','StringArray') + # }}} Index: sm/trunk-jpl/src/m/classes/hydrologysommers.m =================================================================== --- /issm/trunk-jpl/src/m/classes/hydrologysommers.m (revision 23024) +++ (revision ) @@ -1,105 +1,0 @@ -%HYDROLOGYSOMMERS class definition -% -% Usage: -% hydrologysommers=hydrologysommers(); - -classdef hydrologysommers - properties (SetAccess=public) - head = NaN; - gap_height = NaN; - bump_spacing = NaN; - bump_height = NaN; - englacial_input = NaN; - moulin_input = NaN; - reynolds = NaN; - spchead = NaN; - neumannflux = NaN; - relaxation = 0; - storage = 0; - requested_outputs = {}; - end - methods - function self = extrude(self,md) % {{{ - end % }}} - function self = hydrologysommers(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 = {'HydrologyHead','HydrologyGapHeight','EffectivePressure','HydrologyBasalFlux','DegreeOfChannelization'}; - end % }}} - - function self = setdefaultparameters(self) % {{{ - % Set under-relaxation parameter to be 1 (no under-relaxation of nonlinear iteration) - self.relaxation=1; - self.storage=0; - self.requested_outputs={'default'}; - end % }}} - function md = checkconsistency(self,md,solution,analyses) % {{{ - - %Early return - if ~ismember('HydrologySommersAnalysis',analyses) - return; - end - - md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'Inf',1,'timeseries',1); - md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'Inf',1,'timeseries',1); - md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1,'Inf',1); - md = checkfield(md,'fieldname','hydrology.spchead','size',[md.mesh.numberofvertices 1]); - md = checkfield(md,'fieldname','hydrology.relaxation','>=',0); - md = checkfield(md,'fieldname','hydrology.storage','>=',0); - md = checkfield(md,'fieldname','hydrology.requested_outputs','stringrow',1); - end % }}} - function disp(self) % {{{ - disp(sprintf(' hydrologysommers solution parameters:')); - fielddisplay(self,'head','subglacial hydrology water head (m)'); - fielddisplay(self,'gap_height','height of gap separating ice to bed (m)'); - fielddisplay(self,'bump_spacing','characteristic bedrock bump spacing (m)'); - fielddisplay(self,'bump_height','characteristic bedrock bump height (m)'); - fielddisplay(self,'englacial_input','liquid water input from englacial to subglacial system (m/yr)'); - fielddisplay(self,'moulin_input','liquid water input from moulins (at the vertices) to subglacial system (m^3/s)'); - fielddisplay(self,'reynolds','Reynolds'' number'); - fielddisplay(self,'neumannflux','water flux applied along the model boundary (m^2/s)'); - fielddisplay(self,'spchead','water head constraints (NaN means no constraint) (m)'); - fielddisplay(self,'relaxation','under-relaxation coefficient for nonlinear iteration'); - fielddisplay(self,'storage','englacial storage coefficient (void ratio)'); - fielddisplay(self,'requested_outputs','additional outputs requested'); - end % }}} - function marshall(self,prefix,md,fid) % {{{ - - yts=md.constants.yts; - - WriteData(fid,prefix,'name','md.hydrology.model','data',3,'format','Integer'); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','head','format','DoubleMat','mattype',1); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','gap_height','format','DoubleMat','mattype',2); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_spacing','format','DoubleMat','mattype',2); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_height','format','DoubleMat','mattype',2); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','englacial_input','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','reynolds','format','DoubleMat','mattype',2); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','neumannflux','format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','spchead','format','DoubleMat','mattype',1); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','relaxation','format','Double'); - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','storage','format','Double'); - 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: sm/trunk-jpl/src/m/classes/hydrologysommers.py =================================================================== --- /issm/trunk-jpl/src/m/classes/hydrologysommers.py (revision 23024) +++ (revision ) @@ -1,106 +1,0 @@ -from fielddisplay import fielddisplay -from checkfield import checkfield -from WriteData import WriteData - -class hydrologysommers(object): - """ - HYDROLOGYSOMMERS class definition - - Usage: - hydrologysommers=hydrologysommers() - """ - - def __init__(self): # {{{ - self.head = float('NaN') - self.gap_height = float('NaN') - self.bump_spacing = float('NaN') - self.bump_height = float('NaN') - self.englacial_input = float('NaN') - self.moulin_input = float('NaN') - self.reynolds = float('NaN') - self.spchead = float('NaN') - self.neumannflux = float('NaN') - self.relaxation = 0 - self.storage = 0 - self.requested_outputs = [] - - #set defaults - self.setdefaultparameters() - - #}}} - def __repr__(self): # {{{ - - string=' hydrologysommers solution parameters:' - string="%s\n%s"%(string,fielddisplay(self,'head','subglacial hydrology water head (m)')) - string="%s\n%s"%(string,fielddisplay(self,'gap_height','height of gap separating ice to bed (m)')) - string="%s\n%s"%(string,fielddisplay(self,'bump_spacing','characteristic bedrock bump spacing (m)')) - string="%s\n%s"%(string,fielddisplay(self,'bump_height','characteristic bedrock bump height (m)')) - string="%s\n%s"%(string,fielddisplay(self,'englacial_input','liquid water input from englacial to subglacial system (m/yr)')) - string="%s\n%s"%(string,fielddisplay(self,'moulin_input','liquid water input from moulins (at the vertices) to subglacial system (m^3/s)')) - string="%s\n%s"%(string,fielddisplay(self,'reynolds','Reynolds'' number')) - string="%s\n%s"%(string,fielddisplay(self,'neumannflux','water flux applied along the model boundary (m^2/s)')) - string="%s\n%s"%(string,fielddisplay(self,'spchead','water head constraints (NaN means no constraint) (m)')) - string="%s\n%s"%(string,fielddisplay(self,'relaxation','under-relaxation coefficient for nonlinear iteration')) - string="%s\n%s"%(string,fielddisplay(self,'storage','englacial storage coefficient (void ratio)')) - string="%s\n%s"%(string,fielddisplay(self,'requested_outputs','additional outputs requested')) - return string - #}}} - def extrude(self,md): # {{{ - return self - #}}} - def setdefaultparameters(self): # {{{ - # Set under-relaxation parameter to be 1 (no under-relaxation of nonlinear iteration) - self.relaxation=1 - self.storage=0 - self.requested_outputs=['default'] - return self - #}}} - def defaultoutputs(self,md): # {{{ - list = ['HydrologyHead','HydrologyGapHeight','EffectivePressure','HydrologyBasalFlux','DegreeOfChannelization'] - return list - #}}} - - def checkconsistency(self,md,solution,analyses): # {{{ - - #Early return - if 'HydrologySommersAnalysis' not in analyses: - return md - - md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'Inf',1,'timeseries',1) - md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'Inf',1,'timeseries',1) - md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements],'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1,'Inf',1) - md = checkfield(md,'fieldname','hydrology.spchead','size',[md.mesh.numberofvertices]) - md = checkfield(md,'fieldname','hydrology.relaxation','>=',0) - md = checkfield(md,'fieldname','hydrology.storage','>=',0) - md = checkfield(md,'fieldname','hydrology.requested_outputs','stringrow',1) - - return md - # }}} - def marshall(self,prefix,md,fid): # {{{ - yts=md.constants.yts - - WriteData(fid,prefix,'name','md.hydrology.model','data',3,'format','Integer') - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','head','format','DoubleMat','mattype',1) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','gap_height','format','DoubleMat','mattype',2) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_spacing','format','DoubleMat','mattype',2) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','bump_height','format','DoubleMat','mattype',2) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','englacial_input','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','reynolds','format','DoubleMat','mattype',2) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','neumannflux','format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','spchead','format','DoubleMat','mattype',1) - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','relaxation','format','Double') - WriteData(fid,prefix,'object',self,'class','hydrology','fieldname','storage','format','Double') - #process requested outputs - outputs = self.requested_outputs - indices = [i for i, x in enumerate(outputs) if x == 'default'] - if len(indices) > 0: - outputscopy=outputs[0:max(0,indices[0]-1)]+self.defaultoutputs(md)+outputs[indices[0]+1:] - outputs =outputscopy - WriteData(fid,prefix,'data',outputs,'name','md.hydrology.requested_outputs','format','StringArray') - # }}}