.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/jena.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_jena.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_jena.py:


=================================
Sub-hourly ET from Jena, Germany
=================================

.. GENERATED FROM PYTHON SOURCE LINES 6-15

.. code-block:: Python


    if __name__ == '__main__':
        import os
        import site
        # add parent directory to path
        et_dir = os.path.dirname(os.path.dirname(os.path.abspath(__file__)))
        print(et_dir)
        site.addsitedir(et_dir)



.. rst-class:: sphx-glr-script-out

.. code-block:: pytb

    Traceback (most recent call last):
      File "/home/docs/checkouts/readthedocs.org/user_builds/etutil/checkouts/latest/examples/jena.py", line 11, in <module>
        et_dir = os.path.dirname(os.path.dirname(os.path.abspath(__file__)))
    NameError: name '__file__' is not defined




.. GENERATED FROM PYTHON SOURCE LINES 16-50

.. code-block:: Python


    import pandas as pd
    from ETUtil import PenmanMonteith

    _df = pd.read_csv(os.path.join(et_dir, 'data/jen_obj.txt'))
    # https://www.bgc-jena.mpg.de/wetter/Weatherstation.pdf
    # TODO input df contains 'SWDR (W/m**2)', can that be used?
    df = pd.DataFrame()
    df[['rel_hum', 'wind_speed', 'temp', 'tdew',
        'rn']] = _df[['rh (%)', 'wv (m/s)', 'T (degC)', 'Tdew (degC)', 'Rn (W/m**2)']]

    # rn = df['rn']
    # rn = np.where(rn < 0.0, 0.0, rn)
    # df['rn'] = rn

    df.index = pd.to_datetime(_df['index'])
    df.index.freq = pd.infer_freq(df.index)

    units = {'temp': 'Centigrade',
             'rel_hum': 'percent',
             'tdew': 'Fahrenheit',
             'wind_speed': 'MilesPerHour'}

    constants = dict()
    constants['lat_dec_deg'] = 50.9
    constants['altitude'] = 155
    # These values are not accurate
    constants['a_s'] = 0.23
    constants['albedo'] = 0.23
    constants['b_s'] = 0.5
    constants['wind_z'] = 2

    eto_model = PenmanMonteith(df, units=units, constants=constants, verbosity=2)


.. GENERATED FROM PYTHON SOURCE LINES 51-54

.. code-block:: Python


    pet = eto_model()


.. GENERATED FROM PYTHON SOURCE LINES 55-58

.. code-block:: Python


    pet.head()


.. GENERATED FROM PYTHON SOURCE LINES 59-62

.. code-block:: Python


    eto_model.plot_inputs()


.. GENERATED FROM PYTHON SOURCE LINES 63-65

.. code-block:: Python


    eto_model.plot_outputs()


.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 0.001 seconds)


.. _sphx_glr_download_auto_examples_jena.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: jena.ipynb <jena.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: jena.py <jena.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: jena.zip <jena.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_