[Dart-dev] [3965] DART/trunk/models/lorenz_96_2scale: adding test script to exercise state-space diagnostics

[nancy at ucar.edu]

An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090715/4712df39/attachment.html 
-------------- next part --------------
Added: DART/trunk/models/lorenz_96_2scale/matlab/RunAllTests.m
===================================================================
--- DART/trunk/models/lorenz_96_2scale/matlab/RunAllTests.m	                        (rev 0)
+++ DART/trunk/models/lorenz_96_2scale/matlab/RunAllTests.m	2009-07-15 22:30:26 UTC (rev 3965)
@@ -0,0 +1,193 @@
+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_96_2scale';
+ pinfo.truth_time     = [1 365];
+ pinfo.diagn_time     = [1 365];
+ pinfo.var            = 'Y';
+ 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_96_2scale';
+ pinfo.def_var            = 'state';
+ pinfo.num_state_vars     = 396;
+ pinfo.num_ens_members    = 24;
+ pinfo.time_series_length = 365;
+ pinfo.min_state_var      = 1;
+ pinfo.max_state_var      = 396;
+ pinfo.min_ens_mem        = 1;
+ pinfo.max_ens_mem        = 24;
+ pinfo.def_state_vars     = [1 2 3];
+ pinfo.fname              = './Prior_Diag.nc';
+ pinfo.base_var           = 'X';
+ pinfo.base_var_index     = 10;
+ pinfo.base_time          = 125;
+ 
+ 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    = 'Prior_Diag.nc';
+ pinfo.model    = 'Lorenz_96_2scale';
+ pinfo.var1name = 'X';
+ pinfo.var2name = 'Y';
+ pinfo.var3name = 'X';
+ pinfo.var1ind  = 30;
+ pinfo.var2ind  = 300;
+ pinfo.var3ind  = 10;
+ pinfo.ens_mem  = 'ensemble member10';
+ 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_96_2scale';
+ 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            = 'X';
+ pinfo.var_inds       = [1 12 24];
+ 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              = 'Lorenz_96_2scale';
+ pinfo.def_state_vars     = [1 12 24];
+ pinfo.truth_file         = 'True_State.nc';
+ pinfo.diagn_file         = 'Prior_Diag.nc';
+ pinfo.truth_time         = [1 365];
+ pinfo.diagn_time         = [1 365];
+
+ 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_96_2scale';
+ pinfo.base_var        = 'X';
+ pinfo.state_var       = 'Y';
+ pinfo.base_var_index  = 18;
+ pinfo.base_time       = 128;
+ pinfo.state_var_index = 180;
+
+ 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        = 'X';
+ pinfo.state_var       = 'Y';
+ pinfo.base_var_index  = 20;
+ pinfo.base_time       = 30;
+ pinfo.state_var_index = 300;
+
+ PlotJeffCorrel(pinfo)
+


Property changes on: DART/trunk/models/lorenz_96_2scale/matlab/RunAllTests.m
___________________________________________________________________
Name: svn:mime-type
   + text/x-matlab
Name: svn:keywords
   + Date Revision Author HeadURL Id
Name: svn:eol-style
   + native