Changeset 22875

Timestamp:

06/25/18 07:32:44 (7 years ago)

Author:

bdef

Message:

initializing effective pressure on frictioneffectivepressure for coupling

File:

• issm/trunk-jpl/src/m/classes/friction.py (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/friction.py

@@ -38 +38 @@
                 self.q=project3d(md,'vector',self.q,'type','element')
                 #if self.coupling==0: #doesnt work with empty loop, so just skip it?
-                if self.coupling==3:
+                if self.coupling in[3,4]:
                         self.effective_pressure=project3d(md,'vector',self.effective_pressure,'type','node','layer',1)
-                elif self.coupling==4:
-                        raise ValueError('coupling not supported yet')
                 elif self.coupling > 4:
                         raise ValueError('md.friction.coupling larger than 4, not supported yet')
@@ -71 +69 @@
                 WriteData(fid,prefix,'object',self,'fieldname','q','format','DoubleMat','mattype',2)
                 WriteData(fid,prefix,'class','friction','object',self,'fieldname','coupling','format','Integer')
-                if self.coupling==3:
+                if self.coupling in[3,4]:
                         WriteData(fid,prefix,'class','friction','object',self,'fieldname','effective_pressure','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts)
                 elif self.coupling > 4: