2b:null 2c:["$","title",null,{"children":"Price capture in a coupled solar / wind + BESS operation in Ukraine | atsokol | Observable"}] 31:T4d7,Plot.plot((() => { const n = 3; // number of facet columns const keys = Array.from(months.keys()); const index = new Map(keys.map((key, i) => [key,i])); const fx = (key) => index.get(key) % n; const fy = (key) => Math.floor(index.get(key) / n); return { title: "Wind yield variation", subtitle: "Median (P50) capacity factor and P10-P90 range", height: 700, width: 800, y: {insetTop: 10, label: null, tickFormat: d3.format(".0%")}, x: {label: "hour of day"}, fx: {padding: 0.03, label: null, tickFormat: d=> ''}, fy: {label: null, tickFormat: d=> ''}, marks: [ Plot.areaY(windCapFactor, { x: "hour", y1: "capFactor_p10", y2: "capFactor_p90", fx: (d) => fx(d.month), fy: (d) => fy(d.month), curve: "catmull-rom", fill: "#bae4bc", fillOpacity: 0.4 }), Plot.lineY(windCapFactor, { x: "hour", y: "capFactor_p50", fx: (d) => fx(d.month), fy: (d) => fy(d.month), curve: "catmull-rom", stroke: "steelblue", tip: true }), Plot.text(keys, {fx, fy, text: d => months.get(d), frameAnchor: "top-left", dx: 6, dy: 6}), Plot.frame() ] }; })())32:T4da,Plot.plot((() => { const n = 3; // number of facet columns const keys = Array.from(months.keys()); const index = new Map(keys.map((key, i) => [key,i])); const fx = (key) => index.get(key) % n; const fy = (key) => Math.floor(index.get(key) / n); return { title: "Solar yield variation", subtitle: "Median (P50) capacity factor and P10-P90 range", height: 700, width: 800, y: {insetTop: 10, label: null, tickFormat: d3.format(".0%")}, x: {label: "hour of day"}, fx: {padding: 0.03, label: null, tickFormat: d=> ''}, fy: {label: null, tickFormat: d=> ''}, marks: [ Plot.areaY(solarCapFactor, { x: "hour", y1: "capFactor_p10", y2: "capFactor_p90", fx: (d) => fx(d.month), fy: (d) => fy(d.month), curve: "catmull-rom", fill: "#bae4bc", fillOpacity: 0.4 }), Plot.lineY(solarCapFactor, { x: "hour", y: "capFactor_p50", fx: (d) => fx(d.month), fy: (d) => fy(d.month), curve: "catmull-rom", stroke: "steelblue", tip: true }), Plot.text(keys, {fx, fy, text: d => months.get(d), frameAnchor: "top-left", dx: 6, dy: 6}), Plot.frame() ] }; })())33:T4e4,Plot.plot({ title: "Coupled wind farm + battery operation", marginRight: 30, width: 1000, x: {label: null, ticks: "6 hours"}, y: { nice: true, grid: true, label: "MW"}, color: {legend: true}, marks: [ Plot.rectY( aq.from(wind_dispatch.filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7))) .select('time','charge','discharge','yield') .derive({charge: d=> -d.charge}) .fold(['yield','charge','discharge']), { x: "time", y: "value", fill: "key", tip: true } ), Plot.ruleY([0]), Plot.lineY( price_2023.filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7)), { x: "time", y: "price", y2: "price", stroke: "green", strokeWidth: 2, }), Plot.lineY( wind_dispatch .filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7)) .map(d => ({time: d.time, sales: d.sale_direct + d.discharge})) .filter(d => typeof d.sales === "number" && !isNaN(d.sales)) , { x: "time", y: "sales", stroke: "black", strokeWidth: 2, curve: "step" }), ] })34:T4e7,Plot.plot({ title: "Coupled solar farm + battery operation", marginRight: 30, width: 1000, x: {label: null, ticks: "6 hours"}, y: { nice: true, grid: true, label: "MW"}, color: {legend: true}, marks: [ Plot.rectY( aq.from(solar_dispatch.filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7))) .select('time','charge','discharge','yield') .derive({charge: d=> -d.charge}) .fold(['yield','charge','discharge']), { x: "time", y: "value", fill: "key", tip: true } ), Plot.ruleY([0]), Plot.lineY( price_2023.filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7)), { x: "time", y: "price", y2: "price", stroke: "green", strokeWidth: 2, }), Plot.lineY( solar_dispatch .filter(d => d.time >= startDate && d.time <= d3.timeDay.offset(startDate, 7)) .map(d => ({time: d.time, sales: d.sale_direct + d.discharge})) .filter(d => typeof d.sales === "number" && !isNaN(d.sales)) , { x: "time", y: "sales", stroke: "black", strokeWidth: 2, curve: "step" }), ] })35:[{"id":"3faec4cf91484299f3908c8751eb4b9c75ffaae67dd98fb4d1ae053d1e5f08556755f10eea4e208e7a983e4e96d63763e6844686fc6cf046f9c05b4112229821","href":"https://static.observableusercontent.com/files/3faec4cf91484299f3908c8751eb4b9c75ffaae67dd98fb4d1ae053d1e5f08556755f10eea4e208e7a983e4e96d63763e6844686fc6cf046f9c05b4112229821","name":"solar_dispatch_2023.csv","type":"text/csv","lastModified":1734200428309.56,"size":551285},{"id":"929aa9abbb0b54fbdbd9704950c9d172e0c97855edeba212bf72a76fecc537285298383f558b047f1eeadf4756b4e566c33ed0a154827110724e7e03897ed253","href":"https://static.observableusercontent.com/files/929aa9abbb0b54fbdbd9704950c9d172e0c97855edeba212bf72a76fecc537285298383f558b047f1eeadf4756b4e566c33ed0a154827110724e7e03897ed253","name":"wind_dispatch_2023.csv","type":"text/csv","lastModified":1734206454684.644,"size":647858},{"id":"d276aa6decb5698b1d6d6ca458aa485b23e8858073018f221943f429e0977e6810cfcd673f1f7222331e819ea7bea0d74792171621f18777b5e583c6ec76755b","href":"https://static.observableusercontent.com/files/d276aa6decb5698b1d6d6ca458aa485b23e8858073018f221943f429e0977e6810cfcd673f1f7222331e819ea7bea0d74792171621f18777b5e583c6ec76755b","name":"Syvash wind yield 2014-2023@1.csv","type":"text/csv","lastModified":1731962613648.905,"size":4188691},{"id":"ff3de1cb40f1481e8c6c777a99b6d8c96e319dab4f914d7cac2cea378ea3125acff965ad8d12c56b87aac7d63d43c4339e5dd199347bb2ff9bb7b4b6ec2691c3","href":"https://static.observableusercontent.com/files/ff3de1cb40f1481e8c6c777a99b6d8c96e319dab4f914d7cac2cea378ea3125acff965ad8d12c56b87aac7d63d43c4339e5dd199347bb2ff9bb7b4b6ec2691c3","name":"Prohresivka solar yield 2005-2023.csv","type":"text/csv","lastModified":1734191777441.667,"size":4051930}] 2f:{"id":"980d7e1f4e0e0fd0","slug":null,"title":"Price capture in a coupled solar / wind + BESS operation in Ukraine","thumbnail":"f4cd10194f15fa864e8a9dda014b85aeed12e94724f573377993424d6225a7df","access_level":"public","create_time":"$D2024-11-18T20:31:17.050Z","update_time":"$D2025-02-27T13:31:24.400Z","view_time":null,"first_public_time":"$D2025-02-27T13:39:21.592Z","has_unpublished_edits":null,"fork":true,"stars":1,"starred":false,"owner":{"type":"team","id":"67596cf83e89a91c","login":"atsokol","name":"atsokol","avatar_url":"https://avatars-next.observableusercontent.com/377abaada84daeb487e417deeb4779c04a819d8e607bfcf58c388483815ab8e4","tier":"starter_2024"},"version":1001,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"1cc75fc7f25f3d6a","login":"atsokol","name":"Andrii T","avatar_url":"https://avatars-next.observableusercontent.com/e5b1ecc1785c25cd2a558eb586101a317755d10435b6e2122a918d53245d5ab0","tier":null}],"members":[{"type":"individual","id":"1cc75fc7f25f3d6a","login":"atsokol","name":"Andrii T","avatar_url":"https://avatars-next.observableusercontent.com/e5b1ecc1785c25cd2a558eb586101a317755d10435b6e2122a918d53245d5ab0","tier":null,"role":"editor"},{"type":"team","id":"67596cf83e89a91c","login":"atsokol","name":"atsokol","avatar_url":"https://avatars-next.observableusercontent.com/377abaada84daeb487e417deeb4779c04a819d8e607bfcf58c388483815ab8e4","tier":"starter_2024","role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"title":"Price capture in a coupled solar / wind + BESS operation in Ukraine","cells":[{"id":0,"value":"# Price capture in a coupled solar / wind + BESS operation in Ukraine","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":154,"value":"Plot.plot({\n title: \"Wind price capture\",\n subtitle: \"Premium / discount vs average DAM price\",\n marginLeft: 50,\n marginRight: 30,\n width: 800,\n x: {label: null},\n y: { nice: true, domain: [0.7, 1.25], grid: true, label:null, tickFormat: d3.format(\".0%\")},\n marks: [\n Plot.lineY(\n wind_capture, \n { \n x: \"date_month\",\n y: \"capture\",\n stroke: \"steelblue\",\n tip: true,\n curve: \"catmull-rom\"\n }\n ),\n Plot.ruleY([1])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":893,"value":"Plot.plot({\n title: \"Solar price capture\",\n subtitle: \"Premium / discount vs average DAM price\",\n marginLeft: 50,\n marginRight: 30,\n width: 800,\n x: {label: null},\n y: { nice: true, domain: [0.7, 1.25], grid: true, label:null, tickFormat: d3.format(\".0%\")},\n marks: [\n Plot.lineY(\n solar_capture, \n { \n x: \"date_month\",\n y: \"capture\",\n stroke: \"orange\",\n tip: true,\n curve: \"catmull-rom\"\n }\n ),\n Plot.ruleY([1])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":78,"value":"Plot.plot({\n title: \"Wind yield by month and hour of day\",\n subtitle: \"Median (P50) capacity factor by month and time of day\",\n marginLeft: 50,\n marginRight: 30,\n width: 800,\n height: 300,\n x: { line: true, label: \"hour of day\"},\n y: { nice: true, label: null, grid: true, tickFormat: d3.format(\".0%\")},\n color: {interpolate: d3.interpolateRgbBasisClosed([\"#067FF0\", \"#FA7107\"]), \n tickFormat: d => short_months.get(d), \n label: null,\n legend: true},\n marks: [\n Plot.lineY(\n windCapFactor, \n { \n x: \"hour\",\n y: \"capFactor_p50\",\n stroke: \"month\", \n curve: \"catmull-rom\",\n }\n )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":311,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":779,"value":"Plot.plot({\n title: \"Solar yield by month and hour of day\",\n subtitle: \"Median (P50) capacity factor by month and time of day\",\n marginLeft: 50,\n marginRight: 30,\n width: 800,\n height: 300,\n x: { line: true, label: \"hour of day\"},\n y: { nice: true, label: null, grid: true, tickFormat: d3.format(\".0%\")},\n color: {interpolate: d3.interpolateRgbBasisClosed([\"#067FF0\", \"#FA7107\"]), \n tickFormat: d => short_months.get(d), \n label: null,\n legend: true},\n marks: [\n Plot.lineY(\n solarCapFactor, \n { \n x: \"hour\",\n y: \"capFactor_p50\",\n stroke: \"month\", \n curve: \"catmull-rom\",\n }\n )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":789,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":840,"value":"## Battery coupled RES","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":905,"value":"viewof startDate = Inputs.date({label: \"Start date\", value: \"2023-01-01\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":654,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":750,"value":"Plot.plot({\n title: \"Wind price capture with and without battery storage\",\n marginRight: 30,\n width: 1000,\n x: {label: null, ticks: \"month\"},\n y: { nice: true, grid: true, label:null, tickFormat: d3.format(\".0%\")},\n color: {legend: true},\n marks: [\n Plot.lineY(\n wind_battery_capture,\n {\n x: \"date_month\",\n y: \"value\",\n stroke: \"key\",\n strokeWidth: 2,\n tip: true,\n curve: \"catmull-rom\"\n }),\n Plot.ruleY([1]),\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":899,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":944,"value":"Plot.plot({\n title: \"Solar price capture with and without battery storage\",\n marginRight: 30,\n width: 1000,\n x: {label: null, ticks: \"month\"},\n y: { nice: true, grid: true, label:null, tickFormat: d3.format(\".0%\")},\n color: {legend: true},\n marks: [\n Plot.lineY(\n solar_battery_capture,\n {\n x: \"date_month\",\n y: \"value\",\n stroke: \"key\",\n strokeWidth: 2,\n tip: true,\n curve: \"catmull-rom\"\n }),\n Plot.ruleY([1]),\n ]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":332,"value":"months = new Map([\n [1, \"January\"],\n [2, \"February\"],\n [3, \"March\"],\n [4, \"April\"],\n [5, \"May\"],\n [6, \"June\"],\n [7, \"July\"],\n [8, \"August\"],\n [9, \"September\"],\n [10, \"October\"],\n [11, \"November\"],\n [12, \"December\"]\n]);","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":586,"value":"short_months = new Map([\n [1, \"Jan\"],\n [2, \"Feb\"],\n [3, \"Mar\"],\n [4, \"Apr\"],\n [5, \"May\"],\n [6, \"Jun\"],\n [7, \"Jul\"],\n [8, \"Aug\"],\n [9, \"Sep\"],\n [10, \"Oct\"],\n [11, \"Nov\"],\n [12, \"Dec\"]\n]);","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":635,"value":"#### Battery dispatch ","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":732,"value":"wind_battery_capture = {\n const yield_with_battery = wind_dispatch\n .map(d => ({...d, \n yield_wind: d.yield / 400,\n yield_battery: (d.sale_direct + d.discharge) / 400,\n date: d3.utcDay(d.time),\n date_month: d3.utcMonth(d.time)\n }))\n\n\n return aq.from(price_2023)\n .join_left(aq.from(yield_with_battery))\n .derive({\n price_w: d => d.price * d.yield_wind,\n price_wb: d => d.price * d.yield_battery\n })\n .groupby(\"date_month\")\n .rollup({price: d => op.mean(d.price),\n price_w: d => op.sum(d.price_w) / op.sum(d.yield_wind),\n price_wb: d => op.sum(d.price_wb) / op.sum(d.yield_battery)\n })\n .derive({\n capture_wind: d => d.price_w / d.price,\n capture_wind_battery: d => d.price_wb / d.price\n })\n .fold([\"capture_wind\", \"capture_wind_battery\"])\n .objects()\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":946,"value":"solar_battery_capture = {\n const yield_with_battery = solar_dispatch\n .map(d => ({...d, \n yield_solar: d.yield / 400,\n yield_battery: (d.sale_direct + d.discharge) / 400,\n date: d3.utcDay(d.time),\n date_month: d3.utcMonth(d.time)\n }))\n\n\n return aq.from(price_2023)\n .join_left(aq.from(yield_with_battery))\n .derive({\n price_s: d => d.price * d.yield_solar,\n price_sb: d => d.price * d.yield_battery\n })\n .groupby(\"date_month\")\n .rollup({price: d => op.mean(d.price),\n price_s: d => op.sum(d.price_s) / op.sum(d.yield_solar),\n price_sb: d => op.sum(d.price_sb) / op.sum(d.yield_battery)\n })\n .derive({\n capture_solar: d => d.price_s / d.price,\n capture_solar_battery: d => d.price_sb / d.price\n })\n .fold([\"capture_solar\", \"capture_solar_battery\"])\n .objects()\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":638,"value":"wind_dispatch = FileAttachment(\"wind_dispatch_2023.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":896,"value":"solar_dispatch = FileAttachment(\"solar_dispatch_2023.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":486,"value":"#### Wind and solar price capture","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":621,"value":"price_2023 = price_hourly\n .filter(d => d.date.getFullYear() === 2023)\n .map(d => ({time: d.datetime, price: d.price}))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":92,"value":"wind_capture = aq.from(price_hourly.map(d => ({...d, date_month: d3.utcMonth(d.date)})))\n .join(aq.from(wind_yield), ['datetime', 'local_time'])\n .derive({price_w: d => d.price * d.electricity})\n .groupby(\"date_month\")\n .rollup({price: d => op.mean(d.price),\n price_w: d => op.sum(d.price_w) / op.sum(d.electricity)\n })\n .derive({capture: d => d.price_w / d.price })\n .objects()\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":844,"value":"solar_capture = aq.from(price_hourly.map(d => ({...d, date_month: d3.utcMonth(d.date)})))\n .join(aq.from(solar_yield), ['datetime', 'time'])\n .derive({price_w: d => d.price * d.yield})\n .groupby(\"date_month\")\n .rollup({price: d => op.mean(d.price),\n price_w: d => op.sum(d.price_w) / op.sum(d.yield)\n })\n .derive({capture: d => d.price_w / d.price })\n .objects()\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":491,"value":"price_hourly = {\n \n const system = [\"IPS of Ukraine\", \"IPS of Ukraine (synchronized with ENTSO-E systems)\"]\n \n let prices = DAM_price_raw\n .filter(d => system.includes(d.energy_system))\n .map( d => ({...d, datetime: new Date(new Date(d.date).setUTCHours(d.hour-1))}))\n \n let fx_rates = exchange_rates\n .map(d => ({date: d3.utcParse(\"%d.%m.%Y\")(d.exchangedate), rate: d.rate_per_unit}))\n\n const colsToTransform = ['price', 'price_cap'];\n \n return aq.from(prices)\n .join_left(aq.from(fx_rates))\n .objects()\n .map(row => {\n const newRow = { ...row };\n colsToTransform.forEach(col => {newRow[col] = newRow[col] / newRow[\"rate\"]})\n return newRow;\n })\n\n}\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":461,"value":"#### Exchange rate data","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":479,"value":"exchange_rates = fetchp(nbuUrl)\n .then(response => response.json())","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":469,"value":"nbuUrl = {\n const today = new Date()\n const startDate = new Date(\"2019-07-01\")\n\n return \"https://bank.gov.ua/NBU_Exchange/exchange_site?start=\" + \n startDate.toISOString().split('T')[0].replace(/-/g, '') +\n \"&end=\" +\n today.toISOString().split('T')[0].replace(/-/g, '') +\n \"&valcode=\" +\n \"EUR\" +\n \"&sort=exchangedate&order=desc&json\"\n}\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":67,"value":"#### Electricity price data","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":58,"value":"DAM_price_raw = await fetchp(\n `${baseUrl + uiid_dam}/download/?${params}`,\n {\n method: 'GET',\n headers: {'Authorization': `Token e39edd2111e844899c556997772b8bc7`},\n mode : \"cors\"\n }\n)\n.then(response => response.text())\n.then(text => d3.csvParse(text, d3.autoType))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":456,"value":"params = new URLSearchParams({\n format: \"csv\", \n language: \"en\" \n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":453,"value":"uiid_dam = \"5a616fba-fbc9-4073-9532-9161592faca8\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":450,"value":"baseUrl = \"https://map.ua-energy.org/api/v1/resources/\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":63,"value":"#### Wind and solar capacity factors","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":69,"value":"capacity_wind = 1","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":774,"value":"capacity_solar = 1","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":771,"value":"solarCapFactor = {\n let solar_yield_hr_month = solar_yield\n .map(d => ({...d, \n hour: d.time.getHours(), \n month: d.time.getMonth() + 1\n }))\n\n return aq.from(solar_yield_hr_month)\n .params({capacity_solar: capacity_solar})\n .groupby(\"hour\", \"month\")\n .rollup({\n capFactor_p50: d => op.median(d.yield) / capacity_solar,\n capFactor_p10: d => op.quantile(d.yield, 0.1) / capacity_solar,\n capFactor_p90: d => op.quantile(d.yield, 0.9) / capacity_solar,\n })\n .orderby(\"month\", \"hour\")\n .objects()\n \n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":303,"value":"windCapFactor = {\n let wind_yield_hr_month = wind_yield\n .map(d => ({...d, \n hour: d.local_time.getHours(), \n month: d.local_time.getMonth() + 1\n }))\n\n return aq.from(wind_yield_hr_month)\n .params({capacity_wind: capacity_wind})\n .groupby(\"hour\", \"month\")\n .rollup({\n capFactor_p50: d => op.median(d.electricity) / capacity_wind,\n capFactor_p10: d => op.quantile(d.electricity, 0.1) / capacity_wind,\n capFactor_p90: d => op.quantile(d.electricity, 0.9) / capacity_wind,\n })\n .orderby(\"month\", \"hour\")\n .objects()\n \n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":807,"value":"#### Wind and solar yield data","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":19,"value":"wind_yield = FileAttachment(\"Syvash wind yield 2014-2023@1.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":767,"value":"solar_yield = FileAttachment(\"Prohresivka solar yield 2005-2023.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":991,"value":"

Wind yield data taken from here: https://www.re-explorer.org/home

Solar yield data taken from here: https://re.jrc.ec.europa.eu/pvg_tools/en/#HR

","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":46,"value":"import {aq, op} from \"@uwdata/arquero\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":440,"value":"import {fetchp} from '@tomlarkworthy/fetchp'","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 30:"$2f:body"