[MITgcm-support] cube face halo

Chris Hill cnh at mit.edu
Wed Jan 28 17:30:58 EST 2004


 Benny,

 Would it help to use the model to write out the tiled XC, YC etc... with
their overlaps, that's a fairly easy option to activate. Then you would have
the numbers you need. Is Dimitris back yet - he knows how to do this.

Chris

> -----Original Message-----
> From: Benny Cheng [mailto:bcheng at jpl.nasa.gov] 
> Sent: Wednesday, January 28, 2004 1:02 PM
> To: Chris Hill
> Cc: mitgcm-support at mitgcm.org; bnc at pacific.jpl.nasa.gov; 
> adcroft at mit.edu
> Subject: RE: [MITgcm-support] cube face halo
> 
> 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