diff --git a/run_wrf_dh.py b/run_wrf_dh.py
index dff06cf88a0250a9f19faf72b985d0ac887ebe11..b9d4cd02a9e103201d98ff629fa5c2de12410e9d 100644
--- a/run_wrf_dh.py
+++ b/run_wrf_dh.py
@@ -3,7 +3,7 @@ __date__ = '2020-06-29'
 
 import argparse
 # from mlair.data_handler.data_handler_single_station import DataHandlerSingleStation
-from mlair.data_handler.data_handler_wrf_chem import DataHandlerWRF, DataHandlerMainSectWRF
+from mlair.data_handler.data_handler_wrf_chem import DataHandlerWRF, DataHandlerMainSectWRF, DataHandlerMainMinorSectWRF
 from mlair.workflows import DefaultWorkflow
 from mlair.helpers import remove_items
 from mlair.configuration.defaults import DEFAULT_PLOT_LIST
@@ -25,11 +25,14 @@ def load_stations():
 
 def main(parser_args):
     plots = remove_items(DEFAULT_PLOT_LIST, "PlotConditionalQuantiles")
+    plots.append("PlotBootstrapSkillScore")
     workflow = DefaultWorkflow(  # stations=load_stations(),
         # stations=["DEBW087","DEBW013", "DEBW107",  "DEBW076"],
 
         train_model=False, create_new_model=True, network="UBA",
-        evaluate_bootstraps=False,  # plot_list=["PlotCompetitiveSkillScore"],
+        evaluate_bootstraps=True,  number_of_bootstraps=30, create_new_bootstraps=True,  #
+        plot_list=plots,
+        model_name_for_plots='SecModel',
 #         competitors=["test_model", "test_model2"],
 #         competitor_path=os.path.join(os.getcwd(), "data", "comp_test"),
         competitors=["baseline", "sector_baseline"],
@@ -43,9 +46,10 @@ def main(parser_args):
         window_history_size=6,
         # ('Germany', (5.98865807458, 47.3024876979, 15.0169958839, 54.983104153))
         stations=["coords__48_8479__10_0963", "coords__51_8376__14_1417",
-                  "coords__50_7536__7_0827", "coords__51_4070__6_9656",
-                  "coords__49_8421__7_8662", "coords__49_7410__7_1935",
-                  "coords__51_1566__11_8182", "coords__51_4065__6_9660",
+                  # "coords__50_7536__7_0827", "coords__51_4070__6_9656",
+                  # "coords__49_8421__7_8662", "coords__49_7410__7_1935",
+                  # "coords__51_1566__11_8182", "coords__51_4065__6_9660",
+
                   # "coords__50_7536__7_0827", "coords__54_0000__10_000",
 
                   # "coords__50_7333__7_1000", "coords__50_0000__8_0000",
@@ -55,8 +59,9 @@ def main(parser_args):
                   # "coords__45_7777__9_1000", "coords__50_0000__4_0000",
                   ],
 
-        data_handler=DataHandlerWRF,
-        # data_handler=DataHandlerMainSectWRF, #,
+        # data_handler=DataHandlerWRF,
+        data_handler=DataHandlerMainSectWRF, #,
+        # data_handler=DataHandlerMainMinorSectWRF,
         # data_path='/home/felix/Data/WRF-Chem/upload_aura_2021-02-24/2009/',
         # data_path='/home/felix/Data/WRF-Chem/test_cut_nc/',
         # data_path='/home/felix/Data/WRF-Chem/test_cut_nc_joint',
@@ -68,9 +73,9 @@ def main(parser_args):
         external_coords_file='/home/felix/Data/WRF-Chem/test_cut_nc/coords.nc',
         # external_coords_file = "/p/scratch/deepacf/kleinert1/IASS_proc/coords.nc",
 
-        date_format_of_nc_file = "%Y-%m",
-        common_file_starter = "wrfout_d01",
-        time_dim = 'XTIME',
+        date_format_of_nc_file="%Y-%m",
+        common_file_starter="wrfout_d01",
+        time_dim='XTIME',
 
         transformation={
             "T2": {"method": "standardise"},
@@ -90,6 +95,8 @@ def main(parser_args):
         target_var='o3',
         statistics_per_var={'T2': None, 'o3': None, 'wdir10ll': None, 'wspd10ll': None,
                             'no': None, 'no2': None, 'co': None, 'PSFC': None, 'PBLH': None, 'CLDFRA': None, },
+        # separate_vars=["o3", "o3Sect", "o3SectLeft", "o3SectRight"],
+        separate_vars=["o3", "o3Sect"],
         # variables=['T2', 'Q2', 'PBLH', 'U10ll', 'V10ll', 'wdir10ll', 'wspd10ll'],
         # target_var=["T2"],
         # statistics_per_var={'T2': None, 'Q2': None, 'PBLH': None,
@@ -101,33 +108,33 @@ def main(parser_args):
         radius=100,  # km
 
         start='2009-01-01',
-        #end='2009-01-04',
+        end='2009-01-04',
         #end='2009-01-31',
-        end='2009-03-31',
+        # end='2009-03-31',
         
-        #train_start='2009-01-01',
-        #train_end='2009-01-02',
         train_start='2009-01-01',
+        train_end='2009-01-02',
+        # train_start='2009-01-01',
         #train_end='2009-01-15',
-        train_end='2009-02-28',
+        # train_end='2009-02-28',
         
-        #val_start='2009-01-02',
-        #val_end='2009-01-03',
+        val_start='2009-01-02',
+        val_end='2009-01-03',
         ###################################
         #val_start='2009-01-15',
         #val_end='2009-01-22',
         ###################################
-        val_start='2009-03-01',
-        val_end='2009-03-14',
+        # val_start='2009-03-01',
+        # val_end='2009-03-14',
 
-        #test_start='2009-01-03',
-        #test_end='2009-01-04',
+        test_start='2009-01-03',
+        test_end='2009-01-04',
         ###################################
         #test_start='2009-01-22',
         #test_end='2009-01-31',
         ###################################
-        test_start='2009-03-15',
-        test_end='2009-03-31',
+        # test_start='2009-03-15',
+        # test_end='2009-03-31',
         
         sampling='hourly',