[Dart-dev] DART/branches Revision: 12994

[dart at ucar.edu]

thoar at ucar.edu
2019-03-01 15:54:48 -0700 (Fri, 01 Mar 2019)
68
works with obs_diag_output.nc files that have no posterior values.




Modified: DART/branches/recam/diagnostics/matlab/plot_rmse_xxx_evolution.m
===================================================================
--- DART/branches/recam/diagnostics/matlab/plot_rmse_xxx_evolution.m	2019-03-01 21:05:02 UTC (rev 12993)
+++ DART/branches/recam/diagnostics/matlab/plot_rmse_xxx_evolution.m	2019-03-01 22:54:48 UTC (rev 12994)
@@ -73,6 +73,7 @@
 
 default_level = -1;
 default_obsname = 'none';
+default_verbosity = 'yes';
 default_range = [NaN NaN];
 p = inputParser;
 
@@ -82,10 +83,12 @@
     addParameter(p,'obsname',default_obsname, at ischar);
     addParameter(p,'range',default_range, at isnumeric);
     addParameter(p,'level',default_level, at isnumeric);
+    addParameter(p,'verbose',default_verbosity, at ischar);
 else
     addParamValue(p,'obsname',default_obsname, at ischar);
     addParamValue(p,'range',default_range, at isnumeric);
     addParamValue(p,'level',default_level, at isnumeric);
+    addParamValue(p,'verbose',default_verbosity, at ischar);
 end
 p.parse(fname, copy, varargin{:});
 
@@ -101,6 +104,13 @@
     disp(p.Unmatched)
 end
 
+global verbose
+if (strncmpi(p.Results.verbose,'y',1))
+   verbose = 1;
+else
+   verbose = 0;
+end
+
 if (numel(p.Results.range) ~= 2)
     error('range must be an array of length two ... [bottom top]')
 end
@@ -129,8 +139,8 @@
 plotdat.plevel_edges  = local_ncread(fname,'plevel_edges');
 plotdat.hlevel        = local_ncread(fname,'hlevel');
 plotdat.hlevel_edges  = local_ncread(fname,'hlevel_edges');
-[plotdat.ncopies,~]   = nc_dim_info(fname,'copy');
-[plotdat.nregions,~]  = nc_dim_info(fname,'region');
+[plotdat.ncopies, ~]  = nc_dim_info(fname,'copy');
+[plotdat.nregions, ~] = nc_dim_info(fname,'region');
 plotdat.region_names  = strtrim(ncread(fname,'region_names')');
 
 dimensionality        = nc_read_att(fname, '/', 'LocationRank');
@@ -186,17 +196,11 @@
 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');
 
-figuredata = setfigure();
+global figuredata
+figuredata = setfigure('landscape');
 
-prior_green = [  0/255 128/255   0/255];
-poste_blue  = [  0/255   0/255 255/255];
-obs_red     = [215/255  10/255  83/255];
-
-plotdat.ges_color = prior_green;
-plotdat.anl_color = poste_blue;
-plotdat.obs_color = obs_red;
-
 %%---------------------------------------------------------------------
 % Loop around (time-copy-level-region) observation types
 %----------------------------------------------------------------------
@@ -203,19 +207,19 @@
 psfname = cell(plotdat.nvars);
 
 for ivar = 1:plotdat.nvars
-    
+
     % create the variable names of interest.
-    
+
     plotdat.myvarname = plotdat.varnames{ivar};
     plotdat.guessvar  = sprintf('%s_guess',plotdat.varnames{ivar});
     plotdat.analyvar  = sprintf('%s_analy',plotdat.varnames{ivar});
-    
+
     plotdat.trusted   = nc_read_att(fname, plotdat.guessvar, 'TRUSTED');
     if (isempty(plotdat.trusted)), plotdat.trusted = 'NO'; end
-    
+
     % remove any existing postscript file - will simply append each
     % level as another 'page' in the .ps file.
-    
+
     for iregion = 1:plotdat.nregions
         psfname{iregion} = sprintf('%s_rmse_%s_evolution_region%d.ps', ...
             plotdat.varnames{ivar}, plotdat.copystring, iregion);
@@ -222,19 +226,19 @@
         fprintf('Removing %s from the current directory.\n',psfname{iregion})
         system(sprintf('rm %s',psfname{iregion}));
     end
-