pyvale 2025.4.0__py3-none-any.whl

pyvale/examples/analyticdatagen/ex2_1_analyticsensors.py

@@ -0,0 +1,83 @@
+'''
+================================================================================
+Analytic test case data - linear
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+'''
+import matplotlib.pyplot as plt
+import numpy as np
+import pyvale
+
+def main() -> None:
+    (sim_data,_) = pyvale.AnalyticCaseFactory.scalar_linear_2d()
+
+    descriptor = pyvale.SensorDescriptorFactory.temperature_descriptor()
+
+    field_key = 'scalar'
+    t_field = pyvale.FieldScalar(sim_data,
+                                 field_key=field_key,
+                                 spat_dims=2)
+
+    n_sens = (4,1,1)
+    x_lims = (0.0,10.0)
+    y_lims = (0.0,7.5)
+    z_lims = (0.0,0.0)
+    sens_pos = pyvale.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+
+    use_sim_time = False
+    if use_sim_time:
+        sample_times = None
+    else:
+        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
+
+    sensor_data = pyvale.SensorData(positions=sens_pos,
+                                         sample_times=sample_times)
+
+    tc_array = pyvale.SensorArrayPoint(sensor_data,
+                                       t_field,
+                                       descriptor)
+
+    errors_on = {'indep_sys': True,
+                 'rand': True,
+                 'dep_sys': True}
+
+    error_chain = []
+    if errors_on['indep_sys']:
+        error_chain.append(pyvale.ErrSysOffset(offset=-5.0))
+        error_chain.append(pyvale.ErrSysUniform(low=-5.0,
+                                            high=5.0))
+        gen_norm = pyvale.GeneratorNormal(std=1.0)
+
+    if errors_on['rand']:
+        error_chain.append(pyvale.ErrRandNormPercent(std_percent=1.0))
+        error_chain.append(pyvale.ErrRandUnifPercent(low_percent=-1.0,
+                                            high_percent=1.0))
+
+    if errors_on['dep_sys']:
+        error_chain.append(pyvale.ErrSysDigitisation(bits_per_unit=2**8/100))
+        error_chain.append(pyvale.ErrSysSaturation(meas_min=0.0,meas_max=300.0))
+
+    if len(error_chain) > 0:
+        error_integrator = pyvale.ErrIntegrator(error_chain,
+                                                  sensor_data,
+                                                  tc_array.get_measurement_shape())
+        tc_array.set_error_integrator(error_integrator)
+
+    measurements = tc_array.get_measurements()
+
+    pyvale.print_measurements(tc_array,
+                            (0,1), # Sensor 1
+                            (0,1), # Component 1: scalar field = 1 component
+                            (measurements.shape[2]-5,measurements.shape[2]))
+
+    (fig,_) = pyvale.plot_time_traces(tc_array,field_key)
+    plt.show()
+
+    pv_plot = pyvale.plot_point_sensors_on_sim(tc_array,field_key)
+    pv_plot.show(cpos="xy")
+
+if __name__ == '__main__':
+    main()

pyvale/examples/ex1_1_thermal2d.py

@@ -0,0 +1,89 @@
+"""
+================================================================================
+Example: thermocouples on a 2d plate
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+"""
+from pathlib import Path
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: point sensors on a 2D thermal simulation
+    ----------------------------------------------------------------------------
+    - Quick start
+    - Basic sensor array construction using the sensor array factory
+    - Basic visualisation of sensor locations and sensor traces with the pyvale
+      wrapper for pyvista and matplotlib.
+    """
+
+    data_path = pyv.DataSet.thermal_2d_path()
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    field_key = "temperature"
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    n_sens = (3,2,1)
+    x_lims = (0.0,100.0)
+    y_lims = (0.0,50.0)
+    z_lims = (0.0,0.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+    sens_data = pyv.SensorData(positions=sens_pos)
+
+    tc_array = pyv.SensorArrayFactory \
+        .thermocouples_basic_errs(sim_data,
+                                  sens_data,
+                                  field_key,
+                                  spat_dims=2)
+
+    measurements = tc_array.get_measurements()
+    print(f"\nMeasurements for last sensor:\n{measurements[-1,0,:]}\n")
+
+    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_key)
+    # Set this to "interactive" to get an interactive 3D plot of the simulation
+    # and labelled sensor locations, set to "save_fig" to create a vector
+    # graphic using a specified camera position.
+    pv_plot_mode = "interactive"
+
+    if pv_plot_mode == "interactive":
+        pv_plot.camera_position = [(-7.547, 59.753, 134.52),
+                                   (41.916, 25.303, 9.297),
+                                   (0.0810, 0.969, -0.234)]
+        pv_plot.show()
+
+        print(80*"=")
+        print("Camera positions = ")
+        print(pv_plot.camera_position)
+        print(80*"="+"\n")
+
+    if pv_plot_mode == "save_fig":
+        # Determined manually by moving camera and then dumping camera position
+        # to console after window close - see "interactive above"
+        pv_plot.camera_position = [(-7.547, 59.753, 134.52),
+                                   (41.916, 25.303, 9.297),
+                                   (0.0810, 0.969, -0.234)]
+        save_render = Path("src/examples/plate_thermal_2d_sim_view.svg")
+        pv_plot.save_graphic(save_render) # only for .svg .eps .ps .pdf .tex
+        pv_plot.screenshot(save_render.with_suffix(".png"))
+
+    # Set this to "interactive" to get a matplotlib.pyplot with the sensor
+    # traces plotted over time. Set to "save_fig" to save an image of the plot
+    # to file.
+    trace_plot_mode = "interactive"
+
+    (fig,_) = pyv.plot_time_traces(tc_array,field_key)
+
+    if trace_plot_mode == "interactive":
+        plt.show()
+    if trace_plot_mode == "save_fig":
+        save_traces = Path("src/examples/plate_thermal_2d_traces.png")
+        fig.savefig(save_traces, dpi=300, format="png", bbox_inches="tight")
+
+
+if __name__ == "__main__":
+    main()

pyvale/examples/ex1_2_thermal2d.py

@@ -0,0 +1,111 @@
+"""
+================================================================================
+Example: thermocouples on a 2d plate
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+"""
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: Point sensors on a 2D thermal simulation
+    ----------------------------------------------------------------------------
+    - Explanation of the usage of "get_measurements()" and "calc_measurements()"
+
+    NOTES:
+    - A sensor measurement is defined as:
+        measurement = truth + systematic error + random error.
+    - Calling the "get" methods of the sensor array will retrieve the results
+      for the current experiment.
+    - Calling the "calc" methods will generate a new
+      experiment by sampling/calculating the systematic and random errors.
+    """
+    data_path = pyv.DataSet.thermal_2d_path()
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    field_key = list(sim_data.node_vars.keys())[0] # type: ignore
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    n_sens = (4,1,1)
+    x_lims = (0.0,100.0)
+    y_lims = (0.0,50.0)
+    z_lims = (0.0,0.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+    sens_data = pyv.SensorData(positions=sens_pos)
+
+    tc_array = pyv.SensorArrayFactory \
+        .thermocouples_basic_errs(sim_data,
+                                  sens_data,
+                                  field_key,
+                                  spat_dims=2)
+
+
+    measurements = tc_array.get_measurements()
+
+    print("\n"+80*"-")
+    print("For a sensor: measurement = truth + sysematic error + random error")
+    print(f"measurements.shape = {measurements.shape} = "+
+          "(n_sensors,n_field_components,n_timesteps)\n")
+    print("The truth, systematic error and random error arrays have the same "+
+          "shape.")
+
+    print(80*"-")
+    print("Looking at the last 5 time steps (measurements) of sensor 0:")
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+    print(80*"-")
+    print("If we call the \"calc_measurements()\" method then the errors are "+
+          "(re)calculated or sampled.")
+    measurements = tc_array.calc_measurements()
+
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+
+
+    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
+    ax.set_title("Exp 1: called calc_measurements()")
+
+    print(80*"-")
+    print("If we call the \"get_measurements()\" method then the errors are the "+
+          "same:")
+    measurements = tc_array.get_measurements()
+
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+
+    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
+    ax.set_title("Exp 2: called get_measurements()")
+
+    print(80*"-")
+    print("If we call the \"calc_measurements()\" method again we generate/sample"+
+           "new errors:")
+    measurements = tc_array.calc_measurements()
+
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+
+    (_,ax) = pyv.plot_time_traces(tc_array,field_key)
+    ax.set_title("Exp 3: called calc_measurements()")
+
+    print(80*"-")
+
+    plot_on = True
+    if plot_on:
+        plt.show()
+
+
+if __name__ == "__main__":
+    main()

pyvale/examples/ex1_3_thermal2d.py

@@ -0,0 +1,113 @@
+"""
+================================================================================
+Example: thermocouples on a 2d plate
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+"""
+import numpy as np
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: Point sensors on a 2D thermal simulation
+    ----------------------------------------------------------------------------
+    - Full construction of a point sensor array from scratch
+    - Explanation of the different types of error models
+    - There are flags throughout the example allowing the user to toggle on/off
+      parts of the sensor array construction
+    """
+    data_path = pyv.DataSet.thermal_2d_path()
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    use_auto_descriptor = "blank"
+    if use_auto_descriptor == "factory":
+        descriptor = pyv.SensorDescriptorFactory.temperature_descriptor()
+    elif use_auto_descriptor == "manual":
+        descriptor = pyv.SensorDescriptor()
+        descriptor.name = "Temperature"
+        descriptor.symbol = "T"
+        descriptor.units = r"^{\circ}C"
+        descriptor.tag = "TC"
+    else:
+        descriptor = pyv.SensorDescriptor()
+
+    field_key = "temperature" # ("disp_x","disp_y")
+    t_field = pyv.FieldScalar(sim_data,
+                                 field_key=field_key,
+                                 spat_dims=2)
+
+    n_sens = (4,1,1)
+    x_lims = (0.0,100.0)
+    y_lims = (0.0,50.0)
+    z_lims = (0.0,0.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+
+    use_sim_time = False
+    if use_sim_time:
+        sample_times = None
+    else:
+        sample_times = np.linspace(0.0,np.max(sim_data.time),50)
+
+    sensor_data = pyv.SensorData(positions=sens_pos,
+                                   sample_times=sample_times)
+
+    tc_array = pyv.SensorArrayPoint(sensor_data,
+                                       t_field,
+                                       descriptor)
+
+    errors_on = {"indep_sys": True,
+                 "rand": True,
+                 "dep_sys": True}
+
+    error_chain = []
+    if errors_on["indep_sys"]:
+        error_chain.append(pyv.ErrSysOffset(offset=-5.0))
+        error_chain.append(pyv.ErrSysUniform(low=-5.0,
+                                                high=5.0))
+
+    if errors_on["rand"]:
+        error_chain.append(pyv.ErrRandNormPercent(std_percent=1.0))
+        error_chain.append(pyv.ErrRandUnifPercent(low_percent=-1.0,
+                                            high_percent=1.0))
+
+    if errors_on["dep_sys"]:
+        error_chain.append(pyv.ErrSysDigitisation(bits_per_unit=2**8/100))
+        error_chain.append(pyv.ErrSysSaturation(meas_min=0.0,meas_max=300.0))
+
+    if len(error_chain) > 0:
+        error_integrator = pyv.ErrIntegrator(error_chain,
+                                                  sensor_data,
+                                                  tc_array.get_measurement_shape())
+        tc_array.set_error_integrator(error_integrator)
+
+
+    measurements = tc_array.get_measurements()
+
+    print("\n"+80*"-")
+    print("For a sensor: measurement = truth + sysematic error + random error")
+    print(f"measurements.shape = {measurements.shape} = "+
+          "(n_sensors,n_field_components,n_timesteps)\n")
+    print("The truth, systematic error and random error arrays have the same "+
+          "shape.")
+
+    print(80*"-")
+    print("Looking at the last 5 time steps (measurements) of sensor 0:")
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+    print(80*"-")
+
+    pyv.plot_time_traces(tc_array,field_key)
+    plt.show()
+
+
+if __name__ == "__main__":
+    main()

pyvale/examples/ex1_4_thermal2d.py

@@ -0,0 +1,89 @@
+"""
+================================================================================
+Example: thermocouples on a 2d plate
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+"""
+import numpy as np
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: point sensors on a 2D thermal simulation
+    ----------------------------------------------------------------------------
+    - Demonstrates options for controlling plots of points sensor traces using
+      matplotlib
+    """
+    data_path = pyv.DataSet.thermal_2d_path()
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    field_key = list(sim_data.node_vars.keys())[0] # type: ignore
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    n_sens = (4,1,1)
+    x_lims = (0.0,100.0)
+    y_lims = (0.0,50.0)
+    z_lims = (0.0,0.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+
+    sample_times = np.linspace(0.0,np.max(sim_data.time),12)
+
+    sens_data = pyv.SensorData(positions=sens_pos,
+                                  sample_times=sample_times)
+
+    tc_array = pyv.SensorArrayFactory \
+        .thermocouples_basic_errs(sim_data,
+                                  sens_data,
+                                  field_key,
+                                  spat_dims=2)
+
+    err_int = pyv.ErrIntegrator([pyv.ErrSysOffset(offset=-5.0)],
+                                     sens_data,
+                                     tc_array.get_measurement_shape())
+    tc_array.set_error_integrator(err_int)
+
+    measurements = tc_array.get_measurements()
+
+    print(80*"-")
+    print("Looking at the last 5 time steps (measurements) of sensor 0:")
+    pyv.print_measurements(tc_array,
+                              (0,1),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+    print(80*"-")
+
+    trace_props = pyv.TraceOptsSensor()
+
+    trace_props.truth_line = None
+    trace_props.sim_line = None
+    pyv.plot_time_traces(tc_array,field_key,trace_props)
+
+    trace_props.meas_line = "--o"
+    trace_props.truth_line = "-x"
+    trace_props.sim_line = ":+"
+    pyv.plot_time_traces(tc_array,field_key,trace_props)
+
+    trace_props.sensors_to_plot = np.arange(measurements.shape[0]-2
+                                           ,measurements.shape[0])
+    pyv.plot_time_traces(tc_array,field_key,trace_props)
+
+    trace_props.sensors_to_plot = None
+    trace_props.time_min_max = (0.0,100.0)
+    pyv.plot_time_traces(tc_array,field_key,trace_props)
+
+    plt.show()
+
+    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_key)
+    pv_plot.camera_position = [(-7.547, 59.753, 134.52),
+                                   (41.916, 25.303, 9.297),
+                                   (0.0810, 0.969, -0.234)]
+    pv_plot.show()
+
+
+if __name__ == "__main__":
+    main()

pyvale/examples/ex1_5_thermal2d.py

@@ -0,0 +1,105 @@
+'''
+================================================================================
+Example: thermocouples on a 2d plate
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+'''
+import numpy as np
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: point sensors on a 2D thermal simulation
+    ----------------------------------------------------------------------------
+    """
+    data_path = pyv.DataSet.thermal_2d_path()
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    field_key = list(sim_data.node_vars.keys())[0] # type: ignore
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    n_sens = (4,1,1)
+    x_lims = (0.0,100.0)
+    y_lims = (0.0,50.0)
+    z_lims = (0.0,0.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+
+    sample_times = np.linspace(0.0,np.max(sim_data.time),50) # | None
+
+    sensor_data = pyv.SensorData(positions=sens_pos,
+                                    sample_times=sample_times)
+
+    tc_array = pyv.SensorArrayFactory \
+        .thermocouples_no_errs(sim_data,
+                               sensor_data,
+                               field_key,
+                               spat_dims=2)
+
+    #===========================================================================
+    # Examples of full error library
+
+    #---------------------------------------------------------------------------
+    # Standard independent systematic errors
+    err_chain = []
+    err_chain.append(pyv.ErrSysOffset(offset=-1.0))
+    err_chain.append(pyv.ErrSysOffsetPercent(offset_percent=-1.0))
+
+    err_chain.append(pyv.ErrSysUniform(low=-2.0,
+                                        high=2.0))
+    err_chain.append(pyv.ErrSysUniformPercent(low_percent=-2.0,
+                                                    high_percent=2.0))
+
+    err_chain.append(pyv.ErrSysNormal(std=1.0))
+    err_chain.append(pyv.ErrSysNormPercent(std_percent=2.0))
+
+    sys_gen = pyv.GeneratorTriangular(left=-1.0,
+                                         mode=0.0,
+                                         right=1.0)
+    err_chain.append(pyv.ErrSysGenerator(sys_gen))
+
+    #---------------------------------------------------------------------------
+    err_chain.append(pyv.ErrRandNormal(std = 2.0))
+    err_chain.append(pyv.ErrRandNormPercent(std_percent=2.0))
+
+    err_chain.append(pyv.ErrRandUniform(low=-2.0,high=2.0))
+    err_chain.append(pyv.ErrRandUnifPercent(low_percent=-2.0,
+                                               high_percent=2.0))
+
+    rand_gen = pyv.GeneratorTriangular(left=-5.0,
+                                          mode=0.0,
+                                          right=5.0)
+    err_chain.append(pyv.ErrRandGenerator(rand_gen))
+
+    #---------------------------------------------------------------------------
+    err_chain.append(pyv.ErrSysDigitisation(bits_per_unit=2**8/100))
+    err_chain.append(pyv.ErrSysSaturation(meas_min=0.0,meas_max=300.0))
+
+    err_int = pyv.ErrIntegrator(err_chain,
+                                     sensor_data,
+                                     tc_array.get_measurement_shape())
+    tc_array.set_error_integrator(err_int)
+
+
+    #===========================================================================
+
+    measurements = tc_array.calc_measurements()
+    print(80*'-')
+    sens_num = 4
+    print('The last 5 time steps (measurements) of sensor {sens_num}:')
+    pyv.print_measurements(tc_array,
+                              (sens_num-1,sens_num),
+                              (0,1),
+                              (measurements.shape[2]-5,measurements.shape[2]))
+    print(80*'-')
+
+    pyv.plot_time_traces(tc_array,field_key)
+    plt.show()
+
+
+if __name__ == '__main__':
+    main()

pyvale/examples/ex2_1_thermal3d .py

@@ -0,0 +1,87 @@
+'''
+================================================================================
+Example: 3d thermocouples on a monoblock
+
+pyvale: the python validation engine
+License: MIT
+Copyright (C) 2025 The Computer Aided Validation Team
+================================================================================
+'''
+from pathlib import Path
+import matplotlib.pyplot as plt
+import mooseherder as mh
+import pyvale as pyv
+
+
+def main() -> None:
+    """pyvale example: thermocouples on a 3D divertor monoblock heatsink
+    ----------------------------------------------------------------------------
+    """
+    #data_path = pyv.DataSet.thermal_3d_path()
+    data_path = Path.cwd()/"src"/"pyvale"/"simcases"/"case16_out.e"
+    sim_data = mh.ExodusReader(data_path).read_all_sim_data()
+    field_name = 'temperature'
+    # Scale to mm to make 3D visualisation scaling easier
+    sim_data.coords = sim_data.coords*1000.0 # type: ignore
+
+    n_sens = (1,4,1)
+    x_lims = (12.5,12.5)
+    y_lims = (0.0,33.0)
+    z_lims = (0.0,12.0)
+    sens_pos = pyv.create_sensor_pos_array(n_sens,x_lims,y_lims,z_lims)
+
+    sens_data = pyv.SensorData(positions=sens_pos)
+
+    tc_array = pyv.SensorArrayFactory \
+        .thermocouples_basic_errs(sim_data,
+                                  sens_data,
+                                  field_name,
+                                  spat_dims=3)
+
+    measurements = tc_array.get_measurements()
+    print(f'\nMeasurements for sensor at top of block:\n{measurements[-1,0,:]}\n')
+
+    pv_plot = pyv.plot_point_sensors_on_sim(tc_array,field_name)
+
+    # Set this to 'interactive' to get an interactive 3D plot of the simulation
+    # and labelled sensor locations, set to 'save_fig' to create a vector
+    # graphic using a specified camera position.
+    pv_plot_mode = 'interactive'
+
+    if pv_plot_mode == 'interactive':
+        pv_plot.camera_position = [(59.354, 43.428, 69.946),
+                                   (-2.858, 13.189, 4.523),
+                                   (-0.215, 0.948, -0.233)]
+        pv_plot.show()
+
+        print(80*"=")
+        print('Camera positions = ')
+        print(pv_plot.camera_position)
+        print(80*"="+"\n")
+
+    if pv_plot_mode == 'save_fig':
+        # Determined manually by moving camera and then dumping camera position
+        # to console after window close - see 'interactive above'
+        pv_plot.camera_position = [(59.354, 43.428, 69.946),
+                                   (-2.858, 13.189, 4.523),
+                                   (-0.215, 0.948, -0.233)]
+        save_render = Path('src/examples/monoblock_thermal_sim_view.svg')
+        pv_plot.save_graphic(save_render) # only for .svg .eps .ps .pdf .tex
+        pv_plot.screenshot(save_render.with_suffix('.png'))
+
+    # Set this to 'interactive' to get a matplotlib.pyplot with the sensor
+    # traces plotted over time. Set to 'save_fig' to save an image of the plot
+    # to file.
+    trace_plot_mode = 'interactive'
+
+    (fig,_) = pyv.plot_time_traces(tc_array,field_name)
+
+    if trace_plot_mode == 'interactive':
+        plt.show()
+    if trace_plot_mode == 'save_fig':
+        save_traces = Path('src/examples/monoblock_thermal_traces.png')
+        fig.savefig(save_traces, dpi=300, format='png', bbox_inches='tight')
+
+
+if __name__ == '__main__':
+    main()