[MITgcm-support] (Fwd) problems reading theta.bin

[Jean-Michel Campin]

Hi Jane,

I just have a quick look at your email, and it seems that there is
something inconsistent between:
nx=200 (gendata.m) or  delX=200*5.e3 (data)
and 
sNx = 200 with nSx =   2 (SIZE.h) that give you 2*200=400 grid points
in the X direction (more precisely, 2 tiles (nSx=2) of 200 grid points each).

The things that I see in the output confirm my interpretation:

(PID.TID 0000.0001)      nPx =    1 ; /* No. processes in X */
(PID.TID 0000.0001)      nPy =    1 ; /* No. processes in Y */
(PID.TID 0000.0001)      nSx =    2 ; /* No. tiles in X per process */
(PID.TID 0000.0001)      nSy =    1 ; /* No. tiles in Y per process */
(PID.TID 0000.0001)      sNx =  200 ; /* Tile size in X */
(PID.TID 0000.0001)      sNy =   50 ; /* Tile size in Y */
(PID.TID 0000.0001)      OLx =    2 ; /* Tile overlap distance in X */
(PID.TID 0000.0001)      OLy =    2 ; /* Tile overlap distance in Y */
(PID.TID 0000.0001)      nTx =    1 ; /* No. threads in X per process */
(PID.TID 0000.0001)      nTy =    1 ; /* No. threads in Y per process */
(PID.TID 0000.0001)       Nr =   50 ; /* No. levels in the vertical   */
(PID.TID 0000.0001)       nX =  400 ; /* Total domain size in X ( = nPx*nSx*sNx ) */

I don't know if you want 400 or 200 grid points in the X direction,
but I know that the numbers in SIZE.h and data & gendata.m should
match.

For a practical point of view, to run on 1 cpu, you don't need to
use 2 tiles (nSx=2), since it's slightly less efficient that 1 big tile (nSx=1
and sNx= 400 or 200, depending on what you need).

Hope this will be helpful to you.

Jean-Michel