add reference set to CSI branch

README.md

@@ -10,7 +10,7 @@ To compile and run:
 * Run `make` in the test folder to compile
 * Run `./LakeComoSim settings_lakeComo.txt < u_test46.txt` to perform a simulation with a random policy
 
-The simulation is performed using a 1 year trajectory of inflow data simulated by a hydrologic model; the mass balance equation of the lake storage is integrated over time with a hourly time step; the daily release decision is determined by a Gaussian Radial Basis Function with 4 inputs, i.e. sine and cosine of the day of the year, lake level, and a perfect forecast of the cumulated inflows over the next 51 days according to the results of [Denaro et al. 2017](https://www.sciencedirect.com/science/article/pii/S0309170816304651), 1 output, and 5 Gaussian bases, for a total of 46 policy parameters (decision variables). See the settings file in the test folder for further details.
+The simulation is performed using a 1 year trajectory of inflow data simulated by a hydrologic model; the mass balance equation of the lake storage is integrated over time with a hourly time step; the daily release decision is determined by a Gaussian Radial Basis Function with 4 inputs (i.e. sine and cosine of the day of the year, lake level, and a perfect forecast of the cumulated inflows over the next 51 days according to the results of [Denaro et al. 2017](https://www.sciencedirect.com/science/article/pii/S0309170816304651) ), 1 output, and 5 Gaussian bases, for a total of 46 policy parameters (decision variables). See the settings file in the test folder for further details.
 The observed inflows used in the paper are from [Consorzio dell'Adda](http://www.addaconsorzio.it/). 
 
 

refset/optLakeComo_perfectF.reference

@@ -0,0 +1,24 @@
+% Problem = ProblemStub
+% Variables = 29
+% Objectives = 2
+6.923076923076923 2349.077654373516
+7.461538461538462 2336.700087347796
+6.384615384615385 2364.8797093352996
+5.461538461538462 2467.128022279774
+5.923076923076923 2400.7738179730504
+13.692307692307692 2263.7802025393407
+22.692307692307693 2243.622615343896
+8.461538461538462 2315.265073575777
+14.307692307692308 2259.320651687491
+18.153846153846153 2249.8770128635206
+15.923076923076923 2253.6732097665936
+10.846153846153847 2281.346211939603
+9.461538461538462 2297.6071698315955
+12.461538461538462 2271.3175754030626
+7.923076923076923 2329.04941935778
+11.461538461538462 2276.056248292972
+8.923076923076923 2305.7227573078167
+12.846153846153847 2269.1408371340585
+10.384615384615385 2285.4901025333397
+9.923076923076923 2290.4959120505773
+%

refset/script_plotParetoFrontCSI.m

@@ -0,0 +1,27 @@
+% script to create Figure 6 of the paper 
+
+clear
+clc
+
+% load reference sets
+
+load -ascii optLakeComo_perfectF.reference  % perfect inflow forecast cumulated over 51 days 
+load -ascii optLakeComo_base.reference      % baseline with standard policy inputs (doy, level)
+load -ascii optLakeComo_inflowF.reference   % CSI inflow forecast cumulated over 51 days 
+load -ascii optLakeComo_precipF.reference   % CSI seasonal precipitation forecast
+load -ascii optLakeComo_SST.reference       % CSI pre-season observed SST
+
+% Pareto front
+figure; 
+plot( optLakeComo_base(:,end-1), optLakeComo_base(:,end), 'o', ...
+    'Color', [.5 .5 .5], 'MarkerFaceColor', [.5 .5 .5] );
+hold on; plot( optLakeComo_perfectF(:,end-1), optLakeComo_perfectF(:,end), 'ko', ...
+    'MarkerFaceColor', [0 0 0] );
+hold on; plot( optLakeComo_precipF(:,end-1), optLakeComo_precipF(:,end), 'bo', ...
+    'MarkerFaceColor', [0 0 1] );
+hold on; plot( optLakeComo_inflowF(:,end-1), optLakeComo_inflowF(:,end), 'ro', ...
+    'MarkerFaceColor', [1 0 0] );
+hold on; plot( optLakeComo_SST(:,end-1), optLakeComo_SST(:,end), 'go', ...
+    'MarkerFaceColor', [0 1 0] );
+legend('baseline', 'perfect forecast', 'precipitation forecast', 'inflow forecast', 'observed pre-season SST');
+axis([5 15 2200 3200]); xlabel('J^{flood} (flooding days)'); ylabel('J^{supply} (deficit)')