[gradsusr] Maskout two different fields

[Andrew Friedman]

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