[Dart-dev] [4530] DART/trunk/diagnostics/matlab/plot_rank_histogram.m: Added logic to change Y label based on whether or not the outlier observations

nancy at ucar.edu nancy at ucar.edu
Fri Oct 15 12:27:21 MDT 2010


Revision: 4530
Author:   thoar
Date:     2010-10-15 12:27:21 -0600 (Fri, 15 Oct 2010)
Log Message:
-----------
Added logic to change Y label based on whether or not the outlier observations 
were included in the calculation of the rank histogram.

Modified Paths:
--------------
    DART/trunk/diagnostics/matlab/plot_rank_histogram.m

-------------- next part --------------
Modified: DART/trunk/diagnostics/matlab/plot_rank_histogram.m
===================================================================
--- DART/trunk/diagnostics/matlab/plot_rank_histogram.m	2010-10-15 18:13:47 UTC (rev 4529)
+++ DART/trunk/diagnostics/matlab/plot_rank_histogram.m	2010-10-15 18:27:21 UTC (rev 4530)
@@ -98,6 +98,8 @@
 plotdat.latlim1            = nc_attget(fname, nc_global, 'latlim1');
 plotdat.latlim2            = nc_attget(fname, nc_global, 'latlim2');
 plotdat.biasconv           = nc_attget(fname, nc_global, 'bias_convention');
+plotdat.outlierstring      = nc_attget(fname, nc_global, 'outliers_in_histogram');
+plotdat.QCsusedstring      = nc_attget(fname, nc_global, 'DART_QCs_in_histogram');
 
 % Make sure the time index makes sense.
 
@@ -301,8 +303,13 @@
    set(h,'HorizontalAlignment','center')
 
    xlabel({'Observation Rank (among ensemble members)',obsstring})
-   ylabel('count')
 
+   if ( strcmp(plotdat.outlierstring,  'TRUE'))
+      ylabel('count - including outlier observations')
+   else
+      ylabel('count')
+   end
+
    % more annotation ...
 
    if (plotdat.region == 1)


More information about the Dart-dev mailing list