<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Just to put an example behind the email announcement of the Matlab script that</div><div>compares 'two' experments - I'll include the calling sequence and an example plot</div><div>that shows the results of two of Hui's ernesto experiments:</div><div><br></div><div>orient tall,</div><div>wysiwyg;</div><div>files = {'/fs/image/home/hliu/DART/models/wrf/work/ernesto/ctl/obs_diag_output.nc',</div><div><div> '/fs/image/home/hliu/DART/models/wrf/work/ernesto/airsQ/obs_diag_output.nc'};</div><div>titles = {'Control', 'airsQ'};</div><div>varnames = {'RADIOSONDE_U_WIND_COMPONENT', 'RADIOSONDE_TEMPERATURE'};</div><div>qtty = 'rmse'; % rmse, spread, totalspread, bias, etc.</div><div>prpo = 'both'; % prior & posterior</div><div> </div><div>two_experiments_profile(files, titles, varnames, qtty, prpo)</div></div><div><br></div><div>Please make note that the 'region' string (in this case 'Northern Hemisphere')</div><div>comes from whatever was used when 'obs_diag' was run. It this case, I would</div><div>have declared something like: </div><div><br></div><div><div>&obs_diag_nml</div><div> reg_names = 'Ernesto',</div></div><div><br></div><div><img height="945" width="724" apple-width="yes" apple-height="yes" src="cid:68887688-D8E8-4918-AD55-950C74702AD9@scd.ucar.edu"></div><div><br></div><div><br></div><div>Here's the enticement to keep reading. I called it 'two_experiments', but in reality,</div><div>it is N_experiments. As things get more cluttered, you can elect to plot only the</div><div>priors (or posteriors, if you want to incur the Wrath of Jeff). I put a lot of effort</div><div>into improving the reliability and labeling of the two axes (nobs, and 'rmse').</div><div>Each experiment has its own nobs 'possible' and 'used' that are the same color as</div><div>the corresponding line.</div><div><br></div><div>I am in the process of putting a little watermark on the bottom with the filenames</div><div>that are being plotted. At present, I'm not automatically printing these things.</div><div>I might add an option ...</div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="font-family: Helvetica; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Tim Hoar, Associate Scientist</div><div>National Center for Atmospheric Research</div><div><a href="mailto:thoar@ucar.edu">thoar@ucar.edu</a></div><div>303 497 1708</div></div></span> </div><br></body></html>