[Dart-dev] DART/branches Revision: 13016

[dart at ucar.edu]

thoar at ucar.edu
2019-03-14 16:48:48 -0600 (Thu, 14 Mar 2019)
76
progress toward supporting QC 8 if available and no posterior information.




Modified: DART/branches/qc8diags/diagnostics/matlab/plot_evolution.m
===================================================================
--- DART/branches/qc8diags/diagnostics/matlab/plot_evolution.m	2019-03-14 22:16:20 UTC (rev 13015)
+++ DART/branches/qc8diags/diagnostics/matlab/plot_evolution.m	2019-03-14 22:48:48 UTC (rev 13016)
@@ -235,9 +235,7 @@
     logfid = fopen(lgfname,'wt');
     fprintf(logfid,'%s\n',lgfname);
 
-    %% todo FIXME replace with a permute routine to get desired shape
     % get appropriate vertical coordinate variable
-    % regions-levels-copy-time
 
     [dimnames, ~] = nc_var_dims(fname, plotdat.guessvar);
 
@@ -269,8 +267,10 @@
 
     analy_raw = local_ncread(fname, plotdat.analyvar);
     if ( isempty(analy_raw) )
-       analy = guess;
-       analy(:) = NaN;
+        % force analysis to be the same shape as the guess,
+        % and full of NaNs instead of 0s.
+        analy = guess;
+        analy(:) = NaN;
     else
        analy_raw = permute(analy_raw,length(size(analy_raw)):-1:1);
        analy = reshape(analy_raw, plotdat.Nbins,   plotdat.ncopies, ...
@@ -485,7 +485,7 @@
 switch lower(plotdat.copystring)
     case 'bias'
         % plot a zero-bias line
-        zeroline = line(axlims(1:2),[0 0], 'Color',[0 100 0]/255,'Parent',ax1);
+        zeroline = line(axlims(1:2),[0 0], 'Color',[200 200 200]/255,'Parent',ax1);
         set(zeroline,'LineWidth',2.5,'LineStyle','-')
         plotdat.ylabel = sprintf('%s (%s)',plotdat.copystring,plotdat.biasconv);
     otherwise
@@ -550,14 +550,11 @@
 
 nevaluated = sum(plotdat.ges_Nqc1(:) + plotdat.ges_Nqc3(:));
 if (nevaluated > 0)
-   set(get(ax2,'Ylabel'), ...
-       'String','# of obs : o=possible, \ast,\diamondsuit=evaluated', ...
-       'FontSize', figuredata.fontsize)
+   string1 = '# of obs : o=possible; \ast,\diamondsuit=evaluated';
 else
-   set(get(ax2,'Ylabel'), ...
-       'String','# of obs : o=possible, \ast,\diamondsuit=assimilated', ...
-       'FontSize', figuredata.fontsize)
+   string1 = '# of obs : o=possible; \ast,\diamondsuit=assimilated';
 end
+set(get(ax2,'Ylabel'), 'String', string1, 'FontSize', figuredata.fontsize)
 
 
 %=====================================================================

Modified: DART/branches/qc8diags/diagnostics/matlab/plot_rank_histogram.m
===================================================================
--- DART/branches/qc8diags/diagnostics/matlab/plot_rank_histogram.m	2019-03-14 22:16:20 UTC (rev 13015)
+++ DART/branches/qc8diags/diagnostics/matlab/plot_rank_histogram.m	2019-03-14 22:48:48 UTC (rev 13016)
@@ -60,6 +60,10 @@
     error('wrong number of arguments ... ')
 end
 
+%TODO actually implement the varargin ... should be able to specify a 
+%         specific region or level
+
+
 % Make sure the file exists.
 
 if (exist(fname,'file') ~= 2)
@@ -276,7 +280,7 @@
         % plot by region
         
         for iregion = 1:plotdat.nregions
-            figure(iregion); clf; orient(figuredata.orientation); wysiwyg
+            figure(iregion); clf; orient(figuredata.orientation);
             
             plotdat.region   = iregion;
             plotdat.myregion = deblank(plotdat.region_names(iregion,:));

Modified: DART/branches/qc8diags/diagnostics/matlab/plot_rmse_xxx_evolution.m
===================================================================
--- DART/branches/qc8diags/diagnostics/matlab/plot_rmse_xxx_evolution.m	2019-03-14 22:16:20 UTC (rev 13015)
+++ DART/branches/qc8diags/diagnostics/matlab/plot_rmse_xxx_evolution.m	2019-03-14 22:48:48 UTC (rev 13016)
@@ -71,6 +71,7 @@
 %
 % DART $Id$
 
+default_markersize = 8;
 default_level = -1;
 default_obsname = 'none';
 default_verbosity = 'yes';
@@ -84,11 +85,13 @@
     addParameter(p,'range',default_range, at isnumeric);
     addParameter(p,'level',default_level, at isnumeric);
     addParameter(p,'verbose',default_verbosity, at ischar);
+    addParameter(p,'MarkerSize',default_markersize, at isnumeric);
 else
     addParamValue(p,'obsname',default_obsname, at ischar);
     addParamValue(p,'range',default_range, at isnumeric);
     addParamValue(p,'level',default_level, at isnumeric);