<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>To draw the box and whiskers, use 'set gxout bar' and 'set gxout errbar'. Then give the display command two variables, separated by a a semi-colon. Here's the script that I used to draw the first four plots in the ensemble documentation: </div><div><br></div><div><div>'open gfsens.2009010100.ctl'</div><div>'set lon -60'</div><div>'set lat 45'</div><div>'set lev 500'</div><div>'set parea 1 9.7 2.5 6'</div><div>'set t 0.5 65.5'</div><div>'set e 1'</div><div>'clear'</div><div>'set grads off'</div><div>'set tlsupp year'</div><div>'set cmark 0'</div><div>'set vrange 4950 5600'</div><div>'set ylint 100'</div><div>'d hgt'</div><div>'printim edemo1.png '</div><div><br></div><div>* Remaining members</div><div>e=2</div><div>cols='2 3 4 5 6 8 9 10 11 12 13 14 2 3 4 5 6 8 9 10 11 12 13 14 2 3 4 5 6 7 8 9 10 11 12 13 14 '</div><div>while (e<=21)</div><div> 'set e 'e</div><div> c=subwrd(cols,e)</div><div> 'set ccolor 'c</div><div> 'set cmark 0'</div><div> 'd hgt'</div><div> e=e+1</div><div>endwhile</div><div>'printim edemo2.png '</div><div><br></div><div>* all 21 members drawn as shaded boxes</div><div>'clear'</div><div>'set gxout grfill'</div><div>'set e 0.5 21.5'</div><div>'set clevs 5100 5150 5200 5250 5300 5350 5400 5450 '</div><div>'set ccols 14 4 11 5 3 7 12 8 2 '</div><div>'set grads off'</div><div>'set tlsupp year'</div><div>'d hgt'</div><div>'cbarn 0.8 1'</div><div>'printim edemo3.png'</div><div><br></div><div>* Calculate the ensemble mean </div><div>'set e 1'</div><div>'set t 1 last'</div><div>'define ensmean=ave(hgt,e=1,e=21)'</div><div><br></div><div>* Calculate the variance</div><div>diffsq = 'pow(hgt-ensmean,2)'</div><div>variance = 'ave('diffsq',e=1,e=21)'</div><div>'define stddev=sqrt('variance')'</div><div><br></div><div>* Calculate the min/max </div><div>'define ensmin=tloop(min(hgt,e=1,e=21))'</div><div>'define ensmax=tloop(max(hgt,e=1,e=21))'</div><div><br></div><div>* Plot the results</div><div>'clear'</div><div>'set t 0.5 65.5'</div><div>'set vrange 4950 5600'</div><div>'set ylint 100'</div><div>'set tlsupp year'</div><div>'set grads off'</div><div><br></div><div><div>* Draw error bars for min/max</div><div>'set gxout errbar'</div><div>'set ccolor 4'</div><div>'d ensmin;ensmax'</div><div><br></div><div>* Draw bars for +/- standard deviation</div><div>plus = '(ensmean+stddev)'</div><div>minus = '(ensmean-stddev)'</div><div>'set gxout bar'</div><div>'set bargap 50'</div><div>'set baropts outline'</div><div>'set ccolor 3'</div><div>'d 'minus';'plus</div><div><br></div><div>* Draw line for Ensemble mean </div><div>'set gxout line'</div><div>'set cmark 0'</div><div>'set cthick 6'</div><div>'set digsiz 0.05'</div><div>'set ccolor 2'</div><div>'d ensmean'</div><div><br></div><div>'printim edemo4.png'</div><div><br></div></div></div><div><br></div><br><div><div>On Jan 21, 2009, at 7:39 PM, Hsin-Hsing,Chia wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hi Jennifer:<br> Thanks for your effort on the documentation about the ensemble dataset.<br>How to plot the figure 4, which has the blue whiskers in this documentation ?<br><br> Thanks<br><br>Hsin-Hsing<br></div></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><div>--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div><div><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></span></div></span> </div><br></body></html>