integer function RPN_COMM_topo(nxg,minx,maxx,nxl,nxlmax,
     %                   halox,nx0,alongx,fill)

ARGUMENTS

• nxg : number of points of global domain along one axis
• halox : size of halo
• alongx : logical, true if the domain is divided along x
• fill : logical, add one more point to maxx if nxlmax is even (for memory optimization purpose)

OUT
 

minx maxx : dimensions of the local array needed to contain the local tile

nxl : number of local points along axis for this PE

nxlmax : maximum of local points along axis for all PE

nx0 : Position of the first element of the local PE over nxg

DESCRIPTION

For example, a 3x2 PE grid, 23x12 points, halo of 1 along x axis:

+-+-+-+
|3|4|5|
+-+-+-+
|0|1|2|
+-+-+-+

ierr = RPN_COMM_topo(23,minx,maxx,nxl,nxlmax,1,nx0,.true.,.false.)

minx = 0
maxx = 10
nxl = 8 on PE 0,1,3,4, and 7 on PE 2 and 5.
nxlmax = 8
nx0 = 1 on PE 0 and 3, 9 on 1 and 4, and 17 on 2 and 5.

EXAMPLE

SEE ALSO