pyrestoolbox 3.7.3__tar.gz → 3.7.5__tar.gz

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/.gitignore

@@ -46,3 +46,12 @@ BUILD_INSTRUCTIONS.txt
 
 # local virtualenv (389 MB, never committed)
 .venv/
+
+# Internal development artifacts (kept locally, not published)
+reviews/
+benchmarks/
+benchmark_*.py
+zstar_*
+pyresToolbox_rust_acceleration_instructions.md
+pyrestoolbox-mcp-email.txt
+*.zip

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/Cargo.lock

@@ -152,7 +152,7 @@ dependencies = [
 
 [[package]]
 name = "pyrestoolbox-native"
-version = "3.7.0"
+version = "3.7.5"
 dependencies = [
  "pyo3",
  "rug",

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "pyrestoolbox-native"
-version = "3.7.0"
+version = "3.7.5"
 edition = "2021"
 license = "GPL-3.0-or-later"
 description = "Rust acceleration extensions for pyResToolbox"

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pyrestoolbox
-Version: 3.7.3
+Version: 3.7.5
 Classifier: Programming Language :: Python :: 3
 Classifier: Programming Language :: Rust
 Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
@@ -24,13 +24,13 @@ Requires-Python: >=3.8
 Description-Content-Type: text/x-rst; charset=UTF-8
 Project-URL: Homepage, https://github.com/mwburgoyne/pyResToolbox
 
-===================================
+================
 ``pyrestoolbox``
-===================================
+================
 
---------------------------------------------------
+-----------------------------------------------
 A collection of Reservoir Engineering Utilities
---------------------------------------------------
+-----------------------------------------------
 
 This set of functions focuses on those that the author uses often while crafting programming solutions.
 These are the scripts that are often copy/pasted from previous work - sometimes slightly modified - resulting in a trail of slightly different versions over the years. Some attempt has been made here to make this implementation flexible enough such that it can be relied on as-is going forward.
@@ -39,7 +39,7 @@ Note: Version 3.0 consolidates simulation-oriented functions under the simtools
 
 Includes functions to perform calculations including;
 
-- Decline Curve Analysis with Arps and Duong models — fitting, forecasting, EUR, windowed fitting, secondary phase ratio models, and uptime inference
+- Decline Curve Analysis with Arps and Duong models - fitting, forecasting, EUR, windowed fitting, secondary phase ratio models, and uptime inference
 - Material Balance for gas (P/Z with Cole plot and Havlena-Odeh aquifer support) and oil (Havlena-Odeh with drive indices, parameter regression, and tabulated PVT)
 - Inflow Performance Relationships (IPR) for oil and gas wells
 - Vertical Lift Performance (VLP) with four multiphase flow correlations (Hagedorn-Brown, Woldesemayat-Ghajar, Gray, Beggs & Brill)
@@ -47,6 +47,7 @@ Includes functions to perform calculations including;
 - Eclipse VFP table generation (VFPPROD and VFPINJ keywords)
 - PVT Calculations for oil
 - PVT calculation for gas, including up to 100% inerts for CO2, H2S, N2 and H2
+- Gas caloric properties (enthalpy, Cp, Cv, Joule-Thomson coefficient) from the BNS tuned Peng Robinson EOS
 - Gas hydrate formation prediction with thermodynamic inhibitor calculations
 - Return critical parameters for typical components
 - Creation of Black Oil Table information (PVDO, PVDG, PVTO, PVTW keywords)
@@ -55,39 +56,40 @@ Includes functions to perform calculations including;
 - Generation of AQUTAB include file influence functions for use in ECLIPSE
 - Creation of Corey, LET and Jerauld relative permeability tables in Eclipse format, with curve fitting support
 - Calculation of Methane, CO2 and multicomponent gas saturated brine properties (Soreide-Whitson VLE)
+- Multi-salt gas-free brine viscosity: IAPWS-2008 water, an ion-additive salt ratio over any combination of Na, K, Mg, Ca, Sr, Ba, Cl, Br, sulfate and carbonate, and a measured pressure factor
 - Method recommendation engine for selecting appropriate correlations based on fluid composition
 - Sensitivity analysis with parameter sweeps and tornado charts
 
 All public PVT, flow rate, and simulation table functions support both oilfield (psia, deg F, ft) and Eclipse METRIC (barsa, deg C, m) unit systems via an optional ``metric=False`` parameter. See individual module documentation for unit mapping details.
 
 Rust Acceleration (Optional)
------------------------------
+----------------------------
 
 pyResToolbox includes optional Rust-compiled extensions that accelerate computationally intensive algorithms. When the compiled extension is present and loadable, these functions run automatically through Rust with no API changes. When the extension is unavailable, all functions fall back silently to the pure Python implementation.
 
 **Accelerated functions:**
 
-- **Nodal VLP segment loops** — all 8 VLP method functions (4 methods x gas/oil)
-- **Gas Z-factor** — DAK, Hall-Yarborough, and BNS full-pipeline calculations
-- **Gas viscosity** — LGE and LBC correlations
-- **Gas pseudopressure** — Gauss-Legendre quadrature integration
-- **Oil density** — Standing-Witte-McCain-Hill (iterative and above-Pb)
-- **Oil FVF** — McCain density-based method
-- **DCA hyperbolic fitting** — grid search with RANSAC (``fit_decline``, ``fit_decline_cum``)
-- **Material balance** — oil matbal regression objective function
-- **CO2-Brine solubility** — Spycher-Pruess iterative RK-EOS solver
-- **VLE flash** — Soreide-Whitson multi-component Peng-Robinson flash
+- **Nodal VLP segment loops** - all 8 VLP method functions (4 methods x gas/oil)
+- **Gas Z-factor** - DAK, Hall-Yarborough, and BNS full-pipeline calculations
+- **Gas viscosity** - LGE and LBC correlations
+- **Gas pseudopressure** - Gauss-Legendre quadrature integration
+- **Oil density** - Standing-Witte-McCain-Hill (iterative and above-Pb)
+- **Oil FVF** - McCain density-based method
+- **DCA hyperbolic fitting** - grid search with RANSAC (``fit_decline``, ``fit_decline_cum``)
+- **Material balance** - oil matbal regression objective function
+- **CO2-Brine solubility** - Spycher-Pruess iterative RK-EOS solver
+- **VLE flash** - Soreide-Whitson multi-component Peng-Robinson flash
 
 **Behavior:**
 
 - If the Rust extension is not found on disk, pure Python is used with no warning
 - If the extension fails to load (e.g. OS permission restrictions, architecture mismatch), a sentinel file is written to avoid repeated probe attempts on subsequent imports. The sentinel is automatically invalidated when the extension file changes (new build or update)
-- All Rust-accelerated paths use ``try/except`` wrappers — any Rust-side error falls back to the Python implementation transparently
+- All Rust-accelerated paths use ``try/except`` wrappers - any Rust-side error falls back to the Python implementation transparently
 
 **Environment variables:**
 
-- ``PYRESTOOLBOX_NO_RUST=1`` — Force pure Python mode (skip Rust extension entirely)
-- ``PYRESTOOLBOX_RETRY_RUST=1`` — Ignore the sentinel file and retry loading the extension
+- ``PYRESTOOLBOX_NO_RUST=1`` - Force pure Python mode (skip Rust extension entirely)
+- ``PYRESTOOLBOX_RETRY_RUST=1`` - Ignore the sentinel file and retry loading the extension
 
 **Programmatic status check:**
 
@@ -109,7 +111,7 @@ Upgrade previous installations with
 
 
 Module List
-=============
+===========
 
 .. list-table::
    :widths: 30 70
@@ -119,7 +121,7 @@ Module List
    * - `matbal <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/matbal.rst>`_
      - Gas P/Z material balance (OGIP), Cole plot diagnostics, Havlena-Odeh with aquifer influx, Oil Havlena-Odeh material balance (OOIP) with drive indices, Parameter regression with bounds, Tabulated PVT support
    * - `gas <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/gas.rst>`_
-     - Gas Tc & Pc Calculation, Gas Z-Factor Calculation, Gas Viscosity, Gas Viscosity \* Z, Gas Compressibility, Gas Formation Volume Factor, Gas Density, Gas Water of Condensation, Convert P/Z to P, Convert Gas Gradient to SG, Delta Pseudopressure, Gas Condensate FWS SG, Gas Flow Rate Radial, Gas Flow Rate Linear, Gas Hydrate Prediction
+     - Gas Tc & Pc Calculation, Gas Z-Factor Calculation, Gas Viscosity, Gas Viscosity \* Z, Gas Compressibility, Gas Formation Volume Factor, Gas Density, Gas Water of Condensation, Convert P/Z to P, Convert Gas Gradient to SG, Delta Pseudopressure, Gas Condensate FWS SG, Gas Flow Rate Radial, Gas Flow Rate Linear, Gas Hydrate Prediction, Gas Caloric Properties (H, Cp, Cv, Joule-Thomson), Non-Darcy and Partial-Penetration Pseudoskins
    * - `oil <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/oil.rst>`_
      - Oil Density from MW, Oil Critical Properties with Twu, Incremental GOR post Separation, Oil Bubble Point Pressure, Oil GOR at Pb, Oil GOR at P, Oil Compressibility, Oil Density, Oil Formation Volume Factor, Oil Viscosity, Harmonize Pb and Rsb, Estimate soln gas SG from oil, Estimate SG of gas post separator, Calculate weighted average surface gas SG, Oil API to SG, Oil SG to API, Oil Flow Rate Radial, Oil Flow Rate Linear
    * - `nodal <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/nodal.rst>`_
@@ -127,7 +129,7 @@ Module List
    * - `library <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/library.rst>`_
      - Return critical parameters for typical single components
    * - `brine <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/brine.rst>`_
-     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE)
+     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE), multi-salt gas-free brine viscosity
    * - `layer <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/layer.rst>`_
      - Lorenz coefficient from Beta value, Lorenz coefficient from flow fraction, Lorenz coefficient to flow fraction, Lorenz coefficient to permeability array
    * - `simtools <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/simtools.rst>`_
@@ -155,8 +157,8 @@ A simple example below of estimating oil bubble point pressure.
 .. code-block:: python
 
     >>> from pyrestoolbox import oil
-    >>> oil.oil_pbub(api=43, degf=185, rsb=2350, sg_g =0.72, pbmethod ='VALMC')
-    5179.51086900132
+    >>> oil.oil_pbub(api=43, degf=185, rsb=2350, sg_g=0.72, pbmethod='VALMC')
+    5199.2406069808885
 
 Fit a decline curve to production data and estimate ultimate recovery.
 
@@ -168,9 +170,9 @@ Fit a decline curve to production data and estimate ultimate recovery.
     >>> q = 1000 * np.exp(-0.05 * t)
     >>> result = dca.fit_decline(t, q, method='exponential')
     >>> result.qi, result.di
-    (1000.0000000000007, 0.05000000000000006)
+    (999.9999999999953, 0.049999999999999885)
     >>> dca.eur(qi=result.qi, di=result.di, b=0, q_min=10)
-    19799.99999999999
+    19799.999999999953
 
 Estimate gas in place from pressure-production history.
 
@@ -183,9 +185,9 @@ Estimate gas in place from pressure-production history.
     ...     degf=200, sg=0.65
     ... )
     >>> r.ogip
-    87.602774253829
+    87.60264634235122
     >>> r.r_squared
-    0.9734794008096929
+    0.9734793606102145
 
 A set of Gas-Oil relative permeability curves with the LET method
 

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/README.rst

@@ -1,10 +1,10 @@
-===================================
+================
 ``pyrestoolbox``
-===================================
+================
 
---------------------------------------------------
+-----------------------------------------------
 A collection of Reservoir Engineering Utilities
---------------------------------------------------
+-----------------------------------------------
 
 This set of functions focuses on those that the author uses often while crafting programming solutions.
 These are the scripts that are often copy/pasted from previous work - sometimes slightly modified - resulting in a trail of slightly different versions over the years. Some attempt has been made here to make this implementation flexible enough such that it can be relied on as-is going forward.
@@ -13,7 +13,7 @@ Note: Version 3.0 consolidates simulation-oriented functions under the simtools
 
 Includes functions to perform calculations including;
 
-- Decline Curve Analysis with Arps and Duong models — fitting, forecasting, EUR, windowed fitting, secondary phase ratio models, and uptime inference
+- Decline Curve Analysis with Arps and Duong models - fitting, forecasting, EUR, windowed fitting, secondary phase ratio models, and uptime inference
 - Material Balance for gas (P/Z with Cole plot and Havlena-Odeh aquifer support) and oil (Havlena-Odeh with drive indices, parameter regression, and tabulated PVT)
 - Inflow Performance Relationships (IPR) for oil and gas wells
 - Vertical Lift Performance (VLP) with four multiphase flow correlations (Hagedorn-Brown, Woldesemayat-Ghajar, Gray, Beggs & Brill)
@@ -21,6 +21,7 @@ Includes functions to perform calculations including;
 - Eclipse VFP table generation (VFPPROD and VFPINJ keywords)
 - PVT Calculations for oil
 - PVT calculation for gas, including up to 100% inerts for CO2, H2S, N2 and H2
+- Gas caloric properties (enthalpy, Cp, Cv, Joule-Thomson coefficient) from the BNS tuned Peng Robinson EOS
 - Gas hydrate formation prediction with thermodynamic inhibitor calculations
 - Return critical parameters for typical components
 - Creation of Black Oil Table information (PVDO, PVDG, PVTO, PVTW keywords)
@@ -29,39 +30,40 @@ Includes functions to perform calculations including;
 - Generation of AQUTAB include file influence functions for use in ECLIPSE
 - Creation of Corey, LET and Jerauld relative permeability tables in Eclipse format, with curve fitting support
 - Calculation of Methane, CO2 and multicomponent gas saturated brine properties (Soreide-Whitson VLE)
+- Multi-salt gas-free brine viscosity: IAPWS-2008 water, an ion-additive salt ratio over any combination of Na, K, Mg, Ca, Sr, Ba, Cl, Br, sulfate and carbonate, and a measured pressure factor
 - Method recommendation engine for selecting appropriate correlations based on fluid composition
 - Sensitivity analysis with parameter sweeps and tornado charts
 
 All public PVT, flow rate, and simulation table functions support both oilfield (psia, deg F, ft) and Eclipse METRIC (barsa, deg C, m) unit systems via an optional ``metric=False`` parameter. See individual module documentation for unit mapping details.
 
 Rust Acceleration (Optional)
------------------------------
+----------------------------
 
 pyResToolbox includes optional Rust-compiled extensions that accelerate computationally intensive algorithms. When the compiled extension is present and loadable, these functions run automatically through Rust with no API changes. When the extension is unavailable, all functions fall back silently to the pure Python implementation.
 
 **Accelerated functions:**
 
-- **Nodal VLP segment loops** — all 8 VLP method functions (4 methods x gas/oil)
-- **Gas Z-factor** — DAK, Hall-Yarborough, and BNS full-pipeline calculations
-- **Gas viscosity** — LGE and LBC correlations
-- **Gas pseudopressure** — Gauss-Legendre quadrature integration
-- **Oil density** — Standing-Witte-McCain-Hill (iterative and above-Pb)
-- **Oil FVF** — McCain density-based method
-- **DCA hyperbolic fitting** — grid search with RANSAC (``fit_decline``, ``fit_decline_cum``)
-- **Material balance** — oil matbal regression objective function
-- **CO2-Brine solubility** — Spycher-Pruess iterative RK-EOS solver
-- **VLE flash** — Soreide-Whitson multi-component Peng-Robinson flash
+- **Nodal VLP segment loops** - all 8 VLP method functions (4 methods x gas/oil)
+- **Gas Z-factor** - DAK, Hall-Yarborough, and BNS full-pipeline calculations
+- **Gas viscosity** - LGE and LBC correlations
+- **Gas pseudopressure** - Gauss-Legendre quadrature integration
+- **Oil density** - Standing-Witte-McCain-Hill (iterative and above-Pb)
+- **Oil FVF** - McCain density-based method
+- **DCA hyperbolic fitting** - grid search with RANSAC (``fit_decline``, ``fit_decline_cum``)
+- **Material balance** - oil matbal regression objective function
+- **CO2-Brine solubility** - Spycher-Pruess iterative RK-EOS solver
+- **VLE flash** - Soreide-Whitson multi-component Peng-Robinson flash
 
 **Behavior:**
 
 - If the Rust extension is not found on disk, pure Python is used with no warning
 - If the extension fails to load (e.g. OS permission restrictions, architecture mismatch), a sentinel file is written to avoid repeated probe attempts on subsequent imports. The sentinel is automatically invalidated when the extension file changes (new build or update)
-- All Rust-accelerated paths use ``try/except`` wrappers — any Rust-side error falls back to the Python implementation transparently
+- All Rust-accelerated paths use ``try/except`` wrappers - any Rust-side error falls back to the Python implementation transparently
 
 **Environment variables:**
 
-- ``PYRESTOOLBOX_NO_RUST=1`` — Force pure Python mode (skip Rust extension entirely)
-- ``PYRESTOOLBOX_RETRY_RUST=1`` — Ignore the sentinel file and retry loading the extension
+- ``PYRESTOOLBOX_NO_RUST=1`` - Force pure Python mode (skip Rust extension entirely)
+- ``PYRESTOOLBOX_RETRY_RUST=1`` - Ignore the sentinel file and retry loading the extension
 
 **Programmatic status check:**
 
@@ -83,7 +85,7 @@ Upgrade previous installations with
 
 
 Module List
-=============
+===========
 
 .. list-table::
    :widths: 30 70
@@ -93,7 +95,7 @@ Module List
    * - `matbal <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/matbal.rst>`_
      - Gas P/Z material balance (OGIP), Cole plot diagnostics, Havlena-Odeh with aquifer influx, Oil Havlena-Odeh material balance (OOIP) with drive indices, Parameter regression with bounds, Tabulated PVT support
    * - `gas <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/gas.rst>`_
-     - Gas Tc & Pc Calculation, Gas Z-Factor Calculation, Gas Viscosity, Gas Viscosity \* Z, Gas Compressibility, Gas Formation Volume Factor, Gas Density, Gas Water of Condensation, Convert P/Z to P, Convert Gas Gradient to SG, Delta Pseudopressure, Gas Condensate FWS SG, Gas Flow Rate Radial, Gas Flow Rate Linear, Gas Hydrate Prediction
+     - Gas Tc & Pc Calculation, Gas Z-Factor Calculation, Gas Viscosity, Gas Viscosity \* Z, Gas Compressibility, Gas Formation Volume Factor, Gas Density, Gas Water of Condensation, Convert P/Z to P, Convert Gas Gradient to SG, Delta Pseudopressure, Gas Condensate FWS SG, Gas Flow Rate Radial, Gas Flow Rate Linear, Gas Hydrate Prediction, Gas Caloric Properties (H, Cp, Cv, Joule-Thomson), Non-Darcy and Partial-Penetration Pseudoskins
    * - `oil <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/oil.rst>`_
      - Oil Density from MW, Oil Critical Properties with Twu, Incremental GOR post Separation, Oil Bubble Point Pressure, Oil GOR at Pb, Oil GOR at P, Oil Compressibility, Oil Density, Oil Formation Volume Factor, Oil Viscosity, Harmonize Pb and Rsb, Estimate soln gas SG from oil, Estimate SG of gas post separator, Calculate weighted average surface gas SG, Oil API to SG, Oil SG to API, Oil Flow Rate Radial, Oil Flow Rate Linear
    * - `nodal <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/nodal.rst>`_
@@ -101,7 +103,7 @@ Module List
    * - `library <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/library.rst>`_
      - Return critical parameters for typical single components
    * - `brine <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/brine.rst>`_
-     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE)
+     - Calculate suite of brine properties with variable methane, Calculate suite of CO2 saturated brine properties, Multicomponent gas-saturated brine (Soreide-Whitson VLE), multi-salt gas-free brine viscosity
    * - `layer <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/layer.rst>`_
      - Lorenz coefficient from Beta value, Lorenz coefficient from flow fraction, Lorenz coefficient to flow fraction, Lorenz coefficient to permeability array
    * - `simtools <https://github.com/mwburgoyne/pyResToolbox/blob/main/pyrestoolbox/docs/simtools.rst>`_
@@ -129,8 +131,8 @@ A simple example below of estimating oil bubble point pressure.
 .. code-block:: python
 
     >>> from pyrestoolbox import oil
-    >>> oil.oil_pbub(api=43, degf=185, rsb=2350, sg_g =0.72, pbmethod ='VALMC')
-    5179.51086900132
+    >>> oil.oil_pbub(api=43, degf=185, rsb=2350, sg_g=0.72, pbmethod='VALMC')
+    5199.2406069808885
 
 Fit a decline curve to production data and estimate ultimate recovery.
 
@@ -142,9 +144,9 @@ Fit a decline curve to production data and estimate ultimate recovery.
     >>> q = 1000 * np.exp(-0.05 * t)
     >>> result = dca.fit_decline(t, q, method='exponential')
     >>> result.qi, result.di
-    (1000.0000000000007, 0.05000000000000006)
+    (999.9999999999953, 0.049999999999999885)
     >>> dca.eur(qi=result.qi, di=result.di, b=0, q_min=10)
-    19799.99999999999
+    19799.999999999953
 
 Estimate gas in place from pressure-production history.
 
@@ -157,9 +159,9 @@ Estimate gas in place from pressure-production history.
     ...     degf=200, sg=0.65
     ... )
     >>> r.ogip
-    87.602774253829
+    87.60264634235122
     >>> r.r_squared
-    0.9734794008096929
+    0.9734793606102145
 
 A set of Gas-Oil relative permeability curves with the LET method
 

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "maturin"
 
 [project]
 name = "pyrestoolbox"
-version = "3.7.3"
+version = "3.7.5"
 description = "pyResToolbox - A collection of Reservoir Engineering Utilities"
 license = {text = "GPL-3.0-or-later"}
 authors = [{name = "Mark W. Burgoyne", email = "mark.w.burgoyne@gmail.com"}]

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/pyrestoolbox/brine/__init__.py

@@ -1,4 +1,5 @@
 from .brine import *
+from .brine import recommended_method  # which entry point to use, and why
 from . import pr_vphi
 from .pr_vphi import V2_inf as V_phi_pr
 from . import vphi_route

{pyrestoolbox-3.7.3 → pyrestoolbox-3.7.5}/pyrestoolbox/brine/_lib_maoduan_viscosity.py

@@ -1,9 +1,11 @@
 """Mao-Duan (2009) water and NaCl viscosity, the SUPERSEDED route.
 
 Ported from brine_props `salt_viscosity_benchmark.py`. Retained only so
-`viscosity_route(route='mao_duan')` reproduces pre-3.7.3 results; it is NaCl
+`viscosity_route(route='mao_duan')` reproduces pre-3.7.3 results (to within
+0.0094% since the d1/b2 fix below); it is NaCl
 only and carries no pressure dependence on the salt ratio. Coefficients are
-McCain Table 4-14 (Mao & Duan 2009).
+Mao & Duan (2009) Tables 2-3, verified against the primary 2026-08-08; the
+McCain Table 4-14 lineage carried two last-digit defects (d1, b2, <=0.0094%).
 
 Against Kestin's measured NaCl viscosity this ties the default Jones-Dole leg
 (0.443% against 0.415% mean) and loses at reservoir temperature (0.754% against
@@ -14,10 +16,10 @@ import math
 
 from ..plyasunov.iapws_if97 import rho_if97
 
-_MD_D = [0, 2885310.0, -11072.577, -9.0834095, 0.030925651, -0.0000274071,
+_MD_D = [0, 2885317.0, -11072.577, -9.0834095, 0.030925651, -0.0000274071,
          -1928385.1, 5621.6046, 13.82725, -0.047609523, 0.000035545041]
 _MD_A = [-0.21319213, 0.0013651589, -0.0000012191756]
-_MD_B = [0.069161945, -0.00027292263, 0.0000002085244]
+_MD_B = [0.069161945, -0.00027292263, 0.00000020852448]
 _MD_C = [-0.0025988855, 0.0000077989227]