31:T81a,{ //let ys = JSON.parse(JSON.stringify(daily)).sort((a, b) => a.Power - b.Power) let q50 = power_exceed[10].Power //d3.quantile(ys, 0.5, y => y.Power) let q90 = power_exceed[18].Power //d3.quantile(ys, 0.1, y => y.Power) let q95 = power_exceed[19].Power //d3.quantile(ys, 0.05, y => y.Power) let colour = d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[1] var data = [{ x:power_exceed.map(v => v.Exceedance),//(v,i) => 1-(i/ys.length)), y:power_exceed.map(v => v.Power), type: 'scatter', mode: 'lines', hoverinfo: 'none', showlegend:false, line: {color: make_trans(colour,0.8), shape: 'spline'} },{ x: [0.5], y: [q50], text: `P50 ${d3.format('.0f')(q50)}MW ${d3.format('.1f')(q50*24*30/1000)}GWh/mn`, textposition: 'top-right', marker: {color: 'black', size: 10}, showlegend: false, mode: 'markers+text', hoverinfo: 'text' },{ x: [0.9], y: [q90], text: `P90 ${d3.format('.0f')(q90)}MW ${d3.format('.1f')(q90*24*30/1000)}GWh/mn`, textposition: 'bottom-left', marker: {color: 'black', size: 10}, showlegend: false, mode: 'markers+text', hoverinfo: 'text' },{ x: [0.95], y: [q95], text: `P95 ${d3.format('.0f')(q95)}MW ${d3.format('.1f')(q95*24*30/1000)}GWh/mn`, textposition: 'top-right', marker: {color: 'black', size: 10}, showlegend: false, mode: 'markers+text', hoverinfo: 'text' }] var layout = { width: width*.8, height:300, font: font, barmode:'overlay', margin: {l: 65,r: 10,b: 45,t: 40}, title:'Daily Power Exceedance', showlegend: true, legend: {orientation: 'h',x: 0,xanchor: 'left',y: -0.1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', range: [0,1.15], tickformat:'%.0f' }, yaxis: { title: 'MW', zeroline: false, showline: false, tickformat: '.0f', separatethousands:true } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div} ` }32:T64b,{ var sum; function trace(name, field, colour) { return { name: name, x:calmonthly.map(mnth=>mnth.WaterMonth), y:calmonthly.map(mnth=>(mnth[field]/1000)), text:calmonthly.map(mnth=>(d3.format('.0f')(mnth[field]/1000))), type: 'scatter', mode: field=='Energy_P90'?'lines+markers+text':'lines+markers', textposition:'outside', hoverinfo: 'x+y', showlegend:true, line: {color: make_trans(colour,0.8)}, marker: {color: make_trans(colour,1), size: 8}, textposition:'top-right' } } var data = [ trace('P50', 'Energy_median', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[1]), trace('P90 (GWh/day)', 'Energy_P90', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[6]), trace('P95', 'Energy_P95', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[11]) ] var layout = { width: width*.8, height:350, font: font, barmode:'overlay', margin: {l: 65,r: 10,b: 45,t: 40}, title:'Monthly Energy Probability', showlegend: true, legend: {x: 0,xanchor: 'left',y: 1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Oct','Nov','Dec','Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep'], }, yaxis: { title: 'GWh', range:[0,150], zeroline: false, showline: false, separatethousands:true, tickformat: '.0f' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div} ` }33:T636,{ var sum; function trace(name, field, colour) { return { name: name, x:calmonthly.map(mnth=>mnth.WaterMonth), y:calmonthly.map(mnth=>(mnth[field])), text:calmonthly.map(mnth=>(d3.format('.0f')(mnth[field]))), type: 'scatter', mode: field=='Power_P90'?'lines+markers+text':'lines+markers', textposition:'outside', hoverinfo: 'x+y', showlegend:true, line: {color: make_trans(colour,0.8)}, marker: {color: make_trans(colour,1), size: 8}, textposition:'top-right' } } var data = [ trace('P50', 'Power_median', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[1]), trace('P90 (MW)', 'Power_P90', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[6]), trace('P95', 'Power_P95', d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[11]) ] var layout = { width: width*.8, height:350, font: font, barmode:'overlay', margin: {l: 65,r: 10,b: 45,t: 40}, title:'Monthly Power Probability', showlegend: true, legend: {x: 0,xanchor: 'left',y: 1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Oct','Nov','Dec','Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep'], }, yaxis: { title: 'MW', range:[0,190], zeroline: false, showline: false, separatethousands:true, tickformat: '.0f' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div} ` }34:T6bb,{ var data = [ { x:calmonthly.map(mnth=>mnth.WaterMonth), y:calmonthly.map(mnth=>(mnth.Energy_P95/1000)), text:calmonthly.map(mnth=>(d3.format('.0f')(mnth.Energy_P95/1000))), type: 'scatter', mode: 'lines+markers+text', textposition:'outside', hoverinfo: 'x+y', name: '10m Minimum Head', showlegend:true, line: {color: make_trans("#a9a9b4",0.8)}, marker: {color: make_trans("#a9a9b4",1), size: 8}, textposition:'top-right' }, { x:calmonthly_extendhead.map(mnth=>mnth.WaterMonth), y:calmonthly_extendhead.map(mnth=>(mnth.Energy_P95/1000)), text:calmonthly_extendhead.map(mnth=>(d3.format('.0f')(mnth.Energy_P95/1000))), type: 'scatter', mode: 'lines+markers+text', textposition:'outside', hoverinfo: 'x+y', name: '9m Minimum Head', showlegend:true, line: {color: make_trans("#d66000",0.8)}, marker: {color: make_trans("#d66000",1), size: 8}, textposition:'top-right' } ] var layout = { width: width*.8, height:300, font: font, barmode:'overlay', margin: {l: 65,r: 10,b: 45,t: 40}, title:'P95 Monthly Energy', showlegend: true, legend: {x: 0,xanchor: 'left',y: 1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Oct','Nov','Dec','Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep'], }, yaxis: { title: 'GWh', range:[0,110], zeroline: false, showline: false, separatethousands:true, tickformat: '.0f' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div} ` }35:T6a1,{ var data = [ { x:calmonthly.map(mnth=>mnth.WaterMonth), y:calmonthly.map(mnth=>(mnth.Power_P95)), text:calmonthly.map(mnth=>(d3.format('.0f')(mnth.Power_P95))), type: 'scatter', mode: 'lines+markers+text', textposition:'outside', hoverinfo: 'x+y', name: '10m Minimum Head', showlegend:true, line: {color: make_trans("#a9a9b4",0.8)}, marker: {color: make_trans("#a9a9b4",1), size: 8}, textposition:'top-right' }, { x:calmonthly_extendhead.map(mnth=>mnth.WaterMonth), y:calmonthly_extendhead.map(mnth=>(mnth.Power_P95)), text:calmonthly_extendhead.map(mnth=>(d3.format('.0f')(mnth.Power_P95))), type: 'scatter', mode: 'lines+markers+text', textposition:'outside', hoverinfo: 'x+y', name: '9m Minimum Head', showlegend:true, line: {color: make_trans("#d66000",0.8)}, marker: {color: make_trans("#d66000",1), size: 8}, textposition:'top-right' } ] var layout = { width: width*.8, height:300, font: font, barmode:'overlay', margin: {l: 65,r: 10,b: 45,t: 40}, title:'P95 Monthly Power', showlegend: true, legend: {x: 0,xanchor: 'left',y: 1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Oct','Nov','Dec','Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep'], }, yaxis: { title: 'MW', range:[0,150], zeroline: false, showline: false, separatethousands:true, tickformat: '.0f' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div} ` }36:T45b,function histo95(monthNo) { let xs = monthly_p95_excess.filter(v => v.month===monthNo) var data = [{ x: xs.map(v => v.p95_excess/1000), type: 'histogram', hoverinfo: 'x+y', showlegend:false, name: '', histnorm: 'probability', }] var layout = { height: 170, width: width * 0.5, font: font, barmode:'overlay', showlegend: true, legend: {orientation:'v',xanchor: 'left', x:0, y: 1}, margin: {l: 65,r: 10,b: 45,t: 0}, xaxis: { showgrid: true, zeroline: false, ticks:'outside', title: 'GWh', range: [0,80]}, yaxis: { showgrid: false, zeroline: true, tickformat: '%,.2f', ticks:'outside', }}; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

${monthNames[monthNo-1]}
Mean: ${d3.format(',.0f')(d3.mean(xs, v => v.p95_excess)/1000)} GWh
P50: ${d3.format(',.0f')(d3.quantile(xs, 0.5, v => v.p95_excess)/1000)} GWh

` }37:T459,function histo(monthNo) { let xs = monthly_p90_excess.filter(v => v.month===monthNo) var data = [{ x: xs.map(v => v.p90_excess/1000), type: 'histogram', hoverinfo: 'x+y', showlegend:false, name: '', histnorm: 'probability', }] var layout = { height: 170, width: width * 0.5, font: font, barmode:'overlay', showlegend: true, legend: {orientation:'v',xanchor: 'left', x:0, y: 1}, margin: {l: 65,r: 10,b: 45,t: 0}, xaxis: { showgrid: true, zeroline: false, ticks:'outside', title: 'GWh', range: [0,80]}, yaxis: { showgrid: false, zeroline: true, tickformat: '%,.2f', ticks:'outside', }}; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

${monthNames[monthNo-1]}
Mean: ${d3.format(',.0f')(d3.mean(xs, v => v.p90_excess)/1000)} GWh
P50: ${d3.format(',.0f')(d3.quantile(xs, 0.5, v => v.p90_excess)/1000)} GWh

` }38:[] 2f:{"title":"Power and Energy Tranches","cells":[{"id":0,"value":"# Power and Energy Tranches","pinned":false,"hidden":false,"mode":"md"},{"id":112,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":26,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":375,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":393,"value":"html`\n

Impact of Reducing Minimum Head from 10m to 9m

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":396,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":410,"value":"$35","pinned":false,"hidden":false,"mode":"ojs"},{"id":351,"value":"html`\n

Probability of energy in excess of calendar month P90.

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":325,"value":"html`${[...Array(12).keys()].map(m => histo(m+1))}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":419,"value":"html`\n

Probability of energy in excess of calendar month P95.

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":422,"value":"html`${[...Array(12).keys()].map(m => histo95(m+1))}`","pinned":false,"hidden":false,"mode":"ojs"},{"id":424,"value":"$36","pinned":false,"hidden":false,"mode":"ojs"},{"id":275,"value":"$37","pinned":false,"hidden":false,"mode":"ojs"},{"id":261,"value":"monthly_p90_excess = monthly.map(v => {\n return {\n year:v.Year,\n month: v.Month,\n p90_excess: d3.max([v.Energy - calmonthly[waterMonth(v.Month)-1].Energy_P90,0])\n}\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":426,"value":"monthly_p95_excess = monthly.map(v => {\n return {\n year:v.Year,\n month: v.Month,\n p95_excess: d3.max([v.Energy - calmonthly[waterMonth(v.Month)-1].Energy_P95,0])\n}\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":291,"value":"power_exceed = await d3.csv(\n \"https://gis.westernpower.org/data/output_data/\" + folder + model.OutputPrefix + \"_pe_power_exceed.csv?\" +\n Math.floor(Math.random() * 1000),\n d3.autoType\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"monthly = (await d3.csv(\"https://gis.westernpower.org/data/output_data/\" + folder + model.OutputPrefix + \"_pe_monthly.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":false,"hidden":false,"mode":"ojs"},{"id":387,"value":"calmonthly_extendhead=(await d3.csv(\"https://gis.westernpower.org/data/output_data/scenarios/extendhead_pe_calmonthly.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":true,"hidden":false,"mode":"ojs"},{"id":11,"value":"calmonthly=(await d3.csv(\"https://gis.westernpower.org/data/output_data/\" + folder + model.OutputPrefix + \"_pe_calmonthly.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":false,"hidden":false,"mode":"ojs"},{"id":369,"value":"folder = model_name==='Base Case'?'':'scenarios/'","pinned":true,"hidden":false,"mode":"ojs"},{"id":362,"value":"model=models.filter(md=>md.ModelName==model_name)[0]","pinned":true,"hidden":false,"mode":"ojs"},{"id":360,"value":"model_name='Base Case'//","pinned":true,"hidden":false,"mode":"ojs"},{"id":358,"value":"models = (await d3.csv(\"https://gis.westernpower.org/data/input_data/models.csv?\" + Math.floor(Math.random() * 1000),d3.autoType))","pinned":true,"hidden":false,"mode":"ojs"},{"id":9,"value":"import {font} from \"@andyflem/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":253,"value":"function month_days(i) {\n if(i === 4 || i === 6 || i === 9 || i === 11) {\n\t\treturn 30\n\t} else if(i === 2) {\n\t\treturn 28\n\t} else {\n\t\treturn 31\n\t}\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":320,"value":"monthNames = ['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec']","pinned":true,"hidden":false,"mode":"ojs"},{"id":92,"value":"waterMonth = function(month) {\n if (month>=10) {\n return month-9\n } else {\n return month+3\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":60,"value":"make_trans=function(col,opacity) {\n let c=d3.color(col)\n c.opacity=opacity\n return c.formatRgb()\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"d3 = require(\"d3@5\", \"d3-array@2\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":7,"value":"plt = require(\"https://cdn.plot.ly/plotly-latest.min.js\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W38","stdlib":"1","titling":"auto"}