[Dart-dev] DART/releases Revision: 11632
dart at ucar.edu
dart at ucar.edu
Fri May 12 11:55:39 MDT 2017
thoar at ucar.edu
2017-05-12 11:55:38 -0600 (Fri, 12 May 2017)
400
Adding Moha's inflation example scripts that were committed to the
rma_trunk at revision 11630. There are some feature enhancements
to clear the histograms, etc.
In addition, these scripts have the localization graphic and
fixed a bug: the nonlinear 'a' parameter must be positive,
and fixed another bug that did not preserve the true forcing at 8 when
the forcing input for the model was changed.
Modified: DART/releases/Lanai/DART_LAB/matlab/gaussian_product.m
===================================================================
--- DART/releases/Lanai/DART_LAB/matlab/gaussian_product.m 2017-05-11 23:58:18 UTC (rev 11631)
+++ DART/releases/Lanai/DART_LAB/matlab/gaussian_product.m 2017-05-12 17:55:38 UTC (rev 11632)
@@ -10,11 +10,11 @@
% If the parameters of the two gaussians are known, the parameters of the
% resulting gaussian can be calculated.
%
-% See also: oned_model, oned_ensemble, twod_ensemble,
-% run_lorenz_63, run_lorenz_96
+% See also: oned_model.m oned_model_inf.m oned_ensemble.m
+% twod_ensemble.m run_lorenz_63.m run_lorenz_96.m run_lorenz_96_inf.m
-%% DART software - Copyright 2004 - 2016 UCAR. This open source software is
-% provided by UCAR, "as is", without charge, subject to all terms of use at
+%% DART software - Copyright UCAR. This open source software is provided
+% by UCAR, "as is", without charge, subject to all terms of use at
% http://www.image.ucar.edu/DAReS/DART/DART_download
%
% DART $Id$
@@ -190,7 +190,7 @@
'FontUnits', 'normalized', ...
'FontName', atts.fontname, ...
'FontSize', 0.2);
-
+
text_xstart = 0.01;
text_dy = 0.25;
text_y = 0.925 - text_dy - 0.05;
@@ -255,9 +255,9 @@
%This function plots the graph using the inputs in the 4 edit boxes. It
%makes changes to handles, so the function must return an update to
%handles. This is done through the function definition.
-
+
[prior_mean, prior_sd, obs_mean, obs_err_sd, is_err] = g_prod_plot(handles);
-
+
% If there is an error, zero out the posterior text values
% don't try to do posterior computation
if(is_err)
@@ -264,13 +264,13 @@
reset_Posterior();
return;
end
-
+
% Compute the posterior mean, sd and weight
[post_mean, post_sd, weight] = ...
product_of_gaussians(prior_mean, prior_sd, obs_mean, obs_err_sd);
post_handle = plot_gaussian(post_mean, post_sd, 1);
- set(post_handle, 'Color', 'b', 'LineWidth', 2);
-
+ set(post_handle, 'Color', atts.blue, 'LineWidth', 2);
+
%Round post_mean, post_sd and weight to 4 decimal places
post_mean = round(post_mean * 10000);
post_mean = post_mean/10000;
@@ -278,22 +278,22 @@
post_sd = post_sd/10000;
weight = round(weight * 10000);
weight = weight/10000;
-
+
% Print values
str1 = sprintf('Mean = %.4f',post_mean);
set(handles.ui_text_posterior_mean, 'String', str1);
str1 = sprintf('SD = %.4f',post_sd);
set(handles.ui_text_posterior_sd, 'String', str1);
-
+
% Also plot the weighted posterior as dashed
post_handle = plot_gaussian(post_mean, post_sd, weight);
- set(post_handle, 'Color', 'b', 'LineStyle', '--');
+ set(post_handle, 'Color', atts.blue, 'LineStyle', '--');
str1 = sprintf('Weight = %.4f',weight);
set(handles.ui_text_posterior_weight, 'String', str1);
-
+
h = legend('Prior', 'Obs. Likelihood', 'Posterior', 'Weighted Posterior');
set(h, 'box', 'on', 'Location', 'NorthWest')
-
+
end
% These functions plot the graph immediately after the user edits a text box
@@ -330,4 +330,3 @@
% $URL$
% $Revision$
% $Date$
-
Modified: DART/releases/Lanai/DART_LAB/matlab/oned_ensemble.m
===================================================================
--- DART/releases/Lanai/DART_LAB/matlab/oned_ensemble.m 2017-05-11 23:58:18 UTC (rev 11631)
+++ DART/releases/Lanai/DART_LAB/matlab/oned_ensemble.m 2017-05-12 17:55:38 UTC (rev 11632)
@@ -29,11 +29,11 @@
% change the Observation Error SD, lay down an ensemble pretty far away
% from the observation - have fun with it.
%
-% See also: gaussian_product, oned_model, twod_ensemble, run_lorenz_63,
More information about the Dart-dev
mailing list