ansys-pyensight-core 0.8.0__tar.gz → 0.8.1__tar.gz

{ansys_pyensight_core-0.8.0 → ansys_pyensight_core-0.8.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ansys-pyensight-core
-Version: 0.8.0
+Version: 0.8.1
 Summary: A python wrapper for Ansys EnSight
 Author-email: "ANSYS, Inc." <pyansys.core@ansys.com>
 Maintainer-email: "ANSYS, Inc." <pyansys.core@ansys.com>

{ansys_pyensight_core-0.8.0 → ansys_pyensight_core-0.8.1}/pyproject.toml

@@ -6,7 +6,7 @@ build-backend = "flit_core.buildapi"
 
 [project]
 name = "ansys-pyensight-core"
-version = "0.8.0"
+version = "0.8.1"
 description = "A python wrapper for Ansys EnSight"
 readme = "README.rst"
 requires-python = ">=3.9,<4"

{ansys_pyensight_core-0.8.0 → ansys_pyensight_core-0.8.1}/src/ansys/pyensight/core/session.py

@@ -1778,3 +1778,38 @@ class Session:
                 message = f"A newer version of EnSight is required to use this API:{base_msg}"
             raise InvalidEnSightVersion(message)
         return valid
+
+    def find_remote_unused_ports(
+        self,
+        count: int,
+        start: Optional[int] = None,
+        end: Optional[int] = None,
+        avoid: Optional[list[int]] = None,
+    ) -> Optional[List[int]]:
+        """
+        Find "count" unused ports on the host system.  A port is considered
+        unused if it does not respond to a "connect" attempt.  Walk the ports
+        from 'start' to 'end' looking for unused ports and avoiding any ports
+        in the 'avoid' list.  Stop once the desired number of ports have been
+        found.  If an insufficient number of ports were found, return None.
+        An admin user check is used to skip [1-1023].
+
+        Parameters
+        ----------
+        count: int
+            number of unused ports to find
+        start: int
+            the first port to check or None (random start)
+        end: int
+            the last port to check or None (full range check)
+        avoid: list
+            an optional list of ports not to check
+
+        Returns
+        -------
+            the detected ports or None on failure
+        """
+        cmd = "from cei import find_unused_ports\n"
+        cmd += f"ports = find_unused_ports({count}, start={start}, end={end}, avoid={avoid})"
+        self.cmd(cmd, do_eval=False)
+        return self.cmd("ports")

ansys_pyensight_core-0.8.1/src/ansys/pyensight/core/utils/readers.py

@@ -0,0 +1,270 @@
+"""Readers module.
+
+This module contains utilities to do readers specific operations.
+"""
+import os
+import re
+from threading import Thread
+import time
+from types import ModuleType
+from typing import Optional, Tuple, Union
+import uuid
+
+from ansys.pyensight.core import LocalLauncher
+
+try:
+    import ensight
+except ImportError:
+    from ansys.api.pyensight import ensight_api
+
+
+class Readers:
+    """A namespace to access the interfaces of specific readers"""
+
+    def __init__(self, interface: Union["ensight_api.ensight", "ensight"]):
+        self._ensight = interface
+        self._dvs = DVS(self._ensight)
+
+    @property
+    def dvs(self):
+        return self._dvs
+
+
+class DVS:
+    """A namespace to access specific DVS interfaces"""
+
+    def __init__(self, interface: Union["ensight_api.ensight", "ensight"]):
+        self._ensight: Union["ensight_api.ensight", "ensight"] = interface
+        self._dvs_port: Optional[int] = None
+
+    MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT = "stay_at_current"
+    MONITOR_NEW_TIMESTEPS_JUMP_TO_END = "jump_to_end"
+
+    def _launch_live_dvs_242_cmd(
+        self,
+        tmp_name: str,
+        port: int = 0,
+        secret_key: Optional[str] = None,
+        monitor_new_timesteps: str = MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT,
+    ):
+        """Launch a live DVS session from PyEnSight for EnSight 242 and onward"""
+        indent = "    "
+        cmd = "import re\n"
+        cmd += "def dvs_callback():\n"
+        cmd += f'{indent}command_string = f"set_server_port={port}"\n'
+        if secret_key is not None:
+            secret_string = f'{indent}command_string += f"&set_secret_key='
+            secret_string += secret_key + '"\n'
+            cmd += secret_string
+        tmp_name = tmp_name.replace("\\", "\\\\")
+        cmd += f"{indent}reply = ensight.objs.core.CURRENTCASE[0].client_command(command_string)\n"
+        cmd += f"{indent}dvs_port_string = ensight.objs.core.CURRENTCASE[0].client_command('get_server_port')\n"
+        cmd += (
+            f"{indent}dvs_port = int(re.search(r':([0-9]{{4,5}})', dvs_port_string).groups(1)[0])\n"
+        )
+        cmd += f"{indent}with open('{tmp_name}', 'w') as temp_file:\n"
+        cmd += f"{2*indent}temp_file.write(str(dvs_port))\n"
+        cmd += f"{indent}return True\n\n"
+        cmd += "reply = ensight.objs.core.CURRENTCASE[0].client_command_callback(dvs_callback)\n"
+        if monitor_new_timesteps:
+            cmd += f'ensight.solution_time.monitor_for_new_steps("{monitor_new_timesteps}")\n'
+        cmd += 'ensight.part.elt_representation("3D_feature_2D_full")\n'
+        cmd += 'err = ensight.data.replace("notexisting.dvs")\n'
+        return cmd
+
+    def _launch_live_dvs_241_cmd(
+        self,
+        port: int = 0,
+        secret_key: Optional[str] = None,
+        monitor_new_timesteps: str = MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT,
+    ):
+        """Launch a live DVS session from PyEnSight for EnSight 241"""
+        if port == 0:
+            ports = self._ensight._session.find_remote_unused_ports(1)
+            if ports:
+                self._dvs_port = ports[0]
+                port = ports[0]
+        indent = "    "
+        cmd = "import os\n"
+        cmd += f"try:\n{indent}os.remove('remote.dvs')\n"
+        cmd += f"except (FileNotFoundError, OSError):\n{indent}pass\n"
+        cmd += "path = os.path.join(os.getcwd(), 'remote.dvs')\n"
+        cmd += "with open(path, 'w') as dvsfile:\n"
+        cmd += f"{indent}dvsfile.write('#!DVS_CASE 1.0\\n')\n"
+        cmd += f"{indent}dvsfile.write(f'SERVER_PORT_BASE={port}\\n')\n"
+        cmd += f"{indent}dvsfile.write('SERVER_PORT_MULT=1\\n')\n"
+        if secret_key:
+            cmd += f"{indent}dvsfile.write(f'SERVER_SECURITY_SECRET={secret_key}\\n')\n"
+        cmd += f'ensight.solution_time.monitor_for_new_steps("{monitor_new_timesteps}")\n'
+        cmd += 'ensight.part.elt_representation("3D_feature_2D_full")\n'
+        cmd += "ensight.data.replace(path)\n"
+        return cmd
+
+    def _find_dvs_port(self, tmp_name: Optional[str] = None):
+        """Find the dvs port allocated from the input temporary name"""
+        if not tmp_name:
+            raise RuntimeError("Temporary name for dvs port file not available")
+        is_local = isinstance(self._ensight._session._launcher, LocalLauncher)
+        if is_local:
+            with open(tmp_name) as dvs_port_file:
+                self._dvs_port = int(dvs_port_file.read().strip())
+        else:
+            log_content = self._ensight._session._launcher.enshell_log_contents()
+            dvs_port_match = re.search(r"\(0.0.0.0\):([0-9]{4,5})\n", log_content)
+            if dvs_port_match:
+                self._dvs_port = int(str(dvs_port_match.groups(1)[0]).strip())
+
+    @staticmethod
+    def _launch_dvs_callback_in_ensight(
+        port: int,
+        filename: Optional[str],
+        secret_key: Optional[str],
+        monitor_new_timesteps: str = MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT,
+    ) -> None:  # pragma: no cover
+        """Launch a live DVS session in EnSight"""
+        from ceiversion import ensight_suffix  # pylint: disable=import-outside-toplevel
+
+        if int(ensight_suffix) < 242:
+            from cei import find_unused_ports  # pylint: disable=import-outside-toplevel
+
+            if port == 0:
+                ports = find_unused_ports(1)
+                if ports:
+                    port = ports[0]
+            try:
+                os.remove("remote.dvs")
+            except (FileNotFoundError, OSError):
+                pass
+            path = os.path.join(os.getcwd(), "remote.dvs")
+            with open(path, "w") as dvsfile:
+                dvsfile.write("#!DVS_CASE 1.0\n")
+                dvsfile.write(f"SERVER_PORT_BASE={port}\n")
+                dvsfile.write("SERVER_PORT_MULT=1\n")
+                if secret_key is not None:
+                    dvsfile.write(f"SERVER_SECURITY_SECRET={secret_key}\n")
+            if filename:
+                try:
+                    with open(filename, "w") as dvs_port_file:
+                        dvs_port_file.write(str(port))
+                except Exception:
+                    raise RuntimeError(f"Couldn't write allocated DVS port to {filename}")
+            ensight.part.elt_representation("3D_feature_2D_full")
+            ensight.solution_time.monitor_for_new_steps(f"{monitor_new_timesteps}")
+            ensight.data.replace(path)
+        else:
+
+            def dvs_callback():
+                command_string = f"set_server_port={port}"
+                if secret_key is not None:
+                    command_string += f"&set_secret_key={secret_key}"
+                ensight.objs.core.CURRENTCASE[0].client_command(command_string)
+                dvs_port_string = ensight.objs.core.CURRENTCASE[0].client_command("get_server_port")
+                dvs_port_match = re.search(r":([0-9]{4,5})", dvs_port_string)
+                dvs_port = None
+                if dvs_port_match:
+                    dvs_port = int(dvs_port_match.groups(1)[0])
+                if not dvs_port:
+                    raise RuntimeError("DVS couldn't allocate a port")
+                if filename:
+                    try:
+                        with open(filename, "w") as dvs_port_file:
+                            dvs_port_file.write(str(dvs_port))
+                    except Exception:
+                        raise RuntimeError(f"Couldn't write allocated DVS port to {filename}")
+                return True
+
+            ensight.objs.core.CURRENTCASE[0].client_command_callback(dvs_callback)
+            ensight.solution_time.monitor_for_new_steps(f"{monitor_new_timesteps}")
+            ensight.part.elt_representation("3D_feature_2D_full")
+            ensight.data.replace("notexisting.dvs")
+
+    def launch_live_dvs(
+        self,
+        port: int = 0,
+        secret_key: Optional[str] = None,
+        monitor_new_timesteps: str = MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT,
+        start_thread: bool = True,
+        filename: Optional[str] = None,
+    ) -> Tuple[Optional[Thread], Optional[int]]:
+        """To provide an interface to launch an in-situ EnSight DVS session.
+        If in PyEnSight, the function will return a thread which will launch the DVS reader
+        in EnSight, hence the DVS servers, and will also return the port allocated by DVS, to
+        cover the case the port 0 was asked for.
+        If instead the function will be used directly in EnSight via the utils interface, since
+        the reader launch will block the interpreter waiting for new data, the port cannot be returned
+        and cannot be printed up until the first update happens.
+        So, if you need to access the port number in a script and you cannot check the EnSight console,
+        please supply a filename to store the dvs port into.
+
+        Parameters
+        ----------
+        port: int
+            the port number where the first DVS server will be started. In case of a
+            SOS EnSight session, on the following server the DVS servers will be started on the
+            next port, e.g. if the first server starts at port 50055, the second will start
+            at port 50056 and so on
+        secret_key: str
+            an optional secret key to pass in case the DVS clients have been started with a secret key
+            for the underlying gRPC connections. An empty string can be provided if needed
+        monitor_new_timesteps: str
+            set the way EnSight will monitor for new timesteps. Defaults to MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT.
+            The allowed values are MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT
+            and MONITOR_NEW_TIMESTEPS_JUMP_TO_END
+        start_thread: bool
+            True if the thread to be returned needs to be started already. Default is True
+        filename: str
+            An optional path to store the port number in. It will be used only if the utils is being
+            called directly in EnSight.
+
+        Returns
+        -------
+        Thread:
+            a python Thread which holds the dvs load
+        """
+
+        def load_dvs():
+            self._ensight._session.cmd(cmd, do_eval=False)
+
+        if monitor_new_timesteps not in [
+            self.MONITOR_NEW_TIMESTEPS_JUMP_TO_END,
+            self.MONITOR_NEW_TIMESTEPS_STAY_AT_CURRENT,
+        ]:
+            raise RuntimeError(
+                f"{monitor_new_timesteps} value not allowed for an in-situ DVS session"
+            )
+        if not isinstance(self._ensight, ModuleType):  # pragma: no cover
+            self._ensight._session.ensight_version_check("2024 R1")
+        else:
+            self._launch_dvs_callback_in_ensight(
+                port=port,
+                secret_key=secret_key,
+                monitor_new_timesteps=monitor_new_timesteps,
+                filename=filename,
+            )
+            return None, None
+        cmd = ""
+        path = None
+        if int(self._ensight._session.cei_suffix) < 242:
+            cmd = self._launch_live_dvs_241_cmd(
+                port=port, secret_key=secret_key, monitor_new_timesteps=monitor_new_timesteps
+            )
+        else:
+            tmp_name = str(uuid.uuid4())
+            path = os.path.join(self._ensight._session._launcher.session_directory, tmp_name)
+            cmd = self._launch_live_dvs_242_cmd(
+                port=port,
+                secret_key=secret_key,
+                monitor_new_timesteps=monitor_new_timesteps,
+                tmp_name=path,
+            )
+        t = Thread(target=load_dvs)
+        if start_thread:
+            t.start()
+        start = time.time()
+        while not self._dvs_port and time.time() - start < 60:
+            try:
+                self._find_dvs_port(path)
+            except Exception:
+                pass
+            time.sleep(0.5)
+        return t, self._dvs_port

{ansys_pyensight_core-0.8.0 → ansys_pyensight_core-0.8.1}/src/ansys/pyensight/core/utils/variables.py

@@ -486,24 +486,18 @@ class Variables:
         #
         #
         temp_string = "ENS_Force_Net_Tan_ShearForce_X = StatMoment(plist,ENS_Force_Tan_ShearForce_X, 0, Compute_Per_part)"
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
 
         temp_string = "ENS_Force_Net_Tan_ShearForce_Y = StatMoment(plist,ENS_Force_Tan_ShearForce_Y, 0, Compute_Per_part)"
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
 
         temp_string = "ENS_Force_Net_Tan_ShearForce_Z = StatMoment(plist,ENS_Force_Tan_ShearForce_Z, 0, Compute_Per_part)"
         if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
         #
         # get the 3 constant force values XYZ
@@ -520,7 +514,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fx.append(v)
             else:
@@ -531,7 +525,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fy.append(v)
             else:
@@ -542,7 +536,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fz.append(v)
             else:
@@ -553,69 +547,51 @@ class Variables:
         #   they are calc'd to give the user the totals.
         #
         temp_string = "ENS_Force_Total_Net_Tan_ShearForce_X = StatMoment(plist,ENS_Force_Tan_ShearForce_X, 0, Compute_Per_case)"
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
 
         temp_string = "ENS_Force_Total_Net_Tan_ShearForce_Y = StatMoment(plist,ENS_Force_Tan_ShearForce_Y, 0, Compute_Per_case)"
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
 
         temp_string = "ENS_Force_Total_Net_Tan_ShearForce_Z = StatMoment(plist,ENS_Force_Tan_ShearForce_Z, 0, Compute_Per_case)"
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            err_string = "Error, failed to calculate a variable in {}".format(
-                fcn_name
-            )  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
             raise RuntimeError(err_string)  # pragma: no cover
 
         if frame_index > 0 and frame_index < len(self.ensight.objs.core.FRAMES):
             temp_string = "ENS_Force_Net_Tan_ShearForce_R = StatMoment(plist,ENS_Force_Tan_ShearForce_R,0, Compute_Per_part)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
 
             temp_string = "ENS_Force_Net_Tan_ShearForce_T = StatMoment(plist,ENS_Force_Tan_ShearForce_T,0, Compute_Per_part)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
 
             temp_string = "ENS_Force_Net_Tan_ShearForce_A = StatMoment(plist,ENS_Force_Tan_ShearForce_A,0, Compute_Per_part)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
             #
             # Totals
             #
             temp_string = "ENS_Force_Total_Net_Tan_ShearForce_R = StatMoment(plist,ENS_Force_Tan_ShearForce_R,0, Compute_Per_case)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
 
             temp_string = "ENS_Force_Total_Net_Tan_ShearForce_T = StatMoment(plist,ENS_Force_Tan_ShearForce_T,0, Compute_Per_case)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
 
             temp_string = "ENS_Force_Total_Net_Tan_ShearForce_A = StatMoment(plist,ENS_Force_Tan_ShearForce_A,0, Compute_Per_case)"
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                err_string = "Error, failed to calculate a variable in {}".format(
-                    fcn_name
-                )  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
                 raise RuntimeError(err_string)  # pragma: no cover
             #
             # get the 3 constant force values Radial, Theta, Axial
@@ -633,7 +609,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Fr.append(v)
                 else:
@@ -644,7 +620,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Ft.append(v)
                 else:
@@ -655,7 +631,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Fa.append(v)
                 else:
@@ -1027,6 +1003,7 @@ class Variables:
         #
         # Select the part(s) in the list
         # ensight.variables.evaluate("ENS_Force_Net_press_Y = StatMoment(plist,pressure,0,Compute_Per_part)")
+        fcn_name = "sum_press_forces_xyz_rtz"
         self.ensight.utils.parts.select_parts(pobj_list)
         #
         # Calculate the net force X, Y, and Z , 10.2.0(d) now per part constant variable
@@ -1034,20 +1011,23 @@ class Variables:
         force_calc_string = "ENS_Force_Net_press_X = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_X , 0, Compute_Per_part )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         force_calc_string = "ENS_Force_Net_press_Y = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_Y , 0, Compute_Per_part )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         force_calc_string = "ENS_Force_Net_press_Z = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_Z , 0, Compute_Per_part )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         # Calculate the Total force X, Y, and Z , 10.2.0(d) now case constant variable
         #  Totals are a case constants. We don't do anything with these vars
@@ -1056,20 +1036,23 @@ class Variables:
         force_calc_string = "ENS_Force_Total_Net_press_X = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_X , 0, Compute_Per_case )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         force_calc_string = "ENS_Force_Total_Net_press_Y = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_Y , 0, Compute_Per_case )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         force_calc_string = "ENS_Force_Total_Net_press_Z = StatMoment"
         calc_string = "(plist," + "ENS_Force_press_Z , 0, Compute_Per_case )"
         temp_string = force_calc_string + calc_string
-        if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-            return None  # pragma: no cover
+        if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+            err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+            raise RuntimeError(err_string)  # pragma: no cover
         #
         # get a list with a per part force, one for each part, new 10.1.6(b)
         #
@@ -1083,7 +1066,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fx.append(v)
             else:
@@ -1094,7 +1077,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fy.append(v)
             else:
@@ -1105,7 +1088,7 @@ class Variables:
         if val:
             if isinstance(val, list):
                 for v in val:
-                    if not v:  # pragma: no cover
+                    if v is None:  # pragma: no cover
                         return None  # pragma: no cover
                     Fz.append(v)
             else:
@@ -1125,20 +1108,23 @@ class Variables:
             force_calc_string = "ENS_Force_Net_press_R = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_R, 0, Compute_Per_part )"
             temp_string = force_calc_string + calc_string
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             force_calc_string = "ENS_Force_Net_press_T = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_T, 0, Compute_Per_part )"
             temp_string = force_calc_string + calc_string
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             force_calc_string = "ENS_Force_Net_press_A = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_A, 0, Compute_Per_part )"
             temp_string = force_calc_string + calc_string
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             #  Totals are a case constants. We don't do anything with these vars
             #   they are calc'd to give the user the totals.
@@ -1146,20 +1132,23 @@ class Variables:
             force_calc_string = "ENS_Force_Total_Net_press_R = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_R, 0, Compute_Per_case )"
             temp_string = force_calc_string + calc_string
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             force_calc_string = "ENS_Force_Total_Net_press_T = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_T, 0, Compute_Per_case )"
             temp_string = force_calc_string + calc_string
-            if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+            if self._calc_var(pobj_list, temp_string) is False:  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             force_calc_string = "ENS_Force_Total_Net_press_A = StatMoment"
             calc_string = "(plist," + "ENS_Force_press_A, 0, Compute_Per_case )"
             temp_string = force_calc_string + calc_string
             if not self._calc_var(pobj_list, temp_string):  # pragma: no cover
-                return None  # pragma: no cover
+                err_string = f"Error, failed to calculate {temp_string} in {fcn_name}"  # pragma: no cover  # pragma: no cover  # pragma: no cover
+                raise RuntimeError(err_string)  # pragma: no cover
             #
             #   get a list with a per part force, one for each part, new 10.1.6(b)
             #
@@ -1172,7 +1161,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Fr.append(v)
                 else:
@@ -1183,7 +1172,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Ft.append(v)
                 else:
@@ -1194,7 +1183,7 @@ class Variables:
             if val:
                 if isinstance(val, list):
                     for v in val:
-                        if not v:  # pragma: no cover
+                        if v is None:  # pragma: no cover
                             return None  # pragma: no cover
                         Fa.append(v)
                 else: