<div dir="ltr"><div>Anomalies can be defined in many ways. Probably the common one is a climatological anomaly. In that case you&#39;ll need to obtain data containing climatological mean fields. The anomaly in that case is just [field - climatological mean]. I do not know off the top of my head where such fields can be found, but there are websites that provide anomaly graphics, so someone out there knows where you can find those data.<br><br></div>Jeff<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 20, 2016 at 9:50 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 jeff,</p>
<p>I was using this grads server link <span><a href="http://nomads.ncep.noaa.gov:9090/dods/gens/gens20160120/gep_all_00z" target="_blank">http://nomads.ncep.noaa.gov:9090/dods/gens/gens20160120/gep_all_00z</a></span> but ill download the grib files and try it that way. As for the &#39;Anomaly&#39; question I was actually on about GFS data. I would like to produce a 500mb height anomaly like the image attached but dint know how to define this. </p>
<p> </p>
<p>regards,</p>
<p>Stuart.</p>
<p> </p>
<p>On 20-01-2016 15:10, Jeff Duda wrote:</p>
<blockquote type="cite" style="padding-left:5px;border-left:#1010ff 2px solid;margin-left:5px">
<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;</div>
&#39;d tmpprs - ensmean&#39;<br><br></div>
Jeff Duda</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Wed, Jan 20, 2016 at 8:05 AM, Chorley Weather <span>&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"><span style="text-decoration:underline"></span>
<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>
<div>-- <br>
<p>Stuart Markham</p>
<p><img src="http://./?_task=mail&amp;_id=1101846189569faa569e38a&amp;_action=display-attachment&amp;_file=rcmfile30871453304407050656400" alt=""></p>
</div>
</div>
<br>_______________________________________________<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" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br> </blockquote>
</div>
<br><br clear="all"><span class="HOEnZb"><font color="#888888"><br>-- <br>
<div>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms</div>
</font></span></div><span class="HOEnZb"><font color="#888888">
<br>
<pre>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<p> </p>
<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>
</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>