[MITgcm-support] MNC grid boundaries with curvilinear coordinates

Martin Losch Martin.Losch at awi.de
Mon Feb 25 04:09:16 EST 2013


David,

it's great to see python output in this matlab-dominated community (o:

It's not a bug, it's a feature: the default boundary conditions for the MITgcm are double periodic so that the point at i=nx+1 is identical to the point at i=1, and the same for the j-direction. In the mnc-package, the output for variable on u and v-points of each tile is augmented by the first column/row of the neighboring tile, as a service (o: The last tile in x get's the first column of the first tile in x, and the coordinate it necessarily XG(1) = XG(end) - Lx.

If your domain it periodic, that's what you want, right? If not, you can just ignore the last column. In fact, the standard mds output (the *.data/meta files) do not have these extra columns for any of the variables, neither for XG/YG.

Martin

On Feb 22, 2013, at 5:21 PM, David Huard <david.huard at gmail.com> wrote:

> Hi all, 
> 
> I'm using the MITgcm on the Arctic domain and I noticed an issue with the netcdf output grid variables XG and YG. 
> 
> The last row and column of XG and YG are just copies the first row and column, as if it was a pac-man domain. Here is an example: 
> 
> array([[  31.97116089,   32.10241318,   32.23471451, ...,  146.66885376,
>          146.80964661,   31.97116089],
>        [  31.77091026,   31.90175247,   32.033638  , ...,  146.87324524,
>          147.01361084,   31.77091026],
>        [  31.57006645,   31.70048714,   31.83195305, ...,  147.07829285,
>          147.21824646,   31.57006645],
>        ..., 
>        [ -48.62563324,  -48.77060318,  -48.91656113, ..., -129.88247681,
>         -130.03565979,  -48.62563324],
>        [ -48.7310524 ,  -48.87626648,  -49.02246475, ..., -129.77496338,
>         -129.92831421,  -48.7310524 ],
>        [  31.97116089,   32.10241318,   32.23471451, ...,  146.66885376,
>          146.80964661,   31.97116089]])
> 
> 
> The consequence is that the last row and column of the tracer coordinates (XC, YC) do not lie in between the surrounding grid coordinates, which is what I would be expecting. 
> 
> Note that I generated the full grid using the python gluemncbig script. 
> 
> Cheers, 
> 
> David Huard
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list