Python tools for use with GEOS-Chem

From Geos-chem
Revision as of 16:26, 14 June 2019 by Bmy (Talk | contribs) (Core Python packages for Earth Science data)

Jump to: navigation, search

On this page, we provide information about Python-language tools for visualizing or plotting GEOS-Chem data.


The GEOS-Chem Support Team is currently developing GCPy, which is used to create the evaluation plots from GEOS-Chem benchmark simulations. While GCPy is still in development, we encourage interested users to start using GCPy and to report any issues.

Also, if you have your own Python plotting or analysis code, and would like to make it a part of GCPy, please contact the GEOS-Chem Support Team.

Third-party tools

Several third-party Python packages have been developed for GEOS-Chem. If you are interested in using any of these, we recommend that you follow up with the developers directly. The GEOS-Chem Support Team is not responsible for maintaining these tools.

Developer Packages/Description Status
Jiawei Zhuang xESMF: A universal regridder. Leverages ESMF's regridding capabilities to easily regrid data between lon-lat and cubed-sphere grids. Other grids are supported as well. Available on github
Daniel Rothenberg xbpch: Backend for reading bpch output into xarray/dask Available on github
Barron Henderson Several software packages, including:
  1. PseudoNetCDF: a NetCDF like system including visualization (maps, profiles, timeseries, etc)
  2. Process analysis diagnostics: A tool for examining the change in each species due to each process and reaction.
Available on github

Core Python packages for Earth Science data

For the most part, you don't need to use a lot of Python packages to read and plot GEOS-Chem data. For example, GCPy package mostly relies on the following packages, which can be downloaded with either PyPI and/or Conda. You can get

  • xarray: For reading netCDF, HDF5, GRIB, etc. data into a common data structure for later manipulation
  • matplotlib: For general plotting operations
  • cartopy: For plotting data on world maps
  • numpy: For math operations
  • xESMF: For regrididdng between cubed-sphere to lat-lon grids
  • pandas: For spreadsheet-like data analysis

--Bob Yantosca (talk) 16:15, 14 June 2019 (UTC)

Tutorials and resources

If you are new to Python, we strongly recommend that you take our GEOS-Chem Python Tutorial. It will walk you though the steps of creating sample plots from GEOS-Chem output:

You may also find these resources useful:

--Bob Yantosca (talk) 16:21, 14 June 2019 (UTC)