<div dir="ltr"><div><div>Daniel,<br></div>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. <br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 2, 2017 at 6:12 AM, Daniel Nairn <span dir="ltr">&lt;<a href="mailto:dannairn@gmail.com" target="_blank">dannairn@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small">No I don&#39;t think so, what I&#39;ve done is break down the script (as below) again and this now actually looks okay</div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><div class="gmail_default">&#39;define cape = capesfc&#39;</div><div class="gmail_default">&#39;define myhel = HLCY0_1000m&#39;</div><div class="gmail_default">&#39;define myval = maskout(cape,myhel-100)&#39;</div><div class="gmail_default">&#39;define myval2 = const(myval, 0, -u)&#39;</div><div class="gmail_default">&#39;d myval2&#39;</div><div class="gmail_default"><br></div><div class="gmail_default">I&#39;m still unsure as to what I&#39;m doing wrong when its all combined into one line, but at least this appears to be working. Thank you.</div></div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif;font-size:small"><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On 2 October 2017 at 11:52, Andrew Friedman <span dir="ltr">&lt;<a href="mailto:andfried@gmail.com" target="_blank">andfried@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Daniel,<br>
Could it be an issue with contour shading? If you pick a point where you know there is a value, is the output correct?<br>
<span class="m_-7879378552651008287HOEnZb"><font color="#888888">Andrew<br>
</font></span><div class="m_-7879378552651008287HOEnZb"><div class="m_-7879378552651008287h5"><br>
&gt; On Oct 2, 2017, at 9:00 AM, Daniel Nairn &lt;<a href="mailto:dannairn@gmail.com" target="_blank">dannairn@gmail.com</a>&gt; wrote:<br>
&gt;<br>
&gt; Hello Andrew<br>
&gt;<br>
&gt; Yes I simplified the script as suggested at that does look much better, thank you. Not quite sure what I&#39;ve done wrong with the const loops though.<br>
&gt;<br>
&gt; On 1 October 2017 at 21:55, Andrew Friedman &lt;<a href="mailto:andfried@gmail.com" target="_blank">andfried@gmail.com</a>&gt; wrote:<br>
&gt; Hi Daniel,<br>
&gt;<br>
&gt; 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.<br>
&gt;<br>
&gt; Andrew<br>
&gt;<br>
&gt; &gt; On Oct 1, 2017, at 12:59 PM, Daniel Nairn &lt;<a href="mailto:dannairn@gmail.com" target="_blank">dannairn@gmail.com</a>&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; Hello<br>
&gt; &gt;<br>
&gt; &gt; 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<br>
&gt; &gt;<br>
&gt; &gt; &#39;define cape = capesfc(e=2)&#39;<br>
&gt; &gt; &#39;define myhel = HLCY0_1000m(e=2)&#39;<br>
&gt; &gt; &#39;define capehel = const(const(maskout(cape,myhel<wbr>-100),1),0.0,-u)&#39;<br>
&gt; &gt; &#39;d capehel&#39;<br>
&gt; &gt;<br>
&gt; &gt; I&#39;m not convinced this is correct as my final image output does not show any CAPE values where expected. Can anybody help at all?<br>
&gt; &gt;<br>
&gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
<br>
<br>
______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
</div></div></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Jeff Duda<br>Post-doctoral research fellow<br>University of Oklahoma School of Meteorology<br></div></div></div></div></div>
</div>