I finally got it working! I was using a variable (define sst TMPsfc-273.15), then I just tried without the variable and now it draws the entire globe. :D<br><br>Thanks for the help!<br><br><div class="gmail_quote">2009/5/15 Fil <span dir="ltr"><<a href="mailto:fil@meteopt.com">fil@meteopt.com</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,<br><br>I get the same results with grib1 data. I also tried in a 64-bit box and I didn't get any error, but I get the same resulting map. Tried grads 1.9 and opengrads 2 and the result is the same.<br>
<br>Filipe<br>
<br>btw: grib1 is a LOT faster than grib2.<div><div></div><div class="h5"><br><br><div class="gmail_quote">2009/5/15 Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</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;">
<div>Wacky! <div>My earlier sample was grib, but I get the same good results with the grib2 version of this data.... on my mac. When I tried it on my 64-bit linux box, I got an error I've never seen before: </div>
<div><br></div><div><div>ga-> d sst</div><div> jpc_decode return = 0 </div><div>Constant field. Value = 271.33</div><div>ga-> </div><div><br></div><div>That "jpc_decode" message is in the grib2 library in the dec_jpeg2000.c souce file. I have no idea what it means. Or where it got the 271.33 number. I will have to poke at it further. Could be a classic case of operator error. In the meanwhile ... Can you try using the grib1 version of this data instead? Go back to the original FTP site that you started with: <a href="ftp://polar.ncep.noaa.gov/pub/sst/ophi/" target="_blank">ftp://polar.ncep.noaa.gov/pub/sst/ophi/</a></div>
<div><br></div><div>Jennifer</div><div>p.s. Why do these puzzles show up on Friday afternoons?</div><div><br></div></div><div><br></div><div><br><div><div><div></div><div><div>On May 15, 2009, at 2:27 PM, Fil wrote:</div>
<br></div></div><blockquote type="cite"><div><div></div><div>Hi Jennifer,<br><br>With 0.083333 in the ctl it does shows "LON set to 0 360". Maybe it's something I'm doing wrong but even then it won't draw lon < 0º. I'm sending the files I'm using for testing. The grib can be downloaded from:<br>
<br><a href="http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/sst.20090514/rtgssthr_grb_0.083.grib2" target="_blank">http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/sst.20090514/rtgssthr_grb_0.083.grib2</a><br><br>
Thanks.<br><br><div class="gmail_quote"> 2009/5/15 Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</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;">
<div>Fil, <div>I have downloaded a data sample and done some testing. When you open your descriptor file with GrADS, look at the info written to the command window. Does it say "LON set to 0.042 359.944" or "LON set to 0 360"? </div>
<div><div><br></div><div>The answer will depend on whether you have enough significant digits in the grid increment value in the XDEF and YDEF statements in your descriptor file. When I use "0.083333" instead of "0.083", then GrADS calculates that the file does indeed 'wrap' the globe and 'set lon -180 180' works as it should. </div>
<div><br></div><div>Jennifer</div><div><br></div><div><br></div><div>GrADS automatically issues a 'set lon 0 360' command when it opens a file that wraps the globe. It </div><div><br></div></div><div><div><div></div>
<div><br><div><div>On May 15, 2009, at 1:31 PM, Fil wrote:</div><br><blockquote type="cite">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" target="_blank">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> 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><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></blockquote></div><br></div></div><font color="#888888"><div> <span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a></div><div><br>
</div><br></span></span></div></span> </div><br></font></div></div></blockquote></div><br> </div></div><span><<a href="http://grads.gs" target="_blank">grads.gs</a>></span><span><sst.ctl></span><span><sst.png></span></blockquote>
</div><div><br><div> <span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
<span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; border-spacing: 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div>
--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a></div><div><br>
</div><br></span></span></div></span> </div><br></div></div></div></blockquote></div><br>
</div></div></blockquote></div><br>