[MITgcm-devel] grid rotation

[Martin Losch]

Hi Jean-Michel,

> Just a comment: regarding the the grid-rotation calls, them
> could have been put in ini_grid.F (and then used also with
> curvilinear grids, when x & y are lon & lat), no ?

Is the grid rotation relevant to any other but spherical grids? If so,  
then I agree that it should be moved to ini_grid.F.
Another related issue: we could call calcAngles (angleCS, angleSN)  
online also for curvilinear grids, if they are not read  
( if .not.anglesAreSet at the end of ini_curvilinear_grid ); currently  
they are not require except for non-hydrostatic code on a curvilinear  
grid, but it's nice to have them have the correct values as they are  
always part of the write_grid.F output (and consequently they are  
wrong, i.e. 1 and 0, in most curvilinear grid cases).

Martin