<html><body><div style="font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 10pt; color: #000000"><div>Hi</div><div><br></div><div>You need to open 2 files (3 if you want the control to be yet a different curve)</div><div><br></div><div>in my scripts I do it this way</div><div><br></div><div> ///</div><div> 'sdfopen <a href="http://nomads.ncep.noaa.gov:9090/dods/gens/gens'date'/gep_all_'hour'z'" data-mce-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'" data-mce-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' " data-mce-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><a href="http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z' " data-mce-href="http://nomads.ncep.noaa.gov:9090/dods/gfs_0p25/gfs'date'/gfs_0p25_'hour'z' "></a><br data-mce-bogus="1"></div><div> ///</div><div><br></div><div><a></a>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><br></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><br></div><div>Then you can play with the data.</div><div><br></div><div>eg.</div><div><br></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><br></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><br></div><div>Do this for all the members of choice.</div><div><br></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><br></div><div> ///</div><div> 'define mean = (temp2m_mem2 + temp2m_mem3 + ... +temp2m_mem21)/20'</div><div> 'd mean'</div><div> ///</div><div><br></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><br></div><div> ///</div><div> 'define temp2_oper = tmp2m.3-273.15'</div><div> 'd temp2_oper'</div><div> ///</div><div><br></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><br></div><div><img src="cid:c6761b73ac260ba8cbad72b7c7c4c08f36a10f71@zimbra" data-mce-src="cid:c6761b73ac260ba8cbad72b7c7c4c08f36a10f71@zimbra" height="307" width="385"></div><div><br></div><div>Best regards</div><div><br></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;" data-mce-style="color: #000; font-weight: normal; font-style: normal; text-decoration: none; font-family: Helvetica,Arial,sans-serif; font-size: 12pt;"><b>Van: </b>"Chorley Weather" <weatherstu@chorleyweather.com><br><b>Aan: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br><b>Verzonden: </b>Woensdag 10 februari 2016 13:13:50<br><b>Onderwerp: </b>[gradsusr] Ensemble Graphs<br><div><br></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="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG" alt="" src="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG" data-mce-src="http://www.chorleyweather.com/forecast-charts/SampleDatasets/logo%20450x100.PNG"></p></div><br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr</div><div><br></div></div></body></html>