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
168 stars 165 forks source link

[BUG/ISSUE] Number of ND64 tracers gets set to zero in input_mod.F #67

Closed yantosca closed 5 years ago

yantosca commented 5 years ago

Describe the bug The ND64 bpch radiative flux diagnostics are not printed out, even when ND64 is turned on in input.geos.

To Reproduce

  1. Create a geosfp_4x5_standard run dir.
  2. make -j8 all
  3. make -j8 run

Expected behavior The ND64 tracers should be saved to the bpch file, and metadata printed to tracerinfo.dat.

Required information

Input and log files to attach

---------------
Diag    L   Tracers being saved to disk
ND05   72   1 - 20
ND06    1   1 -  4
ND07   72   1 - 15
ND08    1   1 -  4
ND11    1   1 -  5
ND13   72   1 -  1
ND21   72   1 - 60
ND22   72   1 - 78
ND31   73   1 -  1
ND32   72   1 -  1
ND38   72   1 -162
ND39   72   1 -162
ND42   72   1 - 25
ND43   72   1 -  4
ND45   72   1 -162
ND54   72   1 -  1
ND55   72   1 -  3
ND64   72   1  - 0  **NOTE: Should be 1 - 3 ***
ND66   72   1 -  7
ND67   72   1 - 23
ND68   72   1 -  8
yantosca commented 5 years ago

The following commits should now resolve this issue: https://github.com/geoschem/geos-chem/commit/872c3a04ef3a55e4a5480a96a931658a5c79b54d and https://github.com/geoschem/geos-chem/commit/114193f053fe024caa0896f263c1cc61cf35cbc5. These are currently in the bugfix/RRTMG branch, which will be merged into GEOS-Chem 12.7.0.

This now gives plots that look like this, that follow the day-night pattern: test