[yeharav]

Hi all,

I am trying to load OF decomposed case using the python simple of pv4.4.

So far I used the script

X = OpenFOAMReader(FileName='X.foam')
X.MeshRegions = ['internalMesh']
X.CellArrays = ['U']
X.CaseType = 'Decomposed Case'
print(X.TimestepValues)

in OpenFOAM python console.

I see that I have to press the 'apply' button for things to load (and get the decomposed case).

What is the scripting equivalent. I have tried X.ReadPipeline() and X.ReadPipelineInformation(), and X.refresh() to no avail.

Any Ideas?

[tomf]

Hi,

I would suggest to use the trace functionality. This would mean to do everything manually once, while tracing the steps, at the end stop the trace and than modify the python script that paraview gives you to your needs.

Regards,
Tom

[yeharav]

Hi,

Thanks for the suggestion.

The script was actually produced by the trace feature.

My aim is to use the script outside of paraview (pvbatch, for example).
where the code that was produced by the trace feature did not load the time steps of the decomposed case.

When I used it in the paraview, I saw that I have to press the 'apply' button for things to work. However, I will try it again.

[yeharav]

The script is

X OpenFOAMReader(FileName='X.foam',CaseType = 'Decomposed Case')
X.MeshRegions = ['internalMesh']
X.CellArrays = ['U']

X.UpdatePipeline([the time]).

To access the data we use the

import vtk.numpy_interface.dataset_adapter as dsa

slice1 = Slice(X) 
rawdata =  servermanager.Fetch(slice)
data = dsa.WrapDataObject(rawData)
data.PointData ...