harden: reject non-physical (degree-valued) cone and uptilt angles

WindIO stores angles in radians (SI): the IEA references carry cone_angle
0.0698 (4 deg), uptilt 0.1047 (6 deg) and root twist 0.2094 (12 deg), and
the auto-RNA reads them directly, which the IEA-22 deck cross-check
confirms to <0.5 % on CM. A hand-authored file that leaves an angle in
degrees (cone_angle: 4 meaning 4 deg) would otherwise be evaluated at
4 rad = 229 deg and silently corrupt the coned lever and axial offset.
Guard cone_angle and uptilt against a physical +/-45 deg bound so that
mistake fails loudly with a message naming the radians convention.

Author: SMI-Lab-Inha

src/pybmodes/io/windio.py

@@ -732,6 +732,27 @@ def _finite_float(value: Any, what: str) -> float:
     return f
 
 
+def _finite_angle(value: Any, what: str, max_abs: float = float(np.pi / 4)) -> float:
+    """Coerce ``value`` to a finite angle in radians within a physical bound.
+
+    WindIO stores angles in radians (SI). A rotor precone or shaft tilt beyond
+    ~45 deg is non-physical for a horizontal-axis turbine and almost always
+    means degrees were supplied by mistake (e.g. ``cone_angle: 4`` read as
+    4 rad = 229 deg instead of 4 deg), which would silently corrupt the coned
+    lever and axial offset. Reject it with a clear message rather than
+    evaluate the trig at the wrong angle.
+    """
+    ang = _finite_float(value, what)
+    if abs(ang) > max_abs:
+        raise ValueError(
+            f"{what} = {ang:.4f} rad ({np.degrees(ang):.1f} deg) is outside the "
+            f"physical range +/-{np.degrees(max_abs):.0f} deg. WindIO angles are "
+            f"in radians; a value near a small integer usually means degrees "
+            f"were supplied."
+        )
+    return ang
+
+
 def _positive_mass(value: Any, what: str) -> float:
     """Coerce ``value`` to a positive, finite mass (kg), else raise."""
     m = _finite_float(value, what)
@@ -1055,7 +1076,7 @@ def _nac_geom(key: str) -> float:
         return _finite_float(val, f"nacelle.drivetrain.{key}")
 
     overhang = _nac_geom("overhang")
-    uptilt = _nac_geom("uptilt")
+    uptilt = _finite_angle(_nac_geom("uptilt"), "nacelle.drivetrain.uptilt")
     dist_tt_hub = _nac_geom("distance_tt_hub")
 
     # --- Hub: components.<hub>.elastic_properties_mb ---
@@ -1082,7 +1103,7 @@ def _nac_geom(key: str) -> float:
     )
     if hub_r < 0.0:
         raise ValueError(f"hub diameter must be non-negative; got {hub_diam!r}.")
-    cone_ang = _finite_float(hub.get("cone_angle", 0.0), "hub cone_angle")
+    cone_ang = _finite_angle(hub.get("cone_angle", 0.0), "hub cone_angle")
     cone_cos = float(np.cos(cone_ang))
     cone_sin = float(np.sin(cone_ang))
 

tests/test_windio_rna.py

@@ -200,6 +200,30 @@ def _mk(cone: float):
     assert coned.mass == pytest.approx(flat.mass)
 
 
+def test_rna_rejects_degree_valued_cone_angle(tmp_path: pathlib.Path) -> None:
+    """WindIO angles are radians; a degrees-valued cone_angle (e.g. 4 read as
+    4 rad = 229 deg) is non-physical and rejected rather than silently
+    evaluated at the wrong angle (Codex review on #130)."""
+    pytest.importorskip("yaml")
+    from pybmodes.io.windio import read_windio_rna
+
+    o = _base_ontology()
+    o["components"]["hub"]["cone_angle"] = 4.0  # 4 deg mistakenly left un-converted
+    with pytest.raises(ValueError, match="radians"):
+        read_windio_rna(_write(o, tmp_path, "cone_deg.yaml"))
+
+
+def test_rna_rejects_degree_valued_uptilt(tmp_path: pathlib.Path) -> None:
+    """A degrees-valued uptilt is likewise rejected as out of physical range."""
+    pytest.importorskip("yaml")
+    from pybmodes.io.windio import read_windio_rna
+
+    o = _base_ontology()
+    o["components"]["nacelle"]["drivetrain"]["uptilt"] = 6.0  # 6 deg un-converted
+    with pytest.raises(ValueError, match="radians"):
+        read_windio_rna(_write(o, tmp_path, "uptilt_deg.yaml"))
+
+
 @pytest.mark.parametrize("n_bl", [1, 2])
 def test_rna_rejects_one_or_two_bladed_rotor(
     tmp_path: pathlib.Path, n_bl: int