[gradsusr] gradsusr Digest, Vol 48, Issue 43

Tue Feb 25 06:48:33 EST 2014

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 ???

Kunal Bali
Research Scholar
Radio & Atmospheric Science Division
CSIR - National Physical Laboratory
New Delhi - 110012

On Mon, Feb 24, 2014 at 11:58 PM, <gradsusr-request at gradsusr.org> wrote:

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