[MITgcm-devel] Re: TARGET_NEC_SX in ini_curvilinear_grid.F

[Jean-Michel Campin]

Hi Martin,

On Mon, Feb 16, 2009 at 09:42:26AM +0100, Martin Losch wrote:
> Hi Jean-Michel,
> both you suggestion in this email and the fix that you checked in work, 
> so I guess now this part is OK. Thanks.

Good thing to know (since we already remove your fix).

> However, I have another problem (that I almost forgot about). In my 2cpu 
> configuration (nPx=1,nPy=2,nSx=1,nSy=1) without EXCH2 it appears that 
> only tile001.mitgrid is read:
>> sx8::test> grep tile00* STDOUT.*
>> STDOUT.0000:(PID.TID 0000.0001) tile:   1 ; Read from file  
>> tile001.mitgrid
>> STDOUT.0001:(PID.TID 0001.0001) tile:   1 ; Read from file  
>> tile001.mitgrid
> with the obvious consequences for the grid. (np2=1 and nsy=2 does not  
> work either). It looks like the code is set up to read only cubed sphere 
> faces (or domains with tiles in X only). Maybe it's a good idea to catch 
> that case (#ifndef OLD_GRID_IO and if (sny*npy>1) then STOP)?
>
> Martin

We can change it to make it work (with 1 grid-file per tile). 
It's just that we need to agree on how the tile & grid-file name will be, 
since we only have 1 number for both Bi,Bj (<- I mean global tile indices, 
in (nSx*nPx,nSy*nPy)),
and I would like to keep the same file name shape: 
 tile[3_digit_number].mitgrid

The easier would be:
 tile_number = Bi + (Bj-1)*(nSx*nPx)
so that it goes from 1:nSx*nPx for Bj=1 
then       1+nSx*nPx:2*nSx*nPx for Bj=2
and so on.

Would this be usefull ? Or do you prefer a "STOP" ?

Cheers,
Jean-Michel