Hi Carissa,<br><br> From my experience, you are right. It's a small stack problem. The jpeg2000 compression<br>uses the Jasper library which requires a large stack when processing large grids. You can<br>change the stack size by ulimit -s XXXXX. (Only works once per terminal session.)<br>
<br> As for the NOMADS implementation, you have two choices<br><br> (1) increase the stack size on the NOMADS machine. Not sure how this affects the memory footprint.<br> (probably OS dependent) <br> (2) change the sst from jpeg2000 packing into complex packing<br>
wgrib2 IN.grb -set_grib_type c2 -grib_out OUT.grb<br><br> The advantage of (2) is that complex packing is much faster than jpeg2000 to decode.<br><br> Wesley Ebisuzaki<br><br><br><br><div class="gmail_quote">
On Thu, Mar 14, 2013 at 12:25 PM, Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">Hi, Carissa -- <div>I can draw your temperature field on my mac laptop, but get the same error you do on my linux server ("Constant field. Value= 271.35"). There is a non-informative message "jpc_decode return" that comes from the grib2 library. The 'jpc_decode' routine is part of the jasper library, which is required by the grib2 library, but I can't find anything about it in the jasper documentation. In the library's code, there are comments that this routine is the JPEG-2000 decoder. So the problem likely has something to do with the JPEG compression. I have no idea why it would work on AIX and OSX but not linux, but I don't think ulimit or the GrADS version has anything to do with it (I was testing with GrADS 2.0.2, but the grib2 I/O code hasn't changed in a long time). Can you try using wgrib2 or some other tool to change the compression and see if that makes a difference? I don't have any other bright ideas … perhaps Wesley can be more helpful. </div>
<div>--Jennifer</div><div><br></div><div><br></div><div><br><div><div>On Mar 14, 2013, at 11:17 AM, Carissa Klemmer - NOAA Federal wrote:</div><br><blockquote type="cite">Users,<br><br>I work for NCEP on the team that supports the operational NOMADS site. We are in the process of adding a new product to GrADS-OpenDAP but are running into issues. I have included the syntax for the control file used below to make the index file. Once the control and index file are created I attempt to view the data plotted through grads. <br>
Here is the problem. I can not get the data to plot anything but a constant value on our sandbox nomads server or on my desktop. I can however see the plot successfully on our operational supercomputer. I have suspicions of the ulimit stack size, but have not been able to confirm this due to the test machines system constraints. So can anyone else get this data to successfully plot? If so please share your wisdom. Why does the exact same control file and grads version work on one system but not another?<br>
<br>The supercomputer is running grads version 2.0.a7 on AIX with ulimit -s of 4000000.<br>NOMADS sandbox is running 2.0.a3.p3 on LINUX with ulimit -s of 16832<br>I have downloaded both 2.0.2 and 2.0.a7 to my desktop (LINUX) with a ulimit -s 8000. <br>
<br>1. The grib2 file: rtgssthr_grb_0.083.grib2<br>which can be found for download under: <a href="http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/sst.20130313/" target="_blank">http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/sst.20130313/</a><br>
<br>2. The control file (sst.ctl) was created from g2ctl as follows:<br>dset rtgssthr_grb_0.083.grib2<br>index sst.idx<br>undef 9.999E+20<br>title /com/gfs/prod/sst.20130312/rtgssthr_grb_0.083.grib2<br>* produced by g2ctl v0.0.8.2<br>
* command line options: /com/gfs/prod/sst.20130312/rtgssthr_grb_0.083.grib2<br>* griddef=1:0:(4320 x 2160):grid_template=0:winds(N/S): lat-lon grid:(4320 x 2160) units 1e-06 input WE:NS output WE:SN res 48 lat 89.9 58000 to -89.958000 by 0.083000 lon 0.042000 to 359.958000 by 0.083000 #points=9331200:winds(N/S)<br>
<br>dtype grib2<br>ydef 2160 linear -89.958000 0.083<br>xdef 4320 linear 0.042000 0.083000<br>tdef 1 linear 00Z12mar2013 1mo<br>zdef 1 linear 1 1<br>vars 1<br>TMPsfc 0,1,0 0,0,0 ** surface Temperature [K]<br>ENDVARS<br>
<br>3. to create the index file I ran:<br>>gribmap -0 -v -i sst.ctl<br><br>4.To view the output<br>>grads<br>>d tmpsfc<br><br>As even more proof that I somewhat know what I am doing, I can successfully create a control and index file using the same steps as above and view pretty graphics for the half degree SST, rtgssthr_grb_0.5.grib2. <br>
<br>Thanks for any insight. <br clear="all"><br>-- <br><font><span style="color:rgb(102,102,102);font-family:georgia,serif">Carissa Klemmer
</span></font><br style="color:rgb(102,102,102);font-family:georgia,serif"><font><span style="color:rgb(102,102,102);font-family:georgia,serif">NCEP Central Operations
</span></font><br style="color:rgb(102,102,102);font-family:georgia,serif"><font><span style="color:rgb(102,102,102);font-family:georgia,serif">Production Management Branch Dataflow Team
</span></font><br style="color:rgb(102,102,102);font-family:georgia,serif"><font><span style="color:rgb(102,102,102);font-family:georgia,serif"><a href="tel:301-683-3835" value="+13016833835" target="_blank">301-683-3835</a></span></font>
_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div><br><div>
<span style="font-size:12px"><span style="border-spacing:0px 0px;text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;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>
<br></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br>