search

geoschem / gcpy

Python toolkit for GEOS-Chem. Contains basic plotting scripts, plus the suite of GEOS-Chem benchmarking utilities.
https://gcpy.readthedocs.io
Other
50 stars 24 forks source link

Question about GCPY functions #194

Closed TylerElgiar closed 1 year ago

TylerElgiar commented 1 year ago

I have been running the GEOSChem model simulating atmospheric mercury chemistry over the continental united states. I am trying to produce a couple different surface maps for my runs. I am trying to sum together all of the different oxidized species of mercury to produce one map of the summed species. Is this something easily done in GCPY? If so, how do I do this?

yantosca commented 1 year ago

Thanks for writing @TylerElgiar. GCPy does not have any plotting scripts for Hg simulations yet. I know that Ari Feinberg (@arifein) had some Hg benchmarking scripts that may (or may not) do what you want so I would check with him first.

We do plan on bringing the Hg benchmarking scripts into GCPy at some point but this hasn't been done yet.

arifein commented 1 year ago

Hi @TylerElgiar. I have a git repository for the python benchmarking scripts for GC-Hg here: https://github.com/arifein/pythonHgBenchmark. You can see what the output plots look like in the Figures/ directory of that repository. My strategy for the many different oxidized Hg species is to sum them up in a bash script that postprocesses the GC output with cdo/nco commands. I do this also for the oxidized Hg fluxes, to get overall values (which might be useful for you). The postprocessing script for v14 of GC is here: postproc_runs_Hg_v14.sh. Happy to get any issues/updates that you find for these benchmarking scripts.

TylerElgiar commented 1 year ago

Thank you @arifein. Your script was helpful in summing together HgII species for my maps.