Paraview python script for creating data using only CLI(terminal interface only) and saving in csv file

ParaView Support
,
1

Hello All,
I have less knowledge and experience about paraview python scripting.
I have paraview python script, partly generated by paraview and partly created by me (thanks to the ‘macro record’ option creator, I hope the person live long enough).
Script functionality:
1.The script essentially takes the file readable by paraview from openfoam avilable in that directory
2.processes some graphical option changes and adjustments
3.then moves towards calculation of ‘isovolumedisplay’ and from that calculates ‘alpha.water’ using ‘integrateVariable’
4.Shows calculated value in shreadsheet format.
From that I copy the value manually and paste it in Excel cell.
Program:

#!usr/bin/python
#script to automate the process of calculating volume fraction (also called as void fraction) alpha 
#and show graphically 

#### import the simple module from the paraview
from paraview.simple import *
#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()

# create a new 'OpenFOAMReader'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM = OpenFOAMReader(registrationName='3D_Case20_halfGeom_(70)_FlR_0_0009631_Water.OpenFOAM', FileName='E:\\CCS\\3D_cases\\H2O_Cases\\halfGeom\\3D_Case20_halfGeom_(70)_FlR_0_0009631_Water\\3D_Case20_halfGeom_(70)_FlR_0_0009631_Water.OpenFOAM')
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.MeshRegions = ['internalMesh']
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.CellArrays = ['U', 'UMean', 'UPrime2Mean', 'alpha.water', 'alpha.waterMean', 'alpha.waterPrime2Mean', 'p', 'pMean', 'pPrime2Mean', 'p_rgh']

# get animation scene
animationScene1 = GetAnimationScene()

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps()

# Properties modified on a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.LabelSize = '64-bit'

# get active view
renderView1 = GetActiveViewOrCreate('RenderView')

# show data in view
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay = Show(a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM, renderView1, 'UnstructuredGridRepresentation')

# get color transfer function/color map for 'p'
pLUT = GetColorTransferFunction('p')

# get opacity transfer function/opacity map for 'p'
pPWF = GetOpacityTransferFunction('p')

# trace defaults for the display properties.
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.Representation = 'Surface'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ColorArrayName = ['POINTS', 'p']
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.LookupTable = pLUT
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectTCoordArray = 'None'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectNormalArray = 'None'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectTangentArray = 'None'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OSPRayScaleArray = 'p'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OSPRayScaleFunction = 'PiecewiseFunction'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectOrientationVectors = 'U'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleFactor = 0.01
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectScaleArray = 'p'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GlyphType = 'Arrow'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GlyphTableIndexArray = 'p'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GaussianRadius = 0.0005
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SetScaleArray = ['POINTS', 'p']
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleTransferFunction = 'PiecewiseFunction'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityArray = ['POINTS', 'p']
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityTransferFunction = 'PiecewiseFunction'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.DataAxesGrid = 'GridAxesRepresentation'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.PolarAxes = 'PolarAxesRepresentation'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScalarOpacityFunction = pPWF
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScalarOpacityUnitDistance = 0.0018
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityArrayName = ['POINTS', 'p']
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ExtractedBlockIndex = 1

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 606.7, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 606.7, 1.0, 0.5, 0.0]

# reset view to fit data
renderView1.ResetCamera()

# get the material library
materialLibrary1 = GetMaterialLibrary()

# show color bar/color legend
a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SetScalarBarVisibility(renderView1, True)

# update the view to ensure updated data information
renderView1.Update()

# create a new 'Clip'
clip1 = Clip(registrationName='Clip1', Input=a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM)
clip1.ClipType = 'Plane'
clip1.HyperTreeGridClipper = 'Plane'
clip1.Scalars = ['POINTS', 'p']
clip1.Value = 491.24

# init the 'Plane' selected for 'ClipType'
clip1.ClipType.Origin = [0, 0, 0.0065]

# init the 'Plane' selected for 'HyperTreeGridClipper'
clip1.HyperTreeGridClipper.Origin = [0, 0, 0.0065]

# Properties modified on clip1.ClipType
clip1.ClipType.Origin = [0.0, -0.045, 0.0]
clip1.ClipType.Normal = [0.0, 1.0, 0.0]

# show data in view
clip1Display = Show(clip1, renderView1, 'UnstructuredGridRepresentation')

# trace defaults for the display properties.
clip1Display.Representation = 'Surface'
clip1Display.ColorArrayName = ['POINTS', 'p']
clip1Display.LookupTable = pLUT
clip1Display.SelectTCoordArray = 'None'
clip1Display.SelectNormalArray = 'None'
clip1Display.SelectTangentArray = 'None'
clip1Display.OSPRayScaleArray = 'p'
clip1Display.OSPRayScaleFunction = 'PiecewiseFunction'
clip1Display.SelectOrientationVectors = 'U'
clip1Display.ScaleFactor = 0.0024
clip1Display.SelectScaleArray = 'p'
clip1Display.GlyphType = 'Arrow'
clip1Display.GlyphTableIndexArray = 'p'
clip1Display.GaussianRadius = 0.00012
clip1Display.SetScaleArray = ['POINTS', 'p']
clip1Display.ScaleTransferFunction = 'PiecewiseFunction'
clip1Display.OpacityArray = ['POINTS', 'p']
clip1Display.OpacityTransferFunction = 'PiecewiseFunction'
clip1Display.DataAxesGrid = 'GridAxesRepresentation'
clip1Display.PolarAxes = 'PolarAxesRepresentation'
clip1Display.ScalarOpacityFunction = pPWF
clip1Display.ScalarOpacityUnitDistance = 0.00113
clip1Display.OpacityArrayName = ['POINTS', 'p']
clip1Display.ExtractedBlockIndex = 1

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
clip1Display.ScaleTransferFunction.Points = [418.45, 0.0, 0.5, 0.0, 592.33, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
clip1Display.OpacityTransferFunction.Points = [418.45, 0.0, 0.5, 0.0, 592.33, 1.0, 0.5, 0.0]

# hide data in view
Hide(a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM, renderView1)

# show color bar/color legend
clip1Display.SetScalarBarVisibility(renderView1, True)

# update the view to ensure updated data information
renderView1.Update()

# Properties modified on clip1
clip1.Invert = 0

# update the view to ensure updated data information
renderView1.Update()

# create a new 'Clip'
clip2 = Clip(registrationName='Clip2', Input=clip1)
clip2.ClipType = 'Plane'
clip2.HyperTreeGridClipper = 'Plane'
clip2.Scalars = ['POINTS', 'p']
clip2.Value = 491.24

# init the 'Plane' selected for 'ClipType'
clip2.ClipType.Origin = [0, 0.0038, 0.0065]

# init the 'Plane' selected for 'HyperTreeGridClipper'
clip2.HyperTreeGridClipper.Origin = [0, 0.0038, 0.0065]

# Properties modified on clip2.ClipType
clip2.ClipType.Origin = [0.0, 0.045, 0.0]
clip2.ClipType.Normal = [0.0, 1.0, 0.0]

# show data in view
clip2Display = Show(clip2, renderView1, 'UnstructuredGridRepresentation')

# trace defaults for the display properties.
clip2Display.Representation = 'Surface'
clip2Display.ColorArrayName = ['POINTS', 'p']
clip2Display.LookupTable = pLUT
clip2Display.SelectTCoordArray = 'None'
clip2Display.SelectNormalArray = 'None'
clip2Display.SelectTangentArray = 'None'
clip2Display.OSPRayScaleArray = 'p'
clip2Display.OSPRayScaleFunction = 'PiecewiseFunction'
clip2Display.SelectOrientationVectors = 'U'
clip2Display.ScaleFactor = 0.009
clip2Display.SelectScaleArray = 'p'
clip2Display.GlyphType = 'Arrow'
clip2Display.GlyphTableIndexArray = 'p'
clip2Display.GaussianRadius = 0.00045
clip2Display.SetScaleArray = ['POINTS', 'p']
clip2Display.ScaleTransferFunction = 'PiecewiseFunction'
clip2Display.OpacityArray = ['POINTS', 'p']
clip2Display.OpacityTransferFunction = 'PiecewiseFunction'
clip2Display.DataAxesGrid = 'GridAxesRepresentation'
clip2Display.PolarAxes = 'PolarAxesRepresentation'
clip2Display.ScalarOpacityFunction = pPWF
clip2Display.ScalarOpacityUnitDistance = 0.0016
clip2Display.OpacityArrayName = ['POINTS', 'p']
clip2Display.ExtractedBlockIndex = 1

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
clip2Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
clip2Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

# hide data in view
Hide(clip1, renderView1)

# show color bar/color legend
clip2Display.SetScalarBarVisibility(renderView1, True)

# update the view to ensure updated data information
renderView1.Update()

# create a new 'Clip'
clip3 = Clip(registrationName='Clip3', Input=clip2)
clip3.ClipType = 'Plane'
clip3.HyperTreeGridClipper = 'Plane'
clip3.Scalars = ['POINTS', 'p']
clip3.Value = 466.04

# init the 'Plane' selected for 'ClipType'
clip3.ClipType.Origin = [0.0, 0.0, 0.0065]

# init the 'Plane' selected for 'HyperTreeGridClipper'
clip3.HyperTreeGridClipper.Origin = [0.0, 0.0, 0.0065]

# Properties modified on clip3
clip3.Invert = 0

# Properties modified on clip3.ClipType
clip3.ClipType.Origin = [0.0, 0.0, 0.0017]
clip3.ClipType.Normal = [0.0, 0.0, 1.0]

# show data in view
clip3Display = Show(clip3, renderView1, 'UnstructuredGridRepresentation')

# trace defaults for the display properties.
clip3Display.Representation = 'Surface'
clip3Display.ColorArrayName = ['POINTS', 'p']
clip3Display.LookupTable = pLUT
clip3Display.SelectTCoordArray = 'None'
clip3Display.SelectNormalArray = 'None'
clip3Display.SelectTangentArray = 'None'
clip3Display.OSPRayScaleArray = 'p'
clip3Display.OSPRayScaleFunction = 'PiecewiseFunction'
clip3Display.SelectOrientationVectors = 'U'
clip3Display.ScaleFactor = 0.009
clip3Display.SelectScaleArray = 'p'
clip3Display.GlyphType = 'Arrow'
clip3Display.GlyphTableIndexArray = 'p'
clip3Display.GaussianRadius = 0.00045
clip3Display.SetScaleArray = ['POINTS', 'p']
clip3Display.ScaleTransferFunction = 'PiecewiseFunction'
clip3Display.OpacityArray = ['POINTS', 'p']
clip3Display.OpacityTransferFunction = 'PiecewiseFunction'
clip3Display.DataAxesGrid = 'GridAxesRepresentation'
clip3Display.PolarAxes = 'PolarAxesRepresentation'
clip3Display.ScalarOpacityFunction = pPWF
clip3Display.ScalarOpacityUnitDistance = 0.0017
clip3Display.OpacityArrayName = ['POINTS', 'p']
clip3Display.ExtractedBlockIndex = 1

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
clip3Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
clip3Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

# hide data in view
Hide(clip2, renderView1)

# show color bar/color legend
clip3Display.SetScalarBarVisibility(renderView1, True)

# update the view to ensure updated data information
renderView1.Update()

# create a new 'Iso Volume'
isoVolume1 = IsoVolume(registrationName='IsoVolume1', Input=clip3)
isoVolume1.InputScalars = ['POINTS', 'p']
isoVolume1.ThresholdRange = [375.79, 556.30]

# Properties modified on isoVolume1
isoVolume1.InputScalars = ['POINTS', 'alpha.water']
isoVolume1.ThresholdRange = [0.91, 1.0]

# show data in view
isoVolume1Display = Show(isoVolume1, renderView1, 'UnstructuredGridRepresentation')

# trace defaults for the display properties.
isoVolume1Display.Representation = 'Surface'
isoVolume1Display.ColorArrayName = ['POINTS', 'p']
isoVolume1Display.LookupTable = pLUT
isoVolume1Display.SelectTCoordArray = 'None'
isoVolume1Display.SelectNormalArray = 'None'
isoVolume1Display.SelectTangentArray = 'None'
isoVolume1Display.OSPRayScaleArray = 'p'
isoVolume1Display.OSPRayScaleFunction = 'PiecewiseFunction'
isoVolume1Display.SelectOrientationVectors = 'U'
isoVolume1Display.ScaleFactor = 0.009
isoVolume1Display.SelectScaleArray = 'p'
isoVolume1Display.GlyphType = 'Arrow'
isoVolume1Display.GlyphTableIndexArray = 'p'
isoVolume1Display.GaussianRadius = 0.00045
isoVolume1Display.SetScaleArray = ['POINTS', 'p']
isoVolume1Display.ScaleTransferFunction = 'PiecewiseFunction'
isoVolume1Display.OpacityArray = ['POINTS', 'p']
isoVolume1Display.OpacityTransferFunction = 'PiecewiseFunction'
isoVolume1Display.DataAxesGrid = 'GridAxesRepresentation'
isoVolume1Display.PolarAxes = 'PolarAxesRepresentation'
isoVolume1Display.ScalarOpacityFunction = pPWF
isoVolume1Display.ScalarOpacityUnitDistance = 0.0048
isoVolume1Display.OpacityArrayName = ['POINTS', 'p']
isoVolume1Display.ExtractedBlockIndex = 1

# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
isoVolume1Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 554.07, 1.0, 0.5, 0.0]

# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
isoVolume1Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 554.07, 1.0, 0.5, 0.0]

# hide data in view
Hide(clip3, renderView1)

# show color bar/color legend
isoVolume1Display.SetScalarBarVisibility(renderView1, True)

# update the view to ensure updated data information
renderView1.Update()

# set scalar coloring
ColorBy(isoVolume1Display, ('POINTS', 'alpha.water'))

# Hide the scalar bar for this color map if no visible data is colored by it.
HideScalarBarIfNotNeeded(pLUT, renderView1)

# rescale color and/or opacity maps used to include current data range
isoVolume1Display.RescaleTransferFunctionToDataRange(True, False)

# show color bar/color legend
isoVolume1Display.SetScalarBarVisibility(renderView1, True)

# get color transfer function/color map for 'alphawater'
alphawaterLUT = GetColorTransferFunction('alphawater')

# get opacity transfer function/opacity map for 'alphawater'
alphawaterPWF = GetOpacityTransferFunction('alphawater')

# create a new 'Integrate Variables'
integrateVariables1 = IntegrateVariables(registrationName='IntegrateVariables1', Input=isoVolume1)

# Create a new 'SpreadSheet View'
spreadSheetView1 = CreateView('SpreadSheetView')
spreadSheetView1.ColumnToSort = ''
spreadSheetView1.BlockSize = 1024

# show data in view
integrateVariables1Display = Show(integrateVariables1, spreadSheetView1, 'SpreadSheetRepresentation')

# get layout
layout1 = GetLayoutByName("Layout #1")

# add view to a layout so it's visible in UI
AssignViewToLayout(view=spreadSheetView1, layout=layout1, hint=0)

# update the view to ensure updated data information
renderView1.Update()

# update the view to ensure updated data information
spreadSheetView1.Update()

# Properties modified on spreadSheetView1
spreadSheetView1.FieldAssociation = 'Cell Data'

# Properties modified on spreadSheetView1
spreadSheetView1.HiddenColumnLabels = []

# Properties modified on spreadSheetView1
spreadSheetView1.HiddenColumnLabels = ['Cell ID', 'Cell Type', 'UMean', 'UMean_Magnitude', 'UPrime2Mean', 'UPrime2Mean_Magnitude', 'U', 'U_Magnitude', 'pMean', 'pPrime2Mean', 'Block Number']

# set active source
SetActiveSource(integrateVariables1)

SelectIDs(IDs=[-1, 0], FieldType=0, ContainingCells=0)

#================================================================
# addendum: following script captures some of the application
# state to faithfully reproduce the visualization during playback
#================================================================

#--------------------------------
# saving layout sizes for layouts

# layout/tab size in pixels
layout1.SetSize(1079, 835)

#-----------------------------------
# saving camera placements for views

# current camera placement for renderView1
renderView1.CameraPosition = [0.163, -0.067, -0.107]
renderView1.CameraFocalPoint = [0.0, 0.0, 0.0065]
renderView1.CameraViewUp = [0.137, 0.927, -0.35]
renderView1.CameraParallelScale = 0.054

#--------------------------------------------
# uncomment the following to render all views
# RenderAllViews()
# alternatively, if you want to write images, you can use SaveScreenshot(...).

What I want to do now is:
1.Get rid of all the graphical aspects of above mentioned python script and make it CLI basis (I run the script using Terminal and the output should takes place)
2.Later I want to ru nthe script using For loop for multiple ‘recorded time’ values for any case.
3.For each ‘recorded time’ the script should copy the alpha.water value and paste it in excel file.

Hence, with that in mind I changed the script as the following way:
Program_moded:

#!usr/bin/python
#script to automate the process of calculating volume fraction (also called as void fraction) alpha 
#and show graphically

#### import the simple module from the paraview
from paraview.simple import *
I = [20, 30, 40, 60, 70, 90]
J = [0001338, 0004013, 0009631, 0024078, 0033441]

for i in I:
	for j in J:
		# create a new 'OpenFOAMReader'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM = OpenFOAMReader(registrationName='3D_Case20_halfGeom_('I')_FlR_0_'J'_Water.OpenFOAM', FileName='3D_Case20_halfGeom_('I')_FlR_0_'J'_Water.OpenFOAM')
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.MeshRegions = ['internalMesh']
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.CellArrays = ['U', 'UMean', 'UPrime2Mean', 'alpha.water', 'alpha.waterMean', 'alpha.waterPrime2Mean', 'p', 'pMean', 'pPrime2Mean', 'p_rgh']

		# update animation scene based on data timesteps
		animationScene1.UpdateAnimationUsingDataTimeSteps()

		# Properties modified on a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM.LabelSize = '64-bit'

		# get color transfer function/color map for 'p'
		pLUT = GetColorTransferFunction('p')

		# get opacity transfer function/opacity map for 'p'
		pPWF = GetOpacityTransferFunction('p')

		# trace defaults for the display properties.
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.Representation = 'Surface'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ColorArrayName = ['POINTS', 'p']
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.LookupTable = pLUT
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectTCoordArray = 'None'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectNormalArray = 'None'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectTangentArray = 'None'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OSPRayScaleArray = 'p'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OSPRayScaleFunction = 'PiecewiseFunction'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectOrientationVectors = 'U'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleFactor = 0.01
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SelectScaleArray = 'p'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GlyphType = 'Arrow'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GlyphTableIndexArray = 'p'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.GaussianRadius = 0.0005
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.SetScaleArray = ['POINTS', 'p']
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleTransferFunction = 'PiecewiseFunction'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityArray = ['POINTS', 'p']
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityTransferFunction = 'PiecewiseFunction'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.DataAxesGrid = 'GridAxesRepresentation'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.PolarAxes = 'PolarAxesRepresentation'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScalarOpacityFunction = pPWF
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScalarOpacityUnitDistance = 0.0018
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityArrayName = ['POINTS', 'p']
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ExtractedBlockIndex = 1

		# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 606.7, 1.0, 0.5, 0.0]

		# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
		a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAMDisplay.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 606.7, 1.0, 0.5, 0.0]

		# create a new 'Clip'
		clip1 = Clip(registrationName='Clip1', Input=a3D_Case20_halfGeom_70_FlR_0_0009631_WaterOpenFOAM)
		clip1.ClipType = 'Plane'
		clip1.HyperTreeGridClipper = 'Plane'
		clip1.Scalars = ['POINTS', 'p']
		clip1.Value = 491.24

		# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
		clip1Display.ScaleTransferFunction.Points = [418.45, 0.0, 0.5, 0.0, 592.33, 1.0, 0.5, 0.0]

		# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
		clip1Display.OpacityTransferFunction.Points = [418.45, 0.0, 0.5, 0.0, 592.33, 1.0, 0.5, 0.0]

		# create a new 'Clip'
		clip2 = Clip(registrationName='Clip2', Input=clip1)
		clip2.ClipType = 'Plane'
		clip2.HyperTreeGridClipper = 'Plane'
		clip2.Scalars = ['POINTS', 'p']
		clip2.Value = 491.24

		# init the 'Plane' selected for 'ClipType'
		clip2.ClipType.Origin = [0, 0.0038, 0.0065]

		# init the 'Plane' selected for 'HyperTreeGridClipper'
		clip2.HyperTreeGridClipper.Origin = [0, 0.0038, 0.0065]

		# Properties modified on clip2.ClipType
		clip2.ClipType.Origin = [0.0, 0.045, 0.0]
		clip2.ClipType.Normal = [0.0, 1.0, 0.0]

		# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
		clip2Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

		# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
		clip2Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

		# create a new 'Clip'
		clip3 = Clip(registrationName='Clip3', Input=clip2)
		clip3.ClipType = 'Plane'
		clip3.HyperTreeGridClipper = 'Plane'
		clip3.Scalars = ['POINTS', 'p']
		clip3.Value = 466.04

		# init the 'Plane' selected for 'ClipType'
		clip3.ClipType.Origin = [0.0, 0.0, 0.0065]

		# init the 'Plane' selected for 'HyperTreeGridClipper'
		clip3.HyperTreeGridClipper.Origin = [0.0, 0.0, 0.0065]

		# Properties modified on clip3
		clip3.Invert = 0

		# Properties modified on clip3.ClipType
		clip3.ClipType.Origin = [0.0, 0.0, 0.0017]
		clip3.ClipType.Normal = [0.0, 0.0, 1.0]

		# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
		clip3Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

		# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
		clip3Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 556.30, 1.0, 0.5, 0.0]

		# create a new 'Iso Volume'
		isoVolume1 = IsoVolume(registrationName='IsoVolume1', Input=clip3)
		isoVolume1.InputScalars = ['POINTS', 'p']
		isoVolume1.ThresholdRange = [375.79, 556.30]

		# Properties modified on isoVolume1
		isoVolume1.InputScalars = ['POINTS', 'alpha.water']
		isoVolume1.ThresholdRange = [0.91, 1.0]

		# init the 'PiecewiseFunction' selected for 'ScaleTransferFunction'
		isoVolume1Display.ScaleTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 554.07, 1.0, 0.5, 0.0]

		# init the 'PiecewiseFunction' selected for 'OpacityTransferFunction'
		isoVolume1Display.OpacityTransferFunction.Points = [375.79, 0.0, 0.5, 0.0, 554.07, 1.0, 0.5, 0.0]

		# hide data in view
		Hide(clip3, renderView1)

		# show color bar/color legend
		isoVolume1Display.SetScalarBarVisibility(renderView1, True)

		# update the view to ensure updated data information
		renderView1.Update()

		# set scalar coloring
		ColorBy(isoVolume1Display, ('POINTS', 'alpha.water'))

		# Hide the scalar bar for this color map if no visible data is colored by it.
		HideScalarBarIfNotNeeded(pLUT, renderView1)

		# rescale color and/or opacity maps used to include current data range
		isoVolume1Display.RescaleTransferFunctionToDataRange(True, False)

		# show color bar/color legend
		isoVolume1Display.SetScalarBarVisibility(renderView1, True)

		# get color transfer function/color map for 'alphawater'
		alphawaterLUT = GetColorTransferFunction('alphawater')

		# get opacity transfer function/opacity map for 'alphawater'
		alphawaterPWF = GetOpacityTransferFunction('alphawater')

		# create a new 'Integrate Variables'
		integrateVariables1 = IntegrateVariables(registrationName='IntegrateVariables1', Input=isoVolume1)

		# Create a new 'SpreadSheet View'
		spreadSheetView1 = CreateView('SpreadSheetView')
		spreadSheetView1.ColumnToSort = ''
		spreadSheetView1.BlockSize = 1024

		# show data in view
		integrateVariables1Display = Show(integrateVariables1, spreadSheetView1, 'SpreadSheetRepresentation')

		# get layout
		layout1 = GetLayoutByName("Layout #1")

		# add view to a layout so it's visible in UI
		AssignViewToLayout(view=spreadSheetView1, layout=layout1, hint=0)

		# update the view to ensure updated data information
		renderView1.Update()

		# update the view to ensure updated data information
		spreadSheetView1.Update()

		# Properties modified on spreadSheetView1
		spreadSheetView1.FieldAssociation = 'Cell Data'

		# Properties modified on spreadSheetView1
		spreadSheetView1.HiddenColumnLabels = []

		# Properties modified on spreadSheetView1
		spreadSheetView1.HiddenColumnLabels = ['Cell ID', 'Cell Type', 'UMean', 'UMean_Magnitude', 'UPrime2Mean', 'UPrime2Mean_Magnitude', 'U', 'U_Magnitude', 'pMean', 'pPrime2Mean', 'Block Number']

		# set active source
		SetActiveSource(integrateVariables1)

		SelectIDs(IDs=[-1, 0], FieldType=0, ContainingCells=0)

Result: Failure, Paraview sent out error, and crashed

Where did I go wrong?
How can I improve it?
Is there anyone having similar scripting idea?

2

You should inspect the Paraview error, or paste the error message here.

3

This is new code then created:

#### import the simple module from the paraview
from paraview.simple import *

#### disable automatic camera reset on 'Show'
paraview.simple._DisableFirstRenderCameraReset()

# create a new 'OpenFOAMReader'
a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam = OpenFOAMReader(registrationName='3D_Case33_halfGeom_(70)_FlR_0_0009142_0-51xMPZ.foam', FileName='3D_Case33_halfGeom_(70)_FlR_0_0009142_0-51xMPZ.foam')
a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam.LabelSize = '64-bit'
a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam.MeshRegions = ['patch/corrugatedwall']
a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam.CellArrays = ['alpha.mpz', 'alpha.mpzMean', 'alpha.mpzPrime2Mean', 'p_rgh']

# get animation scene
animationScene1 = GetAnimationScene()

# update animation scene based on data timesteps
animationScene1.UpdateAnimationUsingDataTimeSteps()

UpdatePipeline(time=1.0, proxy=a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam)

# Properties modified on animationScene1
animationScene1.AnimationTime = 21.0

# get the time-keeper
timeKeeper1 = GetTimeKeeper()

# create a new 'Threshold'
threshold1 = Threshold(registrationName='Threshold1', Input=a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam)
threshold1.Scalars = ['POINTS', 'alpha.mpz']
threshold1.LowerThreshold = 0.91
threshold1.UpperThreshold = 1.0
threshold1.Invert = 0

UpdatePipeline(time=21.0, proxy=threshold1)

# create a new 'Integrate Variables'
integrateVariables1 = IntegrateVariables(registrationName='IntegrateVariables1', Input=threshold1)
# Properties modified on integrateVariables1
integrateVariables1.DivideCellDataByVolume = 0

UpdatePipeline(time=21.0, proxy=integrateVariables1)

This is the error I’m getting:

ERROR: In /builds/gitlab-kitware-sciviz-ci/build/superbuild/paraview/src/VTK/IO/Geometry/vtkOpenFOAMReader.cxx, line 11307
vtkPOpenFOAMReader (0x7efec0004320): Cannot open directory: ./constant/

ERROR: In /builds/gitlab-kitware-sciviz-ci/build/superbuild/paraview/src/VTK/Common/ExecutionModel/vtkExecutive.cxx, line 753
vtkPVCompositeDataPipeline (0x1824cb40): Algorithm vtkPOpenFOAMReader(0x7efec0004320) returned failure for request: vtkInformation (0x18a49750)
  Debug: Off
  Modified Time: 422486
  Reference Count: 1
  Registered Events: (none)
  Request: REQUEST_INFORMATION
  ALGORITHM_AFTER_FORWARD: 1
  FORWARD_DIRECTION: 0



ERROR: In /builds/gitlab-kitware-sciviz-ci/build/superbuild/paraview/src/VTK/IO/Geometry/vtkOpenFOAMReader.cxx, line 11307
vtkPOpenFOAMReader (0x7efec0004320): Cannot open directory: ./constant/

ERROR: In /builds/gitlab-kitware-sciviz-ci/build/superbuild/paraview/src/VTK/Common/ExecutionModel/vtkExecutive.cxx, line 753
vtkPVCompositeDataPipeline (0x1824cb40): Algorithm vtkPOpenFOAMReader(0x7efec0004320) returned failure for request: vtkInformation (0x18a49750)
  Debug: Off
  Modified Time: 422486
  Reference Count: 1
  Registered Events: (none)
  Request: REQUEST_INFORMATION
  ALGORITHM_AFTER_FORWARD: 1
  FORWARD_DIRECTION: 0
4

I’m using ParaView-5.10.
Please I need help

5

Please share your data.

ParaView-5.10

There is no such version yet, are you reffering to 5.10-RC1 ?

6

Yes its 5.10.0-RC1, I thought if its 5.10.0-RC1, it may be referred as 5.10.
I also found the error issue.
in the line:

a3D_Case33_halfGeom_70_FlR_0_0009142_051xMPZfoam = OpenFOAMReader(registrationName='3D_Case33_halfGeom_(70)_FlR_0_0009142_0-51xMPZ.foam', FileName='3D_Case33_halfGeom_(70)_FlR_0_0009142_0-51xMPZ.foam')

if the filename is not having complete path starting from its root directory, then it spits out error about ./constant. directory.
I was trying to implement the filename as just file name present in the folder.
I believe in python if the .py file is in that specific folder where other files are also present which are called during code execution, then having full path starting from root directory does not make sense. Thus, I removed whole path and placed the file in the respective OpenFOAM simulation folder from where the .foam file is called upon execution.
However, I guess it does not work like that, and it spit out error saying, ‘cant open directory ./constant/’.

I’m still having troubles of capturing ‘alpha.water’ value from spreadsheetview inside paraview and store it in designated excel file.
Then after that I want capture these ‘alpha.water’ values over each time steps however that also I’m not able to figure out.
How to do that?

7

Either Copy paste or save to csv file within the spreadsheet view.

Then after that I want capture these ‘alpha.water’ values over each time steps

Not that trivial to do, share your data.

8

The file is the OpenFOAM simulation case
This is shared link.
And for that the paraview python script is uploaded which shows the ‘alpha.water’ in spreadsheetview of paraview. make sure to add complete path name in [filename=’’], the pathname starting from the directory. Othrwise same error may appear, which is posted in above replies of this post.
interfacialAreaMacro-Watr-60-00241.py (5.0 KB)

9

PlotDataOverTime or PlotSelectionOverTime should do the trick