<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>I have just committed changes to the obs_diag.f90 program that relate</div><div>to the creation of the rank histogram.</div><div><br></div><div>1) I added a namelist variable for obs_diag_nml:</div><div><br></div><div> outliers_in_histogram = .false.</div><div><br></div><div> .false. means that DART QC == 7 observations (those rejected by the outlier</div><div> threshold) will NOT be considered during the creation of the rank histogram.</div><div><br></div><div> .true. means they WILL be considered - and you should expect to see more mass</div><div> on the first and last bins ...</div><div><br></div><div>2) The matlab script 'DART/diagnostics/matlab/plot_rank_histogram.m' will read</div><div> the global attribute "outliers_in_histogram" and change the Y axis label</div><div> to reflect whether or not the outliers were used. If they ARE used, the label</div><div> will read: 'count - including outlier observations'; if they are not used, </div><div> the label will read: 'count'.</div><div><br></div><div>Happy Ranking ...</div><div><br></div><div>Tim</div><br><div><div>On Oct 8, 2010, at 4:54 PM, Tim Hoar wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>I extended obs_diag.f90 to also output rank histogram information if</div><div>there are any 'prior ensemble member xxx' copies in the obs_seq.final files.</div><div>(They exist if you specify input.nml:filter_nml:num_output_obs_members == <ens_size>.)</div><div><br></div><div>After you create the obs_diag_output.nc file, you can run </div><div><br></div><div>DART/diagnostics/matlab/plot_rank_histogram.m</div><div><br></div><div>to view the results. It has the 'expected' help file - available from the matlab</div><div>command line.</div><div><br></div><div>At present, ALL observations with DART QC 0,1,2,3, and 7 are used to create the</div><div>rank histogram. 7 is the outlier observation, so even these are being 'counted'.</div><div>The observation error -as specified in each observation- is used to generate</div><div>the sampling error signal for the ensemble members' expected observations that</div><div>define the bin edges for the rank histogram.</div><div><br></div><div>obs_diag.f90 has also been updated to provide better input guidance for the</div><div>instance when the user has specified a time window that does not intersect with</div><div>the observations in the files specified.</div><div><br></div><div>Please</div><div><br></div><div>svn update DART/diagnostics</div><div><br></div><div>rebuild your obs_diag, and let me know how I can improve it. I have one idea already,</div><div>so stay tuned for a minor tweak to the matlab script on Monday.</div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><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></div></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><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>