.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/timeseries/plot_cmaqpandora.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_timeseries_plot_cmaqpandora.py>`
        to download the full example code

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

.. _sphx_glr_auto_examples_timeseries_plot_cmaqpandora.py:


CMAQ vs Pandora
===============

Timeseries comparison of NO2 columns from CMAQ and Pandora over a small window
over New York/New Jersey for 2018-07-01.

.. GENERATED FROM PYTHON SOURCE LINES 8-35



.. image-sg:: /auto_examples/timeseries/images/sphx_glr_plot_cmaqpandora_001.png
   :alt: plot cmaqpandora
   :srcset: /auto_examples/timeseries/images/sphx_glr_plot_cmaqpandora_001.png
   :class: sphx-glr-single-img





.. code-block:: default


    import pyrsig


    bbox = (-74.5, 40., -73.5, 41)
    cmaqkey = 'cmaq.equates.conus.integrated.NO2_COLUMN'
    datakey = 'pandora.L2_rnvs3p1_8.nitrogen_dioxide_vertical_column_amount'
    # Or use TropOMI or any other NO2 columns data
    # datakey = 'tropomi.offl.no2.nitrogendioxide_tropospheric_column'

    # Get a CMAQ file from RSIG
    api = pyrsig.RsigApi(
        bbox=bbox, bdate='2018-07-01T12', edate='2018-07-01T23:59:59'
    )
    ds = api.to_ioapi(cmaqkey)

    # pair_rsigcmaq will match the CMAQ bbox, bdate, and edate
    df = pyrsig.cmaq.pair_rsigcmaq(ds, 'NO2_COLUMN', datakey)

    pdf = df.groupby(['time']).mean(numeric_only=True)
    z1 = pdf['nitrogen_dioxide_vertical_column_amount']
    z2 = (pdf['CMAQ_NO2_COLUMN'] * 1e15)
    ax = z1.plot(marker='+', linestyle='none', label='Pandora')
    ax = z2.plot(ax=ax, color='green', label='CMAQ')
    ax.legend()
    ax.set(ylim=(0, None), ylabel='NO2 [molecules/cm2]', xlabel='Time [UTC]')
    ax.figure.savefig('cmaq_pandora.png')


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

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


.. _sphx_glr_download_auto_examples_timeseries_plot_cmaqpandora.py:

.. only:: html

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




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

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

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

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


.. only:: html

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

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