search

geoschem / CAM

Community Atmosphere Model
Other
0 stars 0 forks source link

photolysis_init does not function on ne30 compsets #36

Open jimmielin opened 3 days ago

jimmielin commented 3 days ago

because CAM_PIO_OPENFILE is a collective operation and on SE grids the assumption that all CPUs have the same # of chunks does not work.

this works for drydep_mod because CAM_PIO_OPENFILE is called by each CPU and not per chunk, but for photolysis_init it is called by chunk because it populates fields in State_Chm.

A proper fix to make GEOS-Chem aware of root chunk information needs to be made. I will prepare a pull request to update State_Grid for that.

There are likely more bugs associated with the SE grid runs with GEOS-Chem chemistry. Will update later

jimmielin commented 3 days ago

a fix is available: patches to CAM (https://github.com/jimmielin/CAM/commit/cc5642fa61e0bcf0f647546c2b398bc449c4de0c), patches to GEOS-Chem (https://github.com/jimmielin/geos-chem/commit/4176034645808fdbdaa0a8f15dacc0aab1024f9b)

PR to GEOS-Chem science repo: https://github.com/geoschem/geos-chem/pull/2360

includes a new State_Grid%ROOT_ID field to store information about BEGCHUNK.

these are based off 14.3, but should be cherry-picked into 14.4.