<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body style='font-size: 10pt; font-family: Tahoma,Arial,Helvetica,sans-serif'>
<p>Thank you,</p>
<p>That's fine, the errors I am getting relate to the args. Error: Non numeric args to numeric operation. Error line 16. Here is the section of script</p>
<p>function main(args)</p>
<p>* Parse the arguments: date, hour<br />if (args = '')<br />prompt 'Enter forecast date (example, 20110717) --> ' <br />pull date<br />prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> ' <br />pull hour<br />else <br />date = subwrd(args,1)<br />hour = subwrd(args,2)<br />endif</p>
<p><br />'reinit'<br />'sdfopen http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gep_all_'hour'z' *All GEFS members<br />'sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z' *GFS oper</p>
<p>* Get the Time axis info<br />tsize = subwrd(_tdef,2)<br />_t1 = 1 ;* 2nd half of timeseries<br />_t2 = 45<br />tsize = _t2 - _t1 + 1<br />'set t '_t1' '_t2<br />'q dims'<br />times = sublin(result,5)<br />_time1 = subwrd(times,6) <br />_time2 = subwrd(times,8)<br />_tdim = _time1' '_time2</p>
<p>tincr = subwrd(_tdef,5)<br />_tdef = 'tdef 'tsize' linear '_time1' 'tincr</p>
<p>huh = subwrd(_tdef,4)</p>
<p><br />* Get info from the descriptor file<br />'q ctlinfo'<br />_ctl = result<br />_undef = getctl(undef)<br />_tdef = getctl(tdef)<br />_zdef = getctl(zdef)</p>
<p>maps = 3</p>
<p>* Get the Time axis info<br />tsize = subwrd(_tdef,2)<br />_t1 = 1 ;* 2nd half of timeseries<br />_t2 = 45<br />tsize = _t2 - _t1 + 1<br />'set t '_t1' '_t2<br />'q dims'<br />times = sublin(result,5)<br />_time1 = subwrd(times,6) <br />_time2 = subwrd(times,8)<br />_tdim = _time1' '_time2</p>
<p>tincr = subwrd(_tdef,5)<br />_tdef = 'tdef 'tsize' linear '_time1' 'tincr</p>
<p>huh = subwrd(_tdef,4)</p>
<p> </p>
<p>Stuart.</p>
<p> </p>
<p> </p>
<p>On 11-02-2016 10:46, Chorley Weather wrote:</p>
<blockquote type="cite" style="padding-left:5px; border-left:#1010ff 2px solid; margin-left:5px"><!-- html ignored --><!-- head ignored --><!-- meta ignored -->
<p>Hi, </p>
<p>Thanks for the detailed information, I have tried these ways of getting them to run together but getting a few errors, would be be possible to see one of your scripts. I understand they maybe private so perhaps email me direct at weatherstu@chorleyweather.com. </p>
<p>Many thanks</p>
<p>Stuart.</p>
<p>On 11-02-2016 08:40, dzengiz.tafa@telenet.be wrote:</p>
<blockquote style="padding-left: 5px; border-left: #1010ff 2px solid; margin-left: 5px;"><!-- node type 8 --><!-- node type 8 --><!-- node type 8 -->
<div style="font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 10pt; color: #000000;">
<div>Hi</div>
<div> </div>
<div>You need to open 2 files (3 if you want the control to be yet a different curve)</div>
<div> </div>
<div>in my scripts I do it this way</div>
<div> </div>
<div> ///</div>
<div> 'sdfopen <a href="http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gep_all_'hour'z'">http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gep_all_'hour'z'</a> *All GEFS members<br /> 'sdfopen <a href="http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gec00_'hour'z'">http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gec00_'hour'z'</a> * GEFS control<br /> 'sdfopen <a href="http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z'">http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z'</a><a> *GFS oper </a></div>
<div> ///</div>
<div> </div>
<div>Just replace date with a date string like 20160211 and hour like 00 (run of 00Z). If you are not using opendap, just use the 'open' command, followed by the ctl file you created using grib2ctl & gribmap</div>
<div> </div>
<div>Then the trick is to use the 3 files that are opened. They are opened in the order provided so</div>
<div>* The file-connection with all the members is file number 1</div>
<div>* The file-connection with the control is file number 2 &</div>
<div>* The file-connection with the operational forecast is file number 3</div>
<div> </div>
<div>Then you can play with the data.</div>
<div> </div>
<div>eg.</div>
<div> </div>
<div>Set a member value (1 to 20) and when defining a variable, use the grib var name followed by the opendap connection, or file number. In this case 1, because this is the file connection containing all the members</div>
<div> </div>
<div> ///</div>
<div> * Member 2, (member 1 is the control in the opendap)</div>
<div> 'set e 2'</div>
<div> 'define temp2m_mem2 = tmp2m.1-273.15'</div>
<div> 'd temp2m_mem1'</div>
<div> ///</div>
<div> </div>
<div>Do this for all the members of choice.</div>
<div> </div>
<div>For the mean, just calculate a mean value, summing up all the seperate members & deviding them by the amount of members used in the sum.</div>
<div> </div>
<div> ///</div>
<div> 'define mean = (temp2m_mem2 + temp2m_mem3 + ... +temp2m_mem21)/20'</div>
<div> 'd mean'</div>
<div> ///</div>
<div> </div>
<div>If you want to use the Operational run data or the Control data as used in this example, you can use replace the .1 as in the member example with .2 or .3, according to your files being opened in your version of the script. In the example provided this translates into</div>
<div> </div>
<div> ///</div>
<div> 'define temp2_oper = tmp2m.3-273.15'</div>
<div> 'd temp2_oper'</div>
<div> ///</div>
<div> </div>
<div>As an attachment, I provided the output I get using this method & it pretty much resembles what you would like to create. Hope this helps.</div>
<div> </div>
<div><img src="cid:145519296256bc7b829f090873090022@chorleyweather.com" alt="" width="385" height="307" /></div>
<div> </div>
<div>Best regards</div>
<div> </div>
<hr id="zwchr" />
<div style="color: #000; font-weight: normal; font-style: normal; text-decoration: none; font-family: Helvetica,Arial,sans-serif; font-size: 12pt;"><strong>Van: </strong>"Chorley Weather" <weatherstu@chorleyweather.com><br /><strong>Aan: </strong>"GrADS Users Forum" <gradsusr@gradsusr.org><br /><strong>Verzonden: </strong>Woensdag 10 februari 2016 13:13:50<br /><strong>Onderwerp: </strong>[gradsusr] Ensemble Graphs<br />
<div> </div>
<p>Hi all, </p>
<p>I am in the process of developing GEFS ensemble graphs similar to the image, however, I am stuck on how to run the GEFS mean (red line) and also I'm unsure on how to run the GFS (Black Line) along side the GEFS. Can anybody shed some light on this for me please.....Many thanks in advance. </p>
<p> </p>
<p>******Email re sent.... forgot to attach image******</p>
<div>-- <br />
<p>Stuart Markham</p>
<p><img src="cid:145519296256bc7b829f499819784821@chorleyweather.com" alt="" /></p>
</div>
<br />_______________________________________________<br />gradsusr mailing list<br />gradsusr@gradsusr.org<br />http://gradsusr.org/mailman/listinfo/gradsusr</div>
<div> </div>
</div>
<!-- node type 8 --><br />
<pre>_______________________________________________
gradsusr mailing list
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<p> </p>
<div>-- <br />
<p>Stuart Markham</p>
<p><img src="cid:145519296256bc7b829f804697897837@chorleyweather.com" alt="" /></p>
</div>
</blockquote>
<p> </p>
<div>-- <br />
<p>Stuart Markham</p>
<p><img src="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG" alt="" /></p>
</div>
</body></html>