agrif_create_coordinates writes incorrect u,v points to child grid file

[mdunphy]

The agrif_create_coordinates program writes the wrong values for the u and v points when producing a child grid coordinates.nc file. The code assumes a regular lon,lat grid, but that is not the case for many grids, such as ORCA-based grids.

This was probably overlooked because the model doesn't make use of the u,v coordinates (when jphgr_msh=0), rather it loads scale factors from coordinates.nc, and the child-grid scale factors are correct despite the error. Also T and F grid business (bathy, etc) would be unaffected -- so many simulations will not notice this bug. However, there is at least one case where it matters, which is agrif_create_restart that would be slightly affected due to interpolating velocities to the incorrect u,v points.

Christoph Renkl discovered this bug.

For a fix: in agrif_create_coordinates.f90, change

  G1%glamu = G1%glamf
  G1%glamv = G1%glamt

to

  CALL agrif_interp(G0%glamu,G1%glamu,'U')
  CALL agrif_interp(G0%glamv,G1%glamv,'V')

and also change

  G1%gphiu = G1%gphit
  G1%gphiv = G1%gphif

to

  CALL agrif_interp(G0%gphiu,G1%gphiu,'U')
  CALL agrif_interp(G0%gphiv,G1%gphiv,'V')

[jchanut]

In 9632:

Error: Failed to load processor CommitTicketReference

No macro or processor named 'CommitTicketReference' found