search

geoschem / geos-chem

GEOS-Chem "Science Codebase" repository. Contains GEOS-Chem science routines, run directory generation scripts, and interface code. This repository is used as a submodule within the GCClassic and GCHP wrappers, as well as in other modeling contexts (external ESMs).
http://geos-chem.org
Other
167 stars 160 forks source link

EDGARv7 and EDGARv6_SF files are incompatible with GCHP (plus other files) #2125

Closed yantosca closed 6 months ago

yantosca commented 8 months ago

Name and Institution (Required)

Name: Bob Yantosca Institution: Harvard + GCST

Confirm you have reviewed the following documentation

Description of your issue or question

The HEMCO/CH4/v2023-04/EDGARv7 and HEMCO/CH4/v2023-04/EDGARv6_SF files have several incompatibilities for GCHP. Here is the isCoards output:

EDGARv7

===========================================================================
Filename: v7.0_FT2021_CH4_2019_AGS.0.1x0.1.nc
===========================================================================

The following items adhere to the COARDS standard:
---------------------------------------------------------------------------
-> Dimension "time" adheres to standard usage
-> Dimension "lat" adheres to standard usage
-> Dimension "lon" adheres to standard usage
-> Dimension "time" uses only lowercase letters (requirement for GCHP)
-> time(time)
-> time is a single value
-> time[0] = 0 (requirement for GCHP)
-> time:calendar = "standard" 
-> time:long_name = "Time" 
-> time:units = "hours since 2019-01-01 00:00:00" 
-> "time:units" does not contain decimals (requirement for GCHP)
-> Dimension "lat" uses only lowercase letters (requirement for GCHP)
-> lat(lat)
-> lat is monotonically increasing
-> lat:long_name = "latitude" 
-> lat:units = "degrees_north" 
-> Dimension "lon" uses only lowercase letters (requirement for GCHP)
-> lon(lon)
-> lon is monotonically increasing
-> lon:long_name = "longitude" 
-> lon:units = "degrees_east" 
-> emi_ch4(time,lat,lon)
-> emi_ch4 has proper dimension ordering (requirement for GCHP)
-> emi_ch4:long_name = "Emissions of CH4 - " 
-> emi_ch4:units = "kg m-2 s-1" 
-> conventions: "CF-1.0" 
-> history: "Sun Apr 23 16:12:29 2023: ncecat -u time tmp.nc out.nc\n",
-> title: "Annual Mean Emissions of CH4 - " 
-> references: "European Commission, Joint Research Centre (JRC)/Netherlands Environmental Assessment Agency (PBL). Emission Database for Global Atmospheric Research (EDGAR),http://edgar.jrc.ec.europe.eu" 

The following items DO NOT ADHERE to the COARDS standard:
---------------------------------------------------------------------------
Variable "time" is of type "int" (problem for GCHP)
-> = needs to be declared with (seehttp://edgar.jrc.ec.europa.eu/methodology.php#12souforthedefinitionsofthesinglesources")
-> = has improper dimension ordering (problem for GCHP)
-> =:long_name (or =:standard_name) is missing
-> =:units is missing
-> = needs to be declared with (seehttp://edgar.jrc.ec.europa.eu/methodology.php#12souforthedefinitionsofthesinglesources")
-> = has improper dimension ordering (problem for GCHP)
-> =:long_name (or =:standard_name) is missing
-> =:units is missing

The following optional items are RECOMMENDED:
---------------------------------------------------------------------------
-> Consider adding time:axis = "T"
-> Consider adding lat:axis = "Y"
-> Consider adding lon:axis ="X"
-> Consider adding emi_ch4:_FillValue
-> Consider adding emi_ch4:missing_value
-> Consider adding =:_FillValue
-> Consider adding =:missing_value
-> Consider adding =:_FillValue
-> Consider adding =:missing_value
-> Consider adding the "format" global attribute

---------------------------------------------------------------------------
For detailed instructions on how to edit netCDF files, see:
 1. "Work with netCDF files" and
 2. "Prepare COARDS-compliant netCDF files"
at the GEOS-Chem manual (geos-chem.readthedocs.io).
---------------------------------------------------------------------------

EDGARv6_CH4

===========================================================================
Filename: EDGARv6_CH4_MonthlyScaleFactors_AGS.0.1x0.1.nc
===========================================================================

The following items adhere to the COARDS standard:
---------------------------------------------------------------------------
-> Dimension "time" adheres to standard usage
-> Dimension "lat" adheres to standard usage
-> Dimension "lon" adheres to standard usage
-> Dimension "time" uses only lowercase letters (requirement for GCHP)
-> time(time)
-> time is monotonically increasing
-> time:axis = "T" 
-> time:calendar = "standard" 
-> time:long_name = "Time" 
-> time:units = "hours since 1970-01-01 00:00:00" 
-> "time:units" does not contain decimals (requirement for GCHP)
-> Dimension "lon" uses only lowercase letters (requirement for GCHP)
-> lon(lon)
-> lon is monotonically increasing
-> lon:axis = "X" 
-> lon:long_name = "longitude" 
-> lon:units = "degrees_east" 
-> Dimension "lat" uses only lowercase letters (requirement for GCHP)
-> lat(lat)
-> lat is monotonically increasing
-> lat:axis = "Y" 
-> lat:long_name = "latitude" 
-> lat:units = "degrees_north" 
-> sf_ch4(time,lat,lon)
-> sf_ch4 has proper dimension ordering (requirement for GCHP)
-> sf_ch4:long_name = "CH4 monthly scale factor" 
-> sf_ch4:units = "1" 
-> conventions: "CF-1.0" 
-> history: "Sun Apr 23 15:55:21 2023: ncatted -a comment,sf_ch4,d,, EDGARv6_CH4_MonthlyScaleFactors_AGS.0.1x0.1.nc\n",
-> title: "Monthly Mean (Jan) Emissions of CH4 - " 
-> references: "European Commission, Joint Research Centre (JRC)/Netherlands Environmental Assessment Agency (PBL). Emission Database for Global Atmospheric Research (EDGAR),http://edgar.jrc.ec.europe.eu" 

The following items DO NOT ADHERE to the COARDS standard:
---------------------------------------------------------------------------
-> time[0] != 0 (problem for GCHP)

The following optional items are RECOMMENDED:
---------------------------------------------------------------------------
-> Consider adding sf_ch4:_FillValue
-> Consider adding sf_ch4:missing_value
-> Consider adding the "format" global attribute

---------------------------------------------------------------------------
For detailed instructions on how to edit netCDF files, see:
 1. "Work with netCDF files" and
 2. "Prepare COARDS-compliant netCDF files"
at the GEOS-Chem manual (geos-chem.readthedocs.io).
---------------------------------------------------------------------------
yantosca commented 8 months ago

Tagging @msulprizio

yantosca commented 8 months ago

Also found that the HEMCO/CH4/v2023-04/GC_CH4_LOSS/GCC14_72LM.ch4loss.4x5.nc4 file is also not compatible with GCHP:

===========================================================================
Filename: ./HcoDir/CH4/v2023-04/GC_CH4_LOSS/GCC14_72LM.ch4loss.4x5.nc4
===========================================================================

The following items adhere to the COARDS standard:
---------------------------------------------------------------------------
-> Dimension "time" adheres to standard usage
-> Dimension "lev" adheres to standard usage
-> Dimension "lat" adheres to standard usage
-> Dimension "lon" adheres to standard usage
-> Dimension "ilev" adheres to standard usage

The following items DO NOT ADHERE to the COARDS standard:
---------------------------------------------------------------------------
-> Dimension "nb" is non-standard usage,
   which will cause GCHP to fail during file read.
   PLEASE REMOVE THIS DIMENSION AND ALL VARIABLES THAT USE IT!!!
yantosca commented 8 months ago

Update: I have fixed the gmi.clim.CH4.geos5.2x25.nc file with the netCDF Operators ncks command:

$ ncks -x -v lon_bnds,lat_bnds,nb GCC14_72LM.ch4loss.4x5.nc4 tmp.nc
$ mv tmp.nc GCC14_72LM.ch4loss.4x5.nc4

GCHP also does not like the following variables (which do not have time,lev,lat,lon ordering), so we have removed these as well.

$ ncks -x -v hyai,hybi,ilev,AREA,hyam,hybm  GCC14_72LM.ch4loss.4x5.nc4 tmp.nc
$ mv tmp.nc GCC14_72LM.ch4loss.4x5.nc4
yantosca commented 8 months ago

The JPL WetCharts emissions file is missing the emi_ch4:long_name attribute, which is required by GCHP.

===========================================================================
Filename: ./HcoDir/CH4/v2020-09/JPL_WetCharts/HEensemble/JPL_WetCharts_2010-2019.Ensemble_Mean.0.5x0.5.nc
===========================================================================

The following items adhere to the COARDS standard:
---------------------------------------------------------------------------
... etc not shown...

The following items DO NOT ADHERE to the COARDS standard:
---------------------------------------------------------------------------
-> emi_ch4:long_name (or emi_ch4:standard_name) is missing
-> The "title" global attribute is missing

The following optional items are RECOMMENDED:
---------------------------------------------------------------------------
-> Consider adding emi_ch4:_FillValue
-> Consider adding emi_ch4:missing_value
-> Consider adding the "format" global attribute
stale[bot] commented 6 months ago

This issue has been automatically marked as stale because it has not had recent activity. If there are no updates within 7 days it will be closed. You can add the "never stale" tag to prevent the Stale bot from closing this issue.

yantosca commented 6 months ago

This issue is moot since EDGARv8 will be used as the default inventory for the CH4/carbon simulations.