2b:null 2c:["$","title",null,{"children":"Power and Energy Simulation | Western Power | Observable"}] 31:T7f7,{ var data1 = [ { x: ['Power', 'Flow' ,'Head'], y: [ exceed.Power/180*100, exceed.FlowCanal/1100*100, (exceed.HeadTurbine1-head_min)/head_range*100, ], text: [ d3.format(",.0f")(exceed.Power/180*100) + '%', d3.format(",.0f")(exceed.FlowCanal/1100*100) + '%' , d3.format(",.0f")((exceed.HeadTurbine1-head_min)/head_range*100) + '%' ], marker: { color: ['red','#4169E1',colr(4)], }, textposition: 'outside', hoverinfo: 'none', type: 'bar' } ]; var layout1={ margin: { l: 40, r: 10,b: 30, t: 10 }, height:170, width:width*.2, font: font, yaxis: { range: [0, 120], title: '%', } } const div1 = DOM.element('div'); plt.newPlot(div1, data1, layout1,{staticPlot: true}); return html`

Power: ${d3.format(',.0f')(exceed.Power)} MW
Monthly Energy: ${d3.format(',.0f')(exceed.Power*24*30/1000)} GWh
Turbines: ${Math.ceil(exceed.Turbines)}
Head: ${d3.format(',.1f')(exceed.HeadTurbine1)}m
Flow Exceedance: ${d3.format(',.0f')(exceed.ExceedanceRnd*100)}%

River Flow: ${format_flow(exceed.Flow)}
Environmental Flow: ${format_flow(exceed.EWRTotal + exceed.FlowSpill)} (${d3.format(',.0f')((exceed.EWRTotal + exceed.FlowSpill)/exceed.Flow*100)}%)
Canal Flow: ${format_flow(exceed.FlowCanal)}
Left (Main) Fall Flow: ${format_flow(exceed.FlowChannelA)}
Turbine Flow: ${format_flow(exceed.FlowCanal/Math.ceil(exceed.Turbines))}

${div1}

`; }32:T63c,{ const svg = d3.create("svg") .attr("viewBox", [0, 0, s_width, s_height]); const parent=svg.append("g"); const {nodes, links} = sankey(data); parent.append("g") .attr("stroke", "#000") .selectAll("rect") .data(nodes) .join("rect") .attr("x", d => d.x0) .attr("y", d => d.y0) .attr("height", d => d.y1 - d.y0) .attr("width", d => d.x1 - d.x0) .attr("fill", d => d.color) const link = parent.append("g") .attr("fill", "none") .attr("stroke-opacity", 0.5) .selectAll("g") .data(links) .join("g") .style("mix-blend-mode", "multiply"); const gradient = link.append("linearGradient") .attr("id", d => (d.uid = DOM.uid("link")).id) .attr("gradientUnits", "userSpaceOnUse") .attr("x1", d => d.source.x1) .attr("x2", d => d.target.x0); gradient.append("stop") .attr("offset", "0%") .attr("stop-color", d => d.source.color); gradient.append("stop") .attr("offset", "100%") .attr("stop-color", d => d.target.color); link.append("path") .attr("d", d3.sankeyLinkHorizontal()) .attr("stroke", d => d.uid) .attr("stroke-width", d => Math.max(1, d.width)); parent.append("g") .style("font", "10px Courier New, monospace") .selectAll("text") .data(nodes) .join("text") .attr("x", d => d.x0 < width / 3*2 ? d.x1 + 6 : d.x0 - 6) .attr("y", d => (d.y1 + d.y0) / 2) .attr("dy", "0.35em") .attr("text-anchor", d => d.x0 < width / 3*2 ? "start" : "end") .text(d => d.label); return svg.node(); }33:T1a6e,data = { var turbines=Math.ceil(exceed.Turbines) var turbine_flow=exceed.FlowCanal/turbines; var links; var nodes; if(exceed.checked){ nodes = [ {label: "River " + format_flow(exceed.Flow), name: "River",color:'#4169E1'}, {label: "Left Channel " + d3.format('.0%')(exceed.FlowLeftChannel/exceed.Flow), name: "Left Channel",color:'#B0C4DE'}, {label: "Intake " + d3.format('.0%')(exceed.FlowCanal/exceed.Flow),name: "Intake",color:"#CD5C5C"}, {label: "Canal " + d3.format('.0%')(exceed.FlowCanal/exceed.Flow), name: "Canal",color:"#CD5C5C"}, {label: "Turbine 1 " + d3.format('.0%')(turbine_flow/exceed.Flow), name: "Turbine 1",color:"red"}, {label: "Channel A " + d3.format('.0%')(exceed.FlowChannelA/exceed.Flow), name: "Channel A",color:'#0000FF'}, {label: "Channel C " + d3.format('.0%')(exceed.FlowChannelC/exceed.Flow), name: "Channel C",color:'#0000FF'}, {label: "Channel D " + d3.format('.0%')(exceed.FlowChannelD/exceed.Flow), name: "Channel D",color:'#0000FF'}, {label: "Left Fall " + d3.format('.0%')((exceed.FlowChannelA+exceed.FlowChannelC+exceed.FlowChannelD)/exceed.Flow), name: "Left Fall",color:'#0000FF'}, {label: "Right Channel "+ d3.format('.0%')((exceed.Flow-exceed.FlowLeftChannel)/exceed.Flow), name: "Right Channel",color:'#B0C4DE'}, {label: "Channel E " + d3.format('.0%')(exceed.FlowChannelE/exceed.Flow), name: "Channel E",color:'#1E90FF'}, {label: "Channel FG " + d3.format('.0%')(exceed.FlowChannelFG/exceed.Flow), name: "Channel FG",color:'#1E90FF'}, {label: "Right Fall "+ d3.format('.0%')((exceed.FlowChannelE + exceed.FlowChannelFG)/exceed.Flow), name: "Right Fall",color:'#1E90FF'}, {label: "Gorge "+ d3.format('.0%')((exceed.FlowChannelE + exceed.FlowChannelFG +exceed.FlowChannelA + exceed.FlowChannelC + exceed.FlowChannelD)/exceed.Flow), name: "Gorge",color:'#B0C4DE'}, {label: "Zambezi ", name: "Zambezi",color:'#4169E1'} ] links = [ {source: "River", target: "Left Channel",value: exceed.FlowLeftChannel}, {source: "River", target: "Right Channel",value: exceed.Flow-exceed.FlowLeftChannel}, {source: "Left Channel", target: "Intake",value: exceed.FlowCanal}, {source: "Left Channel", target: "Channel A",value: exceed.FlowChannelA}, {source: "Left Channel", target: "Channel C",value: exceed.FlowChannelC}, {source: "Left Channel", target: "Channel D",value: exceed.FlowChannelD}, {source: "Channel A", target: "Left Fall",value: exceed.FlowChannelA}, {source: "Channel C", target: "Left Fall",value: exceed.FlowChannelC}, {source: "Channel D", target: "Left Fall",value: exceed.FlowChannelD}, {source: "Left Fall", target: "Gorge",value: exceed.FlowChannelA+exceed.FlowChannelC+exceed.FlowChannelD}, {source: "Intake", target: "Canal",value: exceed.FlowCanal}, {source: "Canal", target: "Turbine 1",value: turbine_flow}, {source: "Turbine 1", target: "Zambezi",value: turbine_flow}, {source: "Right Channel", target: "Channel E",value: exceed.FlowChannelE}, {source: "Right Channel", target: "Channel FG",value: exceed.FlowChannelFG}, {source: "Channel E", target: "Right Fall",value: exceed.FlowChannelE}, {source: "Channel FG", target: "Right Fall",value: exceed.FlowChannelFG}, {source: "Right Fall", target: "Gorge",value: exceed.FlowChannelE + exceed.FlowChannelFG}, {source: "Gorge", target: "Zambezi",value: exceed.FlowChannelE + exceed.FlowChannelFG +exceed.FlowChannelA + exceed.FlowChannelC + exceed.FlowChannelD} ]; if (exceed.FlowTurbine2>0){ nodes.splice(6,0,{label: "Turbine 2 " + d3.format('.0%')(turbine_flow/exceed.Flow), name: "Turbine 2",color:"red"}) links.splice(13,0,{source: "Canal", target: "Turbine 2",value: turbine_flow}) links.splice(15,0,{source: "Turbine 2", target: "Zambezi",value: turbine_flow}) } if (exceed.FlowTurbine3>0){ nodes.splice(7,0,{label: "Turbine 3 " + d3.format('.0%')(turbine_flow/exceed.Flow),name: "Turbine 3",color:"red"}) links.splice(17,0,{source: "Canal", target: "Turbine 3",value: turbine_flow}) links.splice(19,0,{source: "Turbine 3", target: "Zambezi",value: turbine_flow}) } if (exceed.FlowTurbine4>0){ nodes.splice(8,0,{label: "Turbine 4 " + d3.format('.0%')(turbine_flow/exceed.Flow), name: "Turbine 4",color:"red"}) links.splice(21,0,{source: "Canal", target: "Turbine 4",value: turbine_flow}) links.splice(23,0,{source: "Turbine 4", target: "Zambezi",value: turbine_flow}) } } else { nodes = [ {label: "River",name: "River",color:'#4169E1'}, {label: "Left Channel", name: "Left Channel", color:'#B0C4DE'}, {label: "Channel A", name: "Channel A",color:'#0000FF'}, {label: "Channel C", name: "Channel C",color:'#0000FF'}, {label: "Channel D", name: "Channel D",color:'#0000FF'}, {label: "Left Fall", name: "Left Fall",color:'#0000FF'}, {label: "Right Channel", name: "Right Channel",color:'#B0C4DE'}, {label: "Channel FG", name: "Channel FG",color:'#1E90FF'}, {label: "Channel E", name: "Channel E",color:'#1E90FF'}, {label: "Right Fall", name: "Right Fall",color:'#1E90FF'}, {label: "Gorge", name: "Gorge",color:'#B0C4DE'}, {label: "Zambezi", name: "Zambezi",color:'#4169E1'} ] links = [ {source: "River", target: "Left Channel",value: natural.FlowLeftChannel}, {source: "River", target: "Right Channel",value: natural.FlowRightChannel}, {source: "Left Channel", target: "Channel A",value: natural.FlowChannelA}, {source: "Left Channel", target: "Channel C",value: natural.FlowChannelC}, {source: "Left Channel", target: "Channel D",value: natural.FlowChannelD}, {source: "Channel A", target: "Left Fall",value: natural.FlowChannelA}, {source: "Channel C", target: "Left Fall",value: natural.FlowChannelC}, {source: "Channel D", target: "Left Fall",value: natural.FlowChannelD}, {source: "Left Fall", target: "Gorge",value: natural.FlowChannelA+natural.FlowChannelC+natural.FlowChannelD}, {source: "Right Channel", target: "Channel FG",value: natural.FlowChannelFG}, {source: "Right Channel", target: "Channel E",value: natural.FlowChannelE}, {source: "Channel FG", target: "Right Fall",value: natural.FlowChannelFG}, {source: "Channel E", target: "Right Fall",value: natural.FlowChannelE}, {source: "Right Fall", target: "Gorge",value: natural.FlowChannelE + natural.FlowChannelFG}, {source: "Gorge", target: "Zambezi",value: natural.FlowChannelE + natural.FlowChannelFG +natural.FlowChannelA + natural.FlowChannelC + natural.FlowChannelD} ]; } return {nodes, links}; }34:T541,html ` `35:[{"id":"6cb6e657a7f88d181b73cd865eba2fe6423b5058b0a9f299deeab3f78c93167cd4af140a76feb6b61a9f9548c896bda0c5e7330699b6c757a9e322507d0fb771","href":"https://static.observableusercontent.com/files/6cb6e657a7f88d181b73cd865eba2fe6423b5058b0a9f299deeab3f78c93167cd4af140a76feb6b61a9f9548c896bda0c5e7330699b6c757a9e322507d0fb771","name":"ewr_sections-1.jpg","type":"image/jpeg","lastModified":1592990709501.828,"size":83780}] 2f:{"id":"d7fe980b106768f0","slug":"simulation","title":"Power and Energy Simulation","thumbnail":"7de9d44295d47351cc25ff5da3f18c1efea7c2a2c7043324ba3c97245498248a","access_level":"public","create_time":"$D2020-06-23T09:20:27.844Z","update_time":"$D2025-07-23T12:24:51.681Z","view_time":null,"first_public_time":"$D2020-06-24T09:06:17.711Z","fork":false,"stars":0,"starred":false,"owner":{"type":"team","id":"3aaf4970407a12be","login":"westernpower","name":"Western Power","avatar_url":"https://avatars-next.observableusercontent.com/0c11a3086911e8caa71a7767a24de09ff93d85f36798166fa6c2d0b61107b761","tier":"starter_2024"},"version":455,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"96924893937b5930","login":"andyflem","name":"Andy Fleming","avatar_url":"https://avatars-next.observableusercontent.com/fd3844143a4c29872f0df28a8ebce77bb56cebf2fa74bb669872e4def82ad809","tier":null}],"members":[{"type":"team","id":"3aaf4970407a12be","login":"westernpower","name":"Western Power","avatar_url":"https://avatars-next.observableusercontent.com/0c11a3086911e8caa71a7767a24de09ff93d85f36798166fa6c2d0b61107b761","tier":"starter_2024","role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"title":"Power and Energy Simulation","cells":[{"id":0,"value":"

Power and Energy Simulation

Ngonye Falls Hydroelectric Project

","pinned":false,"hidden":false,"mode":"html","author":"$undefined","summary":"$undefined"},{"id":34,"value":"viewof exceed = {\n const form = html``;\n form.oninput = () => {\n form.querySelector(\"[id=flow_disp]\").textContent = \n format_flow(pefdc.filter(f=>f.ExceedanceRnd*100==(100-form.exceed_select.valueAsNumber))[0].Flow);\n form.value = pefdc.filter(f=>f.ExceedanceRnd*100==(100-form.exceed_select.valueAsNumber))[0];\n form.value.checked=form.chckbox.checked;\n };\n form.oninput();\n return form;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":116,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":20,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":156,"value":"md`## Appendix`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":429,"value":"container = function(captionText, width_over) {\n const div=DOM.element('div')\n div.setAttribute(\"class\",\"border m-2 p-2\")\n div.setAttribute(\"style\",`width: ${width_over?width_over:(width>800?800:width)}px`)\n if (captionText) {\n const caption=DOM.element('figcaption')\n caption.setAttribute(\"style\",\"max-width:100%\")\n caption.append(captionText)\n div.append(caption)\n }\n return div\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":310,"value":"md`### Data`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":420,"value":"calmonths=(await d3.csv(\"https://gis.westernpower.org/data/output_data/\" + model.OutputPrefix + \"_pe_calmonthly.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":24,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":44,"value":"model_name='Base Case'","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":46,"value":"model=models.filter(md=>md.ModelName==model_name)[0]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":42,"value":"models = (await d3.csv(\"https://gis.westernpower.org/data/input_data/models.csv?\" + Math.floor(Math.random() * 1000),d3.autoType))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":352,"value":"text_col={\n if(exceed.checked){\n return 'black'\n } else {\n return '#DCDCDC'\n } \n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":212,"value":"naturals=natural_split.filter(ns=>ns.Flow>=exceed.Flow)[0]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":229,"value":"natural={\n var tmp = {\n FlowChannelA: exceed.Flow * naturals.A,\n FlowChannelC: exceed.Flow * naturals.C,\n FlowChannelD: exceed.Flow * naturals.D, \n FlowChannelE: exceed.Flow * naturals.E,\n FlowChannelFG: exceed.Flow * naturals.FG\n }\n tmp.FlowLeftChannel = tmp.FlowChannelA + tmp.FlowChannelC + tmp.FlowChannelD\n tmp.FlowRightChannel=tmp.FlowChannelE + tmp.FlowChannelFG\n return tmp\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":210,"value":"natural_split=(await d3.csv(\"https://gis.westernpower.org/data/input_data/flowsplit_natural.csv?\" + Math.floor(Math.random() * 1000),d3.autoType))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":40,"value":"pefdc = (await d3.csv(\"https://gis.westernpower.org/data/output_data/\" + model.OutputPrefix + \"_pe_fdc.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":183,"value":"head_range=head_max-head_min","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":171,"value":"head_min=d3.min(pefdc.map(f=>f.HeadTurbine1))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":177,"value":"head_max=d3.max(pefdc.map(f=>f.HeadTurbine1))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":166,"value":"md`### Formatting`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":27,"value":"sankey = {\n const sankey = d3.sankey()\n .nodeId(d => d.name)\n .nodeAlign(d3[\"sankeyCenter\"])\n .nodeWidth(12)\n .nodePadding(4)\n .nodeSort(null)\n .linkSort(null)\n .extent([[1, 5], [s_width - 1, s_height - 5]]);\n return ({nodes, links}) => sankey({\n nodes: nodes.map(d => Object.assign({}, d)),\n links: links.map(d => Object.assign({}, d))\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":54,"value":"function format_flow(flow) {\n return d3.format(',.0f')(flow) + 'm\\u00b3/s';\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":18,"value":"color = {\n const color = d3.scaleOrdinal(d3.schemePastel2);\n return name => color(name.replace(/ .*/, \"\"));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":29,"value":"s_width=900\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":32,"value":"s_height=300","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":256,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":341,"value":"md`### Imports`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":4,"value":"import {commonStyle} from \"@westernpower/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":6,"value":"import {colr} from \"@westernpower/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":8,"value":"import {font} from \"@westernpower/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":10,"value":"import {margin} from \"@westernpower/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":12,"value":"import {plotTitle} from \"@westernpower/commonlibrary\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":14,"value":"d3 = require(\"d3@5\", \"d3-sankey@0.12\",\"d3-hsv@0.1\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":16,"value":"plt = require(\"https://cdn.plot.ly/plotly-latest.min.js\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 30:"$2f:body"

Power and Energy Simulation

Ngonye Falls Hydroelectric Project

River Flow

Simulate Powerstation