[Thread Prev][Thread Next][Index]

Re: [ferret_users] Visualizing GEBCO with Ferret



Hi Tim,

here my final version, how to deal with GEBCO-topography files:

1. Unfortunately the name of the depth field is "z". This name is used by ferret for axes. I did not find a ferret-solution, you have to rename it using nco-tools. It is easy, if nco-tools are installed. Otherwise I recommend to do so, they are very useful.

ncrename -v z,za  your-file

2. Ferret prefers a different way to have axes and data defined and linked with each other. Hence, you have to reformat the data. The data format in the files you may download are for other tools described at the web-pages.

I am not sure about the following!!! I read the manual, but could not find some statement, where grid points are exactly defined. So I downloaded a 1 deg x 1 deg rectangle and got 61 coordinate points for the "1min"-topography and 120 points for the"08" highly resolving data set respectively.
Hence,
- in the 1 deg file the coordinates are cell edges
- in the 0.5 deg file the coordinates are cell mid points

Please inspect your file with ncdump to be sure, that your final grid after reformatting the data has the same size.

Write a script map_to_axes.jnl. For sure there are mor elegant solutions.
use "$1"
define axis/x=`X_RANGE[i=1]`:`X_RANGE[i=2]`/npoints=`dimension[i=1]`/edges /units="degrees_east" xax define axis/y=`Y_RANGE[i=1]`:`Y_RANGE[i=2]`/npoints=`dimension[i=2]`/edges /units="degrees_north" yax !Please remove the edges - specification, if you read highly resolved topography files.!!!!!!!!!!!!!!!!!!!!!
define grid/x=xax/y=yax xygrid
let dummy=0*x[gx=xygrid]+0*y[gy=xygrid]
let depth=reshape(za,dummy)
SET MEMORY/SIZE=500
save/file=zw.nc/clobber depth
can data/all
can var depth
can grid xygrid

! Now flip the y-axis. Read data again and write to a new file.
SET MEMORY/SIZE=500
use/order=X-Y zw.nc
let/unit=m/title="gebco topography" topo=depth
save/file=depth_$1/clobber topo

Now you can use ferret to make plots.
use depth_your-filename
shade topo

I recommend, not to use the full data sets with the shade command. It is huge and the plot - resolution is not fine enough for the details.
Otherwise please read the manual.
Good luck,
Martin

Tim Loomis wrote:

All,
I'm new to Ferret and would like to find someone who has used it to access the gebco_08.nc. I need some high resolution images out of it. Any help would be appreciated. Any tutorials out there?....Tim






[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement