[gradsusr] Maskout two different fields

[Jennifer M Adams]

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