1 
# $Header: /u/gcmpack/MITgcm/utils/cs_grid/README,v 1.3 2005/09/06 17:38:06 jmc Exp $ 
2 
# $Name: $ 
3 

4 
Some tools (matlab scripts for now) specific to the CubedSphere grid : 
5 

6 
split_C_cub.m :: split a 2d,3d global CSfield into 6 faces ; 
7 
apply to cellcentered CSfield ; extend the face array by 
8 
adding 1 row & 1 column (like the exchanges) 
9 
split_Z_cub.m :: split a 2d,3d global CSfield into 6 faces ; 
10 
apply to cellcorner CSfield ; extend the face array by 
11 
adding 1 row & 1 column (like the exchanges) 
12 
split_UV_cub.m :: split the 2 Cgrid components of a vector 
13 
(2d,3d global CSfield) into 6 faces ; extend the face array 
14 
by adding 1 row or 1 column (like the exchanges) 
15 

16 
cubeZ2latlon.m :: interpolate a CSfield to regular Latlon grid ; 
17 
similar to cube2latlon.m but apply to a field that is 
18 
located at the corner point of the Cgrid cell. 
19 
(deal with the 2 "missing corners") 
20 

21 
calc_vort_cs.m :: compute vorticity (vertical component) on CSgrid 
22 
from MITgcm output velocity fields u,c (on Cgrid); 
23 
Use the model discretization of vorticity. 
24 

25 
the following scripts are reading a file that contains the orientation 
26 
of the local direction of the CSgrid relative to EW,NS direction. 
27 
file= proj_cs[]_2uEvN.bin :: cos & sin of rotation angle, 
28 
format:: real*8, bigendian, size=32x6x32 x2(cos+sin) 
29 
for standard cs32 grid, the corresponding file is: 
30 
verification/global_ocean.cs32x15/input/proj_cs32_2uEvN.bin 
31 

32 
rotate_csAg_EN.m :: Rotate cell centered (Agrid) U,V vector components 
33 
of CSgrid to East,North directions. works for 2D or 3D vector. 
34 

35 
rotate_csCg_EN.m :: Rotate MITgcm Cgrid U,V vector components of CSgrid 
36 
to East,North directions. done in 2 steps: 1rst compute the average 
37 
at the cellcenter and then rotate the vector. 
38 

39 
uvLatLon2cube.m :: interpolate a 2 components vector field U,V (on Agrid) 
40 
to the MITgcm CSgrid (Cgrid). 
41 

42 
More in dir bk_line: 
43 
Scripts that generate "brokenlines" (used to compute Meridional Transport 
44 
streamfunction and vertically integrated horizontal streamfunction) 
45 
+ a README file + an example for each one. 