<div dir="ltr"><div><div>Joseph,<br><br></div>I would use copygb.  Copygb has a few types of interpolation<br></div><div>methods and you can define fields that are vector interpolated.<br></div><div>(Affects winds near the poles.)<br></div><div><br></div>Wesley<br><div><br><br><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 25, 2018 at 9:59 PM, Joseph Basconcillo <span dir="ltr"><<a href="mailto:jbasconcillo@pagasa.dost.gov.ph" target="_blank">jbasconcillo@pagasa.dost.gov.ph</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi, I have NCEP-RA2 data that I want to regrid to N768 resolution, which is about 0.17 deg in midlatitude.<br>
<br>
I used:<br>
<br>
'define newvar=re(var, 0.17)'<br>
<br>
And it worked. <br>
<br>
However, the file size is smaller than I expected. The original size is 63MB (2.5 deg) and when I regrid them to 0.17 deg, the file size resulted to about 17.5MB. I have my doubts considering I re-gridded the data to higher resolution. File size should be higher. That was what I expected.<br>
<br>
Then when I used to tloop and area-average to test the data:<br>
<br>
'd tloop(aave(var,global))'<br>
<br>
A prompt appeared that "entire grid are set to missing data". And it kept looping for few minutes. My guess is that the re-gridding did not "regrid" through time and considered only the first timestep. Can you help me on this?<br>
<br>
Also, is there an option for regridding method (e.g. bilinear)? <br>
<br>
Thank you.<br>
<br>
Yours,<br>
<br>
Joseph Basconcillo<br>
______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
</blockquote></div><br></div></div></div></div></div>