[MITgcm-support] cube face halo

Benny Cheng bcheng at jpl.nasa.gov
Wed Jan 28 13:04:21 EST 2004


yes, I'm trying to interpolate u,v from the cube to the regular lat,long
u,v, like 1/4 deg model grid. Looks like I'll need the lat and long
values inside the halo regions to do a proper interpolation of u,v along
the edges of the cube. I suppose I could get them from the STDOUT.*
files, but I need some help in decoding the output there... ,


e.g. for
STDOUT.0000, I see tile 1 exchanges info with tile 180 and tile 211 ,
but where is i=-424:-422  j=513: 423  on tile 1 ? isn't tile 1 from
i=-2:513, j=-2:513 ? 

also, for the u grid, is it the case then that xC,yG on tile 1 matches
exactly with xG,yC on tile 180 ? 



- Benny


(PID.TID 0000.0001)       NEIGHBOUR    1 = TILE    7 Comm = MSG ( PROC
=    7)
(PID.TID 0000.0001)       NEIGHBOUR    2 = TILE  211 Comm = MSG ( PROC
=  211)
(PID.TID 0000.0001)       NEIGHBOUR    3 = TILE    2 Comm = MSG ( PROC
=    2)
(PID.TID 0000.0001)       NEIGHBOUR    4 = TILE  180 Comm = MSG ( PROC
=  180)
(PID.TID 0000.0001) Tile    1(proc =   1) sends points i=  -2:  88, j= 
-2:   0
(PID.TID 0000.0001)                          to points i=  -2:  88, j= 
-2:   0 in tile    7(proc =   7)
(PID.TID 0000.0001) Tile    1(proc =   1) sends points i=  -2:  88,
j=-424:-422
(PID.TID 0000.0001)                          to points i=  -2:  88, j= 
86:  88 in tile  211(proc = 211)
(PID.TID 0000.0001) Tile    1(proc =   1) sends points i=  -2:   0, j= 
-2:  88
(PID.TID 0000.0001)                          to points i=  -2:   0, j= 
-2:  88 in tile    2(proc =   2)
(PID.TID 0000.0001) Tile    1(proc =   1) sends points i=-424:-422, j=
513: 423
(PID.TID 0000.0001)                          to points i=  -2:  88, j= 
86:  88 in tile  180(proc = 180)
(PID.TID 0000.0001) Tile    1(proc =   7)recv to points i=  -1:  87, j= 
86:  88from tile   7(proc =   7)
(PID.TID 0000.0001) Tile    1(proc = 211)recv to points i=  -2:  88, j= 
-2:   0from tile 211(proc = 211)
(PID.TID 0000.0001) Tile    1(proc =   2)recv to points i=  86:  88, j= 
-2:  88from tile   2(proc =   2)
(PID.TID 0000.0001) Tile    1(proc = 180)recv to points i=  -2:   0, j= 
-2:  88from tile 180(proc = 180)

On Wed, 2004-01-28 at 08:32, Chris Hill wrote:
> Benny,
> 
>  Can you give us a bit more of a clue as to what you need to do? 
> 
> Chris
> 
> > -----Original Message-----
> > From: mitgcm-support-bounces at mitgcm.org 
> > [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of 
> > Alistair Adcroft
> > Sent: Wednesday, January 28, 2004 10:50 AM
> > To: mitgcm-support at mitgcm.org
> > Subject: RE: [MITgcm-support] cube face halo
> > 
> > Not sure what you mean - the halo generally contains the data 
> > from the next tile and on the cube the only difference is the 
> > diagaonl corner regions of the halo are not updated. The cube 
> > exchange routine that should be straight forward to read are 
> > in eesupp/src/*cube*.F .
> > 
> > Perhaps I misunderstand the question?
> > 
> > A.
> > --
> > Dr Alistair Adcroft            http://www.mit.edu/~adcroft
> > MIT Climate Modeling Initiative        tel: (617) 253-5938
> > EAPS 54-1523,  77 Massachusetts Ave,  Cambridge,  MA,  USA
> > 
> > -----Original Message-----
> > From: mitgcm-support-bounces at mitgcm.org
> > [mailto:mitgcm-support-bounces at mitgcm.org] On Behalf Of Benny Cheng
> > Sent: Tuesday, January 27, 2004 7:47 PM
> > To: mitgcm-support at mitgcm.org
> > Subject: [MITgcm-support] cube face halo
> > 
> > 
> > 
> > Is there a quick way to get the halo information surrounding 
> > a FACE of the cube for the cube-sphere code without running 
> > the model? basically, need to know how the information along 
> > the edges of the cube are exchanged...
> > 
> >  - Benny
> > 
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> > 
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://dev.mitgcm.org/mailman/listinfo/mitgcm-support
> > 
> 
> 




More information about the MITgcm-support mailing list