[Dart-dev] [3966] DART/trunk/models: adding test script to exercise state-space diagnostics
nancy at ucar.edu
nancy at ucar.edu
Wed Jul 15 16:36:19 MDT 2009
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090715/2d1ef7da/attachment-0001.html
-------------- next part --------------
Added: DART/trunk/models/9var/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/9var/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/9var/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,214 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = './True_State.nc';
+ pinfo.diagn_file = './Prior_Diag.nc';
+ pinfo.model = '9var';
+ pinfo.var = 'state';
+ pinfo.truth_time = [1 1000];
+ pinfo.diagn_time = [1 1000];
+ pinfo.var_inds = [1 2 3 4 5 6 7 8 9];
+
+ close all; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ close all; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ close all; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ close all; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo;
+ pinfo.model = '9var';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 9;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 1000;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 9;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 2 3 4 5 6 7 8 9];
+ pinfo.fname = './Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.base_var_index = 4;
+ pinfo.base_time = 34;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = '9var';
+ pinfo.var1name = 'state';
+ pinfo.var2name = 'state';
+ pinfo.var3name = 'state';
+ pinfo.var1ind = 1;
+ pinfo.var2ind = 2;
+ pinfo.var3ind = 3;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
+%------------------------------------------------------------
+%plot_reg_factor
+%------------------------------------------------------------
+% plot_reg_factor
+
+%------------------------------------------------------------
+%plot_sawtooth
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_sawtooth
+ fprintf('Finished %s pausing, hit any key\n','plot_sawtooth'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = '9var';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 9;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 1000;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 9;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 2 3 4 5 6 7 8 9];
+ pinfo.prior_file = 'Prior_Diag.nc';
+ pinfo.posterior_file = 'Posterior_Diag.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.diagn_time = [-1 -1];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.truth_time = [-1 -1];
+ pinfo.var = 'state';
+ pinfo.var_inds = [1 2 3 4 5 6 7 8 9];
+ pinfo.copyindices = [7 12 17];
+ pinfo.prior_times = [1 1000];
+ pinfo.posterior_times = [1 1000];
+
+ PlotSawtooth(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotSawtooth'); pause
+
+%------------------------------------------------------------
+%plot_smoother_err
+%------------------------------------------------------------
+% plot_smoother_err
+
+%------------------------------------------------------------
+%plot_total_err
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_total_err
+ fprintf('Finished %s pausing, hit any key\n','plot_total_err'); pause
+end
+
+ clear pinfo; clf
+ pinfo.model = '9var';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 9;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 1000;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 9;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 2 3 4 5 6 7 8 9];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 1000];
+ pinfo.diagn_time = [1 1000];
+
+ PlotTotalErr(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotTotalErr'); pause
+
+%------------------------------------------------------------
+%plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_var_var_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_var_var_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.model = '9var';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 4;
+ pinfo.base_time = 235;
+ pinfo.state_var_index = 8;
+
+ PlotVarVarCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotVarVarCorrel'); pause
+
+%------------------------------------------------------------
+%plot_jeff_correl - virtually identical to plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_jeff_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_jeff_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 3;
+ pinfo.base_time = 300;
+ pinfo.state_var_index = 2;
+
+ clf; PlotJeffCorrel(pinfo)
+
Property changes on: DART/trunk/models/9var/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
+ text/x-matlab
Name: svn:keywords
+ Date Revision Author HeadURL Id
Name: svn:eol-style
+ native
Added: DART/trunk/models/bgrid_solo/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/bgrid_solo/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/bgrid_solo/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,239 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 2];
+ pinfo.diagn_time = [1 2];
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.var = 'u';
+ pinfo.level = 2;
+ pinfo.levelindex = 2;
+ pinfo.longitude = 258.00;
+ pinfo.lonindex = 43;
+ pinfo.latitude = 42.00;
+ pinfo.latindex = 22;
+
+ close all; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ clf; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ clf; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ clf; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.num_state_vars = 4;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 2;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.vars = {'ps' 't' 'u' 'v'};
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 't';
+ pinfo.comp_var = 'u';
+ pinfo.base_time = 0;
+ pinfo.base_tmeind = 1;
+ pinfo.base_lvl = 2;
+ pinfo.base_lvlind = 2;
+ pinfo.base_lat = 39.0000;
+ pinfo.base_latind = 22;
+ pinfo.base_lon = 255.0000;
+ pinfo.base_lonind = 43;
+ pinfo.comp_lvl = 2;
+ pinfo.comp_lvlind = 2;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.var1name = 'ps';
+ pinfo.var2name = 't';
+ pinfo.var3name = 'u';
+ pinfo.var1_lvl = 1;
+ pinfo.var2_lvl = 1;
+ pinfo.var3_lvl = 1;
+ pinfo.var1_lvlind = 1;
+ pinfo.var2_lvlind = 1;
+ pinfo.var3_lvlind = 1;
+ pinfo.var1_lat = 39.0000;
+ pinfo.var2_lat = 39.0000;
+ pinfo.var3_lat = 36.0000;
+ pinfo.var1_latind = 22;
+ pinfo.var2_latind = 22;
+ pinfo.var3_latind = 21;
+ pinfo.var1_lon = 255.0000;
+ pinfo.var2_lon = 255.0000;
+ pinfo.var3_lon = 252.0000;
+ pinfo.var1_lonind = 43;
+ pinfo.var2_lonind = 43;
+ pinfo.var3_lonind = 42;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
+%------------------------------------------------------------
+%plot_reg_factor
+%------------------------------------------------------------
+% plot_reg_factor
+
+%------------------------------------------------------------
+%plot_sawtooth
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_sawtooth
+ fprintf('Finished %s pausing, hit any key\n','plot_sawtooth'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.vars = {'ps' 't' 'u' 'v'};
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.prior_file = 'Prior_Diag.nc';
+ pinfo.posterior_file = 'Posterior_Diag.nc';
+ pinfo.var_names = 'u';
+ pinfo.level = 3;
+ pinfo.levelindex = 3;
+ pinfo.latitude = 42.0000;
+ pinfo.longitude = 258.0000;
+ pinfo.latindex = 22;
+ pinfo.lonindex = 43;
+ pinfo.copies = 3;
+ pinfo.copyindices = [7 12 17];
+ pinfo.truth_time = [1 2];
+ pinfo.prior_times = [1 2];
+ pinfo.posterior_times = [1 2];
+
+ PlotSawtooth(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotSawtooth'); pause
+
+%------------------------------------------------------------
+%plot_smoother_err
+%------------------------------------------------------------
+% plot_smoother_err
+
+%------------------------------------------------------------
+%plot_total_err
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_total_err
+ fprintf('Finished %s pausing, hit any key\n','plot_total_err'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.num_state_vars = 4;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 2;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.vars = {'ps' 't' 'u' 'v'};
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 2];
+ pinfo.diagn_time = [1 2];
+
+ PlotTotalErr(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotTotalErr'); pause
+
+%------------------------------------------------------------
+%plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_var_var_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_var_var_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.model = 'FMS_Bgrid';
+ pinfo.base_var = 't';
+ pinfo.comp_var = 'u';
+ pinfo.base_time = 0;
+ pinfo.base_tmeind = 1;
+
+ pinfo.base_lvl = 4;
+ pinfo.base_lat = -39.0000;
+ pinfo.base_lon = 183.0000;
+ pinfo.base_latind = 9;
+ pinfo.base_lvlind = 4;
+ pinfo.base_lonind = 31;
+
+ pinfo.comp_lvl = 2;
+ pinfo.comp_lat = 42.0000;
+ pinfo.comp_lon = 258.0000;
+ pinfo.comp_lvlind = 2;
+ pinfo.comp_latind = 22;
+ pinfo.comp_lonind = 43;
+
+ PlotVarVarCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotVarVarCorrel'); pause
+
+%------------------------------------------------------------
+%plot_jeff_correl - virtually identical to plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_jeff_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_jeff_correl'); pause
+end
+
+ PlotJeffCorrel(pinfo)
+
Property changes on: DART/trunk/models/bgrid_solo/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
+ text/x-matlab
Name: svn:keywords
+ Date Revision Author HeadURL Id
Name: svn:eol-style
+ native
Modified: DART/trunk/models/bgrid_solo/work/input.nml
===================================================================
--- DART/trunk/models/bgrid_solo/work/input.nml 2009-07-15 22:30:26 UTC (rev 3965)
+++ DART/trunk/models/bgrid_solo/work/input.nml 2009-07-15 22:36:18 UTC (rev 3966)
@@ -31,7 +31,7 @@
first_obs_seconds = -1,
last_obs_days = -1,
last_obs_seconds = -1,
- num_output_state_members = 0,
+ num_output_state_members = 20,
num_output_obs_members = 0,
output_interval = 1,
num_groups = 1,
Added: DART/trunk/models/forced_lorenz_96/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/forced_lorenz_96/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/forced_lorenz_96/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,207 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = './True_State.nc';
+ pinfo.diagn_file = './Prior_Diag.nc';
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.var = 'state';
+ pinfo.truth_time = [1 1000];
+ pinfo.diagn_time = [1 1000];
+ pinfo.var_inds = [1 13 27];
+
+ clf; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ clf; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ clf; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ clf; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 80;
+ pinfo.num_model_vars = 40;
+ pinfo.num_force_vars = 40;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 1000;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 80;
+ pinfo.min_model_var = 1;
+ pinfo.max_model_var = 40;
+ pinfo.min_force_var = 1;
+ pinfo.max_force_var = 40;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 13 27];
+ pinfo.def_force_vars = [41 53 67];
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.base_var_index = 30;
+ pinfo.base_time = 300;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.var1name = 'state';
+ pinfo.var2name = 'state';
+ pinfo.var3name = 'state';
+ pinfo.var1ind = 10;
+ pinfo.var2ind = 20;
+ pinfo.var3ind = 30;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
+%------------------------------------------------------------
+%plot_reg_factor
+%------------------------------------------------------------
+% plot_reg_factor
+
+%------------------------------------------------------------
+%plot_sawtooth
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_sawtooth
+ fprintf('Finished %s pausing, hit any key\n','plot_sawtooth'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.def_var = 'state';
+ pinfo.def_state_vars = [1 13 27];
+ pinfo.def_force_vars = [41 53 67];
+ pinfo.prior_file = 'Prior_Diag.nc';
+ pinfo.posterior_file = 'Posterior_Diag.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.diagn_time = [1 -1];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.truth_time = [1 -1];
+ pinfo.var = 'state';
+ pinfo.var_inds = [1 20 40];
+ pinfo.copyindices = [7 12 17];
+
+ PlotSawtooth(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotSawtooth'); pause
+
+%------------------------------------------------------------
+%plot_smoother_err
+%------------------------------------------------------------
+% plot_smoother_err
+
+%------------------------------------------------------------
+%plot_total_err
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_total_err
+ fprintf('Finished %s pausing, hit any key\n','plot_total_err'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.def_var = 'state';
+ pinfo.def_state_vars = [1 13 27];
+ pinfo.def_force_vars = [41 53 67];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 1000];
+ pinfo.diagn_time = [1 1000];
+
+ PlotTotalErr(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotTotalErr'); pause
+
+%------------------------------------------------------------
+%plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_var_var_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_var_var_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.model = 'Forced_Lorenz_96';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 30;
+ pinfo.base_time = 300;
+ pinfo.state_var_index = 60;
+
+ PlotVarVarCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotVarVarCorrel'); pause
+
+%------------------------------------------------------------
+%plot_jeff_correl - virtually identical to plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_jeff_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_jeff_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 20;
+ pinfo.base_time = 300;
+ pinfo.state_var_index = 10;
+
+ PlotJeffCorrel(pinfo)
+
Property changes on: DART/trunk/models/forced_lorenz_96/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
+ text/x-matlab
Name: svn:keywords
+ Date Revision Author HeadURL Id
Name: svn:eol-style
+ native
Added: DART/trunk/models/ikeda/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/ikeda/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/ikeda/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,187 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = './True_State.nc';
+ pinfo.diagn_file = './Prior_Diag.nc';
+ pinfo.model = 'Ikeda';
+ pinfo.var = 'state';
+ pinfo.truth_time = [1 100];
+ pinfo.diagn_time = [1 100];
+ pinfo.var_inds = [1 2];
+
+ clf; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ clf; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ clf; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ clf; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.model = 'Ikeda';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 2;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 100;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 2;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 2];
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.base_var_index = 1;
+ pinfo.base_time = 35;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = 'Ikeda';
+ pinfo.var1name = 'state';
+ pinfo.var2name = 'state';
+ pinfo.var1ind = 1;
+ pinfo.var2ind = 2;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
+%------------------------------------------------------------
+%plot_reg_factor
+%------------------------------------------------------------
+% plot_reg_factor
+
+%------------------------------------------------------------
+%plot_sawtooth
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_sawtooth
+ fprintf('Finished %s pausing, hit any key\n','plot_sawtooth'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'Ikeda';
+ pinfo.def_var = 'state';
+ pinfo.prior_file = 'Prior_Diag.nc';
+ pinfo.posterior_file = 'Posterior_Diag.nc';
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.truth_time = [1 -1];
+ pinfo.var = 'state';
+ pinfo.var_inds = [1 2];
+ pinfo.copyindices = [7 12 17];
+ pinfo.prior_times = [1 100];
+ pinfo.posterior_times = [1 100];
+
+ PlotSawtooth(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotSawtooth'); pause
+
+%------------------------------------------------------------
+%plot_smoother_err
+%------------------------------------------------------------
+% plot_smoother_err
+
+%------------------------------------------------------------
+%plot_total_err
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_total_err
+ fprintf('Finished %s pausing, hit any key\n','plot_total_err'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'Ikeda';
+ pinfo.def_var = 'state';
+ pinfo.def_state_vars = [1 2];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 100];
+ pinfo.diagn_time = [1 100];
+
+ PlotTotalErr(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotTotalErr'); pause
+
+%------------------------------------------------------------
+%plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_var_var_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_var_var_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.model = 'Ikeda';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 2;
+ pinfo.base_time = 40;
+ pinfo.state_var_index = 1;
+
+ PlotVarVarCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotVarVarCorrel'); pause
+
+%------------------------------------------------------------
+%plot_jeff_correl - virtually identical to plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_jeff_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_jeff_correl'); pause
+end
+
+ PlotJeffCorrel(pinfo)
+
Property changes on: DART/trunk/models/ikeda/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
+ text/x-matlab
Name: svn:keywords
+ Date Revision Author HeadURL Id
Name: svn:eol-style
+ native
Added: DART/trunk/models/lorenz_04/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/lorenz_04/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/lorenz_04/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,204 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = './True_State.nc';
+ pinfo.diagn_file = './Prior_Diag.nc';
+ pinfo.model = 'Lorenz_04';
+ pinfo.var = 'state';
+ pinfo.truth_time = [1 100];
+ pinfo.diagn_time = [1 100];
+ pinfo.var_inds = [100 200 300];
+
+ clf; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ clf; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ clf; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ clf; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo; clf
+pinfo.model = 'Lorenz_04';
+pinfo.def_var = 'state';
+pinfo.num_state_vars = 960;
+pinfo.num_ens_members = 24;
+pinfo.time_series_length = 100;
+pinfo.min_state_var = 1;
+pinfo.max_state_var = 960;
+pinfo.min_ens_mem = 1;
+pinfo.max_ens_mem = 24;
+pinfo.def_state_vars = [1 320 640];
+pinfo.fname = 'Prior_Diag.nc';
+pinfo.base_var = 'state';
+pinfo.base_var_index = 100;
+pinfo.base_time = 20;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = 'Lorenz_04';
+ pinfo.var1name = 'state';
+ pinfo.var2name = 'state';
+ pinfo.var3name = 'state';
+ pinfo.var1ind = 100;
+ pinfo.var2ind = 200;
+ pinfo.var3ind = 300;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
+%------------------------------------------------------------
+%plot_reg_factor
+%------------------------------------------------------------
+% plot_reg_factor
+
+%------------------------------------------------------------
+%plot_sawtooth
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_sawtooth
+ fprintf('Finished %s pausing, hit any key\n','plot_sawtooth'); pause
+end
+
+ clear pinfo; close all
+ pinfo.model = 'Lorenz_04';
+ pinfo.def_var = 'state';
+ pinfo.def_state_vars = [1 320 640];
+ pinfo.prior_file = 'Prior_Diag.nc';
+ pinfo.posterior_file = 'Posterior_Diag.nc';
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.truth_time = [1 -1];
+ pinfo.var = 'state';
+ pinfo.var_inds = [1 320 640];
+ pinfo.copyindices = [7 12 17];
+ pinfo.prior_times = [1 100];
+ pinfo.posterior_times = [1 100];
+
+ PlotSawtooth(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotSawtooth'); pause
+
+%------------------------------------------------------------
+%plot_smoother_err
+%------------------------------------------------------------
+% plot_smoother_err
+
+%------------------------------------------------------------
+%plot_total_err
+%------------------------------------------------------------
+if (interactive)
+ clear; close all; plot_total_err
+ fprintf('Finished %s pausing, hit any key\n','plot_total_err'); pause
+end
+
+ clear pinfo; close all
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 960;
+ pinfo.num_ens_members = 24;
+ pinfo.time_series_length = 100;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 960;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 24;
+ pinfo.def_state_vars = [1 320 640];
+ pinfo.truth_file = 'True_State.nc';
+ pinfo.diagn_file = 'Prior_Diag.nc';
+ pinfo.truth_time = [1 100];
+ pinfo.diagn_time = [1 100];
+
+ PlotTotalErr(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotTotalErr'); pause
+
+%------------------------------------------------------------
+%plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_var_var_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_var_var_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.model = 'Lorenz_04';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 200;
+ pinfo.base_time = 30;
+ pinfo.state_var_index = 800;
+
+ PlotVarVarCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotVarVarCorrel'); pause
+
+%------------------------------------------------------------
+%plot_jeff_correl - virtually identical to plot_var_var_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_jeff_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_jeff_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.state_var = 'state';
+ pinfo.base_var_index = 800;
+ pinfo.base_time = 30;
+ pinfo.state_var_index = 200;
+
+ PlotJeffCorrel(pinfo)
+
Property changes on: DART/trunk/models/lorenz_04/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
+ text/x-matlab
Name: svn:keywords
+ Date Revision Author HeadURL Id
Name: svn:eol-style
+ native
Added: DART/trunk/models/lorenz_63/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/lorenz_63/matlab/RunAllTests.m (rev 0)
+++ DART/trunk/models/lorenz_63/matlab/RunAllTests.m 2009-07-15 22:36:18 UTC (rev 3966)
@@ -0,0 +1,196 @@
+function RunAllTests(dummy)
+% RunAllTests.m
+
+%------------------------------------------------------------
+% Data Assimilation Research Testbed -- DART
+% Copyright 2004-2009, Data Assimilation Research Section
+% University Corporation for Atmospheric Research
+% Licensed under the GPL -- www.gpl.org/licenses/gpl.html
+%
+% <next few lines under version control, do not edit>
+% $URL:
+% http://subversion.ucar.edu/DAReS/DART/trunk/matlab/get_state_copy.m $
+% $Id$
+% $Revision$
+% $Date$
+%------------------------------------------------------------
+
+if (nargin() > 0)
+ interactive = 1;
+else
+ interactive = 0;
+end
+
+if (interactive)
+ plot_bins
+ fprintf('Finished %s pausing, hit any key\n','plot_bins'); pause
+ plot_ens_err_spread
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_err_spread'); pause
+ plot_ens_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_time_series'); pause
+ plot_ens_mean_time_series
+ fprintf('Finished %s pausing, hit any key\n','plot_ens_mean_time_series'); pause
+end
+
+ clear pinfo; close all;
+ pinfo.truth_file = './True_State.nc';
+ pinfo.diagn_file = './Prior_Diag.nc';
+ pinfo.model = 'Lorenz_63';
+ pinfo.var = 'state';
+ pinfo.truth_time = [1 100];
+ pinfo.diagn_time = [1 100];
+ pinfo.var_inds = [1 2 3];
+
+ clf; PlotBins(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotBins'); pause
+
+ clf; PlotEnsErrSpread(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsErrSpread'); pause
+
+ clf; PlotEnsTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsTimeSeries'); pause
+
+ clf; PlotEnsMeanTimeSeries(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotEnsMeanTimeSeries'); pause
+
+%------------------------------------------------------------
+%plot_correl
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_correl
+ fprintf('Finished %s pausing, hit any key\n','plot_correl'); pause
+end
+
+ clear pinfo; clf
+ pinfo.model = 'Lorenz_63';
+ pinfo.def_var = 'state';
+ pinfo.num_state_vars = 3;
+ pinfo.num_ens_members = 22;
+ pinfo.time_series_length = 100;
+ pinfo.min_state_var = 1;
+ pinfo.max_state_var = 3;
+ pinfo.min_ens_mem = 1;
+ pinfo.max_ens_mem = 22;
+ pinfo.def_state_vars = [1 2 3];
+ pinfo.fname = './Prior_Diag.nc';
+ pinfo.base_var = 'state';
+ pinfo.base_var_index = 1;
+ pinfo.base_time = 20;
+
+ PlotCorrel(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotCorrel'); pause
+
+%------------------------------------------------------------
+%plot_phase_space
+%------------------------------------------------------------
+if (interactive)
+ clear; clf; plot_phase_space
+ fprintf('Finished %s pausing, hit any key\n','plot_phase_space'); pause
+end
+
+ clear pinfo; clf
+ pinfo.fname = 'True_State.nc';
+ pinfo.model = 'Lorenz_63';
+ pinfo.var1name = 'state';
+ pinfo.var2name = 'state';
+ pinfo.var3name = 'state';
+ pinfo.var1ind = 1;
+ pinfo.var2ind = 2;
+ pinfo.var3ind = 3;
+ pinfo.ens_mem = 'true state';
+ pinfo.ltype = 'k-';
+
+ PlotPhaseSpace(pinfo)
+ fprintf('Finished %s pausing, hit any key\n','PlotPhaseSpace'); pause
+
@@ Diff output truncated at 40000 characters. @@
More information about the Dart-dev
mailing list