E3SM-Project / E3SM-Project/scream

Incompatible domain grid coordinate with bi-grid ne1024pg2

Open
#1,140 6 comments 0 reactions 0 assignees View on GitHub
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.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.