Skip to content

Conversation

yantosca
Copy link
Contributor

@yantosca yantosca commented Sep 5, 2025

Name and Institution (Required)

Name: Bob Yantosca
Institution: Harvard + GCST

Describe the update

This PR adds a GCPy example script that can display the lons & lats for the GEOS-Chem Classic horizontal grids. This script was used to generate the GEOS-Chem Classic horizontal grids documentation on geos-chem.readthedocs.io. This will make it easy to update the ReadTheDocs documentation if other grids are added to GEOS-Chem Classic in the future.

Expected changes

Usage is as follows:

$ python -m gcpy.examples.grids.display_gcclassic_grid_info
Please select a GEOS-Chem Classic horizontal grid:

1.  4.0°   x 5.0°     global

2.  2.0°   x 2.5°     global

3.  0.5°   x 0.625°   global
4.  0.5°   x 0.625°   nested AS
5.  0.5°   x 0.625°   nested EU
6.  0.5°   x 0.625°   nested NA

7   0.25°  x 0.3125°  global
8.  0.25°  x 0.3125°  nested AF
9.  0.25°  x 0.3125°  nested AS
10. 0.25°  x 0.3125°  nested EU
11. 0.25°  x 0.3125°  nested ME
12. 0.25°  x 0.3125°  nested NA
13  0.25°  x 0.3125°  nested OC
14. 0.25°  x 0.3125°  nested SA
15. 0.25°  x 0.3125°  nested RU

16. 0.125° x 0.15625° global
17. 0.125° x 0.15625° nested AF
18. 0.125° x 0.15625° nested AS
19. 0.125° x 0.15625° nested EU
20. 0.125° x 0.15625° nested NA
21. 0.125° x 0.15625° nested SA
    
Enter a selection >>> 6
resolution
0.5x0.625
longitude centers
 -140.00000 -139.37500 -138.75000 -138.12500 -137.50000 -136.87500 -136.25000 -135.62500
 -135.00000 -134.37500 -133.75000 -133.12500 -132.50000 -131.87500 -131.25000 -130.62500
 -130.00000 -129.37500 -128.75000 -128.12500 -127.50000 -126.87500 -126.25000 -125.62500
 -125.00000 -124.37500 -123.75000 -123.12500 -122.50000 -121.87500 -121.25000 -120.62500
 -120.00000 -119.37500 -118.75000 -118.12500 -117.50000 -116.87500 -116.25000 -115.62500
 -115.00000 -114.37500 -113.75000 -113.12500 -112.50000 -111.87500 -111.25000 -110.62500
 -110.00000 -109.37500 -108.75000 -108.12500 -107.50000 -106.87500 -106.25000 -105.62500
 -105.00000 -104.37500 -103.75000 -103.12500 -102.50000 -101.87500 -101.25000 -100.62500
 -100.00000  -99.37500  -98.75000  -98.12500  -97.50000  -96.87500  -96.25000  -95.62500
  -95.00000  -94.37500  -93.75000  -93.12500  -92.50000  -91.87500  -91.25000  -90.62500
  -90.00000  -89.37500  -88.75000  -88.12500  -87.50000  -86.87500  -86.25000  -85.62500
  -85.00000  -84.37500  -83.75000  -83.12500  -82.50000  -81.87500  -81.25000  -80.62500
  -80.00000  -79.37500  -78.75000  -78.12500  -77.50000  -76.87500  -76.25000  -75.62500
  -75.00000  -74.37500  -73.75000  -73.12500  -72.50000  -71.87500  -71.25000  -70.62500
  -70.00000  -69.37500  -68.75000  -68.12500  -67.50000  -66.87500  -66.25000  -65.62500
  -65.00000  -64.37500  -63.75000  -63.12500  -62.50000  -61.87500  -61.25000  -60.62500
  -60.00000  -59.37500  -58.75000  -58.12500  -57.50000  -56.87500  -56.25000  -55.62500
  -55.00000  -54.37500  -53.75000  -53.12500  -52.50000  -51.87500  -51.25000  -50.62500
  -50.00000  -49.37500  -48.75000  -48.12500  -47.50000  -46.87500  -46.25000  -45.62500
  -45.00000  -44.37500  -43.75000  -43.12500  -42.50000  -41.87500  -41.25000  -40.62500
  -40.00000
latitude centers
   10.00000   10.50000   11.00000   11.50000   12.00000   12.50000   13.00000   13.50000
   14.00000   14.50000   15.00000   15.50000   16.00000   16.50000   17.00000   17.50000
   18.00000   18.50000   19.00000   19.50000   20.00000   20.50000   21.00000   21.50000
   22.00000   22.50000   23.00000   23.50000   24.00000   24.50000   25.00000   25.50000
   26.00000   26.50000   27.00000   27.50000   28.00000   28.50000   29.00000   29.50000
   30.00000   30.50000   31.00000   31.50000   32.00000   32.50000   33.00000   33.50000
   34.00000   34.50000   35.00000   35.50000   36.00000   36.50000   37.00000   37.50000
   38.00000   38.50000   39.00000   39.50000   40.00000   40.50000   41.00000   41.50000
   42.00000   42.50000   43.00000   43.50000   44.00000   44.50000   45.00000   45.50000
   46.00000   46.50000   47.00000   47.50000   48.00000   48.50000   49.00000   49.50000
   50.00000   50.50000   51.00000   51.50000   52.00000   52.50000   53.00000   53.50000
   54.00000   54.50000   55.00000   55.50000   56.00000   56.50000   57.00000   57.50000
   58.00000   58.50000   59.00000   59.50000   60.00000   60.50000   61.00000   61.50000
   62.00000   62.50000   63.00000   63.50000   64.00000   64.50000   65.00000   65.50000
   66.00000   66.50000   67.00000   67.50000   68.00000   68.50000   69.00000   69.50000
   70.00000
longitude edges
 -140.31250 -139.68750 -139.06250 -138.43750 -137.81250 -137.18750 -136.56250 -135.93750
 -135.31250 -134.68750 -134.06250 -133.43750 -132.81250 -132.18750 -131.56250 -130.93750
 -130.31250 -129.68750 -129.06250 -128.43750 -127.81250 -127.18750 -126.56250 -125.93750
 -125.31250 -124.68750 -124.06250 -123.43750 -122.81250 -122.18750 -121.56250 -120.93750
 -120.31250 -119.68750 -119.06250 -118.43750 -117.81250 -117.18750 -116.56250 -115.93750
 -115.31250 -114.68750 -114.06250 -113.43750 -112.81250 -112.18750 -111.56250 -110.93750
 -110.31250 -109.68750 -109.06250 -108.43750 -107.81250 -107.18750 -106.56250 -105.93750
 -105.31250 -104.68750 -104.06250 -103.43750 -102.81250 -102.18750 -101.56250 -100.93750
 -100.31250  -99.68750  -99.06250  -98.43750  -97.81250  -97.18750  -96.56250  -95.93750
  -95.31250  -94.68750  -94.06250  -93.43750  -92.81250  -92.18750  -91.56250  -90.93750
  -90.31250  -89.68750  -89.06250  -88.43750  -87.81250  -87.18750  -86.56250  -85.93750
  -85.31250  -84.68750  -84.06250  -83.43750  -82.81250  -82.18750  -81.56250  -80.93750
  -80.31250  -79.68750  -79.06250  -78.43750  -77.81250  -77.18750  -76.56250  -75.93750
  -75.31250  -74.68750  -74.06250  -73.43750  -72.81250  -72.18750  -71.56250  -70.93750
  -70.31250  -69.68750  -69.06250  -68.43750  -67.81250  -67.18750  -66.56250  -65.93750
  -65.31250  -64.68750  -64.06250  -63.43750  -62.81250  -62.18750  -61.56250  -60.93750
  -60.31250  -59.68750  -59.06250  -58.43750  -57.81250  -57.18750  -56.56250  -55.93750
  -55.31250  -54.68750  -54.06250  -53.43750  -52.81250  -52.18750  -51.56250  -50.93750
  -50.31250  -49.68750  -49.06250  -48.43750  -47.81250  -47.18750  -46.56250  -45.93750
  -45.31250  -44.68750  -44.06250  -43.43750  -42.81250  -42.18750  -41.56250  -40.93750
  -40.31250  -39.68750
latitude edges
    9.75000   10.25000   10.75000   11.25000   11.75000   12.25000   12.75000   13.25000
   13.75000   14.25000   14.75000   15.25000   15.75000   16.25000   16.75000   17.25000
   17.75000   18.25000   18.75000   19.25000   19.75000   20.25000   20.75000   21.25000
   21.75000   22.25000   22.75000   23.25000   23.75000   24.25000   24.75000   25.25000
   25.75000   26.25000   26.75000   27.25000   27.75000   28.25000   28.75000   29.25000
   29.75000   30.25000   30.75000   31.25000   31.75000   32.25000   32.75000   33.25000
   33.75000   34.25000   34.75000   35.25000   35.75000   36.25000   36.75000   37.25000
   37.75000   38.25000   38.75000   39.25000   39.75000   40.25000   40.75000   41.25000
   41.75000   42.25000   42.75000   43.25000   43.75000   44.25000   44.75000   45.25000
   45.75000   46.25000   46.75000   47.25000   47.75000   48.25000   48.75000   49.25000
   49.75000   50.25000   50.75000   51.25000   51.75000   52.25000   52.75000   53.25000
   53.75000   54.25000   54.75000   55.25000   55.75000   56.25000   56.75000   57.25000
   57.75000   58.25000   58.75000   59.25000   59.75000   60.25000   60.75000   61.25000
   61.75000   62.25000   62.75000   63.25000   63.75000   64.25000   64.75000   65.25000
   65.75000   66.25000   66.75000   67.25000   67.75000   68.25000   68.75000   69.25000
   69.75000   70.25000

gcpy/examples/grids/display_gcclassic_grid_info.py
- Added example script to display information (lons & lats)
  about one of the GEOS-Chem Classic horizontal grids.

gcpy/examples/__init__.py
- Now import all scripts from gcpy/examples/grids folder

CHANGELOG.md
- Updated accordingly

Signed-off-by: Bob Yantosca <yantosca@seas.harvard.edu>
@yantosca yantosca added this to the 1.7.0 milestone Sep 5, 2025
@yantosca yantosca requested a review from msulprizio September 5, 2025 21:56
@yantosca yantosca self-assigned this Sep 5, 2025
@yantosca yantosca added category: Feature Request New feature or request topic: Example Scripts Related to GCPy example scripts labels Sep 5, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Feature Request New feature or request topic: Example Scripts Related to GCPy example scripts
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant