<div class="gmail_quote">On Sun, Feb 15, 2009 at 12:29 PM, Steven <span dir="ltr"><<a href="mailto:sweiss@iafrica.com">sweiss@iafrica.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
<br>
I'm using a 0.5x0.5 resolution for GFS data. Ideally within a grid box,<br>
there could be a few towns and cities that I would like produce data for<br>
e.g. tmp2m at a fixed lat/lon. Is there anyway to get the resolution finer -<br>
say to 0.25x0.25 in order to get more values?<br></blockquote><div><br></div><div>A simple, brute force way is use an interpolation function such as linterp() or re() to defined your variable on a higher resoution grid, then work from that defined variable. For this to work you have to create a fake.ctl that defines this high resolution grid, e.g.,</div>
<div><br></div><div>dset fake.dat</div><div>options template</div><div>xdef ...</div><div>ydef ...</div><div>zdef ..</div><div>tdef ...</div><div>vars 1</div><div>var 0 0 dymmy</div><div>endvars</div><div><br></div><div>(the data file fake.dat does not have to exist.) Then,</div>
<div><br></div><div>ga-> open model</div><div>ga-> open fake</div><div>ga-> set dfile 2 (now you are in the fine resolution dimension environment)</div><div>ga> define ts_fine = re(ts.1,0.25)</div><div><br></div>
<div>then use ts_fine in your plots. Another possibility is to use PDEF to have this interpolation done on the fly, behind the scenes.</div><div><br></div><div><br></div><div> Arlindo</div><div><br></div></div>-- <br>Arlindo da Silva<br>
<a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>