barronh

import time
import glob
import cmaqsatproc as csp
import matplotlib.pyplot as plt
import matplotlib.colors as mc
import matplotlib.animation as animation
import numpy as np

# Get a list of paths
pat = '/work/ROMO/2021_MP/cmaq54mp/2021hb_MP_cb6r5ae7_21k/12US1/output/CCTM_CGRID_2021hb_MP_cb6r5ae7_21k.12US1_35_2021010?.nc'

barronh

__doc__ = """
# Make CMAQ WRF-Like

     author: Barron H. Henderson
     date: 2020-07-30
     updated: 2024-06-04

# Overview

* Need: Many tools exist for WRFChem that we would like to use with CMAQ.

barronh

barronh

#!/usr/bin/env python
__doc__ = """
Overview
========

Module designed to make backups of scripts (.sh, .csh, .py), makefiles,
configuration files, and README files.

Contents
--------

barronh

__doc__ = """
Offline CMAQ Python WBD approximator
====================================

github.com/USEPA/CMAQ/emis/emis/DUST_EMIS.F routines and functions converted
from Fortran to Python. All modules are replaced by inline code and direct
access to MCIP files. The original DUST_EMIS.F was developed by Foroutan et
al. (https://dx.doi.org/10.1002/2016MS0008230.

To-Do:

barronh

__all__ = ['pair_rsigcmaq']
__doc__ = """
pair_rsigcmaq is designed to pair RSIG data with CMAQ files.

Example:

    from pair_rsigcmaq import pair_rsigcmaq
    import pandas a pd

    outpath = pair_rsigcmaq('CMAQ_ACONC_20160101', 'O3', 'airnow.ozone')

barronh

#!/bin/bash

POSITIONAL_ARGS=()
PREFIX="ENV"
while [[ $# -gt 0 ]]; do
  case $1 in
    -h|--help)
      echo "Usage"
      echo "-----"
      echo "    $0 [--prefix=PREFIX] TEMPLATE_ENV STARTDAY [ENDDAY [SCRIPT]]"

barronh

"""
This example calculates statistics for a single day from the EPA EQUATES CMAQ simulation for 2018-06-14
and creates a taylor diagram.
"""


import pyrsig
import pandas as pd
import pyproj
import taylordiagram as td

barronh

__doc__ = """
# Overview

This example shows how to open High Resolution Rapid Refresh (HRRR) meteorology
"surfaces." These surfaces are slices in vertical space to create single level
maps of variables like temperature, wind speed components, pressure, etc. Each
surface is opened using the Zarr archive described in Amazon's opendata registry
described at https://registry.opendata.aws/noaa-hrrr-pds/.

# Examples

barronh

__all__ = ['getarea']


def getarea(nlat, nlon, late=None, R=6371000.):
    """
    Arguments
    ---------
    nlat : int
        Number of latitude bands
    nlon : int