[Dart-dev] DART/trunk Revision: 11647

[dart at ucar.edu]

gharamti at ucar.edu
2017-05-19 13:59:46 -0600 (Fri, 19 May 2017)
70
Sped-up the plotting. Added the ability to zoom into the polar plot. 



Modified: DART/trunk/DART_LAB/matlab/oned_model_inf.m
===================================================================
--- DART/trunk/DART_LAB/matlab/oned_model_inf.m	2017-05-19 16:04:17 UTC (rev 11646)
+++ DART/trunk/DART_LAB/matlab/oned_model_inf.m	2017-05-19 19:59:46 UTC (rev 11647)
@@ -1110,15 +1110,6 @@
             handles.error  = prior_error;
             handles.spread = prior_spread;
 
-            legend([h_e h_s], 'Error', 'Spread', 'Location', 'NorthEast');
-            set(legend,'FontName', atts.fontname, 'FontSize', atts.fontsize);
-            legend boxon
-
-            axlims    = axis;
-            axlims(3) = 0.0;
-            axis(axlims)
-
-
             % Update the prior rank histogram figure
             axes(handles.h_prior_rank_histogram);
 
@@ -1133,23 +1124,20 @@
 
             % Plot the figure window for this update
             axes(handles.axes);
-            hold off;
+            cla;
 
             % Find the limits of the plot
             % The height of the obs likelihood controls the vertical axis
             y_max = 1 / (sqrt(2 * pi) * handles.obs_error_sd);
-            xmin = -10;
-            xmax =  10;
+            xmin  = -10;
+            xmax  =  10;
 
             % Put on a black axis line using data limits
-
             plot([xmin xmax], [0, 0], 'k', 'Linewidth', 2);
             hold on;
             ens_axis = [xmin xmax -0.2 y_max + 0.02];
-
-            % Turn off the negative labels
-            set(gca, 'YTick', [0 0.1 0.2 0.3 0.4], 'XLim',[ens_axis(1) ens_axis(2)], 'YLim',[ens_axis(3) ens_axis(4)]);
             grid on;
+            
 
             % Plot the prior ensemble members in green
             % Plotting ticks instead of asterisks makes bins clearer
@@ -1180,13 +1168,6 @@
             h = line([base_x + text_width / 8, 0], [-0.08, -0.03]);
             set(h, 'Color', 'k');
 
-            % Label this plot
-            xlabel('State'               ,'FontName', atts.fontname, 'FontSize', atts.fontsize);
-            title('Latest Ensemble Prior','FontName', atts.fontname, 'FontSize', atts.fontsize);
-
-            L = legend([hg_prior hg_truth],'Prior','Truth');
-            set(L, 'FontName', atts.fontname, 'FontSize', atts.fontsize);
-
             % Update the permanent storage of the rank values
             handles.prior_rank(ens_rank) = handles.prior_rank(ens_rank) + 1;
 
@@ -1214,7 +1195,6 @@
 
             % Adjust the horizontal range of the plot windows as needed
             if( mod(handles.time_step, 5) == 0)
-
                 axes(handles.h_state_evolution);
                 axlims    = axis;
                 axlims(1) = handles.time_step - 4;
@@ -1222,7 +1202,6 @@
                 axis(axlims)
 
                 % Want the lower y limit to stay 0 for error spread
-
                 axes(handles.h_err_spread_evolution);
                 axlims    = axis;
                 axlims(1) = handles.time_step - 4;
@@ -1231,7 +1210,6 @@
                 axis(axlims)
 
                 % Want the lower y limit to stay 0 for inflation
-
                 axes(handles.h_inflation_evolution);
                 axlims    = axis;
                 axlims(1) = handles.time_step - 4;
@@ -1300,7 +1278,6 @@
             handles.post_rank(ens_rank) = handles.post_rank(ens_rank) + 1;
 
             %% Plot the segment for the updated error
-
             axes(handles.h_err_spread_evolution);
 
             post_error = calculate_rmse(new_ens, 0.0);
@@ -1320,15 +1297,9 @@
                 [handles.spread, post_spread]);
             set(h, 'Color', atts.red, 'LineWidth', 2.0);
 
-            % Want the lower y limit to stay 0 for error spread
-            axlims    = axis;
-            axlims(3) = 0.0;