| 1 |
In general: the overlap corners are perillious for any field located on tracer points. |
| 2 |
They are never actively used in the model computation so the values there |
| 3 |
can be meaningless in forward mode or in adjoint mode -- only the exchange |
| 4 |
can make a difference there. |
| 5 |
Using a regular lat/lon grid and exch "1", you find the expected value |
| 6 |
(corner of another tile) but this cannot be expected in general. |
| 7 |
In the Cube Sphere topology for example, the overlap corner values simply |
| 8 |
do not correspond to a model grid point!!! The exch2 routines then put a |
| 9 |
default value there or copy a neighboring value. This is true for XC and |
| 10 |
YC after they have been exchanged for example. |
| 11 |
=> we are never going to use the values that are found in the overlap corners, |
| 12 |
either for XC/YC or for T/S |