...
 
Commits (2)
......@@ -170,8 +170,8 @@ Consider an electricity market with two generator types, one with the cost funct
Consider the two-bus power system shown in Figure \ref{twobus}, where the two nodes represent two markets, each with different total demand $D_i$, and one generator at each node producing $P_i$. At node A the demand is $D_A = 2000 \si{\mega\watt}$, whereas at node B the demand is $D_B = 1000 \si{\mega\watt}$. Furthermore, there is a transmission line with a capacity denoted by $F_{AB}$. The marginal cost of production of the generators connected to buses A and B are given respectively by the following expressions:
\begin{align*}
MC_A & = 20 + 0.03 P_A \hspace{1cm}\eur/\si{\mega\watt\hour} \\
MC_B & = 15 + 0.02 P_B \hspace{1cm} \eur/\si{\mega\watt\hour}
MC_A & = 20 + 0.02 P_A \hspace{1cm}\eur/\si{\mega\watt\hour} \\
MC_B & = 15 + 0.03 P_B \hspace{1cm} \eur/\si{\mega\watt\hour}
\end{align*}
Assume that the demands $D_A$ and $D_B$ are constant and insensitive to price, that energy is sold at its marginal cost of production and that there are no limits on the output of the generators.
......
......@@ -85,12 +85,13 @@
"outputs": [],
"source": [
"#make the code as Python 3 compatible as possible\n",
"from __future__ import print_function, division, absolute_import\n",
"#from __future__ import print_function, division, absolute_import\n",
"import pypsa\n",
"import numpy as np\n",
"import pandas as pd\n",
"import os\n",
"import matplotlib.pyplot as plt\n",
"import cartopy.crs as ccrs\n",
"plt.rcParams['figure.dpi']= 400\n",
"plt.style.use('ggplot')\n",
"%matplotlib inline"
......@@ -338,14 +339,14 @@
" n_rows_g = n_graphs_g // n_cols_g\n",
"else:\n",
" n_rows_g = n_graphs_g // n_cols_g + 1\n",
" \n",
"fig, axes = plt.subplots(nrows=n_rows_g, ncols=n_cols_g,\n",
" figsize=(width_factor*n_cols_g,height_factor*n_rows_g),\n",
" subplot_kw={\"projection\":ccrs.PlateCarree()})\n",
"\n",
"width_factor = 4\n",
"height_factor = 4\n",
"\n",
"fig, axes = plt.subplots(nrows=n_rows_g, ncols=n_cols_g,\n",
" figsize=(width_factor*n_cols_g,height_factor*n_rows_g),\n",
" subplot_kw={\"projection\":ccrs.PlateCarree()})\n",
"\n",
"for i,tech in enumerate(techs):\n",
" i_row = i // n_cols_g\n",
" i_col = i % n_cols_g\n",
......@@ -450,7 +451,7 @@
"\n",
"fig, ax = plt.subplots(figsize=(16,8))\n",
"\n",
"p_by_carrier.plot(\n",
"p_by_carrier_gw.plot(\n",
" kind=\"area\",\n",
" ax=ax,\n",
" linewidth=0,\n",
......