[gradsusr] Omit Missing Values

[Brian Walsh]

Hello,

I'm plotting ocean currents and the grib file contains values of -99 for grid points near land. So, when I try to plot the combined current speed/direction (u, v), it has values very high near and over land. I can use basemap to cover the land, but there are grid points near land and over water that are still showing (see attached sample image where red is the extreme/missing values).

I think the const command is what I need, but I'm not sure how to set this up when I want to combine the u, v parameters. I tried:

        'define ucurrent=const(uogrddsl, -99, -u)'
        'define vcurrent=const(vogrddsl, -99, -u)'
        'e=mag('ucurrent','vcurrent')'
        'd e'

But, no luck. Any ideas?

Thanks

Brian


This email and any attached files are intended solely for the use of the individual(s) and copied addressee(s) and may contain privileged and/or confidential information. Any distribution, use or copying by any means of this information is strictly prohibited. If you are not the intended recipient(s), please notify the sender by telephone or return the email and delete the material from your system(s).
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130123/3c64a7eb/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SURFACECURRENTS.012Hour.jpg
Type: image/jpeg
Size: 216853 bytes
Desc: SURFACECURRENTS.012Hour.jpg
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130123/3c64a7eb/attachment-0003.jpg