Generating global avg timeseries from gx3v7 output

[chris_fletcher@utoronto_ca]

Hi, I am trying to generate timeseries of global mean volume-weighted SALT and TEMP from gx3v7 history output. The POP diagnostics script tool_lib/gen_globavgTS_netcdf.csh is hardwired for gx1v5 -- this is because it points to the TCELL volume file for gx1v5 resolution.

So I am wondering if there is an equivalent TCELL volume file for gx3v7 resolution? I assume there must be, since I found these plots online:

http://www.cesm.ucar.edu/experiments/cesm1.0/diagnostics/b40.t31x3.037/pdts/diagts_TEMP.gif​

I also tried to create my own cell volume grid from a pop history file, using dz, HTE and HTN for the height, width and length of each cell, respectively. But when I did a test and used this method to compute the volume from a gx1v5 history file and compared the results to the gx1v5_tcell_volume.nc from NCAR, there were some substantial differences that made me think my method is not exactly correct.

So, would it be possible for either (a) someone to send me the gx3v7 tcell volume file, or (b) someone to post the correct method for calculating tcell volume from history files?

Thanks in advance.

 

[bates]

Hi,

The code you are looking at is not used in the ocean diagnostics package, but was created as an offline tool to generate globally averaged volume weighted temperature and salinity files. We do not have a tcell volume file for the gx3v7 grid, but we do have one for the gx3v5 grid. I have attached that here in case it might work for you. I am not sure what the differences between version 5 and 7 for the gx3 grid, but if I find out, I'll post here again. They are both nominal 3 degree grids.

We also don't have the code that created these tcell volume files, but it should be relatively easy to make one on your own. From the history files, use TAREA*dz to create the cell volume and be sure to mask it using kmt so that you are only using ocean points.

Hope this helps!
Susan
NCAR Ocean Model Science Liaison