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 163 forks source link

The DryDep collection is not activated for GCHP fullchem benchmarks #2139

Closed yantosca closed 5 months ago

yantosca commented 9 months ago

Name and Institution (Required)

Name: Bob Yantosca Institution: Harvard + GCST

Confirm you have reviewed the following documentation

New GEOS-Chem feature or discussion

The DryDep diagnostic is commented out by default in the GCHP fullchem benchmark simulation. I believe this had been done due to avoid an ESMF "Input file too long" error. However, we need the DryDep diagnostic output to plot the Ox budget and drydep velocities.

Potential fix: Activate the DryDep diagnostic, but only archive the fields that are needed.

Tagging @msulprizio @lizziel

yantosca commented 9 months ago

Activating the DryDep collection (with other collections) caused the problem in #2135. I am testing if just activating the DryDep collection will allow GCHP benchmark simulations to run.

yantosca commented 9 months ago

Using a reduced definition for the DryDep collection gets past the error. I think this is due to some leftover species which may have been removed in the GCHP HISTORY.rc file. Using this definition for the DryDep collection works:

```console
#==============================================================================
# %%%%% THE DryDep COLLECTION %%%%%
#
# Dry depositon fluxes and velocities (for all dry-depositing species)
#
# Available for all simuations that have dry-depositing species
#==============================================================================
  DryDeptemplate:        '%y4%m2%d2_%h2%n2z.nc4',
  DryDep.format:          'CFIO',
  DryDep.timestampStart:  .true.
  DryDep.monthly:         0
  DryDep.frequency:       010000
  DryDep.duration:        010000
  DryDep.mode:            'time-averaged'
  DryDep.fields:          'DryDepVel_ACET       ', 'GCHPchem',
                          'DryDepVel_HNO3       ', 'GCHPchem',
                      'DryDepVel_NH3        ', 'GCHPchem',
                          'DryDepVel_NH4        ', 'GCHPchem',
                          'DryDepVel_NIT        ', 'GCHPchem',
                          'DryDepVel_NITs       ', 'GCHPchem',
                          'DryDepVel_O3         ', 'GCHPchem',
                          'DryDepVel_SO4        ', 'GCHPchem',
                          'DryDep_O3            ', 'GCHPchem',
::
lizziel commented 9 months ago

I can try to improve the error handling for this type of problem.

yantosca commented 9 months ago

@lizziel: I've opened draft PR #2144 where I will test this. This can go into 14.4.0.

stale[bot] commented 7 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 5 months ago

PR #2289 has been merged into GEOS-Chem 14.4.0. We can now close out this issue.