Changeset 22859


Ignore:
Timestamp:
06/21/18 06:02:26 (7 years ago)
Author:
bdef
Message:

CHG: modifying naming related to hydrology sub-step

Location:
issm/trunk-jpl/src/m/classes
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/hydrologydc.m

    r22286 r22859  
    55
    66classdef hydrologydc
    7         properties (SetAccess=public) 
     7        properties (SetAccess=public)
    88                water_compressibility    = 0;
    99                isefficientlayer         = 0;
    … …  
    4040  end
    4141        methods
    42     % {{{ function self = extrude(self,md) 
    43                 function self = extrude(self,md) 
     42    % {{{ function self = extrude(self,md)
     43                function self = extrude(self,md)
    4444                        self.spcsediment_head=project3d(md,'vector',self.spcsediment_head,'type','node','layer',1);
    4545                        self.sediment_transmitivity=project3d(md,'vector',self.sediment_transmitivity,'type','node','layer',1);
    … …  
    4949        self.mask_eplactive_node=project3d(md,'vector',self.mask_eplactive_node,'type','node','layer',1);
    5050                        end
    51                 end 
     51                end
    5252    % }}}
    53                 % {{{ function self = hydrologydc(varargin) 
    54                 function self = hydrologydc(varargin) 
     53                % {{{ function self = hydrologydc(varargin)
     54                function self = hydrologydc(varargin)
    5555                        switch nargin
    5656                                case 0
    … …  
    5959                                        error('constructor not supported');
    6060                        end
    61                 end 
     61                end
    6262                % }}}
    6363    function list = defaultoutputs(self,md) % {{{
    64         list = {'SedimentHead','SedimentHeadResidual','EffectivePressure'};
     64        list = {'SedimentHeadHydrostep','SedimentHeadResidual','EffectivePressureHydrostep'};
    6565        if self.isefficientlayer,
    66             list=[list,{'EplHead','HydrologydcMaskEplactiveNode','HydrologydcMaskEplactiveElt','EplHeadSlopeX','EplHeadSlopeY','HydrologydcEplThickness'}];
     66            list=[list,{'EplHeadHydrostep','HydrologydcMaskEplactiveNode','HydrologydcMaskEplactiveElt','EplHeadSlopeX','EplHeadSlopeY','HydrologydcEplThicknessHydrostep'}];
    6767        end
    6868        if self.steps_per_step>1,
    69             list = [list,'EffectivePressureTimeAverage','SedimentHeadTimeAverage'];
     69            list = [list,'EffectivePressure','SedimentHead'];
    7070            if self.isefficientlayer,
    71                 list = [list,'EplHeadTimeAverage','HydrologydcEplThicknessTimeAverage'];
     71                list = [list,'EplHead','HydrologydcEplThickness'];
    7272            end
    7373        end
    74     end % }}}   
     74    end % }}}
    7575
    7676                function self = initialize(self,md) % {{{
    … …  
    8181
    8282                end % }}}
    83                 % {{{ function self = setdefaultparameters(self) 
    84 
    85                 function self = setdefaultparameters(self) 
     83                % {{{ function self = setdefaultparameters(self)
     84
     85                function self = setdefaultparameters(self)
    8686
    8787                        %Parameters from de Fleurian 2014
    … …  
    114114                        self.epl_conductivity         = 8.0e-02;
    115115                        self.eplflip_lock             = 0;
    116                 end 
    117 
    118                 % }}}
    119                 % {{{ function md = checkconsistency(self,md,solution,analyses) 
    120 
    121                 function md = checkconsistency(self,md,solution,analyses) 
     116                end
     117
     118                % }}}
     119                % {{{ function md = checkconsistency(self,md,solution,analyses)
     120
     121                function md = checkconsistency(self,md,solution,analyses)
    122122                        %Early return
    123123                        if ~ismember('HydrologyDCInefficientAnalysis',analyses) & ~ismember('HydrologyDCEfficientAnalysis',analyses),
    … …  
    163163                                if (self.epl_colapse_thickness>self.epl_initial_thickness),
    164164                                        md = checkmessage(md,'Colapsing thickness for EPL larger than initial thickness');
    165                                 end 
    166                         end
    167                 end 
     165                                end
     166                        end
     167                end
    168168
    169169                % }}}
    170170                % {{{ function disp(self)
    171171
    172                 function disp(self) 
     172                function disp(self)
    173173                        disp(sprintf('   hydrology Dual Porous Continuum Equivalent parameters:'));
    174174                        disp(sprintf('   - general parameters'));
    … …  
    220220                        end
    221221
    222                 end 
    223 
    224                 % }}}
    225                 % {{{ function marshall(self,prefix,md,fid) 
    226 
    227                 function marshall(self,prefix,md,fid) 
     222                end
     223
     224                % }}}
     225                % {{{ function marshall(self,prefix,md,fid)
     226
     227                function marshall(self,prefix,md,fid)
    228228                        WriteData(fid,prefix,'name','md.hydrology.model','data',1,'format','Integer');
    229229                        WriteData(fid,prefix,'object',self,'fieldname','water_compressibility','format','Double');
    … …  
    246246
    247247                        WriteData(fid,prefix,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
    248                         WriteData(fid,prefix,'object',self,'fieldname','sediment_compressibility','format','Double');                   
    249                         WriteData(fid,prefix,'object',self,'fieldname','sediment_porosity','format','Double');                 
     248                        WriteData(fid,prefix,'object',self,'fieldname','sediment_compressibility','format','Double');
     249                        WriteData(fid,prefix,'object',self,'fieldname','sediment_porosity','format','Double');
    250250                        WriteData(fid,prefix,'object',self,'fieldname','sediment_thickness','format','Double');
    251                         WriteData(fid,prefix,'object',self,'fieldname','sediment_transmitivity','format','DoubleMat','mattype',1');             
    252 
    253                         if self.isefficientlayer==1,   
    254                                 WriteData(fid,prefix,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);   
     251                        WriteData(fid,prefix,'object',self,'fieldname','sediment_transmitivity','format','DoubleMat','mattype',1');
     252
     253                        if self.isefficientlayer==1,
     254                                WriteData(fid,prefix,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
    255255                                WriteData(fid,prefix,'object',self,'fieldname','mask_eplactive_node','format','DoubleMat','mattype',1);
    256                                 WriteData(fid,prefix,'object',self,'fieldname','epl_compressibility','format','Double');                       
    257                                 WriteData(fid,prefix,'object',self,'fieldname','epl_porosity','format','Double');                       
     256                                WriteData(fid,prefix,'object',self,'fieldname','epl_compressibility','format','Double');
     257                                WriteData(fid,prefix,'object',self,'fieldname','epl_porosity','format','Double');
    258258                                WriteData(fid,prefix,'object',self,'fieldname','epl_initial_thickness','format','Double');
    259259                                WriteData(fid,prefix,'object',self,'fieldname','epl_colapse_thickness','format','Double');
    … …  
    270270      end
    271271      WriteData(fid,prefix,'data',outputs,'name','md.hydrology.requested_outputs','format','StringArray');
    272                 end 
     272                end
    273273
    274274                % }}}

  • issm/trunk-jpl/src/m/classes/hydrologydc.py

    r22286 r22859  
    4545                self.epl_conductivity         = 0
    4646                self.eplflip_lock             = 0
    47                                  
     47
    4848                #set defaults
    4949                self.setdefaultparameters()
    … …  
    5858                string="%s\n%s"%(string,fielddisplay(self,'rel_tol','tolerance of the nonlinear iteration for the transfer between layers [dimensionless]'))
    5959                string="%s\n%s"%(string,fielddisplay(self,'max_iter','maximum number of nonlinear iteration'))
    60                 string="%s\n%s"%(string,fielddisplay(self,'steps_per_step','number of hydrology steps per time step')) 
     60                string="%s\n%s"%(string,fielddisplay(self,'steps_per_step','number of hydrology steps per time step'))
    6161                string="%s\n%s"%(string,fielddisplay(self,'basal_moulin_input','water flux at a given point [m3 s-1]'))
    6262                string="%s\n%s"%(string,fielddisplay(self,'requested_outputs','additional outputs requested'))
    … …  
    6666                string="%s\n\t\t%s"%(string,'2: hydrostatic pressure')
    6767                string="%s\n\t\t%s"%(string,'3: normal stress')
    68        
     68
    6969                if self.sedimentlimit_flag==1:
    7070                        string="%s\n%s"%(string,fielddisplay(self,'sedimentlimit','user defined upper limit for the inefficient layer [m]'))
    … …  
    7373                string="%s\n\t\t%s"%(string,'0: no transfer')
    7474                string="%s\n\t\t%s"%(string,'1: constant leakage factor: leakage_factor')
    75                          
     75
    7676                if self.transfer_flag is 1:
    7777                        string="%s\n%s"%(string,fielddisplay(self,'leakage_factor','user defined leakage factor [m]'))
    … …  
    8080                string="%s\n\t\t%s"%(string,'0: Confined only')
    8181                string="%s\n\t\t%s"%(string,'1: Confined-Unconfined')
    82                          
     82
    8383                string="%s\n%s"%(string,'   - for the sediment layer')
    8484                string="%s\n%s"%(string,fielddisplay(self,'spcsediment_head','sediment water head constraints (NaN means no constraint) [m above MSL]'))
    … …  
    9595                        string="%s\n%s"%(string,fielddisplay(self,'epl_porosity','epl [dimensionless]'))
    9696                        string="%s\n%s"%(string,fielddisplay(self,'epl_max_thickness','epl maximal thickness [m]'))
    97                         string="%s\n%s"%(string,fielddisplay(self,'epl_initial_thickness','epl initial thickness [m]'))                 
     97                        string="%s\n%s"%(string,fielddisplay(self,'epl_initial_thickness','epl initial thickness [m]'))
    9898                        string="%s\n%s"%(string,fielddisplay(self,'epl_colapse_thickness','epl colapsing thickness [m]'))
    9999                        string="%s\n%s"%(string,fielddisplay(self,'epl_thick_comp','epl thickness computation flag'))
    … …  
    111111                return self
    112112        #}}}
    113         def setdefaultparameters(self): #{{{ 
     113        def setdefaultparameters(self): #{{{
    114114                #Parameters from de Fleurian 2014
    115115                self.water_compressibility    = 5.04e-10
    … …  
    140140                self.epl_conductivity         = 8.0e-02
    141141                self.eplflip_lock             = 0
    142                
     142
    143143                return self
    144144        # }}}
    145145
    146146        def defaultoutputs(self,md): # {{{
    147                 list = ['SedimentHead','SedimentHeadResidual','EffectivePressure']
    148                 if self.isefficientlayer==1:
    149                         list.extend(['EplHead','HydrologydcMaskEplactiveNode','HydrologydcMaskEplactiveElt','EplHeadSlopeX','EplHeadSlopeY','HydrologydcEplThickness'])
     147                list = ['SedimentHeadHydrostep','SedimentHeadResidual','EffectivePressureHydrostep']
     148                if self.isefficientlayer==1:
     149                        list.extend(['EplHeadHydrostep','HydrologydcMaskEplactiveNode','HydrologydcMaskEplactiveElt','EplHeadSlopeX','EplHeadSlopeY','HydrologydcEplThicknessHydrostep'])
    150150                if self.steps_per_step>1:
    151                         list.extend(['EffectivePressureTimeAverage','SedimentHeadTimeAverage'])
     151                        list.extend(['EffectivePressure','SedimentHead'])
    152152                        if self.isefficientlayer==1:
    153                                 list.extend(['EplHeadTimeAverage','HydrologydcEplThicknessTimeAverage'])
     153                                list.extend(['EplHead','HydrologydcEplThickness'])
    154154                return list
    155155        #}}}
    … …  
    162162                return self
    163163        # }}}
    164         def checkconsistency(self,md,solution,analyses): #{{{ 
     164        def checkconsistency(self,md,solution,analyses): #{{{
    165165
    166166                #Early return
    … …  
    206206                        md = checkfield(md,'fieldname','hydrology.epl_conductivity','numel',[1],'>',0.)
    207207        # }}}
    208         def marshall(self,prefix,md,fid): #{{{ 
     208        def marshall(self,prefix,md,fid): #{{{
    209209                WriteData(fid,prefix,'name','md.hydrology.model','data',1,'format','Integer')
    210210                WriteData(fid,prefix,'object',self,'fieldname','water_compressibility','format','Double')
    … …  
    227227                WriteData(fid,prefix,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts)
    228228                WriteData(fid,prefix,'object',self,'fieldname','sediment_compressibility','format','Double')
    229                 WriteData(fid,prefix,'object',self,'fieldname','sediment_porosity','format','Double')                   
     229                WriteData(fid,prefix,'object',self,'fieldname','sediment_porosity','format','Double')
    230230                WriteData(fid,prefix,'object',self,'fieldname','sediment_thickness','format','Double')
    231                 WriteData(fid,prefix,'object',self,'fieldname','sediment_transmitivity','format','DoubleMat','mattype',1)               
    232 
    233                 if self.isefficientlayer==1:   
    234                         WriteData(fid,prefix,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts)     
     231                WriteData(fid,prefix,'object',self,'fieldname','sediment_transmitivity','format','DoubleMat','mattype',1)
     232
     233                if self.isefficientlayer==1:
     234                        WriteData(fid,prefix,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts)
    235235                        WriteData(fid,prefix,'object',self,'fieldname','mask_eplactive_node','format','DoubleMat','mattype',1)
    236                         WriteData(fid,prefix,'object',self,'fieldname','epl_compressibility','format','Double')                 
    237                         WriteData(fid,prefix,'object',self,'fieldname','epl_porosity','format','Double')                       
     236                        WriteData(fid,prefix,'object',self,'fieldname','epl_compressibility','format','Double')
     237                        WriteData(fid,prefix,'object',self,'fieldname','epl_porosity','format','Double')
    238238                        WriteData(fid,prefix,'object',self,'fieldname','epl_max_thickness','format','Double')
    239                         WriteData(fid,prefix,'object',self,'fieldname','epl_initial_thickness','format','Double')                       
     239                        WriteData(fid,prefix,'object',self,'fieldname','epl_initial_thickness','format','Double')
    240240                        WriteData(fid,prefix,'object',self,'fieldname','epl_colapse_thickness','format','Double')
    241                         WriteData(fid,prefix,'object',self,'fieldname','epl_thick_comp','format','Integer')                     
     241                        WriteData(fid,prefix,'object',self,'fieldname','epl_thick_comp','format','Integer')
    242242                        WriteData(fid,prefix,'object',self,'fieldname','epl_conductivity','format','Double')
    243243                        WriteData(fid,prefix,'object',self,'fieldname','eplflip_lock','format','Integer')
Note: See TracChangeset for help on using the changeset viewer.