<div dir="ltr">how could i plot the ASCII file in Grads or How could i convert ASCII( with 3 column lat, lon and values)  to .nc file ???<div><br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr">

<div>Kunal Bali<br></div><div>Research Scholar <br></div><div>Radio &amp; Atmospheric Science Division <br></div><div>CSIR - National Physical Laboratory<br></div><div>New Delhi - 110012<br><br></div><div><br></div><div>
<br>
</div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div>


<br><br><div class="gmail_quote">On Mon, Feb 24, 2014 at 11:58 PM,  <span dir="ltr">&lt;<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Send gradsusr mailing list submissions to<br>
        <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of gradsusr digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Wind-chill parameter GrADS - (L.B.)<br>
   2. Scatter plot, buffer overflow (Joowan Kim)<br>
   3. Re: Scatter plot, buffer overflow (Jennifer Adams)<br>
   4. Re: help with fndlvl function? (Freezing level script)<br>
      (Edward Pollock)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 24 Feb 2014 12:52:08 -0500<br>
From: &quot;L.B.&quot; &lt;<a href="mailto:bcbass2989@gmail.com">bcbass2989@gmail.com</a>&gt;<br>
Subject: Re: [gradsusr] Wind-chill parameter GrADS -<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:CAM1sYasmnA7TDdHaHAzqh1tn-4828YDjcY-wt1D9d2OR8PGzJA@mail.gmail.com">CAM1sYasmnA7TDdHaHAzqh1tn-4828YDjcY-wt1D9d2OR8PGzJA@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Leopoldo,<br>
<br>
Based on the first link you mentioned, that formula looks similar to the<br>
old windchill formula. The newer one which you used in your two<br>
calculations are correct. I should mention that the windchill formula is<br>
meant for temperatures of 50F(10C) or less. This is the reason why your<br>
windchill values are only slightly lower than your actual temperature<br>
values.<br>
<br>
*William (L.B.) LaForce*<br>
Meteorologist<br>
Central Michigan University &#39;13<br>
<a href="http://www.tornadoinsanity.com" target="_blank">www.tornadoinsanity.com</a><br>
<br>
<br>
On Mon, Feb 24, 2014 at 8:43 AM, Leopoldo Alvarez &lt;<a href="mailto:lap.bolsa@gmail.com">lap.bolsa@gmail.com</a>&gt;wrote:<br>
<br>
&gt;  Hi Willians,<br>
&gt;<br>
&gt;  Look at this webpage-&gt;<br>
&gt; <a href="http://www.lenntech.com/calculators/wind/wind-chill.htm" target="_blank">http://www.lenntech.com/calculators/wind/wind-chill.htm</a><br>
&gt; 30 km/h (wind)<br>
&gt; 18?C (temperature)<br>
&gt; 12.1?C (windchill)<br>
&gt;<br>
&gt; - - - - - - - - - - -<br>
&gt;<br>
&gt;  In our examples ....<br>
&gt; I use your formulation...<br>
&gt;  - [1] I use your formulation corrected (convert KPH to MPH and K to F to<br>
&gt; input into the equation)<br>
&gt;<br>
&gt; &#39;far=((9/5)*t)-459.67&#39;<br>
&gt; &#39;mph=tg1*2.236936&#39;<br>
&gt; &#39;define chill =<br>
&gt; 35.74+0.6215*far-35.75*pow(mph,0.16)+0.4275*far*pow(mph,0.16)&#39;<br>
&gt; &#39;chill2=(chill-32)/1.8&#39;<br>
&gt; &#39;d chill2&#39;<br>
&gt;<br>
&gt;  Date / wind (MPH) / temp (?F) / windchill (?F)<br>
&gt; 2014-FEB-19-13Z  17.1496   66.0072   64.8952<br>
&gt; 2014-FEB-19-14Z  19.4474   66.7105   65.5746<br>
&gt; 2014-FEB-19-15Z  19.4474   66.2803   65.0114<br>
&gt;<br>
&gt;<br>
&gt;  - [2] This is the other proposal that you have shown me ...<br>
&gt;<br>
&gt; &#39;define chill =<br>
&gt; 13.12+0.6215*(t-273.15)-11.37*pow((tg1*3.6),0.16)+0.3965*(t-273.15)*pow((tg1*3.6),0.16)&#39;<br>
&gt;<br>
&gt;  Date / wind (KPH) / temp (?C) / windchill (?C)<br>
&gt; 2014-FEB-19-13Z   27.5996  18.8929  18.2662<br>
&gt; 2014-FEB-19-14Z   31.2976  19.2836  18.6438<br>
&gt; 2014-FEB-19-15Z   31.2976  19.0446  18.3308<br>
&gt;<br>
&gt; I still think that with a wind of 30 km/h and a temperature of 18.8 ?C,<br>
&gt; the windchill can not be a &#39;lower tenths&#39;..<br>
&gt;  I&#39;m confused, and I have no security ...<br>
&gt; <a href="http://www.easysurf.cc/cnver16.htm" target="_blank">http://www.easysurf.cc/cnver16.htm</a><br>
&gt; <a href="http://es.easycalculation.com/weather/wind-chill.php" target="_blank">http://es.easycalculation.com/weather/wind-chill.php</a><br>
&gt;<br>
&gt; Best,<br>
&gt;<br>
&gt; Leopoldo Alvarez<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140224/baa60ba3/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140224/baa60ba3/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 24 Feb 2014 11:12:33 -0700<br>
From: Joowan Kim &lt;<a href="mailto:joowan@ucar.edu">joowan@ucar.edu</a>&gt;<br>
Subject: [gradsusr] Scatter plot, buffer overflow<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID:<br>
        &lt;CAE7bxZSTFe_u12Nhk2ggi+RXhLQkZWhZgJ-37_JC6U4Tb1=<a href="mailto:UYA@mail.gmail.com">UYA@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Dear Grads users<br>
<br>
I&#39;m plotting a simple scatter plot using grads (v2.0.a9) x86_64, like<br>
<br>
open ../data/model.ctl<br>
set gxout scatter<br>
set lev 500<br>
d u; v<br>
<br>
But if I try it, it dies showing errors like below<br>
Is there any one who experience this issue, and having a solution?<br>
<br>
Thanks!<br>
<br>
*** buffer overflow detected ***: grads terminated<br>
======= Backtrace: =========<br>
/lib64/libc.so.6(__fortify_fail+0x37)[0x30d6f02827]<br>
/lib64/libc.so.6[0x30d6f00710]<br>
/lib64/libc.so.6[0x30d6effe0b]<br>
/lib64/libc.so.6(__snprintf_chk+0x7a)[0x30d6effcda]<br>
grads(gascat+0x250)[0x45d070]<br>
grads(gaplot+0x3fc)[0x45e27c]<br>
grads(gadspl+0x576)[0x4a68b6]<br>
grads(gacmd+0x1607)[0x4a5217]<br>
grads(gaexec+0x294)[0x4a6da4]<br>
grads(gacmd+0xea0)[0x4a4ab0]<br>
grads(main+0x5c7)[0x419477]<br>
/lib64/libc.so.6(__libc_start_main+0xfd)[0x30d6e1ed1d]<br>
<br>
<br>
<br>
--<br>
Joowan Kim<br>
Atmospheric Chemistry Division<br>
National Center for Atmospheric Research<br>
PO Box 3000  Boulder, CO 80307<br>
e-mail : <a href="mailto:joowan@ucar.edu">joowan@ucar.edu</a><br>
         <a href="mailto:joowan.k@gmail.com">joowan.k@gmail.com</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140224/f0093d47/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140224/f0093d47/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 24 Feb 2014 13:28:08 -0500<br>
From: Jennifer Adams &lt;<a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a>&gt;<br>
Subject: Re: [gradsusr] Scatter plot, buffer overflow<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:EE3724F2-4EC7-4503-BA17-D26A7D37AE1A@cola.iges.org">EE3724F2-4EC7-4503-BA17-D26A7D37AE1A@cola.iges.org</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Another user had this problem, and I was unable to reproduce it. Are you using a COLA-provided build of GrADS?<br>
--Jennifer<br>
<br>
On Feb 24, 2014, at 1:12 PM, Joowan Kim wrote:<br>
<br>
&gt; Dear Grads users<br>
&gt;<br>
&gt; I&#39;m plotting a simple scatter plot using grads (v2.0.a9) x86_64, like<br>
&gt;<br>
&gt; open ../data/model.ctl<br>
&gt; set gxout scatter<br>
&gt; set lev 500<br>
&gt; d u; v<br>
&gt;<br>
&gt; But if I try it, it dies showing errors like below<br>
&gt; Is there any one who experience this issue, and having a solution?<br>
&gt;<br>
&gt; Thanks!<br>
&gt;<br>
&gt; *** buffer overflow detected ***: grads terminated<br>
&gt; ======= Backtrace: =========<br>
&gt; /lib64/libc.so.6(__fortify_fail+0x37)[0x30d6f02827]<br>
&gt; /lib64/libc.so.6[0x30d6f00710]<br>
&gt; /lib64/libc.so.6[0x30d6effe0b]<br>
&gt; /lib64/libc.so.6(__snprintf_chk+0x7a)[0x30d6effcda]<br>
&gt; grads(gascat+0x250)[0x45d070]<br>
&gt; grads(gaplot+0x3fc)[0x45e27c]<br>
&gt; grads(gadspl+0x576)[0x4a68b6]<br>
&gt; grads(gacmd+0x1607)[0x4a5217]<br>
&gt; grads(gaexec+0x294)[0x4a6da4]<br>
&gt; grads(gacmd+0xea0)[0x4a4ab0]<br>
&gt; grads(main+0x5c7)[0x419477]<br>
&gt; /lib64/libc.so.6(__libc_start_main+0xfd)[0x30d6e1ed1d]<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Joowan Kim<br>
&gt; Atmospheric Chemistry Division<br>
&gt; National Center for Atmospheric Research<br>
&gt; PO Box 3000  Boulder, CO 80307<br>
&gt; e-mail : <a href="mailto:joowan@ucar.edu">joowan@ucar.edu</a><br>
&gt;          <a href="mailto:joowan.k@gmail.com">joowan.k@gmail.com</a><br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
--<br>
Jennifer M. Adams<br>
Center for Ocean-Land-Atmosphere Studies (COLA)<br>
111 Research Hall, Mail Stop 2B3<br>
George Mason University<br>
4400 University Drive<br>
Fairfax, VA 22030<br>
<br>
<br>
<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140224/57077992/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140224/57077992/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Mon, 24 Feb 2014 10:28:35 -0800<br>
From: Edward Pollock &lt;<a href="mailto:edwardp@ualberta.ca">edwardp@ualberta.ca</a>&gt;<br>
Subject: Re: [gradsusr] help with fndlvl function? (Freezing level<br>
        script)<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:CAALjUfhkwy1zxgjxLzCU4Fb3R8vZPSf1E3kFr2PXvYvJrAxFqQ@mail.gmail.com">CAALjUfhkwy1zxgjxLzCU4Fb3R8vZPSf1E3kFr2PXvYvJrAxFqQ@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Hi all,<br>
<br>
There have been a few enquiries recently about plotting freezing levels.<br>
<br>
For data on sigma levels with sigma level height given, one can adapt Bob<br>
Hart&#39;s script, <a href="http://isen.gs" target="_blank">isen.gs</a>, to do this. (See the script library on the grads<br>
website.)<br>
<br>
That script can also be adapted to plot any variable on any surface, such<br>
as height on potential vorticity, height on potential temperature, wind on<br>
potential temperature, or wind on potential vorticity etc.<br>
<br>
(Ideally, you should use do this calculation in height coordinates with a<br>
polynomial interpolation or Fourier series, but I think a linear<br>
interpolation between sigma levels works well if you have enough sigma<br>
levels in your data file.)<br>
<br>
I&#39;ve added my version of the script below. (It&#39;s set up to work with fwrite<br>
and several output files and times at the moment. The gen*ctl files are<br>
generic non-dimensional files which are needed for fwrite purposes. Feel<br>
free to modify the script as you please.)<br>
<br>
Cheers,<br>
<br>
Ted Pollock<br>
<br>
**********************<a href="http://frzlvl.gs" target="_blank">frzlvl.gs</a>*************************<br>
<br>
function frzlvl(args)<br>
<br>
dom=1<br>
while(dom&lt;=3)<br>
<br>
&#39;set gxout fwrite&#39;<br>
&#39;set fwrite FRZ.d&#39;dom&#39;.dat&#39;<br>
<br>
tt=1<br>
while(tt&lt;=48)<br>
<br>
&#39;open gen.d&#39;dom&#39;.3d.ctl&#39;<br>
&#39;open gen.d&#39;dom&#39;.2d.ctl&#39;<br>
<br>
&#39;set t &#39;tt<br>
&#39;set z 1 23&#39;<br>
&#39;define height=z+ter.2&#39;<br>
<br>
&#39;set z 1 23&#39;<br>
&#39;define temp=t&#39;<br>
&#39;set z 2 23&#39;<br>
&#39;define tempp=t(z-1)&#39;<br>
&#39;set z 1 22&#39;<br>
&#39;define tempm=t(z+1)&#39;<br>
<br>
&#39;define<br>
tabove=0.5*maskout(maskout(temp,temp-273.15),temp(z-1)-temp(z+1))+0.5*ma<br>
skout(maskout(temp,273.15-tempm),temp(z-1)-temp(z+1))&#39;<br>
&#39;define<br>
tbelow=0.5*maskout(maskout(temp,tempp-273.15),temp(z-1)-temp(z+1))+0.5*m<br>
askout(maskout(temp,273.15-temp),temp(z-1)-temp(z+1))&#39;<br>
<br>
&#39;define fabove=tabove*0+height&#39;<br>
&#39;define fbelow=tbelow*0+height&#39;<br>
<br>
&#39;set z 1&#39;<br>
<br>
&#39;define fabove=max(maskout(fabove,temp(z+0)-temp(z+1)),z=1,z=23)&#39;<br>
&#39;define fbelow=max(maskout(fbelow,temp(z+0)-temp(z+1)),z=1,z=23)&#39;<br>
&#39;define tabove=max(maskout(tabove,temp(z+0)-temp(z+1)),z=1,z=23)&#39;<br>
&#39;define tbelow=max(maskout(tbelow,temp(z+0)-temp(z+1)),z=1,z=23)&#39;<br>
<br>
&#39;define slope=(fabove-fbelow)/(tabove-tbelow)&#39;<br>
&#39;define b=fbelow-slope*tbelow&#39;<br>
&#39;define interp=slope*273.15+b&#39;<br>
<br>
&#39;d<br>
const(maskout(interp,interp),0,-u)+const(const(interp,0),1,-u)*(ter.2+((t-273<br>
.15)/0.0075))&#39;<br>
&#39;d ter.2&#39;<br>
<br>
<br>
&#39;close 2&#39;<br>
&#39;close 1&#39;<br>
<br>
tt=tt+1<br>
endwhile<br>
<br>
&#39;disable fwrite&#39;<br>
<br>
dom=dom+1<br>
endwhile<br>
<br>
&#39;quit&#39;<br>
<br>
return<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140224/322a86e0/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140224/322a86e0/attachment.html</a><br>
<br>
------------------------------<br>
<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>
<br>
End of gradsusr Digest, Vol 48, Issue 43<br>
****************************************<br>
</blockquote></div><br></div>