[gradsusr] Omit Missing Values

[Jeff Duda]

Brian,
Use the maskout function.
'define ucurrent = maskout(uogrdds,98-uogrdds)'
'define vcurrent = maskout(vogrdds,98-vogrdds)'
'e = mag(ucurrent,vcurrent)'
'd e'

I think that should work.  Let me know if it doesn't.

Jeff Duda

On Wed, Jan 23, 2013 at 11:32 AM, Brian Walsh <brianw at provair.com> wrote:

>  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).
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130123/6917355e/attachment-0003.html