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

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

.. _sphx_glr_auto_examples_cmaqready_example.py:


Preparing AQMBC Output for CMAQ
===============================

CMAQ expects BCON files to have 25 instantaneous outputs that cover from 00Z
on the simulation day to 00Z on the next day. `aqmbc` is run on individual files,
which can be a single hour or a few hours. Therefore, often the `aqmbc` outputs
need to be combined.

This example uses single hour outputs from GEOS-CF to make CMAQ-ready BCON
inputs. That is what `aqmbc.cmaq.cmaqready` is for. To do its jobs, `cmaqready`
takes a `date` to be made, the `inpaths` to make it from, and the outpath to
save it to.

The example below creates two days of output from 18 input files. The
inputs cover 2023-04-30T21:30Z to 2023-05-03T00:30Z in 3-hour intervals.
To create the inputs, modify the geoscf_example.py script to run for more
hours. The outputs each cover 25 hours from 00Z to 00Z+1day. To run the
example, you must have already created the inputs. After the run is complete,
you will have the outputs.

Required Inputs:
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-04-30T2130Z_12US1_BCON.nc
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-05-01T0030Z_12US1_BCON.nc
 - ... 14 others ...
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-05-02T2130Z_12US1_BCON.nc
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-05-03T0030Z_12US1_BCON.nc

Produced Outputs:
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-05-01_12US1_BCON_hourly.nc
 - geoscf_mcx_tavg_1hr_g1440x721_v36_2023-05-02_12US1_BCON_hourly.nc

.. GENERATED FROM PYTHON SOURCE LINES 35-37

Import libraries
----------------

.. GENERATED FROM PYTHON SOURCE LINES 37-41

.. code-block:: default


    from pandas import date_range, to_timedelta
    import aqmbc


.. GENERATED FROM PYTHON SOURCE LINES 42-44

User Inputs
-----------

.. GENERATED FROM PYTHON SOURCE LINES 44-56

.. code-block:: default


    # Use the domain name for your domain
    gdnam = '12US1'

    # Define the date range of interest
    outdates = date_range('2023-05-01', '2023-05-02', freq='1d')

    # Define a strftime templates for paths from aqmbc (intmpl) that need to be
    # concatenated and interpolated to make CMAQ-ready files (outtmpl)
    intmpl = f'geoscf_mcx_tavg_1hr_g1440x721_v36_%Y-%m-%dT%H%MZ_{gdnam}_BCON.nc'
    outtmpl = f'geoscf_mcx_tavg_1hr_g1440x721_v36_%Y-%m-%d_{gdnam}_BCON_hourly.nc'


.. GENERATED FROM PYTHON SOURCE LINES 57-59

Daily Processing
----------------

.. GENERATED FROM PYTHON SOURCE LINES 59-69

.. code-block:: default


    # For each output date, concatenate input files and interpolate hourly output.
    for date in outdates:
        # Define output path
        outpath = date.strftime(outtmpl)
        # Get paths for date-1day 21:30 ... date+1day 00:30 by 3-hours
        indates = date_range(date, freq='3h', periods=10) + to_timedelta('-2.5h')
        inpaths = indates.strftime(intmpl)
        # Make 1 output with 25h from 10 single-hour inputs
        aqmbc.cmaq.cmaqready(date, inpaths, outpath=outpath)


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

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


.. _sphx_glr_download_auto_examples_cmaqready_example.py:

.. only:: html

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




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

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

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

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


.. only:: html

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

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