[gradsusr] Maskout two different fields

[Daniel Nairn]

Jennifer - perfect, thank you! Not sure how I missed the extra const but
that fixed it.

Jeff - thanks for the reply.

On 2 October 2017 at 17:14, Jennifer M Adams <jadams21 at gmu.edu> wrote:

> Hi, Daniel —
> Your original expression:
>
> > > > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)’
>
> has an extra const() in there that changes the valid cape values to 1
> before changing the undef values to 0.0. It should be:
>
>   'define capehel = const(maskout(cape,myhel-100),0.0,-u)’
>
> But an even simpler way to do this is to use a logical expression and the
> if() function:
> 'define capehel = if(hlcy>=100,cape,0.0)'
>
> —Jennifer
>
> > On Oct 2, 2017, at 7:12 AM, Daniel Nairn <dannairn at gmail.com> wrote:
> >
> > No I don't think so, what I've done is break down the script (as below)
> again and this now actually looks okay
> >
> > 'define cape = capesfc'
> > 'define myhel = HLCY0_1000m'
> > 'define myval = maskout(cape,myhel-100)'
> > 'define myval2 = const(myval, 0, -u)'
> > 'd myval2'
> >
> > I'm still unsure as to what I'm doing wrong when its all combined into
> one line, but at least this appears to be working. Thank you.
> >
> >
> >
> > On 2 October 2017 at 11:52, Andrew Friedman <andfried at gmail.com> wrote:
> > Hi Daniel,
> > Could it be an issue with contour shading? If you pick a point where you
> know there is a value, is the output correct?
> > Andrew
> >
> > > On Oct 2, 2017, at 9:00 AM, Daniel Nairn <dannairn at gmail.com> wrote:
> > >
> > > Hello Andrew
> > >
> > > Yes I simplified the script as suggested at that does look much
> better, thank you. Not quite sure what I've done wrong with the const loops
> though.
> > >
> > > On 1 October 2017 at 21:55, Andrew Friedman <andfried at gmail.com>
> wrote:
> > > Hi Daniel,
> > >
> > > What if you simplify it by just looking at the inner loop:
> maskout(cape,myhel-100)? Does the output look correct? Then you can build
> on it by adding the successive const loops.
> > >
> > > Andrew
> > >
> > > > On Oct 1, 2017, at 12:59 PM, Daniel Nairn <dannairn at gmail.com>
> wrote:
> > > >
> > > > Hello
> > > >
> > > > I am trying to plot a chart (from one ensemble member) that displays
> SBCAPE where 0-1km helicity is greater than 100. This is how my current
> script looks
> > > >
> > > > 'define cape = capesfc(e=2)'
> > > > 'define myhel = HLCY0_1000m(e=2)'
> > > > 'define capehel = const(const(maskout(cape,myhel-100),1),0.0,-u)'
> > > > 'd capehel'
> > > >
> > > > I'm not convinced this is correct as my final image output does not
> show any CAPE values where expected. Can anybody help at all?
> > > >
> > > > _______________________________________________
> > > > 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
> >
> >
> > _______________________________________________
> > 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
>
> --
> 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/20171002/fe67961e/attachment.html