scopeinpython.thermal
=======================

Scalar-per-timestep thermal/aerodynamic building blocks used by the
SCOPE energy-balance loop (``ebal.R``, not yet ported). Direct ports of
``SCOPEinR/R/Monin_ObuKhov.R``, ``SCOPEinR/R/resistances.R``, and
``SCOPEinR/R/heatfluxes.R`` -- the first pieces of the thermal
energy-balance chain to be ported.

.. note::
   ``get_resistances``'s ``rac``/``rws`` outputs reproduce a real R
   quirk: they use the *stability-uncorrected* eddy diffusivity (the
   local ``Kh`` variable in ``SCOPEinR::get.resistances``), not the
   stability-corrected value returned as the function's own ``Kh``
   output field. Confirmed by direct reading of the R source (the local
   ``Kh`` is never reassigned after the corrected value is written into
   ``resist_out[['Kh']]``) -- reproduced exactly rather than "fixed",
   since there's no independent way to tell whether this was intentional.

.. automodule:: scopeinpython.thermal
   :members:
   :undoc-members:
   :show-inheritance: