[Thread Prev][Thread Next][Index]

Re: How to mask data over ocean


reply in haste, repent in leisure....

I made several mistakes in the ferret command to replace the water (or below sea level) points in my original reply...

The corrections commands should be something like:
yes? if rose[y=25S:34S,x=135E:140E] le 0 then temp=100 else rose
^^ ^^^^^^^^^

Sorry for any confusion.


On Saturday, September 29, 2001, at 02:08 PM, Mark Verschell wrote:


I can't promise a nifty way, but at least this is a way that will work, IF you are working with relatively small parts of the globe...

First, create a subset of the etopo5 data set over the region you're interested in (from your example 10S:45S, 110E:155E)

Then, apply corrections to the data set
yes? if rose[y=25S:34S,x=135E:140E] eq 0 then temp=100

and do this for the other regions, shading the mask to make sure you make no errors, changing the variable names each time of course

(You can zoom in and really look at it closely to make your region selections more carefully, just use explicit regions in the plot command rather than using set region otherwise you mess up your new land mask)

Then save this new mask and use that to mask your rainfall data.

Hope this helps,


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement