You are right, I swapped the results. It seems that for some reason grads won't draw longitude < 0º with 'set lon' starting with a negative number. This won't be a problem (drawing first lon > 0º and then lon < 0º) except for the robinson projection which has to start with -180.<br>
<br>Adding more digits unfortunately didn't work.<br><br>Thank you very much for your help.<br><br><div class="gmail_quote">2009/5/14 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;">Sorry that didn't fix all cases for you. It worked for me running grads<br>
v1.9b4 and v2.0.a3 on IBM AIX.<br>
<br>
I'm confused by your result with:<br>
'set lon -180 180'<br>
you only see -180º to 0º. And with:<br>
'set lon 180 540'<br>
you only see 0º to 180º. Do you have those results swapped in the text<br>
of your msg?<br>
<br>
If you don't get advice from anyone else, I could only suggest you add a<br>
couple more significant digits (3's in this case) to your increment values.<br>
<br>
Good luck.<br>
Diane<br>
<br>
<br>
Fil wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
Hi, Diane<br>
<br>
Thanks a lot for your help, your solution works perfectly for polar<br>
stereographic projections. But on other projections it still doesn't<br>
show anything from -180º to 0º longitude. If I try to plot just for<br>
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<br>
only draws from 0º to 180º. With lon -360 0 it won't draw anything. And<br>
with 0 360 it will draw the entire globe, but with the problem (for me)<br>
of getting centered in lon 180º which cuts Europe in half.<br>
<br>
Thanks.<br>
<br>
2009/5/13 Diane Stokes <<a href="mailto:Diane.Stokes@noaa.gov" target="_blank">Diane.Stokes@noaa.gov</a><br></div>
<mailto:<a href="mailto:Diane.Stokes@noaa.gov" target="_blank">Diane.Stokes@noaa.gov</a>>><div><div></div><div class="h5"><br>
<br>
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>
<br>
Diane<br>
<br>
<br>
<br>
Filipe wrote:<br>
<br>
Hello all,<br>
<br>
I'm trying to plot sst data from<br>
<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<br>
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º,<br>
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<br>
using<br>
latest grib2ctl.pl, wgrib v1.8.0.12x and grads 2.05a (I also<br>
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>
<br>
<br>
<br>
<br>
--<br>
Fil - MeteoPT.com<br>
</div></div></blockquote>
</blockquote></div><br>