geoschem / integrated_methane_inversion

Integrated Methane Inversion workflow repository.
https://imi.readthedocs.org
MIT License
29 stars 24 forks source link

[FEATURE REQUEST] Expand IMI to other species such as CO2 #61

Open msulprizio opened 2 years ago

msulprizio commented 2 years ago

As part of JPL's Earth Information System (EIS) there are plans to expand the IMI to include CO2 and additional observation operators.

msulprizio commented 1 year ago

I don't think this should be too difficult to do. To get started, you would need to:

  1. Add an option for species to the config.yml file (e.g. Species: CH4 as the default, but options would also include CO2)
  2. Add new code to these lines in template.sh to create a template run directory for CO2. For example:
    # Commands to feed to createRunDir.sh
    # Create a GEOS-FP 0.25x0.3125 nested NA CH4 run directory by default
    # (Grid and meteorology fields will be replaced below by the settings
    #  in config.yml)
    if ( $Species == "CH4" ); then
       cmd="3\n2\n4\n4\n2\n${RunDirs}\n${runDir}\nn\n"
    elif ( $Species == "CO2" ); then
       cmd="4\n2\n4\n4\n2\n${RunDirs}\n${runDir}\nn\n"
    else
       printf "ERROR: Species ${Species} is not currently supported by the IMI"
       exit 1
    fi

The values for cmd corresponds to the entries passed to createRunDir.sh in GCClassic.

  1. Modify the scripts in the src directory to replace CH4-specific code with options for other species. I recommend grepping the entire directory for CH4 to see where changes need to be made to allow for other species.

Also tagging @laestrada @djvaron @akurganskiy @hannahnesser for input.

djvaron commented 1 year ago

One more thing we'll need is a CO2 observation operator to replace the TROPOMI_operator.py for methane. The operator is called in jacobian.py.

github-actions[bot] commented 8 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 issue from closing this issue.

msulprizio commented 8 months ago

@hannahnesser is currently working on this

msulprizio commented 2 months ago

See draft PR #238