Changeset 28030 for issm/trunk-jpl/src/m/classes/frictioncoulomb.py
- Timestamp:
- 12/19/23 21:01:48 (15 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/frictioncoulomb.py
r27458 r28030 7 7 8 8 9 class frictioncoulomb( object):10 """ FRICTIONCOULOMBclass definition9 class frictioncoulomb(): 10 """frictioncoulomb class definition 11 11 12 12 Usage: … … 14 14 """ 15 15 16 def __init__(self ): # {{{16 def __init__(self, *args): # {{{ 17 17 self.coefficient = np.nan 18 18 self.coefficientcoulomb = np.nan … … 23 23 self.effective_pressure_limit = 0 24 24 25 self.setdefaultparameters() 25 nargs = len(args) 26 27 if nargs == 0: 28 self.setdefaultparameters() 29 else: 30 raise Exception('constructor not supported') 26 31 # }}} 27 32 … … 77 82 78 83 def marshall(self, prefix, md, fid): # {{{ 84 yts = md.constants.yts 85 79 86 WriteData(fid, prefix, 'name', 'md.friction.law', 'data', 7, 'format', 'Integer') 80 WriteData(fid, prefix, 'object', self, 'fieldname', 'coefficient', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', md.constants.yts)87 WriteData(fid, prefix, 'object', self, 'fieldname', 'coefficient', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', yts) 81 88 WriteData(fid, prefix, 'object', self, 'fieldname', 'coefficientcoulomb', 'format', 'DoubleMat', 'mattype', 1) 82 89 WriteData(fid, prefix, 'object', self, 'fieldname', 'p', 'format', 'DoubleMat', 'mattype', 2) … … 85 92 WriteData(fid, prefix, 'object', self, 'class', 'friction', 'fieldname', 'effective_pressure_limit', 'format', 'Double') 86 93 if self.coupling == 1: 87 WriteData(fid, prefix, 'class', 'friction', 'object', self, 'fieldname', 'effective_pressure', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', md.constants.yts)94 WriteData(fid, prefix, 'class', 'friction', 'object', self, 'fieldname', 'effective_pressure', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', yts) 88 95 elif self.coupling == 2: 89 96 raise ValueError('not implemented yet')
Note:
See TracChangeset
for help on using the changeset viewer.