Example of cube22 plots and matlab script to read cs510 output and generate plots.
Each of the cube runs produces a set of binary outputs. These are archived on
lou.nas.nasa.gov at NASA Advanced Supercomputing Division. The full list of output files is listed
here. The description below
illustrates basic reading in and plotting of binary outputs from
the CS510 cube sphere ECCO2 experiments using matlab. A more extensive
collection of matlab plotting programs can be found here under the MITgcm web site.
The images below show fields from the cube22 simulation.
These images were generated using matlab from the following cube22 output files
- ETAN.0000337464.data. This file holds the sea-surface height
in meters (identified by the file prefix ETAN)
at timestep 337464 (identified by the file suffix 0000337464.
- UVEL.0000322128.data . This file holds the local X-direction
component of velocity in m/s (identified by the file prefix UVEL)
at timestep 322128 (identified by the file suffix 000322128).
- XC.data. This file holds the longitudinal coordinate in degrees east/west
(identified by the file prefix XC) of the cell centers of the CS510
mesh.
- YC.data. This file holds the latitudinal coordinate in degress north/south
(identified by the file prefix YC) of the cell centers of the CS510
mesh.
The matlab program used to make the plots can be seen here. The input files
read by the program can be downloaded from here (some files are quite large). The matlab program produces the plots shown below. The program shows how to read in CS510
simulation outputs, how to parse the simulation timestep into a calendar date and time and
how to plot the outputs to show the orientation of the cube mesh faces.
- Sea-surface height plot
.
- Surface local x-direction velocity plot
.
- Longitudinal coordinate of cell center plot
.
- Latitudinal coordinate of cell center plot
.
All ECCO2 model and data products are freely available.
Questions and requests to Chris Hill
or Dimitris Menemenlis.