[Dart-dev] DART/branches Revision: 12992

[dart at ucar.edu]

thoar at ucar.edu
2019-03-01 11:37:12 -0700 (Fri, 01 Mar 2019)
263
The 'standard' DART colorscheme is now implemented.

obs counts are red, prior info is green, posterior info blue
cleanly supports observation diagnostics that do not have posterior info.
Doesn't really do anything with many of the DART QC values, including #8.




Modified: DART/branches/recam/diagnostics/matlab/plot_profile.m
===================================================================
--- DART/branches/recam/diagnostics/matlab/plot_profile.m	2019-02-28 00:23:40 UTC (rev 12991)
+++ DART/branches/recam/diagnostics/matlab/plot_profile.m	2019-03-01 18:37:12 UTC (rev 12992)
@@ -3,8 +3,8 @@
 % Part of the observation-space diagnostics routines.
 %
 % 'obs_diag' produces a netcdf file containing the diagnostics.
-% obs_diag condenses the obs_seq.final information into summaries for a few specified
-% regions - on a level-by-level basis.
+% 'obs_diag' condenses the obs_seq.final information into summaries for a few
+% specified regions - on a level-by-level basis.
 %
 % The number of observations possible reflects only those observations
 % that have incoming QC values of interest. Any observation with a DART
@@ -45,7 +45,11 @@
 %
 % fname   = 'obs_diag_output.nc';
 % copy    = 'totalspread';
-% plotdat = plot_profile(fname, copy, 'obsname', 'RADIOSONDE_TEMPERATURE');
+% obsname = 'RADIOSONDE_TEMPERATURE';
+% plotdat = plot_profile(fname, copy, 'obsname', obsname);
+%
+% EXAMPLE 3: You can tell it to be quiet, too.
+% plotdat = plot_profile(fname, copy, 'verbose', 'no');
 
 %% DART software - Copyright UCAR. This open source software is provided
 % by UCAR, "as is", without charge, subject to all terms of use at
@@ -58,6 +62,7 @@
 %---------------------------------------------------------------------
 
 default_obsname = 'none';
+default_verbosity = 'yes';
 p = inputParser;
 
 addRequired(p,'fname', at ischar);
@@ -64,8 +69,10 @@
 addRequired(p,'copy', at ischar);
 if (exist('inputParser/addParameter','file') == 2)
     addParameter(p,'obsname',default_obsname, at ischar);
+    addParameter(p,'verbose',default_verbosity, at ischar);
 else
     addParamValue(p,'obsname',default_obsname, at ischar);
+    addParamValue(p,'verbose',default_verbosity, at ischar);
 end
 
 p.parse(fname, copy, varargin{:});
@@ -80,6 +87,13 @@
     disp(p.Unmatched)
 end
 
+global verbose
+if (strncmpi(p.Results.verbose,'y',1))
+   verbose = 1;
+else
+   verbose = 0;
+end
+
 if (exist(fname,'file') ~= 2)
     error('file/fname <%s> does not exist',fname)
 end
@@ -135,11 +149,20 @@
 
 plotdat.nvars         = length(plotdat.varnames);
 plotdat.copyindex     = get_copy_index(fname,copy);
+plotdat.rmseindex     = get_copy_index(fname,'rmse');
 plotdat.Npossindex    = get_copy_index(fname,'Nposs');
 plotdat.Nusedindex    = get_copy_index(fname,'Nused');
+plotdat.NQC0index     = get_copy_index(fname,'N_DARTqc_0');
+plotdat.NQC1index     = get_copy_index(fname,'N_DARTqc_1');
+plotdat.NQC2index     = get_copy_index(fname,'N_DARTqc_2');
+plotdat.NQC3index     = get_copy_index(fname,'N_DARTqc_3');
+plotdat.NQC4index     = get_copy_index(fname,'N_DARTqc_4');
 plotdat.NQC5index     = get_copy_index(fname,'N_DARTqc_5');
 plotdat.NQC6index     = get_copy_index(fname,'N_DARTqc_6');
+plotdat.NQC7index     = get_copy_index(fname,'N_DARTqc_7');
+plotdat.NQC8index     = get_copy_index(fname,'N_DARTqc_8','fatal','no');
 
+global figuredata
 figuredata = setfigure();
 
 %%---------------------------------------------------------------------
@@ -158,24 +181,24 @@
 end
 
 for ivar = varlist
-    
+
     % create the variable names of interest.
-    
+
     plotdat.myvarname = plotdat.varnames{ivar};
     plotdat.guessvar  = sprintf('%s_VPguess',plotdat.varnames{ivar});
     plotdat.analyvar  = sprintf('%s_VPanaly',plotdat.varnames{ivar});
-    
+
     plotdat.trusted   = nc_read_att(fname, plotdat.guessvar, 'TRUSTED');
     if (isempty(plotdat.trusted)), plotdat.trusted = 'NO'; end
-    
+
     % get appropriate vertical coordinate variable