api581-engine 0.1.4__py3-none-any.whl

api581/__init__.py

@@ -0,0 +1,80 @@
+"""API 581 calculation engine. No implicit standard-compliance certification."""
+
+from .core import VERSION as __version__, Context, DomainError, REGISTRY, evaluate
+from . import pof, release, effects, risk, thinning, cracking, level1, consequences, level1_summary, cost_tables, scc_mechanisms, external_corrosion, external_cracking, fatigue, htha, brittle, temper, bundle, prd, prd_consequences, prd_history, level1_type1, level1_toxic, pool_fire, inspection_effectiveness
+from .result import CalculationResult
+from . import fireball
+from . import blast
+from . import physical_explosion
+from . import toxic_level2
+from . import level2_summary
+from . import evaporation
+from . import source_partition
+from . import ignition
+from . import two_phase_events
+from . import tank
+from . import tank_environment
+from . import tank_financial
+from . import tank_source
+from . import tank_integrated
+from . import tank_shell
+from . import slabx_adapter
+from . import spatial
+from . import jet_handoff
+from . import toxic_criteria
+from . import native_toxic
+from . import soil_corrosion
+from . import tank_corrosion
+from . import tank_thinning
+from . import clad_thinning
+from . import thinning_conditions
+from . import lining
+from . import cooling_water
+from . import cooling_selection
+from . import co2_corrosion
+from . import co2_screening
+from . import oxidation
+from . import acid_sour_water
+from . import alkaline_sour_water
+from . import hydrochloric_acid
+from . import hydrofluoric_acid
+from . import amine_corrosion
+from . import sulfuric_acid
+from . import h2sh2_corrosion
+from . import sulfidic_corrosion
+from . import flash_solver
+from . import saturation
+from . import flash_pool
+from . import area_contribution
+from . import groundwater
+from . import thinning_selection
+from . import thinning_history
+from . import caustic_screening
+from . import amine_screening
+from . import ssc_screening
+from . import hic_screened
+from . import acscc_screening
+from . import pascc_screening
+from . import clscc_screening
+from .workflow import Step, Workflow, run_workflow
+from . import digital_twin
+from . import api520
+from . import schema
+from . import evidence
+from .evidence import ValidationEvidence
+from . import standard
+from . import portfolio
+from . import representative_fluids
+from . import corrosion_rate
+from . import thinning_type
+from .catalog import describe, list_calculations, coverage, coverage_summary, clause_evidence, review_queue, calculation_evidence, validation_queue, ScopeBoundary, scope_boundaries, audit_coverage
+
+__all__ = ["api", "Context", "DomainError", "CalculationResult", "evaluate",
+           "Step", "Workflow", "run_workflow", "describe", "list_calculations", "coverage", "coverage_summary", "clause_evidence", "review_queue", "calculation_evidence", "validation_queue", "ScopeBoundary", "scope_boundaries", "audit_coverage", "digital_twin", "schema", "evidence", "ValidationEvidence", "standard", "__version__"]
+
+
+def __getattr__(name):
+    if name == "api":
+        from importlib import import_module
+        return import_module(".api", __name__)
+    raise AttributeError(name)

api581/_amine_tables.py

@@ -0,0 +1,105 @@
+"""Checked 2016 amine source rows; generated by scripts/extract_amine_source_fixture.py."""
+TABLES = {'mea_dea_us': [[1.0, 3.0, 1.0, 3.0, 3.0, 10.0, 5.0, 15.0, 10.0, 25.0, 15.0, 40.0],
+                [2.0, 6.0, 2.0, 6.0, 6.0, 20.0, 15.0, 40.0, 20.0, 45.0, 30.0, 80.0],
+                [5.0, 10.0, 5.0, 15.0, 15.0, 40.0, 30.0, 60.0, 40.0, 90.0, 60.0, 120.0],
+                [1.0, 3.0, 2.0, 6.0, 5.0, 15.0, 10.0, 30.0, 15.0, 45.0, 20.0, 60.0],
+                [2.0, 6.0, 4.0, 12.0, 10.0, 30.0, 20.0, 60.0, 30.0, 90.0, 40.0, 80.0],
+                [5.0, 15.0, 8.0, 25.0, 20.0, 60.0, 40.0, 80.0, 60.0, 120.0, 120.0, 150.0],
+                [2.0, 6.0, 3.0, 9.0, 7.0, 20.0, 10.0, 30.0, 20.0, 60.0, 25.0, 75.0],
+                [4.0, 10.0, 6.0, 20.0, 15.0, 40.0, 20.0, 50.0, 40.0, 80.0, 50.0, 100.0],
+                [8.0, 25.0, 15.0, 45.0, 30.0, 60.0, 40.0, 80.0, 80.0, 120.0, 100.0, 150.0],
+                [2.0, 6.0, 4.0, 10.0, 7.0, 20.0, 15.0, 40.0, 25.0, 70.0, 30.0, 80.0],
+                [4.0, 10.0, 8.0, 25.0, 15.0, 45.0, 30.0, 60.0, 50.0, 100.0, 100.0, 150.0],
+                [8.0, 25.0, 15.0, 40.0, 35.0, 70.0, 60.0, 100.0, 100.0, 140.0, 150.0, 180.0],
+                [3.0, 9.0, 5.0, 15.0, 10.0, 30.0, 15.0, 45.0, 35.0, 70.0, 45.0, 100.0],
+                [6.0, 15.0, 10.0, 30.0, 20.0, 60.0, 45.0, 90.0, 70.0, 130.0, 90.0, 150.0],
+                [10.0, 30.0, 20.0, 40.0, 40.0, 80.0, 90.0, 120.0, 120.0, 150.0, 150.0, 180.0],
+                [3.0, 9.0, 7.0, 20.0, 10.0, 30.0, 25.0, 75.0, 40.0, 100.0, 50.0, 120.0],
+                [6.0, 20.0, 15.0, 45.0, 20.0, 60.0, 50.0, 100.0, 80.0, 140.0, 100.0, 150.0],
+                [10.0, 30.0, 30.0, 60.0, 45.0, 90.0, 100.0, 150.0, 140.0, 180.0, 160.0, 200.0],
+                [4.0, 10.0, 9.0, 30.0, 15.0, 40.0, 30.0, 100.0, 50.0, 120.0, 60.0, 150.0],
+                [8.0, 15.0, 20.0, 40.0, 30.0, 60.0, 60.0, 100.0, 90.0, 140.0, 100.0, 150.0],
+                [15.0, 35.0, 40.0, 80.0, 60.0, 100.0, 100.0, 150.0, 140.0, 180.0, 160.0, 200.0],
+                [5.0, 15.0, 10.0, 30.0, 20.0, 60.0, 40.0, 100.0, 60.0, 120.0, 70.0, 150.0],
+                [10.0, 30.0, 20.0, 60.0, 40.0, 80.0, 70.0, 120.0, 100.0, 150.0, 120.0, 150.0],
+                [20.0, 45.0, 40.0, 80.0, 60.0, 100.0, 100.0, 150.0, 150.0, 180.0, 170.0, 220.0]],
+ 'mea_dea_si': [[0.03, 0.08, 0.03, 0.08, 0.08, 0.25, 0.13, 0.38, 0.25, 0.64, 0.38, 1.02],
+                [0.05, 0.15, 0.05, 0.15, 0.15, 0.51, 0.38, 1.02, 0.51, 1.14, 0.76, 2.03],
+                [0.13, 0.25, 0.13, 0.38, 0.38, 1.02, 0.76, 1.52, 1.02, 2.29, 1.52, 3.05],
+                [0.03, 0.08, 0.05, 0.15, 0.13, 0.38, 0.25, 0.76, 0.38, 1.14, 0.51, 1.52],
+                [0.05, 0.15, 0.1, 0.3, 0.25, 0.76, 0.51, 1.52, 0.76, 2.29, 1.02, 2.03],
+                [0.13, 0.38, 0.2, 0.64, 0.51, 1.52, 1.02, 2.03, 1.52, 3.05, 3.05, 3.81],
+                [0.05, 0.15, 0.08, 0.23, 0.18, 0.51, 0.25, 0.76, 0.51, 1.52, 0.64, 1.91],
+                [0.1, 0.25, 0.15, 0.51, 0.38, 1.02, 0.51, 1.27, 1.02, 2.03, 1.27, 2.54],
+                [0.2, 0.64, 0.38, 1.14, 0.76, 1.52, 1.02, 2.03, 2.03, 3.05, 2.54, 3.81],
+                [0.05, 0.15, 0.1, 0.25, 0.18, 0.51, 0.38, 1.02, 0.64, 1.78, 0.76, 2.03],
+                [0.1, 0.25, 0.2, 0.64, 0.38, 1.14, 0.76, 1.52, 1.27, 2.54, 2.54, 3.81],
+                [0.2, 0.64, 0.38, 1.02, 0.89, 1.78, 1.52, 2.54, 2.54, 3.56, 3.81, 4.57],
+                [0.08, 0.23, 0.13, 0.38, 0.25, 0.76, 0.38, 1.14, 0.89, 1.78, 1.14, 2.54],
+                [0.15, 0.38, 0.25, 0.76, 0.51, 1.52, 1.14, 2.29, 1.78, 3.3, 2.29, 3.81],
+                [0.25, 0.76, 0.51, 1.02, 1.02, 2.03, 2.29, 3.05, 3.05, 3.81, 3.81, 4.57],
+                [0.08, 0.23, 0.18, 0.51, 0.25, 0.76, 0.64, 1.91, 1.02, 2.54, 1.27, 3.05],
+                [0.15, 0.51, 0.38, 1.14, 0.51, 1.52, 1.27, 2.54, 2.03, 3.56, 2.54, 3.81],
+                [0.25, 0.76, 0.76, 1.52, 1.14, 2.29, 2.54, 3.81, 3.56, 4.57, 4.06, 5.08],
+                [0.1, 0.25, 0.23, 0.76, 0.38, 1.02, 0.76, 2.54, 1.27, 3.05, 1.52, 3.81],
+                [0.2, 0.38, 0.51, 1.02, 0.76, 1.52, 1.52, 2.54, 2.29, 3.56, 2.54, 3.81],
+                [0.38, 0.89, 1.02, 2.03, 1.52, 2.54, 2.54, 3.81, 3.56, 4.57, 4.06, 5.08],
+                [0.13, 0.38, 0.25, 0.76, 0.51, 1.52, 1.02, 2.54, 1.52, 3.05, 1.78, 3.81],
+                [0.25, 0.76, 0.51, 1.52, 1.02, 2.03, 1.78, 3.05, 2.54, 3.81, 3.05, 3.81],
+                [0.51, 1.14, 1.02, 2.03, 1.52, 2.54, 2.54, 3.81, 3.81, 4.57, 4.32, 5.59]],
+ 'mdea_us': [[1.0, 3.0, 1.0, 3.0, 3.0, 10.0, 5.0, 15.0, 10.0, 25.0, 15.0, 40.0],
+             [2.0, 6.0, 2.0, 6.0, 6.0, 20.0, 15.0, 40.0, 20.0, 45.0, 30.0, 80.0],
+             [5.0, 10.0, 5.0, 15.0, 15.0, 40.0, 30.0, 60.0, 40.0, 90.0, 60.0, 120.0],
+             [1.0, 3.0, 2.0, 6.0, 5.0, 15.0, 10.0, 30.0, 15.0, 45.0, 20.0, 60.0],
+             [2.0, 6.0, 4.0, 12.0, 10.0, 30.0, 20.0, 60.0, 30.0, 90.0, 40.0, 80.0],
+             [5.0, 15.0, 8.0, 25.0, 20.0, 60.0, 40.0, 80.0, 60.0, 120.0, 120.0, 150.0],
+             [2.0, 6.0, 3.0, 9.0, 7.0, 20.0, 10.0, 30.0, 20.0, 60.0, 25.0, 75.0],
+             [4.0, 10.0, 6.0, 20.0, 15.0, 40.0, 20.0, 50.0, 40.0, 80.0, 50.0, 100.0],
+             [8.0, 25.0, 15.0, 45.0, 30.0, 60.0, 40.0, 80.0, 80.0, 120.0, 100.0, 150.0],
+             [2.0, 6.0, 4.0, 10.0, 7.0, 20.0, 15.0, 40.0, 25.0, 70.0, 30.0, 80.0],
+             [4.0, 10.0, 8.0, 25.0, 15.0, 45.0, 30.0, 60.0, 50.0, 100.0, 100.0, 150.0],
+             [8.0, 25.0, 15.0, 40.0, 35.0, 70.0, 60.0, 100.0, 100.0, 140.0, 150.0, 180.0],
+             [3.0, 9.0, 5.0, 15.0, 10.0, 30.0, 15.0, 45.0, 35.0, 70.0, 45.0, 100.0],
+             [6.0, 15.0, 10.0, 30.0, 20.0, 60.0, 45.0, 90.0, 70.0, 130.0, 90.0, 150.0],
+             [10.0, 30.0, 20.0, 40.0, 40.0, 80.0, 90.0, 120.0, 120.0, 150.0, 150.0, 180.0],
+             [3.0, 9.0, 7.0, 20.0, 10.0, 30.0, 25.0, 75.0, 40.0, 100.0, 50.0, 120.0],
+             [6.0, 20.0, 15.0, 45.0, 20.0, 60.0, 50.0, 100.0, 80.0, 140.0, 100.0, 150.0],
+             [10.0, 30.0, 30.0, 60.0, 45.0, 90.0, 100.0, 150.0, 140.0, 180.0, 160.0, 200.0],
+             [4.0, 10.0, 9.0, 30.0, 15.0, 40.0, 30.0, 100.0, 50.0, 120.0, 60.0, 150.0],
+             [8.0, 15.0, 20.0, 40.0, 30.0, 60.0, 60.0, 100.0, 90.0, 140.0, 100.0, 150.0],
+             [15.0, 35.0, 40.0, 80.0, 60.0, 100.0, 100.0, 150.0, 140.0, 180.0, 160.0, 200.0],
+             [5.0, 15.0, 10.0, 30.0, 20.0, 60.0, 40.0, 100.0, 60.0, 120.0, 70.0, 150.0],
+             [10.0, 30.0, 20.0, 60.0, 40.0, 80.0, 70.0, 120.0, 100.0, 150.0, 120.0, 150.0],
+             [20.0, 45.0, 40.0, 80.0, 60.0, 100.0, 100.0, 150.0, 150.0, 180.0, 170.0, 220.0]],
+ 'mdea_si': [[0.03, 0.08, 0.03, 0.08, 0.08, 0.25, 0.13, 0.38, 0.25, 0.64, 0.38, 1.02],
+             [0.05, 0.15, 0.05, 0.15, 0.15, 0.51, 0.38, 1.02, 0.51, 1.14, 0.76, 2.03],
+             [0.13, 0.25, 0.13, 0.38, 0.38, 1.02, 0.76, 1.52, 1.02, 2.29, 1.52, 3.05],
+             [0.03, 0.08, 0.05, 0.15, 0.13, 0.38, 0.25, 0.76, 0.38, 1.14, 0.51, 1.52],
+             [0.05, 0.15, 0.1, 0.3, 0.25, 0.76, 0.51, 1.52, 0.76, 2.29, 1.02, 2.03],
+             [0.13, 0.38, 0.2, 0.64, 0.51, 1.52, 1.02, 2.03, 1.52, 3.05, 3.05, 3.81],
+             [0.05, 0.15, 0.08, 0.23, 0.18, 0.51, 0.25, 0.76, 0.51, 1.52, 0.64, 1.91],
+             [0.1, 0.25, 0.15, 0.51, 0.38, 1.02, 0.51, 1.27, 1.02, 2.03, 1.27, 2.54],
+             [0.2, 0.64, 0.38, 1.14, 0.76, 1.52, 1.02, 2.03, 2.03, 3.05, 2.54, 3.81],
+             [0.05, 0.15, 0.1, 0.25, 0.18, 0.51, 0.38, 1.02, 0.64, 1.78, 0.76, 2.03],
+             [0.1, 0.25, 0.2, 0.64, 0.38, 1.14, 0.76, 1.52, 1.27, 2.54, 2.54, 3.81],
+             [0.2, 0.64, 0.38, 1.02, 0.89, 1.78, 1.52, 2.54, 2.54, 3.56, 3.81, 4.57],
+             [0.08, 0.23, 0.13, 0.38, 0.25, 0.76, 0.38, 1.14, 0.89, 1.78, 1.14, 2.54],
+             [0.15, 0.38, 0.25, 0.76, 0.51, 1.52, 1.14, 2.29, 1.78, 3.3, 2.29, 3.81],
+             [0.25, 0.76, 0.51, 1.02, 1.02, 2.03, 2.29, 3.05, 3.05, 3.81, 3.81, 4.57],
+             [0.08, 0.23, 0.18, 0.51, 0.25, 0.76, 0.64, 1.91, 1.02, 2.54, 1.27, 3.05],
+             [0.15, 0.51, 0.38, 1.14, 0.51, 1.52, 1.27, 2.54, 2.03, 3.56, 2.54, 3.81],
+             [0.25, 0.76, 0.76, 1.52, 1.14, 2.29, 2.54, 3.81, 3.56, 4.57, 4.06, 5.08],
+             [0.1, 0.25, 0.23, 0.76, 0.38, 1.02, 0.76, 2.54, 1.27, 3.05, 1.52, 3.81],
+             [0.2, 0.38, 0.51, 1.02, 0.76, 1.52, 1.52, 2.54, 2.29, 3.56, 2.54, 3.81],
+             [0.38, 0.89, 1.02, 2.03, 1.52, 2.54, 2.54, 3.81, 3.56, 4.57, 4.06, 5.08],
+             [0.13, 0.38, 0.25, 0.76, 0.51, 1.52, 1.02, 2.54, 1.52, 3.05, 1.78, 3.81],
+             [0.25, 0.76, 0.51, 1.52, 1.02, 2.03, 1.78, 3.05, 2.54, 3.81, 3.05, 3.81],
+             [0.51, 1.14, 1.02, 2.03, 1.52, 2.54, 2.54, 3.81, 3.81, 4.57, 4.32, 5.59]],
+ 'stainless': [[0.1, 0.03, 1.0],
+               [0.15, 0.03, 1.0],
+               [0.25, 0.03, 1.0],
+               [0.35, 0.05, 2.0],
+               [0.45, 0.05, 2.0],
+               [0.55, 0.08, 3.0],
+               [0.65, 0.1, 4.0],
+               [0.7, 0.13, 5.0]]}

api581/_h2sh2_tables.py

@@ -0,0 +1,161 @@
+"""2016 H2S/H2 source rates; generated by scripts/extract_h2sh2_source_fixture.py."""
+TABLES = {'carbon_low_alloy_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 8.0, 10.0, 14.0, 18.0],
+                                           [1.0, 1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0, 14.0, 20.0, 26.0, 34.0],
+                                           [1.0, 1.0, 1.0, 2.0, 4.0, 6.0, 8.0, 12.0, 16.0, 22.0, 29.0, 37.0],
+                                           [1.0, 2.0, 3.0, 4.0, 7.0, 11.0, 16.0, 22.0, 31.0, 41.0, 55.0, 71.0],
+                                           [1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 11.0, 15.0, 21.0, 29.0, 38.0, 50.0],
+                                           [1.0, 2.0, 4.0, 6.0, 9.0, 14.0, 21.0, 29.0, 41.0, 55.0, 73.0, 94.0],
+                                           [1.0, 2.0, 3.0, 5.0, 9.0, 13.0, 19.0, 27.0, 38.0, 51.0, 67.0, 87.0],
+                                           [2.0, 4.0, 6.0, 10.0, 16.0, 25.0, 36.0, 51.0, 71.0, 96.0, 130.0, 170.0],
+                                           [1.0, 2.0, 4.0, 7.0, 10.0, 16.0, 23.0, 33.0, 46.0, 62.0, 82.0, 110.0],
+                                           [2.0, 4.0, 8.0, 13.0, 20.0, 30.0, 44.0, 63.0, 87.0, 120.0, 160.0, 200.0],
+                                           [2.0, 3.0, 6.0, 10.0, 15.0, 23.0, 34.0, 48.0, 66.0, 90.0, 120.0, 150.0],
+                                           [3.0, 6.0, 11.0, 18.0, 29.0, 44.0, 64.0, 91.0, 130.0, 170.0, 230.0, 300.0],
+                                           [2.0, 4.0, 7.0, 11.0, 17.0, 26.0, 38.0, 54.0, 75.0, 100.0, 130.0, 170.0],
+                                           [4.0, 7.0, 12.0, 21.0, 32.0, 49.0, 72.0, 100.0, 140.0, 190.0, 250.0, 330.0],
+                                           [3.0, 5.0, 8.0, 13.0, 21.0, 32.0, 47.0, 67.0, 93.0, 130.0, 170.0, 220.0],
+                                           [5.0, 9.0, 15.0, 26.0, 40.0, 61.0, 89.0, 130.0, 180.0, 240.0, 310.0, 410.0]],
+ 'carbon_low_alloy_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.2, 0.25, 0.36, 0.46],
+                                 [0.03, 0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.25, 0.36, 0.51, 0.66, 0.86],
+                                 [0.03, 0.03, 0.03, 0.05, 0.1, 0.15, 0.2, 0.3, 0.41, 0.56, 0.74, 0.94],
+                                 [0.03, 0.05, 0.08, 0.1, 0.18, 0.28, 0.41, 0.56, 0.79, 1.04, 1.4, 1.8],
+                                 [0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.28, 0.38, 0.53, 0.74, 0.97, 1.27],
+                                 [0.03, 0.05, 0.1, 0.15, 0.23, 0.36, 0.53, 0.74, 1.04, 1.4, 1.85, 2.39],
+                                 [0.03, 0.05, 0.08, 0.13, 0.23, 0.33, 0.48, 0.69, 0.97, 1.3, 1.7, 2.21],
+                                 [0.05, 0.1, 0.15, 0.25, 0.41, 0.64, 0.91, 1.3, 1.8, 2.44, 3.3, 4.32],
+                                 [0.03, 0.05, 0.1, 0.18, 0.25, 0.41, 0.58, 0.84, 1.17, 1.57, 2.08, 2.79],
+                                 [0.05, 0.1, 0.2, 0.33, 0.51, 0.76, 1.12, 1.6, 2.21, 3.05, 4.06, 5.08],
+                                 [0.05, 0.08, 0.15, 0.25, 0.38, 0.58, 0.86, 1.22, 1.68, 2.29, 3.05, 3.81],
+                                 [0.08, 0.15, 0.28, 0.46, 0.74, 1.12, 1.63, 2.31, 3.3, 4.32, 5.84, 7.62],
+                                 [0.05, 0.1, 0.18, 0.28, 0.43, 0.66, 0.97, 1.37, 1.91, 2.54, 3.3, 4.32],
+                                 [0.1, 0.18, 0.3, 0.53, 0.81, 1.24, 1.83, 2.54, 3.56, 4.83, 6.35, 8.38],
+                                 [0.08, 0.13, 0.2, 0.33, 0.53, 0.81, 1.19, 1.7, 2.36, 3.3, 4.32, 5.59],
+                                 [0.13, 0.23, 0.38, 0.66, 1.02, 1.55, 2.26, 3.3, 4.57, 6.1, 7.87, 10.41]],
+ 'cr5_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 8.0, 11.0, 14.0],
+                              [1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 8.0, 12.0, 16.0, 21.0, 27.0],
+                              [1.0, 1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 9.0, 13.0, 18.0, 23.0, 30.0],
+                              [1.0, 1.0, 2.0, 4.0, 6.0, 9.0, 13.0, 18.0, 25.0, 33.0, 44.0, 57.0],
+                              [1.0, 1.0, 2.0, 2.0, 4.0, 6.0, 9.0, 12.0, 17.0, 23.0, 31.0, 40.0],
+                              [1.0, 2.0, 3.0, 5.0, 7.0, 11.0, 17.0, 24.0, 33.0, 44.0, 58.0, 76.0],
+                              [1.0, 2.0, 3.0, 4.0, 7.0, 10.0, 15.0, 22.0, 30.0, 41.0, 54.0, 70.0],
+                              [2.0, 3.0, 5.0, 8.0, 13.0, 20.0, 29.0, 41.0, 57.0, 77.0, 100.0, 130.0],
+                              [1.0, 2.0, 3.0, 5.0, 8.0, 13.0, 19.0, 27.0, 37.0, 50.0, 66.0, 85.0],
+                              [2.0, 4.0, 6.0, 10.0, 16.0, 24.0, 36.0, 51.0, 70.0, 94.0, 130.0, 160.0],
+                              [1.0, 3.0, 5.0, 8.0, 12.0, 19.0, 27.0, 39.0, 53.0, 72.0, 95.0, 120.0],
+                              [3.0, 5.0, 9.0, 15.0, 23.0, 35.0, 52.0, 73.0, 100.0, 140.0, 180.0, 240.0],
+                              [2.0, 3.0, 5.0, 9.0, 14.0, 21.0, 31.0, 44.0, 60.0, 81.0, 110.0, 140.0],
+                              [3.0, 6.0, 10.0, 17.0, 26.0, 40.0, 58.0, 82.0, 110.0, 150.0, 200.0, 270.0],
+                              [2.0, 4.0, 7.0, 11.0, 17.0, 26.0, 38.0, 54.0, 75.0, 100.0, 130.0, 170.0],
+                              [4.0, 7.0, 12.0, 21.0, 32.0, 49.0, 72.0, 100.0, 140.0, 190.0, 250.0, 330.0]],
+ 'cr5_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.2, 0.28, 0.36],
+                    [0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.2, 0.3, 0.41, 0.53, 0.69],
+                    [0.03, 0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.23, 0.33, 0.46, 0.58, 0.76],
+                    [0.03, 0.03, 0.05, 0.1, 0.15, 0.23, 0.33, 0.46, 0.64, 0.84, 1.12, 1.45],
+                    [0.03, 0.03, 0.05, 0.05, 0.1, 0.15, 0.23, 0.3, 0.43, 0.58, 0.79, 1.02],
+                    [0.03, 0.05, 0.08, 0.13, 0.18, 0.28, 0.43, 0.61, 0.84, 1.12, 1.47, 1.93],
+                    [0.03, 0.05, 0.08, 0.1, 0.18, 0.25, 0.38, 0.56, 0.76, 1.04, 1.37, 1.78],
+                    [0.05, 0.08, 0.13, 0.2, 0.33, 0.51, 0.74, 1.04, 1.45, 1.96, 2.54, 3.3],
+                    [0.03, 0.05, 0.08, 0.13, 0.2, 0.33, 0.48, 0.69, 0.94, 1.27, 1.68, 2.16],
+                    [0.05, 0.1, 0.15, 0.25, 0.41, 0.61, 0.91, 1.3, 1.78, 2.39, 3.3, 4.06],
+                    [0.03, 0.08, 0.13, 0.2, 0.3, 0.48, 0.69, 0.99, 1.35, 1.83, 2.41, 3.05],
+                    [0.08, 0.13, 0.23, 0.38, 0.58, 0.89, 1.32, 1.85, 2.54, 3.56, 4.57, 6.1],
+                    [0.05, 0.08, 0.13, 0.23, 0.36, 0.53, 0.79, 1.12, 1.52, 2.06, 2.79, 3.56],
+                    [0.08, 0.15, 0.25, 0.43, 0.66, 1.02, 1.47, 2.08, 2.79, 3.81, 5.08, 6.86],
+                    [0.05, 0.1, 0.18, 0.28, 0.43, 0.66, 0.97, 1.37, 1.91, 2.54, 3.3, 4.32],
+                    [0.1, 0.18, 0.3, 0.53, 0.81, 1.24, 1.83, 2.54, 3.56, 4.83, 6.35, 8.38]],
+ 'cr7_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 8.0, 10.0, 13.0],
+                              [1.0, 1.0, 1.0, 2.0, 2.0, 4.0, 5.0, 8.0, 11.0, 14.0, 19.0, 25.0],
+                              [1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 9.0, 12.0, 16.0, 21.0, 28.0],
+                              [1.0, 1.0, 2.0, 3.0, 5.0, 8.0, 11.0, 16.0, 23.0, 30.0, 40.0, 52.0],
+                              [1.0, 1.0, 1.0, 2.0, 4.0, 5.0, 8.0, 11.0, 16.0, 21.0, 28.0, 37.0],
+                              [1.0, 2.0, 3.0, 4.0, 7.0, 10.0, 15.0, 22.0, 30.0, 40.0, 53.0, 69.0],
+                              [1.0, 1.0, 2.0, 4.0, 6.0, 10.0, 14.0, 20.0, 28.0, 37.0, 49.0, 64.0],
+                              [1.0, 3.0, 5.0, 8.0, 12.0, 18.0, 27.0, 38.0, 52.0, 71.0, 94.0, 120.0],
+                              [1.0, 2.0, 3.0, 5.0, 8.0, 12.0, 17.0, 24.0, 34.0, 46.0, 60.0, 78.0],
+                              [2.0, 3.0, 6.0, 9.0, 15.0, 22.0, 33.0, 46.0, 64.0, 86.0, 110.0, 150.0],
+                              [1.0, 2.0, 4.0, 7.0, 11.0, 17.0, 25.0, 35.0, 49.0, 66.0, 87.0, 110.0],
+                              [3.0, 5.0, 8.0, 13.0, 21.0, 32.0, 47.0, 67.0, 93.0, 130.0, 170.0, 220.0],
+                              [2.0, 3.0, 5.0, 8.0, 13.0, 19.0, 28.0, 40.0, 55.0, 74.0, 98.0, 130.0],
+                              [3.0, 5.0, 9.0, 15.0, 24.0, 36.0, 53.0, 76.0, 100.0, 140.0, 190.0, 240.0],
+                              [2.0, 3.0, 6.0, 10.0, 16.0, 24.0, 35.0, 49.0, 68.0, 92.0, 120.0, 160.0],
+                              [4.0, 7.0, 11.0, 19.0, 30.0, 45.0, 66.0, 94.0, 130.0, 180.0, 230.0, 300.0]],
+ 'cr7_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.2, 0.25, 0.33],
+                    [0.03, 0.03, 0.03, 0.05, 0.05, 0.1, 0.13, 0.2, 0.28, 0.36, 0.48, 0.64],
+                    [0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.23, 0.3, 0.41, 0.53, 0.71],
+                    [0.03, 0.03, 0.05, 0.08, 0.13, 0.2, 0.28, 0.41, 0.58, 0.76, 1.02, 1.32],
+                    [0.03, 0.03, 0.03, 0.05, 0.1, 0.13, 0.2, 0.28, 0.41, 0.53, 0.71, 0.94],
+                    [0.03, 0.05, 0.08, 0.1, 0.18, 0.25, 0.38, 0.56, 0.76, 1.02, 1.35, 1.75],
+                    [0.03, 0.03, 0.05, 0.1, 0.15, 0.25, 0.36, 0.51, 0.71, 0.94, 1.24, 1.63],
+                    [0.03, 0.08, 0.13, 0.2, 0.3, 0.46, 0.69, 0.97, 1.32, 1.8, 2.39, 3.05],
+                    [0.03, 0.05, 0.08, 0.13, 0.2, 0.3, 0.43, 0.61, 0.86, 1.17, 1.52, 1.98],
+                    [0.05, 0.08, 0.15, 0.23, 0.38, 0.56, 0.84, 1.17, 1.63, 2.18, 2.79, 3.81],
+                    [0.03, 0.05, 0.1, 0.18, 0.28, 0.43, 0.64, 0.89, 1.24, 1.68, 2.21, 2.79],
+                    [0.08, 0.13, 0.2, 0.33, 0.53, 0.81, 1.19, 1.7, 2.36, 3.3, 4.32, 5.59],
+                    [0.05, 0.08, 0.13, 0.2, 0.33, 0.48, 0.71, 1.02, 1.4, 1.88, 2.49, 3.3],
+                    [0.08, 0.13, 0.23, 0.38, 0.61, 0.91, 1.35, 1.93, 2.54, 3.56, 4.83, 6.1],
+                    [0.05, 0.08, 0.15, 0.25, 0.41, 0.61, 0.89, 1.24, 1.73, 2.34, 3.05, 4.06],
+                    [0.1, 0.18, 0.28, 0.48, 0.76, 1.14, 1.68, 2.39, 3.3, 4.57, 5.84, 7.62]],
+ 'cr9_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 7.0, 9.0, 12.0],
+                              [1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0, 13.0, 17.0, 23.0],
+                              [1.0, 1.0, 1.0, 2.0, 2.0, 4.0, 6.0, 8.0, 11.0, 15.0, 19.0, 25.0],
+                              [1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 11.0, 15.0, 21.0, 28.0, 37.0, 48.0],
+                              [1.0, 1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0, 14.0, 20.0, 26.0, 34.0],
+                              [1.0, 1.0, 2.0, 4.0, 6.0, 10.0, 14.0, 20.0, 27.0, 37.0, 49.0, 64.0],
+                              [1.0, 1.0, 2.0, 4.0, 6.0, 9.0, 13.0, 18.0, 25.0, 34.0, 45.0, 59.0],
+                              [1.0, 2.0, 4.0, 7.0, 11.0, 17.0, 24.0, 35.0, 48.0, 65.0, 86.0, 110.0],
+                              [1.0, 2.0, 3.0, 4.0, 7.0, 11.0, 16.0, 22.0, 31.0, 42.0, 55.0, 72.0],
+                              [2.0, 3.0, 5.0, 9.0, 13.0, 20.0, 30.0, 42.0, 59.0, 79.0, 110.0, 140.0],
+                              [1.0, 2.0, 4.0, 7.0, 10.0, 16.0, 23.0, 32.0, 45.0, 61.0, 80.0, 100.0],
+                              [2.0, 4.0, 7.0, 12.0, 19.0, 30.0, 43.0, 61.0, 85.0, 120.0, 150.0, 200.0],
+                              [1.0, 3.0, 4.0, 7.0, 12.0, 18.0, 26.0, 37.0, 51.0, 68.0, 90.0, 120.0],
+                              [3.0, 5.0, 8.0, 14.0, 22.0, 33.0, 49.0, 69.0, 96.0, 130.0, 170.0, 220.0],
+                              [2.0, 3.0, 6.0, 9.0, 14.0, 22.0, 32.0, 45.0, 63.0, 85.0, 110.0, 150.0],
+                              [3.0, 6.0, 10.0, 17.0, 27.0, 41.0, 60.0, 86.0, 120.0, 160.0, 210.0, 280.0]],
+ 'cr9_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.13, 0.18, 0.23, 0.3],
+                    [0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.25, 0.33, 0.43, 0.58],
+                    [0.03, 0.03, 0.03, 0.05, 0.05, 0.1, 0.15, 0.2, 0.28, 0.38, 0.48, 0.64],
+                    [0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.28, 0.38, 0.53, 0.71, 0.94, 1.22],
+                    [0.03, 0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.25, 0.36, 0.51, 0.66, 0.86],
+                    [0.03, 0.03, 0.05, 0.1, 0.15, 0.25, 0.36, 0.51, 0.69, 0.94, 1.24, 1.63],
+                    [0.03, 0.03, 0.05, 0.1, 0.15, 0.23, 0.33, 0.46, 0.64, 0.86, 1.14, 1.5],
+                    [0.03, 0.05, 0.1, 0.18, 0.28, 0.43, 0.61, 0.89, 1.22, 1.65, 2.18, 2.79],
+                    [0.03, 0.05, 0.08, 0.1, 0.18, 0.28, 0.41, 0.56, 0.79, 1.07, 1.4, 1.83],
+                    [0.05, 0.08, 0.13, 0.23, 0.33, 0.51, 0.76, 1.07, 1.5, 2.01, 2.79, 3.56],
+                    [0.03, 0.05, 0.1, 0.18, 0.25, 0.41, 0.58, 0.81, 1.14, 1.55, 2.03, 2.54],
+                    [0.05, 0.1, 0.18, 0.3, 0.48, 0.76, 1.09, 1.55, 2.16, 3.05, 3.81, 5.08],
+                    [0.03, 0.08, 0.1, 0.18, 0.3, 0.46, 0.66, 0.94, 1.3, 1.73, 2.29, 3.05],
+                    [0.08, 0.13, 0.2, 0.36, 0.56, 0.84, 1.24, 1.75, 2.44, 3.3, 4.32, 5.59],
+                    [0.05, 0.08, 0.15, 0.23, 0.36, 0.56, 0.81, 1.14, 1.6, 2.16, 2.79, 3.81],
+                    [0.08, 0.15, 0.25, 0.43, 0.69, 1.04, 1.52, 2.18, 3.05, 4.06, 5.33, 7.11]],
+ 'cr12_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 9.0, 11.0, 14.0],
+                               [1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 8.0, 11.0, 14.0, 18.0],
+                               [1.0, 1.0, 1.0, 2.0, 2.0, 4.0, 5.0, 7.0, 9.0, 12.0, 15.0, 19.0],
+                               [1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 9.0, 12.0, 15.0, 19.0, 25.0],
+                               [1.0, 1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0, 13.0, 17.0, 22.0, 27.0],
+                               [1.0, 1.0, 2.0, 3.0, 4.0, 6.0, 9.0, 12.0, 16.0, 21.0, 27.0, 34.0],
+                               [1.0, 1.0, 2.0, 3.0, 5.0, 7.0, 10.0, 13.0, 18.0, 23.0, 30.0, 38.0],
+                               [1.0, 2.0, 3.0, 4.0, 7.0, 10.0, 13.0, 18.0, 25.0, 32.0, 42.0, 53.0]],
+ 'cr12_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.13, 0.15, 0.23, 0.28, 0.36],
+                     [0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.2, 0.28, 0.36, 0.46],
+                     [0.03, 0.03, 0.03, 0.05, 0.05, 0.1, 0.13, 0.18, 0.23, 0.3, 0.38, 0.48],
+                     [0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.23, 0.3, 0.38, 0.48, 0.64],
+                     [0.03, 0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.25, 0.33, 0.43, 0.56, 0.69],
+                     [0.03, 0.03, 0.05, 0.08, 0.1, 0.15, 0.23, 0.3, 0.41, 0.53, 0.69, 0.86],
+                     [0.03, 0.03, 0.05, 0.08, 0.13, 0.18, 0.25, 0.33, 0.46, 0.58, 0.76, 0.97],
+                     [0.03, 0.05, 0.08, 0.1, 0.18, 0.25, 0.33, 0.46, 0.64, 0.81, 1.07, 1.35]],
+ 'stainless_printed_us_converted': [[1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 2.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 2.0, 3.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 2.0, 3.0, 3.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 3.0, 4.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 3.0, 4.0, 5.0, 6.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0],
+                                    [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 2.0, 2.0, 4.0, 5.0, 7.0, 9.0]],
+ 'stainless_printed_si': [[0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.05],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.05, 0.08],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.05, 0.08, 0.08],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.08, 0.1],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.13],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.08, 0.1, 0.13, 0.15],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.08, 0.1, 0.13, 0.15],
+                          [0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.05, 0.05, 0.1, 0.13, 0.18, 0.23]]}

api581/_slabx_worker.py

@@ -0,0 +1,174 @@
+"""Internal JSON worker for process-isolated optional native engines."""
+import contextlib
+from decimal import Decimal
+import hashlib
+import importlib
+import importlib.metadata
+import io
+import json
+import math
+from pathlib import Path
+import sys
+from .core import DomainError, ensure, finite_tree
+from .slabx_adapter import NativePool, NativePoolFootprint, NativeJet, NativeJetFootprint
+
+
+def _fingerprint(package):
+    root = Path(importlib.import_module(package).__file__).parent
+    digest = hashlib.sha256()
+    for path in sorted(root.rglob("*.py")):
+        digest.update(path.relative_to(root).as_posix().encode())
+        digest.update(path.read_bytes())
+    return digest.hexdigest()
+
+
+def run_native(x):
+    required = {"slabx":"1.0.8"}
+    if x.engine == "slabx-lh2":
+        required["slabx-lh2"] = "0.2.0"
+    versions = {}
+    for name,expected in required.items():
+        try:
+            versions[name] = importlib.metadata.version(name)
+        except importlib.metadata.PackageNotFoundError as exc:
+            raise DomainError("engine_not_configured",f"Optional package {name} not installed") from exc
+        ensure(versions[name] == expected,f"Native mapping requires {name}=={expected}; installed {versions[name]}","unsupported_engine_version")
+    from slabx.core.plume import run_dispersion
+    from slabx.core.source import EvaporatingPool, HorizontalJet
+    from slabx.post.concentration import concentration_field
+    from slabx.submodels.atmosphere import Atmosphere
+    from slabx.thermo.base import Substance, LegacyThermo, water_backend
+    s,a = x.substance,x.atmosphere
+    # Canonical decimal SI-prefix conversion avoids a one-ULP difference from
+    # a native kg/mol decimal literal changing the adaptive integration grid.
+    molecular_weight = float(Decimal(str(s.molecular_weight_kg_kmol))/Decimal(1000))
+    substance = Substance(name=s.fluid_id,mw=molecular_weight,
+        cp_vapour=s.vapor_heat_capacity_j_kg_k,cp_liquid=s.liquid_heat_capacity_j_kg_k,
+        dh_vap=s.latent_heat_j_kg,T_boil=s.normal_boiling_point_k,rho_liquid=s.liquid_density_kg_m3)
+    atmosphere = Atmosphere(u_ref=a.wind_speed_m_s,z_ref=a.wind_reference_height_m,T=a.temperature_k,
+        rh=a.relative_humidity_fraction*100,z0=a.roughness_length_m,stability=a.stability,p=a.pressure_pa_abs)
+    jet = isinstance(x,NativeJet)
+    buoyant_vertical = False
+    native_exit_velocity = None
+    if jet:
+        from slabx.core.vertical_jet import VerticalJet
+        source_type = HorizontalJet if x.orientation == "horizontal_downwind" else VerticalJet
+        source = source_type(substance=substance,rate=x.mass_source_rate_kg_s,area=x.jet_exit_area_m2,
+            duration=x.release_duration_s,liquid_fraction=x.jet_liquid_mass_fraction,T_source=x.source_temperature_k,height=x.release_height_m)
+        native_exit_velocity = source.exit_velocity(atmosphere)
+        if x.orientation == "vertical_up" and source.rho_source(atmosphere) <= atmosphere.rho:
+            # SLAB User's Guide Eq. 46a-46d: use an effective horizontal
+            # source at the Briggs buoyant rise height and neutralize source
+            # temperature/cp for downstream dispersion.
+            g = 9.80665
+            rho_s = source.rho_source(atmosphere)
+            ws = source.w_source(atmosphere)
+            bs = source.half_width_geometric
+            ua = atmosphere.wind_speed(max(x.release_height_m, 1.0))
+            ustar = max(atmosphere.u_star, 1e-9)
+            fb = (4.0 / math.pi) * g * (atmosphere.rho - rho_s) * ws * bs * bs / atmosphere.rho
+            h_rise = 0.0
+            for _ in range(30):
+                next_h = 1.2 * (fb / (ua * ustar * ustar)) ** (3.0 / 5.0) * (x.release_height_m + h_rise) ** (2.0 / 5.0)
+                h_rise = next_h
+            neutral_temperature = (molecular_weight / atmosphere.mw_moist) * atmosphere.T
+            neutral_cp = (atmosphere.mw_moist / molecular_weight) * atmosphere.cp_moist
+            neutral = Substance(name=f"{s.fluid_id}-neutralized", mw=molecular_weight,
+                cp_vapour=neutral_cp,cp_liquid=neutral_cp,dh_vap=s.latent_heat_j_kg,
+                T_boil=s.normal_boiling_point_k,rho_liquid=s.liquid_density_kg_m3)
+            source = HorizontalJet(substance=neutral,rate=x.mass_source_rate_kg_s,area=x.jet_exit_area_m2,
+                duration=x.release_duration_s,liquid_fraction=0.0,T_source=neutral_temperature,
+                height=x.release_height_m + h_rise)
+            buoyant_vertical = True
+    else:
+        source = EvaporatingPool(substance=substance,rate=x.vapor_source_rate_kg_s,area=x.source_area_m2,duration=x.release_duration_s)
+    water = water_backend()
+    correction = contextlib.nullcontext()
+    if x.engine == "slabx-lh2":
+        from slabx_lh2.water_ice import with_sublimation
+        from slabx_lh2.plume_width import plume_width_coupling
+        water = with_sublimation(water,p_ambient=a.pressure_pa_abs)
+        correction = plume_width_coupling(active=x.lh2_plume_width_coupling)
+    with correction:
+        trajectory,used = run_dispersion(source,atmosphere,LegacyThermo(source.substance),water,
+            x_max=x.plume_distance_limit_m,t_max=x.puff_time_limit_s,n_puff_steps=x.puff_steps)
+        field = concentration_field(trajectory,atmosphere,z=x.receptor_height_m,t_avg=x.averaging_time_s,
+            t_release=x.release_duration_s,meander_closure=x.meander_closure)
+    distance = [float(v) for v in field.x]
+    concentration = [float(v) for v in field.peak]
+    ensure(len(distance) == len(concentration) and 1 < len(distance) <= 100000,"Invalid native grid","engine_result_invalid")
+    ensure(all(math.isfinite(v) for v in distance) and all(math.isfinite(v) and 0 <= v <= 1 for v in concentration),
+        "Nonphysical native concentration/grid","engine_result_invalid")
+    ensure(all(b >= a for a,b in zip(distance,distance[1:])),"Nonmonotonic native grid","engine_result_invalid")
+    endpoint = float(field.distance_to(x.threshold_volume_fraction))
+    reached = any(v >= x.threshold_volume_fraction for v in concentration)
+    endpoint = endpoint if reached and math.isfinite(endpoint) else None
+    censored = concentration[-1] >= x.threshold_volume_fraction
+    premise,ratio = "not_evaluated",None
+    if x.engine == "slabx-lh2":
+        from slabx_lh2.diagnostics import premise_summary
+        diagnostic = premise_summary(trajectory)
+        ratio = diagnostic["max"]
+        ensure(math.isfinite(ratio),"LH2 premise diagnostic not finite","engine_result_invalid")
+        premise = "violated" if diagnostic["violated"] else "within_premise"
+    status = "premise_violated" if premise == "violated" else "boundary_censored" if censored else "interpolated" if reached else "not_reached"
+    output = {"engine":x.engine,"engine_versions":versions,"engine_source_fingerprints":{
+        name:_fingerprint(name.replace("-","_")) for name in versions},
+        "downwind_distance_m":distance,"centreline_volume_fraction":concentration,
+        "native_grid_maximum_downwind_distance_m":max(distance),
+        "native_threshold_distance_m":endpoint,"threshold_distance_m":endpoint if status == "interpolated" else None,
+        "threshold_status":status,"premise_status":premise,"maximum_vertical_to_horizontal_velocity_ratio":ratio,
+        "used_source_area_m2":float(used.area),"effective_source_half_width_m":float(used.half_width),
+        "averaging_time_s":x.averaging_time_s,"receptor_height_m":x.receptor_height_m,
+        "spatial_output_scope":"centreline_only_no_footprint_area"}
+    finite_tree(output)
+    if jet:
+        output.pop("effective_source_half_width_m")
+        output.update(orientation=x.orientation,jet_exit_area_m2=float(x.jet_exit_area_m2),
+            geometric_exit_half_width_m=float(x.jet_exit_area_m2 ** 0.5 / 2.0),native_exit_velocity_m_s=float(native_exit_velocity if buoyant_vertical else used.exit_velocity(atmosphere)),
+            requested_source_temperature_k=x.source_temperature_k,used_source_temperature_k=float(used.T),
+            source_temperature_clamped=float(used.T) != x.source_temperature_k,
+            native_jet_liquid_mass_fraction=float(used.liquid_fraction),native_initial_downwind_distance_m=distance[0],
+            native_near_field_status="briggs_buoyant_effective_horizontal_source" if buoyant_vertical else "not_reconstructed_before_native_handoff")
+    if isinstance(x,(NativePoolFootprint,NativeJetFootprint)):
+        import numpy as np
+        from .spatial import ConcentrationGrid, threshold_grid_area
+        ny = 2*x.crosswind_points_per_side-1
+        ensure(ny*len(distance) <= 1000000,"Native planar grid exceeds1000000 nodes","engine_result_invalid")
+        y = np.linspace(-x.crosswind_half_span_m,x.crosswind_half_span_m,ny)
+        values = np.asarray(field.at_y(y),dtype=float)
+        ensure(values.shape == (ny,len(distance)),"Native crosswind field shape mismatch","engine_result_invalid")
+        ensure(np.all(np.isfinite(values)) and np.all((values >= 0)&(values <= 1)),"Nonphysical native spatial concentrations","engine_result_invalid")
+        grid = ConcentrationGrid(downwind_coordinates_m=distance,crosswind_coordinates_m=y.tolist(),
+            volume_fraction_by_crosswind_then_downwind=values.tolist(),threshold_volume_fraction=x.threshold_volume_fraction,
+            sampling_basis="Native ConcentrationField.at_y at common z/t_avg and source; main-diagonal triangular interpolation")
+        footprint = threshold_grid_area(grid)
+        if premise == "violated":
+            footprint["threshold_area_m2"] = None
+            footprint["footprint_status"] = "premise_violated"
+        output.update(footprint)
+        output.update(crosswind_distance_m=y.tolist(),volume_fraction_by_crosswind_then_downwind=values.tolist(),
+            spatial_output_scope="planar_sampled_threshold_field")
+    return output
+
+
+def main():
+    try:
+        payload = json.loads(sys.stdin.read())
+        operation = payload.pop("_operation","centreline")
+        models = {"centreline":NativePool,"footprint":NativePoolFootprint,"jet_centreline":NativeJet,"jet_footprint":NativeJetFootprint}
+        ensure(operation in models,"Unknown native worker operation","engine_protocol_failure")
+        request = models[operation].model_validate(payload)
+        diagnostics = io.StringIO()
+        with contextlib.redirect_stdout(diagnostics), contextlib.redirect_stderr(diagnostics):
+            result = run_native(request)
+        result["native_messages"] = diagnostics.getvalue()
+        print(json.dumps(result,allow_nan=False))
+    except Exception as exc:
+        code = exc.code if isinstance(exc,DomainError) else "engine_execution_failed"
+        print(json.dumps({"code":code,"error":str(exc)},ensure_ascii=True))
+        raise SystemExit(1)
+
+
+if __name__ == "__main__":
+    main()