[Dart-dev] DART/trunk Revision: 11205

[dart at ucar.edu]

thoar at ucar.edu
2017-03-06 20:05:35 -0700 (Mon, 06 Mar 2017)
270
Making the trunk version of link_obs have the same x axis range for both axes on Figure 2.
Making the rma_trunk correctly put the continental overlay first so the draw order allows for
the scatterplot to be superimposed on top of the continents when viewed from above.




Modified: DART/branches/rma_trunk/diagnostics/matlab/link_obs.m
===================================================================
--- DART/branches/rma_trunk/diagnostics/matlab/link_obs.m	2017-03-07 02:53:05 UTC (rev 11204)
+++ DART/branches/rma_trunk/diagnostics/matlab/link_obs.m	2017-03-07 03:05:35 UTC (rev 11205)
@@ -99,6 +99,13 @@
 obs.region(5)     = min(obs.z); % use observation Z to specify vertical region
 obs.region(6)     = max(obs.z);
 
+% Some observations are all at same z value. Cannot use same value for
+% upper and lower axis bounds.
+if (obs.region(5) == obs.region(6))
+    obs.region(5) = obs.region(5) - 0.5;
+    obs.region(6) = obs.region(6) + 0.5;
+end
+
 %% Now pack the data in the same fashion as the cell array of column labels.
 
 obs.lonindex  = 1;

Modified: DART/branches/rma_trunk/diagnostics/matlab/linked_observations.m
===================================================================
--- DART/branches/rma_trunk/diagnostics/matlab/linked_observations.m	2017-03-07 02:53:05 UTC (rev 11204)
+++ DART/branches/rma_trunk/diagnostics/matlab/linked_observations.m	2017-03-07 03:05:35 UTC (rev 11205)
@@ -48,6 +48,8 @@
 view(fig1ax1,[-37.5 30]);
 grid(fig1ax1,'on');
 
+myworldmap(obs); hold on;
+
 xstring = sprintf('obsmat(:,%d)',obs.lonindex);
 ystring = sprintf('obsmat(:,%d)',obs.latindex);
 zstring = sprintf('obsmat(:,%d)',obs.zindex  );
@@ -64,8 +66,6 @@
 
 set(fig1ax1,'FontSize',18);
 
-myworldmap(obs);
-
 xlabel(obs.colnames{obs.lonindex});
 ylabel(obs.colnames{obs.latindex});
 zlabel(obs.colnames{obs.zindex});

Modified: DART/trunk/diagnostics/matlab/linked_observations.m
===================================================================
--- DART/trunk/diagnostics/matlab/linked_observations.m	2017-03-07 02:53:05 UTC (rev 11204)
+++ DART/trunk/diagnostics/matlab/linked_observations.m	2017-03-07 03:05:35 UTC (rev 11205)
@@ -152,31 +152,7 @@
 
 figure3 = figure(3); clf(figure3); orient tall; wysiwyg
 
-%% Create axes for ObsVal vs. DART QC scatterplot
-fig3ax1 = axes('Parent',figure3,'Position',[0.15 0.675 0.7 0.25]);
-box(fig3ax1,'on');
 
-xstring = sprintf('obsmat(:,%d)',obs.obsindex);
-ystring = sprintf('obsmat(:,%d)',obs.qcindex);
-
-h1 = scatter(obsmat(:,obs.obsindex),obsmat(:,obs.qcindex), ...
-    'Parent',fig3ax1, ...
-    'DisplayName','obs vs qc', ...
-    'XDataSource',xstring, ...
-    'YDataSource',ystring);
-
-set(fig3ax1,'FontSize',14);
-
-xlabel(obs.colnames{obs.obsindex});
-ylabel(obs.colnames{obs.qcindex});
-h = title(obs.ObsTypeString);
-set(h,'Interpreter','none');
-axis([-Inf Inf 0 8])
-grid(fig3ax1,'on');
-
-fprintf('QC summary follows:\n')
-LabelQC(obs.colnames{obs.qcindex}, obs.qc)
-
 %% Create axes for observation vs ensemble
 % This figure is most useful when all the 'bad' obs have been
 % replaced by Matlab's NAN so as not to blow the scale.
@@ -204,6 +180,8 @@
 axis(axlims)
 line([min(axis) max(axis)],[min(axis) max(axis)],'LineWidth',1.5,'Color','k')
 grid(fig3ax2,'on');
+xmin = axlims(1);
+xmax = axlims(2);
 
 if (sum(isfinite(get(h2,'YData')))) == 0
     Print_Empty_Banner(obs.colnames{obs.copyindex});
@@ -212,8 +190,32 @@
 refreshdata
 linkdata on
 
+%% Create axes for ObsVal vs. DART QC scatterplot
+fig3ax1 = axes('Parent',figure3,'Position',[0.15 0.675 0.7 0.25]);
+box(fig3ax1,'on');
 
+xstring = sprintf('obsmat(:,%d)',obs.obsindex);
+ystring = sprintf('obsmat(:,%d)',obs.qcindex);
 
+h1 = scatter(obsmat(:,obs.obsindex),obsmat(:,obs.qcindex), ...
+    'Parent',fig3ax1, ...
+    'DisplayName','obs vs qc', ...