script1

from paraview import simple
from paraview.simple import Show,Render
from paraview.simple import XMLStructuredGridReader
from paraview.simple import IsoVolume
from paraview.simple import GetActiveView,SaveScreenshot
import os

simple._DisableFirstRenderCameraReset()

# List of all files for analysis
os.chdir('/home/metal-machine/para/bin/')
files_list = ['phase_space_vars_0.vts',
              'phase_space_vars_500.vts',
              'phase_space_vars_1000.vts']

# reading vts files from list
view_list = []
iso_list = []

for i in files_list:
    vts_reader = XMLStructuredGridReader(FileName=i)

# Applying IsoVolume filter

    isoVolume1 = IsoVolume(Input=vts_reader)
    isoVolume1.InputScalars = ['POINTS', 'ion_pdf_n_c0']

# Setting-UP threshold values

    isoVolume1.ThresholdRange = [3e-05, 1.0]
    myview = GetActiveView()
    view_list.append(myview)
    iso_list.append(isoVolume1)

Show(iso_list[0],view_list[0])

# Rendering to png
SaveScreenshot("currentview.png" , view = view_list[0] )