All, I've recently needed to generate some ocean basin masks and thought I'd share the result. Starting with the dataset submitted to the ferret users group by Jim Orr http://ferret.pmel.noaa.gov/Ferret/Mail_Archives/fu_2004/msg00323.htmlI modified it slightly and now offer up the following dataset: http://ferret.pmel.noaa.gov/cgi-bin/nph-dods/data/Basins.ncThis dataset has the following variables: ATLANTIC Atlantic Ocean Mask 1:360 1:180 ... ... PACIFIC Pacific Ocean Mask 1:360 1:180 ... ... ARCTIC Arctic Ocean Mask 1:360 1:180 ... ... MARGINAL_SEAS Marginal Seas Mask 1:360 1:180 ... ... INDIAN Indian Ocean Mask 1:360 1:180 ... ...which are defined as 1 in the named ocean basin and missing everywhere else. Here's an example plotting time series of SST for the entire Indian Ocean: use "http://ferret.pmel.noaa.gov/cgi-bin/nph-dods/data/Basins.nc" use "http://ferret.pmel.noaa.gov/cgi-bin/nph-dods/data/coads_sst.nc" let indian_sst = sst * indian[d=1,gxy=sst[d=2]] set region/t="15-jan-1960":"15-jan-1990" plot indian_sst[x=@ave,y=@ave]Enjoy! -- Jon |