<div dir="ltr"><div><div><div><div>From my experience with GEFS data, each ensemble member is contained in a separate file. So I hope your control file uses templating and has an EDEF entry. If it doesn&#39;t, or you don&#39;t have all the data files, then that explains why your &#39;define ensmean=ave(t850,e=1,e=21)&#39; command is just giving you the control member (because it&#39;s the only one present).<br><br></div>I assume by &quot;anomaly&quot; you mean the difference between a given ensemble member and the ensemble mean. In that case, just plot the difference field between an ensemble member and the ensemble mean. For instance,<br><br></div>&#39;set e 2&#39;<br></div>&#39;d tmpprs - ensmean&#39;<br><br></div>Jeff Duda<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 20, 2016 at 8:05 AM, Chorley Weather <span dir="ltr">&lt;<a href="mailto:weatherstu@chorleyweather.com" target="_blank">weatherstu@chorleyweather.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>
<div style="font-size:10pt;font-family:Tahoma,Arial,Helvetica,sans-serif">
<p>Hi all,</p>
<p> </p>
<p>I have a couple of questions, firstly I am trying get the mean run in a spaghetti graph of the GEFS but I seem to keep getting the control run. Here is the script:</p>
<p>&#39;set t 1 65&#39;<br>&#39;set e 1&#39;<br>&#39;set csmooth on&#39;<br>&#39;set cmark 0&#39;<br>&#39;set lat 53&#39;;&#39;set lon -2.23&#39;<br>&#39;set vrange -30 30&#39;<br>&#39;define t850=tmpprs(lev=850)-273.16&#39;<br>&#39;define ensmean=ave(t850,e=1,e=21)&#39;<br>&#39;set gxout line&#39;<br>&#39;set ccolor 2&#39;<br>&#39;set cthick 8&#39;<br>&#39;d ensmean&#39;</p>
<p> </p>
<p>Second question, using the GFS how do I produce anomaly charts with this script, im not sure how to define it.....</p>
<p>&#39;set t &#39;i<br>&#39;set grads off&#39;<br>&#39;set csmooth on&#39;<br>&#39;define geopotdam = hgtprs(lev=500)/10&#39;<br>&#39;color -gxout shaded -30 30 1 -kind blue-&gt;white-&gt;red-&gt; -xcbar 10.1 10.3 0.7 7.74 -fs 10&#39;</p>
<p> </p>
<p>Regards,</p><span class="HOEnZb"><font color="#888888">
<div>-- <br>
<p>Stuart Markham</p>
<p><img src="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG" alt=""></p>
</div>
</font></span></div>
<br>_______________________________________________<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/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br></div>
</div>