<div dir="ltr"><div><span style="color:rgb(0,0,0)"><font face="times new roman, serif">Hello, </font></span></div><div><pre><font face="times new roman, serif"><i style="color:rgb(0,0,0)">I would like to skip a </i><i><font color="#000000">number of stations data when plot it with grads.I find the message below describing a trick used in COLA, i am wondering if this method still valid </font></i></font></pre><pre><font face="times new roman, serif"><i><font color="#000000">or is there another </font></i></font><i style="font-family:'times new roman',serif"><font color="#000000">one.</font></i></pre><pre><i><font color="#000000" face="times new roman, serif">Thanks.</font></i></pre><pre><i><font color="#000000" face="times new roman, serif">Rabah </font></i></pre></div><b style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><div><b style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><br></b></div><div><b style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><br></b></div>Jennifer Adams</b><span style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"> </span><a href="mailto:gradsusr%40gradsusr.org?Subject=Re:%20%5Bgradsusr%5D%20Skip%20stations%20when%20plotting%20station%20data&In-Reply-To=%3C0AADE852-4D26-4343-8C14-D4A89670CA21%40cola.iges.org%3E" title="[gradsusr] Skip stations when plotting station data" style="font-family:'Times New Roman';font-size:medium">jma at cola.iges.org </a><br style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><i style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium">Thu Jan 20 10:16:09 EST 2011</i><span style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"></span><p style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"></p><ul style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><li>Previous message: <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/029182.html">[gradsusr] Skip stations when plotting station data</a></li><li>Next message: <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/029187.html">[gradsusr] How to set time attribute in sdfwrite</a></li><li><b>Messages sorted by:</b> <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/date.html#29192">[ date ]</a> <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/thread.html#29192">[ thread ]</a> <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/subject.html#29192">[ subject ]</a> <a href="http://gradsusr.org/pipermail/gradsusr/2011-January/author.html#29192">[ author ]</a></li></ul><hr style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><pre style="color:rgb(0,0,0)">Brian Doty has solved this problem with a C program that creates a
'filt' variable for the surface metars. The code does the following:
1) Assigns priorities to each report that are based on a combination
of the number of variables available and the intensity of precipitation
2) Filters out observations at a particular radius: first pass based
on priority, subsequent passes based on proximity to other obs.
3) All the obs are then written out with an extra data variable,
called 'filt', that you can use when drawing. The filt variable ranges
between 1 and 6.
For example,
'open <a href="http://monsoondata.org:9090/dods/stn/metar/past48'">http://monsoondata.org:9090/dods/stn/metar/past48'</a>
'set gxout stnmark'
'set lat 25 50'
'set lon -125 -65'
'd maskout(ts,filt-4)' ;* draws fewer obs
'd maskout(ts,filt-2); ;* draws more obs
We use Brian's code operationally at COLA (for the real time metars
behind our GDS) but cannot support it beyond that. Please write to me
offline if you would like a copy.
--Jennifer</pre><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><br></div><br></div></div>
</div>