GEOS-ESM / GEOS-ESM/GEOSgcm_GridComp

Out of bounds error when running stochastic physics

Open
#314 0 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

bug
Dominant language
Fortran
Stars
10
Forks
13
Avg merge
17h 35m
Merged PRs (30d)
14

Description

While testing something for issue #311 I tried running the stochastic physics in a build that was built with the debug options. I set SKEB: 1 and SPPT: 1 in the AGMC.rc which seems to be the flags to turn this on. When running with debugging the code crashed like so with an out of bounds error:

GEOSgcm.x          000000000198A11C  interp_h_                 369  interp_h.F
GEOSgcm.x          0000000001861428  gauss2grid_                80  gauss2grd.f90
GEOSgcm.x          0000000001858DF9  stoch_module_mp_g         311  stoch_module.F90
GEOSgcm.x          000000000184E44C  stoch_module_mp_s         140  stoch_module.F90
GEOSgcm.x          00000000017B6DC2  geos_physicsgridc        2650  GEOS_PhysicsGridComp.F90

So something is not right. When running with the release build it does not crash.

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Reproduce the debug build with SKEB: 1 and SPPT: 1 in AGMC.rc, then inspect the reported locations: interp_h.F line 369, gauss2grd.f90 line 80, and stoch_module.F90 lines 140 and 311. Trace the out-of-bounds access and verify that stochastic physics completes without the debug crash.

Written by the indexing model from the issue text.

Assessment

Tech stack
fortran
Domain
hpc
Issue type
Bug
Difficulty
3/5
Estimated time
1-2 days
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
35/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.