tecio-python 0.2.0__py3-none-any.whl → 0.2.2__py3-none-any.whl

tecio/cli/tecdump.py

@@ -2,10 +2,10 @@ r"""Dump the full contents of a Tecplot data file to standard output.
 
 Inspecting the contents of a binary Tecplot file (variable names, zone metadata,
 auxiliary data, and raw nodal or element values) ordinarily requires either opening the
-file in Tecplot itself or writing a dedicated script against the TecIO API.  For a quick
+file in Tecplot itself or writing a dedicated script against the TecIO API. For a quick
 sanity check during a simulation or post-processing workflow, ``tecdump`` serializes
 every record in the file to plain text, providing a complete and human-readable
-representation of the file contents directly in the terminal.  The output is structured
+representation of the file contents directly in the terminal. The output is structured
 hierarchically: file header and auxiliary data first, followed by per-zone headers and,
 optionally, the underlying variable and connectivity arrays.
 
@@ -27,7 +27,7 @@ arrays can be suppressed entirely to focus on structural metadata.
 
 :Options:
     ``--ignore-zones``
-        Print only the file-level header and auxiliary data, then exit.  Zone records
+        Print only the file-level header and auxiliary data, then exit. Zone records
         and variable arrays are suppressed entirely. Useful for confirming file-level
         metadata such as solver, strand IDs, or global auxiliary entries without the
         overhead of parsing zone data.
@@ -47,7 +47,7 @@ arrays can be suppressed entirely to focus on structural metadata.
 
     ``-maxvals INT``
         Maximum number of values to print from any variable or connectivity array before
-        truncating with a summary line.  Defaults to a small number suitable for a quick
+        truncating with a summary line. Defaults to a small number suitable for a quick
         preview; set to a large value (e.g. ``1000000``) to print an array in full.
 
 :Returns:

tecio/cli/tecextract.py

@@ -1,12 +1,12 @@
 r"""Extract a subset of zones and/or variables from a Tecplot data file.
 
 Large Tecplot files produced by CFD solvers frequently contain many zones and variables,
-only a fraction of which are relevant to a given analysis.  Loading the full file into
+only a fraction of which are relevant to a given analysis. Loading the full file into
 Tecplot or writing a dedicated extraction script is impractical when the goal is simply
 to isolate a surface zone or a handful of flow variables for downstream processing.
 ``tecextract`` addresses this by writing a new Tecplot file containing only the
 requested subset, with the output format determined by the file extension of the output
-path.  This makes it straightforward to reduce file size, change format, or prepare a
+path. This makes it straightforward to reduce file size, change format, or prepare a
 subset for use with another tool in a single command.
 
 :Usage:
@@ -158,7 +158,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         metavar="LIST",
         help=(
             "Comma-separated list of 1-based zone indices to extract "
-            "(e.g. -zones 1,3,5).  Default is all zones."
+            "(e.g. -zones 1,3,5). Default is all zones."
         ),
     )
     parser.add_argument(
@@ -168,7 +168,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         metavar="LIST",
         help=(
             "Comma-separated list of 1-based variable indices to extract "
-            "(e.g. -variables 1,2,5).  Default is all variables."
+            "(e.g. -variables 1,2,5). Default is all variables."
         ),
     )
     parser.add_argument(
@@ -178,8 +178,8 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="PATH",
         help=(
-            "Output file path.  The extension controls the output format "
-            "(.szplt, .plt, .dat).  Defaults to <stem>_extract<ext> in "
+            "Output file path. The extension controls the output format "
+            "(.szplt, .plt, .dat). Defaults to <stem>_extract<ext> in "
             "the same directory as the input."
         ),
     )

tecio/cli/tecfix.py

@@ -156,7 +156,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="PATH",
         help=(
-            "Explicit output file path.  Defaults to <stem>_fixed<ext> "
+            "Explicit output file path. Defaults to <stem>_fixed<ext> "
             "in the same directory as the input file."
         ),
     )
@@ -196,7 +196,7 @@ def _classify_array(arr: npt.NDArray) -> str | None:
 
     """
     # Some readers (e.g. DAT) may return an empty or None array for passive
-    # or otherwise unavailable variables.  Treat these as clean -> there is
+    # or otherwise unavailable variables. Treat these as clean -> there is
     # nothing to inspect and numpy operations like isnan would error or
     # return misleading results on a zero-size array.
     if arr is None or arr.size == 0:
@@ -235,7 +235,7 @@ def _process_zone(
 
     For each variable:
     - If the variable is already passive or shared, it is forwarded exactly
-      as-is.  Sharing references are 1-based zone indices that are identical
+      as-is. Sharing references are 1-based zone indices that are identical
       in the output file because tecfix writes every zone in order without
       skipping any.
     - If the variable is a floating-point type and contains NaN or Inf it is
@@ -478,7 +478,7 @@ def main(argv: Sequence[str] | None = None) -> int:
                     } or None
 
                     # Sanitize solution_time: if it is NaN the zone metadata
-                    # is corrupt.  Reset both solution_time and strand_id to 0
+                    # is corrupt. Reset both solution_time and strand_id to 0
                     # so Tecplot treats the zone as stationary and the writer
                     # does not propagate the bad value downstream.
                     sol_time = zone.solution_time

tecio/cli/tecmerge.py

@@ -3,8 +3,8 @@ r"""Merge zones from multiple Tecplot data files into a single output file.
 Post-processing workflows commonly produce results distributed across multiple files.
 ``tecmerge`` collects all zones from an arbitrary number of input files into a single
 output, reconciling variable lists across sources by taking their union and writing any
-variable absent from a given source as passive.  Input files may be specified explicitly
-or via a quoted glob pattern and may be any mix of supported formats.  When merging
+variable absent from a given source as passive. Input files may be specified explicitly
+or via a quoted glob pattern and may be any mix of supported formats. When merging
 time-step sequences, solution times and strand IDs can be assigned automatically from a
 start time and either a fixed interval or an end time.
 
@@ -44,7 +44,7 @@ start time and either a fixed interval or an end time.
         Solution time of the first input file. Used with ``--assign-time-strands``.
 
     ``-delta VALUE``
-        Constant time increment between successive input files.  Mutually exclusive with
+        Constant time increment between successive input files. Mutually exclusive with
         ``-end``.
 
     ``-end VALUE``
@@ -136,8 +136,8 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         nargs="+",
         metavar="FILE",
         help=(
-            "Input Tecplot files.  Glob patterns are expanded (quote the "
-            "pattern to prevent shell expansion).  Files are merged in the "
+            "Input Tecplot files. Glob patterns are expanded (quote the "
+            "pattern to prevent shell expansion). Files are merged in the "
             "order given / matched."
         ),
     )
@@ -150,7 +150,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         required=True,
         metavar="PATH",
         help=(
-            "Output file path.  The extension controls the output format "
+            "Output file path. The extension controls the output format "
             "(.szplt, .plt, .dat)."
         ),
     )
@@ -166,7 +166,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="STRING",
         help=(
-            "Dataset title for the output file.  Defaults to the title of "
+            "Dataset title for the output file. Defaults to the title of "
             "the first input file."
         ),
     )
@@ -180,7 +180,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         dest="assign_ts",
         help=(
             "Assign evenly-spaced solution times and a strand ID to all "
-            "zones.  Each input file corresponds to one time step.  "
+            "zones. Each input file corresponds to one time step.  "
             "Requires -start and one of -delta or -end."
         ),
     )
@@ -192,7 +192,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         help="Solution time of the first input file.",
     )
     # -delta and -end are mutually exclusive -- either specifies the spacing,
-    # the other specifies the endpoint.  Providing both over-constrains the
+    # the other specifies the endpoint. Providing both over-constrains the
     # problem and is therefore disallowed at the parser level.
     step_group = ts.add_mutually_exclusive_group()
     step_group.add_argument(
@@ -208,7 +208,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="VALUE",
         help=(
-            "Solution time of the last input file.  The step is computed as "
+            "Solution time of the last input file. The step is computed as "
             "(end - start) / (N - 1) where N is the number of input files."
         ),
     )
@@ -217,7 +217,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         type=int,
         default=1,
         metavar="ID",
-        help="Strand ID to assign to all zones.  Default: 1.",
+        help="Strand ID to assign to all zones. Default: 1.",
     )
 
     return parser.parse_args(argv)
@@ -310,7 +310,7 @@ def _write_zone(
 ) -> None:
     """Write one zone to *writer* using the reconciled variable list.
 
-    Variables present in *zone* are copied verbatim.  Variables absent from
+    Variables present in *zone* are copied verbatim. Variables absent from
     *zone* (``None`` entries in *local_index_map*) are written as passive.
 
     Args:

tecio/cli/teconvert.py

@@ -1,7 +1,7 @@
 r"""Convert between Tecplot data file formats.
 
 The three Tecplot file formats, SZL binary (``.szplt``), PLT binary (``.plt``), and
-ASCII DAT (``.dat``), are not always interchangeable in practice.  Tecplot itself
+ASCII DAT (``.dat``), are not always interchangeable in practice. Tecplot itself
 stores all data internally in binary form and provides license-free command-line
 utilities (``preplot``, ``tec360 -convert``) for converting from lower-level formats
 upward to SZL, but no supported path exists for the reverse without a active license.
@@ -182,7 +182,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="PATH",
         help=(
-            "Explicit output file path.  Defaults to the input file stem with "
+            "Explicit output file path. Defaults to the input file stem with "
             "the new extension in the same directory as the input."
         ),
     )
@@ -198,7 +198,7 @@ def _copy_zones(reader: Any, writer: Any) -> None:
     """Stream all zones from *reader* into the open *writer*.
 
     Each zone is reproduced at its original data type and value location.
-    Connectivity (node maps) is copied verbatim for FE zones.  Zone-level
+    Connectivity (node maps) is copied verbatim for FE zones. Zone-level
     auxiliary data is forwarded as well.
 
     Args:

tecio/cli/tecscale.py

@@ -3,7 +3,7 @@ r"""Scale and/or offset variable arrays in a Tecplot data file.
 CFD solvers and experimental data sources frequently produce output in differing unit
 systems, and converting between them (e.g.pressure in Pa to kPa, temperature in Kelvin
 to Celsius or lengths in meters to feet) is a routine step before visualisation or
-comparison.  Performing these transformations would otherwise require loading the file
+comparison. Performing these transformations would otherwise require loading the file
 in Tecplot or writing a dedicated script. ``tecscale`` applies the linear transformation
 
 .. math::
@@ -11,7 +11,7 @@ in Tecplot or writing a dedicated script. ``tecscale`` applies the linear transf
     v' = v \cdot s + b
 
 to every value in a selected variable array, where :math:`s` is the scale factor and
-:math:`b` is the additive offset.  The transformation can be restricted to a single zone
+:math:`b` is the additive offset. The transformation can be restricted to a single zone
 and the output format is controlled by the output file extension, allowing unit
 conversion and format conversion to be performed in a single step.
 
@@ -144,14 +144,14 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         type=float,
         default=1.0,
         metavar="FLOAT",
-        help="Multiplicative scale factor.  Default: 1.0.",
+        help="Multiplicative scale factor. Default: 1.0.",
     )
     parser.add_argument(
         "-offset",
         type=float,
         default=0.0,
         metavar="FLOAT",
-        help="Additive offset applied after scaling.  Default: 0.0.",
+        help="Additive offset applied after scaling. Default: 0.0.",
     )
     parser.add_argument(
         "-zone",
@@ -159,7 +159,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="INDEX",
         help=(
-            "1-based zone index to apply the transformation to.  Default is all zones."
+            "1-based zone index to apply the transformation to. Default is all zones."
         ),
     )
     parser.add_argument(
@@ -169,7 +169,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=None,
         metavar="PATH",
         help=(
-            "Output file path.  The extension controls the output format.  "
+            "Output file path. The extension controls the output format.  "
             "Defaults to <stem>_scaled<ext> in the same directory."
         ),
     )
@@ -225,7 +225,7 @@ def _resolve_variable(spec: str, var_names: list[str]) -> int | None:
             return i
 
     print(
-        f"Error: variable '{spec}' not found.  Available: {var_names}",
+        f"Error: variable '{spec}' not found. Available: {var_names}",
         file=sys.stderr,
     )
     return None

tecio/cli/tecslice.py

@@ -1,7 +1,7 @@
 r"""Slice a Tecplot data file along IJK indices and/or solution time.
 
 Structured CFD datasets frequently contain more spatial resolution or time steps than
-are necessary for a given analysis.  Reducing the data to a relevant subset (like every
+are necessary for a given analysis. Reducing the data to a relevant subset (like every
 other grid point, a specific index range along one axis, or a window of solution times)
 would otherwise require loading the full file in Tecplot or writing a dedicated script.
 ``tecslice`` applies these reductions directly to the file using a compact
@@ -212,7 +212,7 @@ def _parse_ijk_slice(value: str) -> IjkSliceSpec:
         )
     except ValueError as exc:
         raise argparse.ArgumentTypeError(
-            f"Invalid IJK slice {value!r}: {exc}.  All components must be integers."
+            f"Invalid IJK slice {value!r}: {exc}. All components must be integers."
         ) from exc
 
 
@@ -313,7 +313,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         required=True,
         metavar="PATH",
         help=(
-            "Output file path.  The extension controls the output format "
+            "Output file path. The extension controls the output format "
             "(.szplt, .plt, .dat)."
         ),
     )
@@ -408,10 +408,10 @@ def _ijk_spec_to_slice(spec: IjkSliceSpec | None) -> slice:
     # 1-based inclusive end -> Python exclusive stop.
     #
     # Positive step: 1-based end N -> py_stop N.
-    #   arr[:N] (0-based) gives elements 0..N-1 = 1-based 1..N.  Correct.
+    #   arr[:N] (0-based) gives elements 0..N-1 = 1-based 1..N. Correct.
     #
     # Negative step: 1-based end N -> py_stop N-1.
-    #   arr[s:N-1:-1] includes 0-based index N-1 = 1-based N.  Correct.
+    #   arr[s:N-1:-1] includes 0-based index N-1 = 1-based N. Correct.
     #   Special case: 1-based end=1 -> 0-based index 0 -> py_stop=None
     #   (arr[s:0:-1] would miss index 0; arr[s:None:-1] includes it).
     if spec.end is not None:
@@ -445,7 +445,7 @@ def _build_time_filter(
 ) -> set[int]:
     """Return the set of 0-based zone indices to include in the output.
 
-    Zones with strand ID == 0 are always included.  For each eligible
+    Zones with strand ID == 0 are always included. For each eligible
     strand the zones are sorted by solution time, the ``[start, end]``
     window is applied (inclusive), then the ``skip`` stride is applied to
     the windowed list.
@@ -514,7 +514,7 @@ def _build_protected_set(
 
     A zone is "protected" if it is excluded by the time filter but is
     referenced as a variable-sharing source by at least one zone that will
-    be written.  Protected zones are written with all non-shared variables
+    be written. Protected zones are written with all non-shared variables
     set passive and their strand ID forced to 0 so they act as inert grid
     anchors without contributing to the time series.
 
@@ -554,7 +554,7 @@ def _collect_zone_arrays(
 
     Returns:
         ``(data, locs, passive_vars, var_sharing)`` where each list has
-        length ``num_vars``.  Passive and shared entries have empty arrays.
+        length ``num_vars``. Passive and shared entries have empty arrays.
 
     """
     data: list[np.ndarray] = []
@@ -653,9 +653,9 @@ def _write_zone_protected(
     """Write a zone as a protected grid anchor.
 
     All variables that are not already passive or shared are set passive so
-    that no solution data is carried into the output.  The strand ID is
+    that no solution data is carried into the output. The strand ID is
     forced to 0 so the zone is treated as stationary and does not appear in
-    any time-series animation.  Sharing references are preserved verbatim so
+    any time-series animation. Sharing references are preserved verbatim so
     that zones which share coordinates from this zone continue to work.
 
     Args:
@@ -754,7 +754,7 @@ def _slice_and_write_ordered(
 
         if var.value_location == ValueLocation.CELL_CENTERED:
             # Cell array has shape (I-1, J-1, K-1).
-            # Cell c sits between nodes c and c+1.  For a nodal selection
+            # Cell c sits between nodes c and c+1. For a nodal selection
             # with stride s, the fully-bounded cells are exactly the
             # lower-node indices of each consecutive selected-node pair:
             #   positive step -> selected_nodes[:-1]  (each is the lower node)
@@ -837,7 +837,7 @@ def main(argv: Sequence[str] | None = None) -> int:
 
     if not do_ijk and not do_time:
         print(
-            "Warning: no slice flags provided.  Output will be a verbatim copy.",
+            "Warning: no slice flags provided. Output will be a verbatim copy.",
             file=sys.stderr,
         )
 

tecio/cli/tecstats.py

@@ -26,7 +26,7 @@ corrective operation.
         are reported.
 
     ``-variable INDEX``
-        Restrict output to the variable at the given one-based index.  If omitted, all
+        Restrict output to the variable at the given one-based index. If omitted, all
         variables are reported.
 
     ``-csv``
@@ -154,14 +154,14 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         type=int,
         default=None,
         metavar="INDEX",
-        help="1-based zone index to report.  Default is all zones.",
+        help="1-based zone index to report. Default is all zones.",
     )
     parser.add_argument(
         "-variable",
         type=int,
         default=None,
         metavar="INDEX",
-        help="1-based variable index to report.  Default is all variables.",
+        help="1-based variable index to report. Default is all variables.",
     )
     parser.add_argument(
         "-csv",
@@ -169,7 +169,7 @@ def _parse_args(argv: Sequence[str] | None = None) -> argparse.Namespace:
         default=False,
         dest="write_csv",
         help=(
-            "Write statistics to a CSV file.  The filename is derived "
+            "Write statistics to a CSV file. The filename is derived "
             "automatically from the input file stem with a _stats suffix "
             "always appended, preceded by optional _zone_N and _var_N "
             "segments when -zone or -variable are active "

tecio/dat/__init__.py

@@ -6,7 +6,7 @@ files, but are useful for small datasets, debugging, and interoperability
 with other tools.
 
 This module provides :class:`Read` for parsing and :class:`Write` for
-producing Tecplot 360 ASCII data files (``.dat`` / ``.tec``).  Both
+producing Tecplot 360 ASCII data files (``.dat`` / ``.tec``). Both
 classes mirror the interfaces of :class:`szl.Read` / :class:`szl.Write`
 and :class:`plt.Read` / :class:`plt.Write` so that downstream code can
 switch between file formats by changing only the file extension passed to

tecio/dat/_read.py

@@ -3,7 +3,7 @@ r"""Tecplot ASCII DAT file reader API.
 Supported ``DATAPACKING`` modes:
     * ``BLOCK`` — one contiguous value block per variable (Tecplot default).
     * ``POINT`` — one row of all variable values per node, followed by a separate
-      row-per-cell section for any cell-centred variables.  This is the format most
+      row-per-cell section for any cell-centred variables. This is the format most
       commonly produced by third-party exporters and tools that treat the file like a
       CSV with a header.
 """
@@ -213,9 +213,9 @@ def _next_token_is_value(text: str, i: int) -> bool:
     """Return ``True`` if the token at *i* is a value, not the next key.
 
     Supports the legacy whitespace-separated ``KEY VALUE`` form in zone headers
-    (e.g. ``STRANDID 2``).  A token is treated as the *next key* — and therefore not a
+    (e.g. ``STRANDID 2``). A token is treated as the *next key* — and therefore not a
     value for the preceding key — when it is immediately followed (after optional
-    whitespace) by ``=``.  A quoted token is always a value.
+    whitespace) by ``=``. A quoted token is always a value.
 
     Example:
         >>> _next_token_is_value("STRANDID 2", 9)
@@ -270,7 +270,7 @@ def _kv_split(text: str) -> dict[str, str]:
             i += 1
             continue
 
-        # Determine the separator between key and value.  The modern form is
+        # Determine the separator between key and value. The modern form is
         # ``KEY = VALUE``; some legacy files (and ``preplot``) also accept a bare
         # whitespace separator, e.g. ``STRANDID 2``.
         while i < n and text[i] in " \t":
@@ -362,7 +362,7 @@ def _parse_legacy_element_type(text: str) -> ZoneType:
     """Parse a legacy ``ET=`` (element type) value into a :class:`ZoneType`.
 
     Recognises the canonical Tecplot element names (``TRIANGLE``,
-    ``QUADRILATERAL``, ``TETRAHEDRON``, ``BRICK``, ``LINESEG``).  A prefix-based
+    ``QUADRILATERAL``, ``TETRAHEDRON``, ``BRICK``, ``LINESEG``). A prefix-based
     fallback provides a little tolerance for the minor spelling variants some
     third-party exporters emit, keeping the reader as permissive as ``preplot``.
 
@@ -473,7 +473,7 @@ def _parse_varsharelist(text: str) -> dict[int, int]:
         "([1-3,5]=2)"     → {0: 2, 1: 2, 2: 2, 4: 2}
 
     Note that the bracketed group may contain either a single integer or a
-    hyphen-separated range.  Mixed forms such as ``[1-3,5]`` are also produced by some
+    hyphen-separated range. Mixed forms such as ``[1-3,5]`` are also produced by some
     writers.
 
     Args:
@@ -751,7 +751,7 @@ class ReadVariable:
         """Return :class:`DataType` inferred from the stored NumPy dtype.
 
         A shared variable reports the source zone's actual dtype, since its data array
-        is resolved from that zone (see :attr:`shared_zone`).  Only a passive variable
+        is resolved from that zone (see :attr:`shared_zone`). Only a passive variable
         (which has no data array anywhere in the file) returns :attr:`DataType.FLOAT` as
         a placeholder.
         """
@@ -846,7 +846,7 @@ class ReadZone:
     Attributes:
         datapacking: :class:`~tecio.libtecio.DataPacking` member reflecting the
             ``DATAPACKING`` keyword found in the zone header (``BLOCK`` or
-            ``POINT``).  The data arrays are identical either way; this attribute
+            ``POINT``). The data arrays are identical either way; this attribute
             records how the values were laid out on disk.
         shared_connectivity: 1-based source zone index if this zone's connectivity is
             shared via ``CONNECTIVITYSHAREZONE``, else ``None``.  :attr:`node_map`
@@ -921,7 +921,7 @@ class ReadZone:
     def nodes_per_cell(self) -> int:
         """Nodes per cell based on zone type.
 
-        Uses the module-level ``_NODES_PER_ELEM`` table for simple FE types.  For
+        Uses the module-level ``_NODES_PER_ELEM`` table for simple FE types. For
         ORDERED zones the count is inferred from the number of active dimensions (1-D →
         2, 2-D → 4, 3-D → 8).
 
@@ -954,7 +954,7 @@ class ReadZone:
     ) -> npt.NDArray | None | tuple[npt.NDArray | None, ...]:
         """Return variable data array(s) for this zone.
 
-        A single key (0-based index or exact name) returns one array.  A list of exact
+        A single key (0-based index or exact name) returns one array. A list of exact
         names returns a tuple of arrays in the order given, suitable for unpacking::
 
             p = zone.get_array("p")
@@ -986,7 +986,7 @@ class ReadZone:
 class Read:
     """Read a Tecplot ASCII DAT file into memory.
 
-    The entire file is parsed on construction.  All data is then available through the
+    The entire file is parsed on construction. All data is then available through the
     same attributes and methods as :class:`szl.Read`.
 
     Example:
@@ -1011,8 +1011,12 @@ class Read:
         self._zones: list[ReadZone] = []
         self._zone_list: ZoneList[ReadZone] | None = None
         self._auxdata: ReadAuxData = ReadAuxData()
-        # Index 0 is a None placeholder so that 1-based indexing works directly.
+        # Index 0 is a None placeholder so that 1-based indexing works directly
         self._var_auxdata: list[ReadAuxData | None] = [None]
+        # Raw VARAUXDATA lines seen before the first zone, buffered by
+        # _parse_file_header() and applied once _var_auxdata is allocated with the
+        # correct length (num_vars isn't known until the header finishes parsing)
+        self._deferred_var_aux_lines: list[str] = []
         self._parse()
 
     def __repr__(self) -> str:
@@ -1045,7 +1049,7 @@ class Read:
         """Context manager protocol for read-only file interfaces.
 
         Read classes hold no open file handle between accesses, so there is nothing to
-        release on exit.  Provided for API consistency with the Write classes and to
+        release on exit. Provided for API consistency with the Write classes and to
         support the ``with tecio.open(...) as r:`` pattern.
         """
         pass
@@ -1148,6 +1152,12 @@ class Read:
         # Build per-variable aux data slots now that num_vars is known.
         self._var_auxdata = [None] + [ReadAuxData() for _ in range(self.num_vars)]
 
+        # Now that _var_auxdata exists, apply any VARAUXDATA lines that appeared before
+        # the first zone
+        for raw in self._deferred_var_aux_lines:
+            _apply_varauxdata(raw, self._var_auxdata)
+        self._deferred_var_aux_lines.clear()
+
         while tokens.has_more():
             line = tokens.peek_stripped()
             upper = line.upper()
@@ -1229,9 +1239,11 @@ class Read:
                 if name:
                     self._auxdata._data[name] = value
 
-            # VARAUXDATA lines before the first ZONE are deferred — they need the
-            # variable list, which may not yet be complete.  They are processed in the
-            # main _parse() loop after the header finishes.
+            elif upper_key.startswith("VARAUXDATA"):
+                # Can't process this without num_vars, and therefore _var_auxdata, isn't
+                # known until this header finishes parsing. Buffer the raw line;
+                # _parse() applies these once _var_auxdata is allocated.
+                self._deferred_var_aux_lines.append(line)
 
     def _parse_zone(self, tokens: _LineBuffer) -> None:
         """Parse one ZONE block (header + data blocks + connectivity).
@@ -1288,7 +1300,7 @@ class Read:
         # (legacy).  ``ET`` (element type) only ever appears on finite-element zones —
         # ordered/structured data has no elements — so it merely names the element
         # *shape* once a zone is already known to be FE, and is ignored on ordered
-        # zones.  Modern keywords win when present.
+        # zones. Modern keywords win when present.
         legacy_is_fe: bool | None = None
         legacy_packing: DataPacking | None = None
         if "F" in kv:
@@ -1298,11 +1310,11 @@ class Read:
             zt_raw = kv["ZONETYPE"].rstrip(",").strip().lower()
             zone_type = _STR_TO_ZONETYPE.get(zt_raw, ZoneType.ORDERED)
         elif legacy_is_fe is False:
-            # F=POINT/BLOCK: ordered zone.  A stray ET (if any) does not apply.
+            # F=POINT/BLOCK: ordered zone. A stray ET (if any) does not apply.
             zone_type = ZoneType.ORDERED
         elif legacy_is_fe or "ET" in kv:
             # Finite-element zone: F=FEPOINT/FEBLOCK, or an ET keyword with no F
-            # (some exporters omit F).  The element shape comes from ET, which is
+            # (some exporters omit F). The element shape comes from ET, which is
             # then required.
             if "ET" not in kv:
                 raise ValueError(
@@ -1572,9 +1584,9 @@ class Read:
 
         The spec writes two interleaved sections:
 
-            * **Nodal section** — ``num_nodes`` rows, one per node.  Each row contains
+            * **Nodal section** — ``num_nodes`` rows, one per node. Each row contains
               the values of all active nodal variables in dataset order.
-            * **Cell-centered section** — ``num_cells`` rows, one per element.  Each row
+            * **Cell-centered section** — ``num_cells`` rows, one per element. Each row
               contains the values of all active CC variables in dataset order.
 
         When all variables are nodal the cell-centered section is empty and is skipped