original_docstring: 
          This method plots mean WSS (POISEUILLE) for a single mesh 
          
docstring: This method plots mean WSS (POISEUILLE) for a single mesh
docstring_tokens: This method plots mean WSS ( POISEUILLE ) for a single mesh
code: def PlotPWSS(self, meshid):
          '''
          This method plots mean WSS (POISEUILLE) for a single mesh 
          '''
          try:
              from matplotlib.pyplot import plot, xlabel, ylabel, title, legend, savefig, close
          except ImportError:
              MatPlotLibError()
          meshid = str(meshid)
          for element in self.NetworkMesh.Elements:
              if element.Id == meshid:
                  dofs = element.GetPoiseuilleDofs()
                  p1 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):]
                  p2 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):]
                  if mean(p1)>=mean(p2):                        
                      Flow = (p1-p2)/element.R
                  else:
                      Flow = (p2-p1)/element.R
                  
                  Wss = ((4.0*element.eta)/pi) * (Flow/(mean(element.Radius))**3)
                  print element.Name, "Wss(mean) = ", mean(Wss), "Pa\n", "Wss(peak) = ", max(Wss), "Pa"
                  
          plot(self.t, Wss, 'r-',linewidth = 3, label = 'wss')   #red line
          xlabel('Time ($s$)')
          ylabel('Wss ($Pa$)')
          title ('Wall Shear Stress')    
          legend()
          savefig(self.w_images + meshid +'_Pwss.png')
          close()
code_tokens: def PlotPWSS ( self , meshid ) : try : from matplotlib . pyplot import plot , xlabel , ylabel , title , legend , savefig , close except ImportError : MatPlotLibError ( ) meshid = str ( meshid ) for element in self . NetworkMesh . Elements : if element . Id == meshid : dofs = element . GetPoiseuilleDofs ( ) p1 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] p2 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] if mean ( p1 ) >= mean ( p2 ) : Flow = ( p1 - p2 ) / element . R else : Flow = ( p2 - p1 ) / element . R Wss = ( ( 4.0 * element . eta ) / pi ) * ( Flow / ( mean ( element . Radius ) ) ** 3 ) print element . Name , "Wss(mean) = " , mean ( Wss ) , "Pa\n" , "Wss(peak) = " , max ( Wss ) , "Pa" plot ( self . t , Wss , 'r-' , linewidth = 3 , label = 'wss' ) xlabel ( 'Time ($s$)' ) ylabel ( 'Wss ($Pa$)' ) title ( 'Wall Shear Stress' ) legend ( ) savefig ( self . w_images + meshid + '_Pwss.png' ) close ( )
short_docstring: This method plots mean WSS (POISEUILLE) for a single mesh
short_docstring_tokens: This method plots mean WSS ( POISEUILLE ) for a single mesh
comment: '''
          This method plots mean WSS (POISEUILLE) for a single mesh 
          '''$SEP$#red line