-
Notifications
You must be signed in to change notification settings - Fork 27
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Output Emissions collection with lev positive up #115
Comments
This feature request is somewhat related to: |
If feature request GEOS-ESM/MAPL#914 is implemented as proposed then no work would actually need to be done for this feature request beyond changing |
Specifying lev positive up or down in history files is now possible in MAPL 2.8.3. When this will go into GCHP is still TBD. |
Hi @lizziel --Just checking in to make sure I've understood the relevant threads to this correctly (followed issues #112 and #10 here): if I'm using v 13.2.1, GEOSChem.Emissions collection output will have the surface at lev = 72, and toa at lev = 1; whereas all other collections (eg: SpeciesConc, AerosolMass, etc) have surface at lev = 1 and toa at lev = 72? |
Yes, that is correct! The restart file is also flipped relative to the other diagnostics. |
This update is being moved to the 14.1 release. |
I wonder if we should push this to 15.0 since it will break GCPy and maybe is more appropriate for a major release. |
GCHP has always output the emissions diagnostic collection with levels inverted relative to all other diagnostic collections. Most emissions are 2D so this is not relevant for most emission outputs. However, for 3D emissions the convention difference is a pain for post-processing. It also can be a source of error if users are not aware of the difference from other diagnostics.
The lev positive down convention for HEMCO diagnostics originated in its use in GEOS which always stores and outputs variables with that convention. Changing it to lev positive up should be a straightforward change. The existing code for lev positive down should be kept but put into a MODEL_GEOS cpp ifdef.
The text was updated successfully, but these errors were encountered: