<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 & 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"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></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 'help' 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 "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today'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: "L.B." <<a href="mailto:bcbass2989@gmail.com">bcbass2989@gmail.com</a>><br>
Subject: Re: [gradsusr] Wind-chill parameter GrADS -<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:CAM1sYasmnA7TDdHaHAzqh1tn-4828YDjcY-wt1D9d2OR8PGzJA@mail.gmail.com">CAM1sYasmnA7TDdHaHAzqh1tn-4828YDjcY-wt1D9d2OR8PGzJA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<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 '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 <<a href="mailto:lap.bolsa@gmail.com">lap.bolsa@gmail.com</a>>wrote:<br>
<br>
> Hi Willians,<br>
><br>
> Look at this webpage-><br>
> <a href="http://www.lenntech.com/calculators/wind/wind-chill.htm" target="_blank">http://www.lenntech.com/calculators/wind/wind-chill.htm</a><br>
> 30 km/h (wind)<br>
> 18?C (temperature)<br>
> 12.1?C (windchill)<br>
><br>
> - - - - - - - - - - -<br>
><br>
> In our examples ....<br>
> I use your formulation...<br>
> - [1] I use your formulation corrected (convert KPH to MPH and K to F to<br>
> input into the equation)<br>
><br>
> 'far=((9/5)*t)-459.67'<br>
> 'mph=tg1*2.236936'<br>
> 'define chill =<br>
> 35.74+0.6215*far-35.75*pow(mph,0.16)+0.4275*far*pow(mph,0.16)'<br>
> 'chill2=(chill-32)/1.8'<br>
> 'd chill2'<br>
><br>
> Date / wind (MPH) / temp (?F) / windchill (?F)<br>
> 2014-FEB-19-13Z 17.1496 66.0072 64.8952<br>
> 2014-FEB-19-14Z 19.4474 66.7105 65.5746<br>
> 2014-FEB-19-15Z 19.4474 66.2803 65.0114<br>
><br>
><br>
> - [2] This is the other proposal that you have shown me ...<br>
><br>
> 'define chill =<br>
> 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)'<br>
><br>
> Date / wind (KPH) / temp (?C) / windchill (?C)<br>
> 2014-FEB-19-13Z 27.5996 18.8929 18.2662<br>
> 2014-FEB-19-14Z 31.2976 19.2836 18.6438<br>
> 2014-FEB-19-15Z 31.2976 19.0446 18.3308<br>
><br>
> I still think that with a wind of 30 km/h and a temperature of 18.8 ?C,<br>
> the windchill can not be a 'lower tenths'..<br>
> I'm confused, and I have no security ...<br>
> <a href="http://www.easysurf.cc/cnver16.htm" target="_blank">http://www.easysurf.cc/cnver16.htm</a><br>
> <a href="http://es.easycalculation.com/weather/wind-chill.php" target="_blank">http://es.easycalculation.com/weather/wind-chill.php</a><br>
><br>
> Best,<br>
><br>
> Leopoldo Alvarez<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>
-------------- 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 <<a href="mailto:joowan@ucar.edu">joowan@ucar.edu</a>><br>
Subject: [gradsusr] Scatter plot, buffer overflow<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Message-ID:<br>
<CAE7bxZSTFe_u12Nhk2ggi+RXhLQkZWhZgJ-37_JC6U4Tb1=<a href="mailto:UYA@mail.gmail.com">UYA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear Grads users<br>
<br>
I'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 <<a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a>><br>
Subject: Re: [gradsusr] Scatter plot, buffer overflow<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID: <<a href="mailto:EE3724F2-4EC7-4503-BA17-D26A7D37AE1A@cola.iges.org">EE3724F2-4EC7-4503-BA17-D26A7D37AE1A@cola.iges.org</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<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>
> Dear Grads users<br>
><br>
> I'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>
> _______________________________________________<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>
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 <<a href="mailto:edwardp@ualberta.ca">edwardp@ualberta.ca</a>><br>
Subject: Re: [gradsusr] help with fndlvl function? (Freezing level<br>
script)<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID:<br>
<<a href="mailto:CAALjUfhkwy1zxgjxLzCU4Fb3R8vZPSf1E3kFr2PXvYvJrAxFqQ@mail.gmail.com">CAALjUfhkwy1zxgjxLzCU4Fb3R8vZPSf1E3kFr2PXvYvJrAxFqQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<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'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've added my version of the script below. (It'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<=3)<br>
<br>
'set gxout fwrite'<br>
'set fwrite FRZ.d'dom'.dat'<br>
<br>
tt=1<br>
while(tt<=48)<br>
<br>
'open gen.d'dom'.3d.ctl'<br>
'open gen.d'dom'.2d.ctl'<br>
<br>
'set t 'tt<br>
'set z 1 23'<br>
'define height=z+ter.2'<br>
<br>
'set z 1 23'<br>
'define temp=t'<br>
'set z 2 23'<br>
'define tempp=t(z-1)'<br>
'set z 1 22'<br>
'define tempm=t(z+1)'<br>
<br>
'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))'<br>
'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))'<br>
<br>
'define fabove=tabove*0+height'<br>
'define fbelow=tbelow*0+height'<br>
<br>
'set z 1'<br>
<br>
'define fabove=max(maskout(fabove,temp(z+0)-temp(z+1)),z=1,z=23)'<br>
'define fbelow=max(maskout(fbelow,temp(z+0)-temp(z+1)),z=1,z=23)'<br>
'define tabove=max(maskout(tabove,temp(z+0)-temp(z+1)),z=1,z=23)'<br>
'define tbelow=max(maskout(tbelow,temp(z+0)-temp(z+1)),z=1,z=23)'<br>
<br>
'define slope=(fabove-fbelow)/(tabove-tbelow)'<br>
'define b=fbelow-slope*tbelow'<br>
'define interp=slope*273.15+b'<br>
<br>
'd<br>
const(maskout(interp,interp),0,-u)+const(const(interp,0),1,-u)*(ter.2+((t-273<br>
.15)/0.0075))'<br>
'd ter.2'<br>
<br>
<br>
'close 2'<br>
'close 1'<br>
<br>
tt=tt+1<br>
endwhile<br>
<br>
'disable fwrite'<br>
<br>
dom=dom+1<br>
endwhile<br>
<br>
'quit'<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>