Back to home page

MITgcm

 
 

    


Warning, /utils/matlab/cs_grid/README is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit 4ec37fd8 on 2023-10-03 22:00:32 UTC
d09e20f1d3 Jean*0001 Matlab scripts (written by Alistair) specific to Cubed-Shere grid:
                0002  [initially in utils/matlab and moved here Sep. 15, 2005]
                0003 
                0004 Interpolate to Lat-Lon grid:
4ec37fd829 Jean*0005  cube2latlon.m
                0006  cube2latlon_fast.m
                0007  cube2latlon_preprocess.m
d09e20f1d3 Jean*0008  uvcube2latlon.m
4ec37fd829 Jean*0009  uvcube2latlon_fast.m
9d7ae97fe2 Jean*0010  uvcube2latlongrid.m  (faster version, from Guillaume Maze)
d09e20f1d3 Jean*0011 
                0012 make graphics of cs-grid fields:
4ec37fd829 Jean*0013  merccube.m
                0014  displaytiles.m
                0015  drawedges.m
                0016  fancycube.m
                0017  plotcube.m
d09e20f1d3 Jean*0018 
                0019 others:
4ec37fd829 Jean*0020  permutetiles.m
d09e20f1d3 Jean*0021  tiles.m
                0022 
                0023 --------------------------------------------------------------------------
                0024 Some tools specific to the Cubed-Sphere grid :
896b20e91e Jean*0025 
                0026 split_C_cub.m   :: split a 2d,3d global CS-field into 6 faces ;
4ec37fd829 Jean*0027         apply to cell-centered CS-field ; extend the face array by
896b20e91e Jean*0028         adding 1 row & 1 column (like the exchanges)
                0029 split_Z_cub.m   :: split a 2d,3d global CS-field into 6 faces ;
4ec37fd829 Jean*0030         apply to cell-corner CS-field ; extend the face array by
896b20e91e Jean*0031         adding 1 row & 1 column (like the exchanges)
4ec37fd829 Jean*0032 split_UV_cub.m  :: split the 2 C-grid components of a vector
                0033         (2d,3d global CS-field) into 6 faces ; extend the face array
896b20e91e Jean*0034         by adding 1 row or 1 column (like the exchanges)
                0035 
                0036 cubeZ2latlon.m  :: interpolate a CS-field to regular Lat-lon grid ;
                0037         similar to cube2latlon.m but apply to a field that is
                0038         located at the corner point of the C-grid cell.
                0039         (deal with the 2 "missing corners")
                0040 
                0041 calc_vort_cs.m  :: compute vorticity (vertical component) on CS-grid
4ec37fd829 Jean*0042         from MITgcm output velocity fields u,c (on C-grid);
                0043         Use the model discretization of vorticity.
896b20e91e Jean*0044 
                0045 the following scripts are reading a file that contains the orientation
4ec37fd829 Jean*0046         of the local direction of the CS-grid relative to E-W,N-S direction.
                0047   file= proj_cs[]_2uEvN.bin :: cos & sin of rotation angle,
896b20e91e Jean*0048                       format:: real*8, big-endian, size=32x6x32 x2(cos+sin)
                0049         for standard cs32 grid, the corresponding file is:
4ec37fd829 Jean*0050            verification/global_ocean.cs32x15/input/proj_cs32_2uEvN.bin
896b20e91e Jean*0051 
4ec37fd829 Jean*0052 rotate_uv2uvEN.m :: Rotate CS-grid vector U,V components to East,North directions.
                0053         Input at grid-cell center (A-grid) or grid-cell edges (C-grid);
                0054         works for 2D or 3D vector.
896b20e91e Jean*0055 
                0056 uvLatLon2cube.m :: interpolate a 2 components vector field U,V (on A-grid)
                0057         to the MITgcm CS-grid (C-grid).
b75ad9d212 Dimi*0058 
4ec37fd829 Jean*0059 mk_psiLine_CS.m :: define path to integrate horizontally horizontal transport
                0060         for computing Barotropic Stream-Function, starting from the N.Pole and
                0061         going southward.  => write matlab file: psiLine_N2S_cs32.mat
                0062 use_psiLine.m :: an example that uses "psiLine_N2S_cs32.mat" to compute & plot
                0063         Barotropic Transport Stream-function.
                0064 
                0065 mk_isoLat_bkl.m :: define "broken-lines" as close to latitude-target vector
                0066         "yLat", for computing Meridional Transport stream-function.
                0067         => write matlab file: isoLat_${ydim}_cs32.mat (with ydim = length(yLat))
                0068 use_isoLat_bkl.m :: an example that uses "isoLat_59_cs32.mat" to compute & plot
                0069         Meridional Transport Stream-function.
                0070 
                0071 %--------------------------------------------------
896b20e91e Jean*0072 More in dir bk_line:
                0073  Scripts that generate "broken-lines" (used to compute Meridional Transport
4ec37fd829 Jean*0074   stream-function) as close as great-circle and some ploting script for CS-grid.
                0075  see bk_line/README file there
896b20e91e Jean*0076 
b75ad9d212 Dimi*0077 And in dir read_cs:
                0078  Scripts for reading and writing cube-sphere binary output in various ways.
                0079  Look at read_cs/read_cs.txt for example usage.