Unified approach to defining number of axial sublayers in PyARC

Currently, the round functions is used to calculate the number of axial layers in a subregion, which is followed by a max statement to ensure that a minimum number of 1 axial layer is defined per subregion. Instead, the math.ceil function should be used to define the number of axial layers. This behavior is found in PyARCUserObject.py, PyPROTEUSNodalObject.py, and PyGriffinConnect.py. To ensure consistency between all of these definitions, a single function should be defined that is called by all of these files.