<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"> Dear Chuck,<br><br> I've just tested and your tip helped me. The script worked like a charm after including the 'c events' command.<br><br> Thank you very much!<br><br> Best,<br><br> Thiago.<br><br>--- On <b>Fri, 15/5/09, Charles Seman <i><Charles.Seman@noaa.gov></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Charles Seman <Charles.Seman@noaa.gov><br>Subject: Re: Converting binary data to ascii<br>To: thi_veloso@yahoo.com.br<br>Date: Friday, 15 May, 2009, 1:38 PM<br><br><div class="plainMail">Dear Thiago,<br><br>Perhaps it may help if you add a 'c events' in your script as follows:<br> 'c events'<br> x=x+1<br> endwhile<br>
y=y+1<br> endwhile<br> t=t+1<br> endwhile<br><br>In a local script here, a test in which the 'c events' command was commented out (code segment show below) generated a "Segmentation fault" message:<br>* 'c events'<br> xx = xx+1<br>endwhile<br> yy = yy+1<br>endwhile<br> zz = zz+1<br>endwhile<br> tt = tt+1<br>endwhile<br><br>I hope this helps,<br>Chuck<br><br>Thiago Veloso wrote:<br>> Dear all,<br>> I have a relatively long (6-year) series of sattelite-based hourly precipitation estimates, which is in binary format. Considering I need to manipulate this data on external programs, I'm trying to extract it to ascii files.<br>> Currently I´m using GrADS v1.9b4 running on a ordinary 32-bit desktop Ubuntu 9.04. Some years ago a colleague (now retired) developed a script which used to work in out previous datasets. But it won't work with these specific
precipitation estimates files (maybe due to the enormous volume). GrADS suddenly crashes with an "segmentation fault" error message.<br>> <br>> Would anybody gently help me to identify what is causing the error? Following is the script code...<br>> <br>> Thanks in advance.<br>> <br>> 'reinit'; 'c';<br>> prompt 'Variable: '; pull var<br>> prompt 'Latitude (lat1 lat2): '; pull lat<br>> prompt 'Longitude (lon1 lon2): '; pull lon<br>> prompt 'Tempo (t1 t2): '; pull t<br>> 'open /usr/local/grads-1.9b4/lib/CMORPH_8km-30-minute.ctl'<br>> ' set lat 'lat; 'set lon 'lon; 'set t 't;<br>> 'q dims'; x=sublin(result,2); y=sublin(result,3); t=sublin(result,5);<br>> x1=subwrd(x,11); y1=subwrd(y,11); t1=subwrd(t,11); x2=subwrd(x,13); y2=subwrd(y,13);
t2=subwrd(t,13);<br>> if(x2=''); x1=subwrd(x, 9); x2=x1; endif if(y2=''); y1=subwrd(y, 9); y2=y1; endif if(t2=''); t1=subwrd(t, 9); t2=t1; endif<br>> x=x1; xmax=x2; y=y1; ymax=y2; t=t1; tmax=t2;<br>> lixo=write(var'.txt', 'lon lat lev time 'r)<br>> while (t<=tmax)<br>> 'set t 't; y=y1;<br>> while (y<=ymax)<br>> 'set y 'y; x=x1;<br>> while(x<=xmax)<br>> 'set x 'x; 'set t 't; r1=xyzt(); 'd 'var; r=subwrd(result, 4);<br>> lixo=write(var'.txt', r1%r,append);<br>> x=x+1<br>>
endwhile<br>> y=y+1<br>> endwhile<br>> t=t+1<br>> endwhile<br>> say ' End of Process!'<br>> <br>> function xyzt()<br>> 'q dims'; x=sublin(result,2); y=sublin(result,3); z=sublin(result,4); t=sublin(result,5);<br>> x=subwrd(x,6); y=subwrd(y,6); z=subwrd(z,6); t=subwrd(t,6);<br>> return(x' 'y' 'z' 't' ')<br>> function countwrd(x)<br>> i=0<br>> while (subwrd(x,i+1)!='')<br>> i=i+1;<br>> endwhile<br>> return(i)<br>> function openfile(arq)<br>> 'open 'arq;
r=sublin(result,2);<br>> if subwrd(r,2)='Error'<br>> 'sdfopen 'arq;<br>> endif<br>> return<br>> <br>> <br><br>-- <br>Please note that <a ymailto="mailto:Charles.Seman@noaa.gov" href="/mc/compose?to=Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>email address, not <a ymailto="mailto:cjs@gfdl.noaa.gov" href="/mc/compose?to=cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>.<br><br>********************************************************************<br>Charles Seman <a ymailto="mailto:Charles.Seman@noaa.gov" href="/mc/compose?to=Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>U.S. Department of Commerce / NOAA / OAR<br>Geophysical Fluid Dynamics Laboratory voice:
(609) 452-6547<br>201 Forrestal Road fax: (609) 987-5063<br>Princeton, NJ 08540-6649 <a href="http://www.gfdl.noaa.gov/%7Ecjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>********************************************************************<br><br>"The contents of this message are mine personally and do not necessarily<br>reflect any position of the Government or NOAA."<br><br></div></blockquote></td></tr></table><br>