[gradsusr] Maskout two different fields

[Jeff Duda]

Daniel,
Why are you applying the const() function in this problem? From your code,
myval2 should be a field of all zeros except for in the masked area. I
would think you would want to plot myval only. If myval looks wrong then
either your cape or helicity fields are full of bad data. I would display
each separately to confirm if that is a problem.

Jeff Duda

On Mon, Oct 2, 2017 at 6: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
>
>


-- 
Jeff Duda
Post-doctoral research fellow
University of Oklahoma School of Meteorology
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171002/abd94d37/attachment-0001.html