E3SM-Project / E3SM-Project/scream
Incompatible domain grid coordinate with bi-grid ne1024pg2
- Dominant language
- No language data
- Stars
- 79
- Forks
- 54
- PR merge metrics
- No merged PRs in 30d
Description
A run using grid ne1024pg2_oRRS19to6v3 (available in branch wlin/atm/ne1024pg2_oRRS18to6v3,[ PR #1131 ](https://github.com/E3SM-Project/scream/pull/1131) failed with
```
ERROR: (seq_domain_check_grid) incompatible domain grid coordinates
12273: seq_domain_check_grid - n:482 d1: 231.582371 d2: 231.582371 diff: 0.00000000000142 eps: 0.00000000000100
```
The grids seen by atm and lnd domains have differences from 12th decimal point in lat/lon. Below is backtracing of the code flow from where it was aborted.
```
shr_abort_mod_mp_ 114 shr_abort_mod.F90
seq_domain_mct_mp 693 seq_domain_mct.F90
seq_domain_mct_mp 378 seq_domain_mct.F90
cime_comp_mod_mp_ 1989 cime_comp_mod.F90
```
The differences arise from the two calls in the following if-block
```
if (atm_present .and. lnd_present .and. samegrid_al) then
if (iamroot) write(logunit,F00) ' --- checking atm/land domains ---'
call seq_domain_check_grid(atmdom_a%data, lnddom_a%data, 'lat' , eps=eps_axgrid, mpicom=mpicom_cplid, mask=maskl)
call seq_domain_check_grid(atmdom_a%data, lnddom_a%data, 'lon' , eps=eps_axgrid, mpicom=mpicom_cplid, mask=maskl)
endif
```
It is likely because different versions of grid and/or map files are used during runtime or for generating certain input files. It would be helpful if we know which files are involved to feed the domain/grid data to the above atmdom_a%data and lnddom_a%data.
Contributor guide
No contributing guide indexed for this repository
Assessment
This issue has not been assessed yet.