Hi, Diane<br><br>Thanks a lot for your help, your solution works perfectly for polar stereographic projections. But on other projections it still doesn't show anything from -180º to 0º longitude. If I try to plot just for North America, I get an "Entire Grid Undefined" message.<br>
<br>With 'set lon -180 180' it only draws from -180º to 0º. With lon 180 540 only draws from 0º to 180º. With lon -360 0 it won't draw anything. And with 0 360 it will draw the entire globe, but with the problem (for me) of getting centered in lon 180º which cuts Europe in half.<br>
<br>Thanks.<br><br><div class="gmail_quote">2009/5/13 Diane Stokes <span dir="ltr"><<a href="mailto:Diane.Stokes@noaa.gov">Diane.Stokes@noaa.gov</a>></span><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi, Filipe.<br>
<br>
I recently ran into the same problem with a different dataset on the<br>
same grid.<br>
<br>
The spatial resolution the data you are viewing is 1/12 of a degree.<br>
The increment in the xdef and ydef statements of your control file is<br>
0.083. (That's how the increment is stored in the grib file, so it's<br>
not a result of anything you did wrong).<br>
<br>
That increment value is not sufficient to wrap the data completely<br>
around the globe. You have 4320 points in your xdef:<br>
4321 x 0.083 = 358.643<br>
<br>
If you zoom in on the north pole, you'll see the data displayed does not<br>
extend far enough north either.<br>
<br>
If you add a few significant digits to the increment values in your xdef<br>
and ydef statements, eg:<br>
<br>
xdef 4320 linear 0.042000 0.083333<br>
ydef 2160 linear -89.958000 0.083333<br>
<br>
you'll get a more accurate mapping of the data. The values above will<br>
also bring the left and right edge of the grid sufficiently close that<br>
GrADS will know that the data does wrap the globe. You'd then be able<br>
to display the field when setting longitude values less than 0 or<br>
greater than 360.<br><font color="#888888">
<br>
Diane</font><div><div></div><div class="h5"><br>
<br>
<br>
Filipe wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello all,<br>
<br>
I'm trying to plot sst data from <a href="ftp://polar.ncep.noaa.gov/pub/sst/ophi/" target="_blank">ftp://polar.ncep.noaa.gov/pub/sst/ophi/</a> ,<br>
but for some reason I get a white strip at longitude 0º when using nps<br>
projection, as you can see from this pic:<br>
<br>
<a href="http://img8.imageshack.us/img8/2815/userlist2.png" target="_blank">http://img8.imageshack.us/img8/2815/userlist2.png</a><br>
<br>
When using any other projection, grads only draws from lon > 0º, as you can<br>
see from these pics:<br>
<br>
<a href="http://img8.imageshack.us/img8/5487/userlist1.png" target="_blank">http://img8.imageshack.us/img8/5487/userlist1.png</a><br>
<a href="http://img8.imageshack.us/img8/7201/userlist3.png" target="_blank">http://img8.imageshack.us/img8/7201/userlist3.png</a><br>
<br>
I don't know what I'm doing wrong, hope someone can help me. I'm using<br>
latest grib2ctl.pl, wgrib v1.8.0.12x and grads 2.05a (I also tested v1.9 and<br>
opengrads 2.05). I used the script <a href="http://pastebin.com/f3dc65253" target="_blank">http://pastebin.com/f3dc65253</a><br>
<br>
Best regards.<br>
</blockquote>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Fil - MeteoPT.com<br>