[gradsusr] smoothed shadings in GRADS

Tue Jan 2 01:09:06 EST 2018

Dear Jennifer,

Thank you for this. I got it working already based on your suggestion.

Sincerely,

*Lyndon Mark P. Olaguera*
PhD Student
Laboratory of Climatology
Department of Geography
Faculty of Urban Environmental Science
Minami-Osawa Campus
Tokyo Metropolitan University
Hachioji Tokyo, Japan

On Thu, Dec 28, 2017 at 9:50 PM, Jennifer M Adams <jadams21 at gmu.edu> wrote:

> Yes, the missing grid points generated by the maskout() function are the
> issue here.
> > 'define sig=maskout(ttest1,ttest1-3.055)’
>  You can use ‘gxout grfill’ to show all the valid grid box values for the
> variable ’sig', or you can skip the maskout() call and use ‘gxout shaded’
> but omit drawing any contours less than 3.055 by using the color -1 in ’set
> ccols':
> > ‘set clevs 3.055’
> > ‘set ccols -1 15’
> > ‘d ttest1'
> Which would draw a gray shaded contour everywhere that ttest1 is >= 3.055
> —Jennifer
>
> > On Dec 28, 2017, at 3:00 AM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
> >
> > Hi Sir Andrew,
> >
> > The zero is actually important because this is a difference plot (see
> attached file). I am only using a single color to shade both negative and
> positive values.
> >
> > I'm not sure how to do your suggestion correctly.
> >
> > I use the const() function and set the missing values to a higher
> positive value (e.g. 1000) but the significant areas disappear.
> > Is this correct?
> >
> > 'define c=const(ave(b,lat=0,lat=18),1000,-u)'
> >
> >
> > Here's my script:
> >
> >
> > 'sdfopen ../omega_P45P50_1979-1993_filt.nc'
> > 'sdfopen ../omega_P45P50_1994-2008_filt.nc'
> > 'sdfopen ../../../winds/uwnd/yearly/diff_P45P50_uwnd_9408minus7993.nc'
> >
> > 'set dfile 2'
> > 'define n2=13'
> > 'set lev 1000 100'
> > 'set t 2'
> > 'define x2 = ave(omega,t=2,t=14)'
> > 'define s2 = sqrt(ave(pow(omega-x2,2),t=2,t=14)*(n2/(n2-1)))'
> >
> > 'set dfile 1'
> > 'define n1=13'
> > 'set lev 1000 100'
> > 'set t 2'
> > 'define x1 = ave(omega,t=2,t=14)'
> > 'define s1 = sqrt(ave(pow(omega-x1,2),t=2,t=14)*(n1/(n1-1)))'
> > 'define denom = sqrt((pow(s1,2)/n1)+(pow(s2,2)/n2))'
> > 'define num = x2-x1'
> > 'define ttest = num/denom'
> > 'define ttest1 = abs(ttest)'
> >
> > 'define sig=maskout(ttest1,ttest1-3.055)'
> >
> >
> > 'set lon 40 180'
> > 'define b=maskout(num,sig)'
> > 'set lev 1000 100'
> > 'set lon 40 180'
> > 'define c=ave(b,lat=0,lat=18)'
> > 'set lev 1000 100'
> > 'set lon 40 180'
> > 'define d=ave(num,lat=0,lat=18)'
> >
> > 'set display color white'
> > 'c'
> > 'set parea 1.5 10.5 1.5 8'
> > 'set lon 40 180'
> > 'set grads off'
> > 'set grid off'
> > 'set xlopts 1 5 0.2'
> > 'set ylopts 1 5 0.2'
> > 'set xlint 20'
> > 'set font 1'
> >
> > 'set gxout shaded'
> > 'set lev 1000 100'
> > 'set y 1'
> > 'set zlog on'
> > 'color -8 8 1 -kind gray->gray'
> > 'set csmooth on'
> > 'd -c*100'
> >
> > 'set dfile 3'
> > 'set gxout stream'
> > 'set strmden 7 0.001 0.1 2'
> > 'set cthick 5'
> > 'set lon 40 180'
> > 'set t 1'
> > 'set y 1'
> > 'define u=ave(uwnd,lat=0,lat=18)'
> > 'd u;-d*100'
> >
> > 'draw ylab Pressure(mb)'
> > 'draw title Difference in Zonal Circulation(EQ-18N)'
> > 'printim sig_omega_9408_minus_7993.png'
> > 'c'
> > 'quit'
> >
> >
> > Many thanks for the help.
> >
> > Sincerely,
> >
> >
> > Lyndon Mark P. Olaguera
> > PhD Student
> > Laboratory of Climatology
> > Department of Geography
> > Faculty of Urban Environmental Science
> > Minami-Osawa Campus
> > Tokyo Metropolitan University
> > Hachioji Tokyo, Japan
> >
> >
> > On Thu, Dec 28, 2017 at 4:12 PM, Andrew Friedman <andfried at gmail.com>
> wrote:
> > This is just a guess, but maybe it is because you can’t smooth across
> missing values?
> > What if you set your missing values to zero? Does that change anything?
> >
> > > On Dec 27, 2017, at 10:18 PM, Jeff Duda <jeffduda319 at gmail.com> wrote:
> > >
> > > I'm not sure what to say then. Perhaps the smth9() function only works
> on horizontal fields (that's all I've ever used it on)?
> > >
> > > Perhaps something very clever will be needed. Or you might have to
> perform the smoothing in an alternate data set outside of grads. Or you may
> be out of luck on this option. I'm not sure. Hopefully someone else can
> chime in with better advice.
> > >
> > > Good luck!
> > >
> > > Jeff
> > >
> > > On Wed, Dec 27, 2017 at 11:03 PM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
> > > Dear Sir Jeff,
> > >
> > > I also tried your suggestion but the result does not change.
> > >
> > > 'd smth9(smth9(smth9(smth9(smth9(smth9(smth9(smth9(-c*100))))))))'
> > >
> > > I am only shading the significant areas via the maskout() function.
> Maybe, this is the reason why Im not getting a smooth shading?
> > >
> > > Sincerely,
> > >
> > > Lyndon Mark P. Olaguera
> > > PhD Student
> > > Laboratory of Climatology
> > > Department of Geography
> > > Faculty of Urban Environmental Science
> > > Minami-Osawa Campus
> > > Tokyo Metropolitan University
> > > Hachioji Tokyo, Japan
> > >
> > >
> > > On Thu, Dec 28, 2017 at 2:44 PM, Jeff Duda <jeffduda319 at gmail.com>
> wrote:
> > > use the smth9() function. Nest several iterations of it to increase
> the smoothing level.
> > >
> > > Jeff Duda
> > >
> > > On Wed, Dec 27, 2017 at 9:41 PM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
> > > Hi Bill,
> > >
> > > I tried this but the result is still the same.
> > >
> > > Sincerely,
> > >
> > > Lyndon Mark P. Olaguera
> > > PhD Student
> > > Laboratory of Climatology
> > > Department of Geography
> > > Faculty of Urban Environmental Science
> > > Minami-Osawa Campus
> > > Tokyo Metropolitan University
> > > Hachioji Tokyo, Japan
> > >
> > >
> > > On Wed, Dec 27, 2017 at 9:16 PM, Bill Reilly <
> bill_reilly at compuserve.com> wrote:
> > > Before displaying the shading use the following:
> > >
> > > 'set csmooth on'
> > > --
> > >
> > > Bill Reilly
> > >
> > >
> > > (+34) 686-110-748
> > >
> > >
> > >
> > > bill at passageweather.com
> > > www.passageweather.com
> > >
> > >
> > >
> > > bill at brwebworks.com
> > > www.brwebworks.com
> > >
> > >
> > >
> > > On 27/12/2017 05:48, Lyndon Mark Olaguera wrote:
> > >> Dear GrADS experts,
> > >>
> > >> I am trying to plot a difference plot.In particular, a
> height-longitude cross section attached in this email. I am shading only
> the significant areas.
> > >>
> > >> Is there a way to smoothen the shading in the plot?Currently, the
> shading appears like "blocks".
> > >> Here's a part of my script.
> > >>
> > >> 'define sig=maskout(ttest1,ttest1-3.055)'
> > >> 'set lon 40 180'
> > >> 'define b=maskout(num,sig)'
> > >> 'set lev 1000 100'
> > >> 'set lon 40 180'
> > >> 'define c=ave(b,lat=0,lat=18)'
> > >>
> > >> 'set display color white'
> > >> 'c'
> > >> 'set parea 1.5 10.5 1.5 8'
> > >> 'set lon 40 180'
> > >> 'set grads off'
> > >> 'set grid off'
> > >> 'set xlopts 1 5 0.2'
> > >> 'set ylopts 1 5 0.2'
> > >> 'set xlint 20'
> > >> 'set font 1'
> > >>
> > >> 'set gxout shaded'
> > >> 'set lev 1000 100'
> > >> 'set y 1'
> > >> 'set zlog on'
> > >> 'color -8 8 1 -gxout shaded -kind gray->gray'
> > >> 'd -c*100'
> > >>
> > >> I'll appreciate any help.
> > >>
> > >> Sincerely,
> > >>
> > >>
> > >> Lyndon Mark P. Olaguera
> > >> PhD Student
> > >> Laboratory of Climatology
> > >> Department of Geography
> > >> Faculty of Urban Environmental Science
> > >> Minami-Osawa Campus
> > >> Tokyo Metropolitan University
> > >> Hachioji Tokyo, Japan
> > >>
> > >>
> > >>
> > >> ______________________________
> > >> _________________
> > >> gradsusr mailing list
> > >>
> > >> gradsusr at gradsusr.org
> > >> http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > >
> > > --
> > > Jeff Duda
> > > Post-doctoral research fellow
> > > University of Oklahoma School of Meteorology
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > >
> > > --
> > > Jeff Duda
> > > Post-doctoral research fellow
> > > University of Oklahoma School of Meteorology
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> > <sig_omega_9408_minus_7993.png>__________________________
> _____________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://secure-web.cisco.com/16TSWt--HQj5FmzwFVhGfzhNh2nTzXoL_
> qZNFM0hcdMNyPnRVe7jkDrLmoBOjAy-ah-kdFBkGrtlXj_fVTPYeGaEz8cZsL9B--
> lXtiaJrMc2xC0DmhUH9cdodpN17SJ_iKfcPq3JDM-qLlKWDoIT6nf5k7lPRjMh3ZF8TeIKF
> bMDtNjC3pRSzZB8gDnstLQixnJWmowQHj5BGa4PHJsmn2yrxzhYA7Aa4PFUt
> BRz7HT7iyWqC8Xoz4qKTNSfYw5a4XKMQjQyeMNgQIpf8ze9B0DBXfoZZ3Yv9
> gsDdXsT1P4MO6uzEJLOFDvBbH35DaHMYJujSLT_TIFPzhkekqwQdsXMPjSPOqATLarj2o
> W7fVgcqYnhLwZKulCsCwMUlABHEu4uTZ3ULukpYoQ_DIg6ChByynR8WEyvA_
> znOYmuKa3IrpB5m1iKtJWy_nVyW94JQ/http%3A%2F%2Fgradsusr.org%2Fmailman%
> 2Flistinfo%2Fgradsusr
>
> --
> Jennifer Miletta Adams
> Center for Ocean-Land-Atmosphere Studies (COLA)
> George Mason University
>
>
>
>
> _______________________________________________
> 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/20180102/f4b3c512/attachment-0001.html 