[MITgcm-devel] transport through passages (CS-grid)

[Dimitris Menemenlis]

Alan, I have just checked in MITgcm/utils/matlab/cs_grid/read_cs
some routines that will help speed up reading of CS510 output.  Description and 
examples are in MITgcm/utils/matlab/cs_grid/read_cs/read_cs.txt

In particular, read_cs_face.m can be used to read a single face of the cube, 
e.g., the Arctic face=3, and read_cs_index.m can be used to read specific 
cube-sphere indices, e.g., the ab_IJuv indices from JM's bk_line/grt_circ_bkl.m
routine.  This is a lot faster than reading the whole CS510 field and then 
extracting the face or the indices that you need.

D.

> Hi Alan,
> 
> I checked in the matlab scripts that should allow to compute the
> transport (for CS-grid run) between 2 points, using broken lines.
> 
> under CVS, in MITgcm/utils/matlab/cs_grid/bk_line
> 
> grt_circ_bkl.m write a .mat file (at the top, need to change 
> the coordinate of the 2 points, xa0,ya0 xb0,yb0, and set kwr=1 
> to write this .mat file)
> 
> and if you change krd=2 to krd=3 in use_bk_line.m,
> you will get a short example how to use this broken-line.
> 
> I tried it on cs-32 only. If something is not clear or not
> working on cs-510 grid, just let me know.
> 
> See you,
> Jean-Michel