new figs

Fig_src/gnuplot/bin/databymultiplefile.bash

@@ -0,0 +1,11 @@
+#!/bin/bash
+
+# get the measure values 
+# for each filename listed in input after 1-th argument
+MEASURE="$1"
+shift
+FILENAMES="$*"
+
+for FILENAME in $FILENAMES; do
+    ../bin/data.awk $MEASURE $FILENAME
+done

Fig_src/gnuplot/bin/databymultiplefileYtl.bash

@@ -0,0 +1,11 @@
+#!/bin/bash
+
+# get the measure values 
+# for each filename listed in input after 1-th argument
+MEASURE="$1"
+shift
+FILENAMES="$*"
+
+for FILENAME in $FILENAMES; do
+    ../bin/dataYtl.awk $MEASURE $FILENAME
+done

Fig_src/gnuplot/bin/gnuplot_figs.gp

@@ -869,6 +869,716 @@ if( table ) {
 }
 
 
+reset
+
+# Delay
+fig="41"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="1"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChangesVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
+reset
+
+# Delay
+fig="42"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="32"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChangesVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
+reset
+
+# Delay
+fig="42bis"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="320"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChangesVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
+reset
+
+# Delay
+fig="45"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="1"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChanges40_kFLPreqVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
+reset
+
+# Delay
+fig="46bis"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="320"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChanges40_kFLPreqVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
+exit
+
+reset
+
+# Delay
+fig="46bis"
+measuretype="P"
+exps="1 2 3 4 5 6"
+batchnum="320"
+expbatch="Exp1-6" . "_N" . batchnum
+study="RangeStudy-". expbatch
+datafile(i)=sprintf("../RangeStudy/Results_case33bw_2times_FailBranchExpon_71_43_45_41_CTRL_OffsetNumChanges40_kFLPreqVar_Exp%d_N%s_s31415_bradipo_results.csv", i, batchnum)
+#datafile2="../RangeStudy/Results_case33bw_2times_CTRL_Exp1-5_N" . batchnum . "_s31415_bradipo_results.csv"
+title="Expected Unavailability for Day (min/day) vs Energy Consumption"
+datafiles = ""
+do for [i=1:words(exps)] { datafiles = datafiles . " " . datafile(i) }
+
+scenarios="sim_vardelta_T00_T0w1_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w1_DPT0diff10_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff1_DPw10_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw1_POLICY_0 sim_vardelta_T00_T0w10_DPT0diff10_DPw10_POLICY_0" # space separed list of scenarios to compare
+measurets="'\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=1, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=1, T^{dew}_a=10\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=1\\)' '\\(\\bar{T}_0=0, T^a_0=10, \\Delta{T}^{dew}_0=10, T^{dew}_a=10\\)'"
+serieinds="1" # list of the indexes corresponding to the scenarios to plot
+m=words(scenarios) # scenarios
+lstyles1="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+lstyles2="1 2 5 6 3 4 7 8" # lt and dt style for each scenario
+
+# check input params
+if(words(measurets)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of scenario labels (".words(measurets).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+# check input params
+if(words(serieinds)>words(scenarios)) {
+				      print "Ops! Fig. #".fig.": Number of series index (".words(serieinds).") greater than number of scenarios (".words(scenarios).")"
+				      exit
+}
+
+odata(i)=datapath(word(scenarios,i))
+serieind(i)=word(serieinds,ston(i))
+lstyle(i)=ston(word(lstyles,ston(i)))
+
+# styles
+set pointsize 1.5
+set style data lines
+lstyle(i)=ston(word(lstyles1,ston(i)))
+set for [i=1:m] style line i lt lstyle(i) lw 1 pt lstyle(i)
+lstyle2(i)=ston(word(lstyles2,ston(i)))
+set for [i=1:m] style line i+m lt lstyle2(i) dt 2 lw 1 pt lstyle2(i)
+#set style line 1 lt lstyle(1) lw 1 pt 1
+#set style line 2 lt lstyle(2) dt lstyle(2) lw 1 pt 
+
+measure="Vt_PD_50"
+measurelabel="\\(P\\)"
+measurey=measurelabel . " (MW)"
+measure2="YUsedEnergy"
+measurelabel2="\\(\\mu_E\\)"
+measurey2=measurelabel2 . " (kWh)"
+measuret(i)=sprintf("%s, %s", measurelabel, word(measurets,ston(i)))
+measuret2(i)=sprintf("%s, %s", measurelabel2, word(measurets,ston(i)))
+
+tabfile=fig."-measure-".measuretype."-".study.".tab"
+tikzfilesa=fig."-measure-".measuretype."-".study."-tikzsa.tex"
+tikzfilesanoleg=fig."-measure-".measuretype."-".study."-tikzsa-noleg.tex"
+tikzfilesalogy=fig."-measure-".measuretype."-".study."-tikzsa-logy.tex"
+tikzfile=fig."-measure-".measuretype."-".study."-tikz.tex"
+
+unset title
+set xlabel "PD (MW)"
+set ylabel measurey offset -1
+#set y2label measurey2
+#set format x "%2.1f"
+#set format y "%3.2f"
+set xtics nomirror
+#set xtics add ("0.2" 0.2)
+set ytics nomirror
+#set y2tics
+set xrange [*:*]
+set yrange [*:*]
+#set y2range [0:2200]
+#set y2tics add ("2200" 2200)
+set border 3
+#set key bottom left
+#set key at graph -0.2, graph 1 spacing 1.8 maxrows 12
+set key at graph 0.37, graph 1 spacing 1.2
+unset logscale x
+unset logscale y
+
+set term tikz
+set output tikzfile
+
+#set label "$V_{B48}$" at 180,0.83
+#set arrow from 204,0.828 to 245,0.85
+
+# plot for each scenario
+plot "< ../bin/databymultiplefileYtl.bash Yt_AllDemandAbs " . datafiles using 4:xticlabels(sprintf("%.0f", $4)) with line lt 1 dashtype 1 title "Total \\(PD\\)", \
+     "< ../bin/databymultiplefileYtl.bash Yt_AllUnsuppliedDemandAbs " . datafiles using 4 with line lt 2 dashtype 2 title "Total \\(UD\\)"
+
+if( tikzsa ) {
+    set term tikz standalone
+    set output tikzfilesa
+    replot
+    unset output
+    cmd="pdflatex  -interaction=nonstopmode " . tikzfilesa
+    system(cmd)
+    }
+
+if( aqua ) {
+    set output
+    set term aqua fig
+    replot
+}
+
+if( table ) {
+    set format x "% g"
+    set format y "%g"
+    set table tabfile
+    replot
+    unset table
+}
+
+
 exit
 
 

Fig_src/gnuplot/figs/11-measure-P-RangeStudy-Exp1_5_N1-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:02 2023
+%% Tue Nov 21 21:24:11 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/11-measure-P-RangeStudy-Exp1_5_N1-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:02 2023
+%% Tue Nov 21 21:24:11 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/12-measure-V-RangeStudy-Exp1_5_N1-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:04 2023
+%% Tue Nov 21 21:24:13 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/12-measure-V-RangeStudy-Exp1_5_N1-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:04 2023
+%% Tue Nov 21 21:24:13 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/41-measure-P-RangeStudy-Exp1-6_N1-tikz.tex

@@ -0,0 +1,63 @@
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:19 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$580000$};
+\draw[gp path] (1.872,1.813)--(2.052,1.813);
+\node[gp node right] at (1.688,1.813) {$600000$};
+\draw[gp path] (1.872,2.642)--(2.052,2.642);
+\node[gp node right] at (1.688,2.642) {$620000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.299)--(2.052,4.299);
+\node[gp node right] at (1.688,4.299) {$660000$};
+\draw[gp path] (1.872,5.127)--(2.052,5.127);
+\node[gp node right] at (1.688,5.127) {$680000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$700000$};
+\draw[gp path] (1.872,6.784)--(2.052,6.784);
+\node[gp node right] at (1.688,6.784) {$720000$};
+\draw[gp path] (1.872,7.613)--(2.052,7.613);
+\node[gp node right] at (1.688,7.613) {$740000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$760000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {691496};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {705288};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {707774};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {729065};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {737724};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {745608};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,5.603)--(3.850,6.175)--(5.829,6.278)--(7.807,7.160)--(9.786,7.518)%
+  --(11.764,7.845);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.714)--(3.850,2.401)--(5.829,2.517)--(7.807,3.457)--(9.786,3.847)%
+  --(11.764,4.323);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables

Fig_src/gnuplot/figs/41-measure-P-RangeStudy-Exp1-6_N1-tikzsa.tex

@@ -0,0 +1,88 @@
+\documentclass[10pt]{article}
+\usepackage[T1]{fontenc}
+\usepackage{textcomp}
+
+\usepackage{ifluatex,ifxetex}
+\ifluatex
+  \usepackage{fontspec}
+  \setmainfont{Free Serif}
+\else\ifxetex
+  \usepackage{fontspec}
+  \setmainfont{FreeSerif}
+\else
+  \usepackage{ucs}
+  \SetUnicodeOption{mathletters}
+  \usepackage[utf8x]{inputenc}
+\fi\fi
+
+\usepackage{gnuplot-lua-tikz}
+\pagestyle{empty}
+\usepackage[active,tightpage]{preview}
+\PreviewEnvironment{tikzpicture}
+\setlength\PreviewBorder{\gpbboxborder}
+\begin{document}
+
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:19 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$580000$};
+\draw[gp path] (1.872,1.813)--(2.052,1.813);
+\node[gp node right] at (1.688,1.813) {$600000$};
+\draw[gp path] (1.872,2.642)--(2.052,2.642);
+\node[gp node right] at (1.688,2.642) {$620000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.299)--(2.052,4.299);
+\node[gp node right] at (1.688,4.299) {$660000$};
+\draw[gp path] (1.872,5.127)--(2.052,5.127);
+\node[gp node right] at (1.688,5.127) {$680000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$700000$};
+\draw[gp path] (1.872,6.784)--(2.052,6.784);
+\node[gp node right] at (1.688,6.784) {$720000$};
+\draw[gp path] (1.872,7.613)--(2.052,7.613);
+\node[gp node right] at (1.688,7.613) {$740000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$760000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {691496};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {705288};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {707774};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {729065};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {737724};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {745608};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,5.603)--(3.850,6.175)--(5.829,6.278)--(7.807,7.160)--(9.786,7.518)%
+  --(11.764,7.845);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.714)--(3.850,2.401)--(5.829,2.517)--(7.807,3.457)--(9.786,3.847)%
+  --(11.764,4.323);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables
+\end{document}

Fig_src/gnuplot/figs/41-measure-P-RangeStudy-Exp1-6_N1.tab

@@ -0,0 +1,22 @@
+
+# Curve 0 of 2, 6 points
+# Curve title: "Total \(PD\)"
+# x y type
+ 0 691496  i
+ 1 705288  i
+ 2 707774  i
+ 3 729065  i
+ 4 737724  i
+ 5 745608  i
+
+
+# Curve 1 of 2, 6 points
+# Curve title: "Total \(UD\)"
+# x y type
+ 0 597602  i
+ 1 614195  i
+ 2 616979  i
+ 3 639678  i
+ 4 649099  i
+ 5 660574  i
+

Fig_src/gnuplot/figs/42-measure-P-RangeStudy-Exp1-6_N32-tikz.tex

@@ -0,0 +1,57 @@
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:20 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$620000$};
+\draw[gp path] (1.872,2.228)--(2.052,2.228);
+\node[gp node right] at (1.688,2.228) {$630000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$650000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$660000$};
+\draw[gp path] (1.872,7.198)--(2.052,7.198);
+\node[gp node right] at (1.688,7.198) {$670000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$680000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {647206};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {673037};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {650855};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {655832};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {658620};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {658522};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,4.366)--(3.850,7.576)--(5.829,4.819)--(7.807,5.438)--(9.786,5.784)%
+  --(11.764,5.772);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.193)--(3.850,1.371)--(5.829,1.658)--(7.807,2.041)--(9.786,2.413)%
+  --(11.764,2.633);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables

Fig_src/gnuplot/figs/42-measure-P-RangeStudy-Exp1-6_N32-tikzsa.tex

@@ -0,0 +1,82 @@
+\documentclass[10pt]{article}
+\usepackage[T1]{fontenc}
+\usepackage{textcomp}
+
+\usepackage{ifluatex,ifxetex}
+\ifluatex
+  \usepackage{fontspec}
+  \setmainfont{Free Serif}
+\else\ifxetex
+  \usepackage{fontspec}
+  \setmainfont{FreeSerif}
+\else
+  \usepackage{ucs}
+  \SetUnicodeOption{mathletters}
+  \usepackage[utf8x]{inputenc}
+\fi\fi
+
+\usepackage{gnuplot-lua-tikz}
+\pagestyle{empty}
+\usepackage[active,tightpage]{preview}
+\PreviewEnvironment{tikzpicture}
+\setlength\PreviewBorder{\gpbboxborder}
+\begin{document}
+
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:20 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$620000$};
+\draw[gp path] (1.872,2.228)--(2.052,2.228);
+\node[gp node right] at (1.688,2.228) {$630000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$650000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$660000$};
+\draw[gp path] (1.872,7.198)--(2.052,7.198);
+\node[gp node right] at (1.688,7.198) {$670000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$680000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {647206};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {673037};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {650855};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {655832};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {658620};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {658522};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,4.366)--(3.850,7.576)--(5.829,4.819)--(7.807,5.438)--(9.786,5.784)%
+  --(11.764,5.772);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.193)--(3.850,1.371)--(5.829,1.658)--(7.807,2.041)--(9.786,2.413)%
+  --(11.764,2.633);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables
+\end{document}

Fig_src/gnuplot/figs/42-measure-P-RangeStudy-Exp1-6_N32.tab

@@ -0,0 +1,22 @@
+
+# Curve 0 of 2, 6 points
+# Curve title: "Total \(PD\)"
+# x y type
+ 0 647206  i
+ 1 673037  i
+ 2 650855  i
+ 3 655832  i
+ 4 658620  i
+ 5 658522  i
+
+
+# Curve 1 of 2, 6 points
+# Curve title: "Total \(UD\)"
+# x y type
+ 0 621673  i
+ 1 623103  i
+ 2 625419  i
+ 3 628497  i
+ 4 631490  i
+ 5 633266  i
+

Fig_src/gnuplot/figs/42bis-measure-P-RangeStudy-Exp1-6_N320-tikz.tex

@@ -0,0 +1,57 @@
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:22 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$630000$};
+\draw[gp path] (1.872,2.228)--(2.052,2.228);
+\node[gp node right] at (1.688,2.228) {$635000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$645000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$650000$};
+\draw[gp path] (1.872,7.198)--(2.052,7.198);
+\node[gp node right] at (1.688,7.198) {$655000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$660000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {655994};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {655994};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {655994};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {655994};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {655994};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {655994};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,7.445)--(3.850,7.445)--(5.829,7.445)--(7.807,7.445)--(9.786,7.445)%
+  --(11.764,7.445);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.600)--(3.850,1.600)--(5.829,1.600)--(7.807,1.600)--(9.786,1.600)%
+  --(11.764,1.600);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables

Fig_src/gnuplot/figs/42bis-measure-P-RangeStudy-Exp1-6_N320-tikzsa.tex

@@ -0,0 +1,82 @@
+\documentclass[10pt]{article}
+\usepackage[T1]{fontenc}
+\usepackage{textcomp}
+
+\usepackage{ifluatex,ifxetex}
+\ifluatex
+  \usepackage{fontspec}
+  \setmainfont{Free Serif}
+\else\ifxetex
+  \usepackage{fontspec}
+  \setmainfont{FreeSerif}
+\else
+  \usepackage{ucs}
+  \SetUnicodeOption{mathletters}
+  \usepackage[utf8x]{inputenc}
+\fi\fi
+
+\usepackage{gnuplot-lua-tikz}
+\pagestyle{empty}
+\usepackage[active,tightpage]{preview}
+\PreviewEnvironment{tikzpicture}
+\setlength\PreviewBorder{\gpbboxborder}
+\begin{document}
+
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:22 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$630000$};
+\draw[gp path] (1.872,2.228)--(2.052,2.228);
+\node[gp node right] at (1.688,2.228) {$635000$};
+\draw[gp path] (1.872,3.470)--(2.052,3.470);
+\node[gp node right] at (1.688,3.470) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$645000$};
+\draw[gp path] (1.872,5.956)--(2.052,5.956);
+\node[gp node right] at (1.688,5.956) {$650000$};
+\draw[gp path] (1.872,7.198)--(2.052,7.198);
+\node[gp node right] at (1.688,7.198) {$655000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$660000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {655994};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {655994};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {655994};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {655994};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {655994};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {655994};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,7.445)--(3.850,7.445)--(5.829,7.445)--(7.807,7.445)--(9.786,7.445)%
+  --(11.764,7.445);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.600)--(3.850,1.600)--(5.829,1.600)--(7.807,1.600)--(9.786,1.600)%
+  --(11.764,1.600);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables
+\end{document}

Fig_src/gnuplot/figs/42bis-measure-P-RangeStudy-Exp1-6_N320.tab

@@ -0,0 +1,22 @@
+
+# Curve 0 of 2, 6 points
+# Curve title: "Total \(PD\)"
+# x y type
+ 0 655994  i
+ 1 655994  i
+ 2 655994  i
+ 3 655994  i
+ 4 655994  i
+ 5 655994  i
+
+
+# Curve 1 of 2, 6 points
+# Curve title: "Total \(UD\)"
+# x y type
+ 0 632473  i
+ 1 632473  i
+ 2 632473  i
+ 3 632473  i
+ 4 632473  i
+ 5 632473  i
+

Fig_src/gnuplot/figs/45-measure-P-RangeStudy-Exp1-6_N1-tikz.tex

@@ -0,0 +1,61 @@
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:23 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$560000$};
+\draw[gp path] (1.872,1.917)--(2.052,1.917);
+\node[gp node right] at (1.688,1.917) {$580000$};
+\draw[gp path] (1.872,2.849)--(2.052,2.849);
+\node[gp node right] at (1.688,2.849) {$600000$};
+\draw[gp path] (1.872,3.781)--(2.052,3.781);
+\node[gp node right] at (1.688,3.781) {$620000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$640000$};
+\draw[gp path] (1.872,5.645)--(2.052,5.645);
+\node[gp node right] at (1.688,5.645) {$660000$};
+\draw[gp path] (1.872,6.577)--(2.052,6.577);
+\node[gp node right] at (1.688,6.577) {$680000$};
+\draw[gp path] (1.872,7.509)--(2.052,7.509);
+\node[gp node right] at (1.688,7.509) {$700000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$720000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {672207};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {681429};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {692377};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {702462};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {709982};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {717515};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,6.214)--(3.850,6.644)--(5.829,7.154)--(7.807,7.624)--(9.786,7.974)%
+  --(11.764,8.325);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.851)--(3.850,2.321)--(5.829,2.862)--(7.807,3.360)--(9.786,3.714)%
+  --(11.764,4.078);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables

Fig_src/gnuplot/figs/45-measure-P-RangeStudy-Exp1-6_N1-tikzsa.tex

@@ -0,0 +1,86 @@
+\documentclass[10pt]{article}
+\usepackage[T1]{fontenc}
+\usepackage{textcomp}
+
+\usepackage{ifluatex,ifxetex}
+\ifluatex
+  \usepackage{fontspec}
+  \setmainfont{Free Serif}
+\else\ifxetex
+  \usepackage{fontspec}
+  \setmainfont{FreeSerif}
+\else
+  \usepackage{ucs}
+  \SetUnicodeOption{mathletters}
+  \usepackage[utf8x]{inputenc}
+\fi\fi
+
+\usepackage{gnuplot-lua-tikz}
+\pagestyle{empty}
+\usepackage[active,tightpage]{preview}
+\PreviewEnvironment{tikzpicture}
+\setlength\PreviewBorder{\gpbboxborder}
+\begin{document}
+
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:23 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$560000$};
+\draw[gp path] (1.872,1.917)--(2.052,1.917);
+\node[gp node right] at (1.688,1.917) {$580000$};
+\draw[gp path] (1.872,2.849)--(2.052,2.849);
+\node[gp node right] at (1.688,2.849) {$600000$};
+\draw[gp path] (1.872,3.781)--(2.052,3.781);
+\node[gp node right] at (1.688,3.781) {$620000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$640000$};
+\draw[gp path] (1.872,5.645)--(2.052,5.645);
+\node[gp node right] at (1.688,5.645) {$660000$};
+\draw[gp path] (1.872,6.577)--(2.052,6.577);
+\node[gp node right] at (1.688,6.577) {$680000$};
+\draw[gp path] (1.872,7.509)--(2.052,7.509);
+\node[gp node right] at (1.688,7.509) {$700000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$720000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {672207};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {681429};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {692377};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {702462};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {709982};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {717515};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,6.214)--(3.850,6.644)--(5.829,7.154)--(7.807,7.624)--(9.786,7.974)%
+  --(11.764,8.325);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.851)--(3.850,2.321)--(5.829,2.862)--(7.807,3.360)--(9.786,3.714)%
+  --(11.764,4.078);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables
+\end{document}

Fig_src/gnuplot/figs/45-measure-P-RangeStudy-Exp1-6_N1.tab

@@ -0,0 +1,22 @@
+
+# Curve 0 of 2, 6 points
+# Curve title: "Total \(PD\)"
+# x y type
+ 0 672207  i
+ 1 681429  i
+ 2 692377  i
+ 3 702462  i
+ 4 709982  i
+ 5 717515  i
+
+
+# Curve 1 of 2, 6 points
+# Curve title: "Total \(UD\)"
+# x y type
+ 0 578578  i
+ 1 588679  i
+ 2 600273  i
+ 3 610961  i
+ 4 618565  i
+ 5 626374  i
+

Fig_src/gnuplot/figs/46bis-measure-P-RangeStudy-Exp1-6_N320-tikz.tex

@@ -0,0 +1,61 @@
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:25 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$580000$};
+\draw[gp path] (1.872,1.917)--(2.052,1.917);
+\node[gp node right] at (1.688,1.917) {$600000$};
+\draw[gp path] (1.872,2.849)--(2.052,2.849);
+\node[gp node right] at (1.688,2.849) {$620000$};
+\draw[gp path] (1.872,3.781)--(2.052,3.781);
+\node[gp node right] at (1.688,3.781) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$660000$};
+\draw[gp path] (1.872,5.645)--(2.052,5.645);
+\node[gp node right] at (1.688,5.645) {$680000$};
+\draw[gp path] (1.872,6.577)--(2.052,6.577);
+\node[gp node right] at (1.688,6.577) {$700000$};
+\draw[gp path] (1.872,7.509)--(2.052,7.509);
+\node[gp node right] at (1.688,7.509) {$720000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$740000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {672592};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {682807};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {693019};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {703233};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {713644};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {723876};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,5.300)--(3.850,5.776)--(5.829,6.252)--(7.807,6.728)--(9.786,7.213)%
+  --(11.764,7.690);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.844)--(3.850,1.792)--(5.829,2.321)--(7.807,2.824)--(9.786,6.116)%
+  --(11.764,6.593);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables

Fig_src/gnuplot/figs/46bis-measure-P-RangeStudy-Exp1-6_N320-tikzsa.tex

@@ -0,0 +1,86 @@
+\documentclass[10pt]{article}
+\usepackage[T1]{fontenc}
+\usepackage{textcomp}
+
+\usepackage{ifluatex,ifxetex}
+\ifluatex
+  \usepackage{fontspec}
+  \setmainfont{Free Serif}
+\else\ifxetex
+  \usepackage{fontspec}
+  \setmainfont{FreeSerif}
+\else
+  \usepackage{ucs}
+  \SetUnicodeOption{mathletters}
+  \usepackage[utf8x]{inputenc}
+\fi\fi
+
+\usepackage{gnuplot-lua-tikz}
+\pagestyle{empty}
+\usepackage[active,tightpage]{preview}
+\PreviewEnvironment{tikzpicture}
+\setlength\PreviewBorder{\gpbboxborder}
+\begin{document}
+
+\begin{tikzpicture}[gnuplot]
+%% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
+%% Tue Nov 21 21:24:25 2023
+\path (0.000,0.000) rectangle (12.500,8.750);
+\gpcolor{color=gp lt color border}
+\gpsetlinetype{gp lt border}
+\gpsetdashtype{gp dt solid}
+\gpsetlinewidth{1.00}
+\draw[gp path] (1.872,0.985)--(2.052,0.985);
+\node[gp node right] at (1.688,0.985) {$580000$};
+\draw[gp path] (1.872,1.917)--(2.052,1.917);
+\node[gp node right] at (1.688,1.917) {$600000$};
+\draw[gp path] (1.872,2.849)--(2.052,2.849);
+\node[gp node right] at (1.688,2.849) {$620000$};
+\draw[gp path] (1.872,3.781)--(2.052,3.781);
+\node[gp node right] at (1.688,3.781) {$640000$};
+\draw[gp path] (1.872,4.713)--(2.052,4.713);
+\node[gp node right] at (1.688,4.713) {$660000$};
+\draw[gp path] (1.872,5.645)--(2.052,5.645);
+\node[gp node right] at (1.688,5.645) {$680000$};
+\draw[gp path] (1.872,6.577)--(2.052,6.577);
+\node[gp node right] at (1.688,6.577) {$700000$};
+\draw[gp path] (1.872,7.509)--(2.052,7.509);
+\node[gp node right] at (1.688,7.509) {$720000$};
+\draw[gp path] (1.872,8.441)--(2.052,8.441);
+\node[gp node right] at (1.688,8.441) {$740000$};
+\draw[gp path] (1.872,0.985)--(1.872,1.165);
+\node[gp node center] at (1.872,0.677) {672592};
+\draw[gp path] (3.850,0.985)--(3.850,1.165);
+\node[gp node center] at (3.850,0.677) {682807};
+\draw[gp path] (5.829,0.985)--(5.829,1.165);
+\node[gp node center] at (5.829,0.677) {693019};
+\draw[gp path] (7.807,0.985)--(7.807,1.165);
+\node[gp node center] at (7.807,0.677) {703233};
+\draw[gp path] (9.786,0.985)--(9.786,1.165);
+\node[gp node center] at (9.786,0.677) {713644};
+\draw[gp path] (11.764,0.985)--(11.764,1.165);
+\node[gp node center] at (11.764,0.677) {723876};
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+\node[gp node center,rotate=-270] at (0.108,4.713) {\(P\) (MW)};
+\node[gp node center] at (6.818,0.215) {PD (MW)};
+\node[gp node right] at (4.248,8.256) {Total \(PD\)};
+\gpcolor{rgb color={0.580,0.000,0.827}}
+\gpsetdashtype{gp dt 1}
+\draw[gp path] (4.432,8.256)--(5.348,8.256);
+\draw[gp path] (1.872,5.300)--(3.850,5.776)--(5.829,6.252)--(7.807,6.728)--(9.786,7.213)%
+  --(11.764,7.690);
+\gpcolor{color=gp lt color border}
+\node[gp node right] at (4.248,7.887) {Total \(UD\)};
+\gpcolor{rgb color={0.000,0.620,0.451}}
+\gpsetdashtype{gp dt 2}
+\draw[gp path] (4.432,7.887)--(5.348,7.887);
+\draw[gp path] (1.872,1.844)--(3.850,1.792)--(5.829,2.321)--(7.807,2.824)--(9.786,6.116)%
+  --(11.764,6.593);
+\gpcolor{color=gp lt color border}
+\gpsetdashtype{gp dt solid}
+\draw[gp path] (1.872,8.441)--(1.872,0.985)--(11.764,0.985);
+%% coordinates of the plot area
+\gpdefrectangularnode{gp plot 1}{\pgfpoint{1.872cm}{0.985cm}}{\pgfpoint{11.764cm}{8.441cm}}
+\end{tikzpicture}
+%% gnuplot variables
+\end{document}

Fig_src/gnuplot/figs/46bis-measure-P-RangeStudy-Exp1-6_N320.tab

@@ -0,0 +1,22 @@
+
+# Curve 0 of 2, 6 points
+# Curve title: "Total \(PD\)"
+# x y type
+ 0 672592  i
+ 1 682807  i
+ 2 693019  i
+ 3 703233  i
+ 4 713644  i
+ 5 723876  i
+
+
+# Curve 1 of 2, 6 points
+# Curve title: "Total \(UD\)"
+# x y type
+ 0 598424  i
+ 1 597316  i
+ 2 608660  i
+ 3 619458  i
+ 4 690103  i
+ 5 700336  i
+

Fig_src/gnuplot/figs/51-measure-P-RangeStudy-Exp1_5_N32-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:05 2023
+%% Tue Nov 21 21:24:14 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/51-measure-P-RangeStudy-Exp1_5_N32-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:05 2023
+%% Tue Nov 21 21:24:14 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/52-measure-V-RangeStudy-Exp1_5_N32-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:06 2023
+%% Tue Nov 21 21:24:15 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/52-measure-V-RangeStudy-Exp1_5_N32-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:06 2023
+%% Tue Nov 21 21:24:15 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/53-measure-P-RangeStudy-Exp1_5_N1-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:08 2023
+%% Tue Nov 21 21:24:16 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/53-measure-P-RangeStudy-Exp1_5_N1-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:08 2023
+%% Tue Nov 21 21:24:16 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/54-measure-V-RangeStudy-Exp1_5_N1-tikz.tex

@@ -1,6 +1,6 @@
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:09 2023
+%% Tue Nov 21 21:24:18 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}

Fig_src/gnuplot/figs/54-measure-V-RangeStudy-Exp1_5_N1-tikzsa.tex

@@ -24,7 +24,7 @@
 
 \begin{tikzpicture}[gnuplot]
 %% generated with GNUPLOT 5.4p8 (Lua 5.3; terminal rev. Jun 2020, script rev. 118)
-%% Fri Nov 10 11:44:09 2023
+%% Tue Nov 21 21:24:18 2023
 \path (0.000,0.000) rectangle (12.500,8.750);
 \gpcolor{color=gp lt color border}
 \gpsetlinetype{gp lt border}