Skip to content

awiesm3 links: correct feom→atm neighbour count and add R→RnfA link - #56

Merged
JanStreffing merged 2 commits into
masterfrom
feat/tco639-next-rmp
Jul 23, 2026
Merged

awiesm3 links: correct feom→atm neighbour count and add R→RnfA link#56
JanStreffing merged 2 commits into
masterfrom
feat/tco639-next-rmp

Conversation

@JanStreffing

Copy link
Copy Markdown
Owner

Follow-up to #55 (which merged the TCO639/NEXT config + weight-gen tooling). These two commits landed on the branch afterwards and carry the actual coupling-correctness fixes, validated against the awiesm3 N56 namcouple + its rmp files:

  • feom→atm GAUSWGT uses 9 neighbours, not 25. The GAUSWGT U map (SST/ice/currents back to the atmosphere) is 1 9 2 in the awiesm3 namcouple. The regenerated rmp_feom_to_A640_GAUSWGT_9.nc (930,508,808 B) matches N56's rmp_feom_to_A640_GAUSWGT.nc (930,820,852 B) to within 0.03%.
  • Added the missing R<nn>→RnfA runoff link (GAUSWGT D 1 25 0.1) — runoff/calving from the atmosphere land grid onto the runoff-mapper grid. awiesm3_feom_links now returns the full 5-link set.
  • The neighbour count is resolution-dependent: 25 is kept for the coarse A096/TCO95 atmosphere (its namcouple value), 9 for higher-resolution atmospheres (A640/TCO639).

Verified end-to-end: full 5-file rmp set generated (cdf2), reordered natural→dist_9600, and the reordered dst_address confirmed a valid bijective relabel with weights byte-preserved.

Match the awiesm3 (N56) namcouple: the feom->atm direction (SST/ice/
currents) uses GAUSWGT U with 9 neighbours (1 9 2), not 25; and the
runoff/calving atm->mapper link (R<nn> -> RnfA, GAUSWGT D 1 25 0.1)
was missing. awiesm3_feom_links now returns the full 5-link set.
Add submit_weightgen_..._regen.sh to regenerate just the two changed
links onto an existing weight set.
The 9-neighbour default (90b1eec) is right for high-res atmospheres (TCO639) but
broke TCO95/A096, whose deployed namcouple uses GAUSWGT U ... 25. A regen with 9
emitted rmp_feom_to_A096_GAUSWGT_9.nc, so OASIS fell back to a stale full-mesh
_25 -> MCT column-count mismatch (208296 vs 211567) at the mesh increment.

Key feom->atm on atm_grid: A096 -> 25 neighbours, else -> 9 (default).
@JanStreffing
JanStreffing merged commit 790b8a0 into master Jul 23, 2026
6 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant