...
Found the root cause of the “speckling” problem in the ne256 trigrid setup reported here: /wiki/spaces/NGDNA/pages/3597565953 This problem is due to the land->atmosphere state map being inconsistent with the land franction on the atmosphere grid. In particular, the map does not provide data in a few grids points where the land fraction is .30. At such a point, the atmosphere field is 30% land + 70% ocean, and since land will send no data, this results in temperature at a single point being consistently 30% too low. This is reported in the metric for “land frac consistency error” given in the page for the lnd2atm state map.
Uses the python “holoviews” package which can directly plot unstructured grid data (by shading each polygon cell) manytimes faster than matplotlib’s ‘polycollection’ approach.
The E3SM v3 trigrid setups have a small (but larger than expected) discrepancy between the land fraction computed from ocn->atm (computed from ocn2atm map) and ocn->lnd->atm , suggesting the use of (computed from domain.lnd + lnd2atm map). The metric for this is the “land frac consistency error” given in the page for the lnd2atm flux map. This error suggests that a mixture of different tools to compute ave maps and domain files, or an inconsistency in the TR aave tool. (The older ne256 trigrid setup has much smaller errors for this metric and it’s flux maps may have been computed with nco or esmf).
Grid | Report | Notes | ||||
---|---|---|---|---|---|---|
ne30pg2_r05_IcoswISC30E3r5 |
| E3SM v3 LR | ||||
northamericax4v1pg2_r025_IcoswISC30E3r5 |
| E3SM v3 NA-RRM | ||||
ne120pg2_r025_RRSwISC6to18E3r4 |
| E3SM v3 HR | ||||
ne256pg2_r0125_oRRS18to6v3 |
| Identifies problem with LND->ATM state map responsible for speckling issues | ||||
ne1024pg2_oRRS18to6v3 |
| Showcase holoviews ability to plot unstructured UHR data. Requires 34GB of memory and takes about 20min on an Anvil node to produce all the plots. | ||||
...