changed gitignore

Author: Abhat-93

.gitignore

@@ -3,3 +3,4 @@
 *.mat
 *.png
 *.fig
+*.pdf

src/+datools/+examples/+sandu/l63_experiments.m

@@ -4,11 +4,11 @@
 % Use this to run Lorenz 63 experiments
 modelname = 'Lorenz63';
 % uncomment the filter you want to run
-%filtername = 'EnKF';
+filtername = 'EnKF';
 %filtername = 'ETKF';
 %filtername = 'ETPF';
 %filtername = 'SIR';
-filtername = 'RHF';
+%filtername = 'RHF';
 
 % oservation variance
 variance = 1;

src/+datools/+examples/+sandu/runexperiments.m

@@ -40,6 +40,8 @@ function runexperiments(user)
         filtertype = 'Particle';
     case 'RHF'
         filtertype = 'Ensemble';
+    case 'SIS_EnKF'
+        filtertype = 'Particle';
 end
 fprintf('Filtername = %s, Observation Variance = %.2f, Runs = %d, spinups = %d\n', ...
     filtername, variance, steps, spinup);
@@ -213,6 +215,16 @@ function runexperiments(user)
                     'Parallel', false, ...
                     'RankHistogram', histvar, ...
                     'Rejuvenation', rejuvenation);
+            case 'SIS_EnKF'
+                filter = datools.statistical.ensemble.(filtername)(model, ...
+                    'Observation', observation, ...
+                    'NumEnsemble', ensN, ...
+                    'ModelError', modelerror, ...
+                    'EnsembleGenerator', ensembleGenerator, ...
+                    'Inflation', inflation, ...
+                    'Parallel', false, ...
+                    'RankHistogram', histvar,...
+                    'Rejuvenation', rejuvenation);
             case 'RHF'
                 filter = datools.statistical.ensemble.(filtername)(model, ...
                     'Observation', observation, ...

src/+datools/+examples/+sandu/runnerscript.m

@@ -17,17 +17,20 @@
 % filtername = 'RHF';
 
 %% User inputs
-user.modelname = 'Lorenz96';
+user.modelname = 'Lorenz63';
 
-user.filtername = 'EnKF';
+user.filtername = 'SIS_EnKF';
 
-user.ensNs = [25, 50, 75, 100];
+%user.ensNs = [25, 50, 75, 100];
+user.ensNs = [1000];
 
-user.infs = [1.01, 1.02, 1.05];
+%user.infs = [1.01, 1.02, 1.05];
+user.infs = [1.01];
 
 user.variance = 1; % observation variance
 
-user.rejs = round(2*logspace(-2, -1, 4), 2);
+%user.rejs = round(2*logspace(-2, -1, 4), 2);
+user.rejs = [0.1];
 
 user.spinups = 2;