Hi Ray,
One other note. When I simply run an OPF on my networks fine, and the buses correspond to what I have defined. So it really seems to be a problem with the load change function, which for some reason is not working. Thanks. ________________________________ From: Tayeb Meridji Sent: Wednesday, April 10, 2019 8:36 AM To: MATPOWER discussion forum Subject: Re: Problems creating a new case in Matpower - 14 Generator Model Hi Ray, No, bus numbers in the base case i have created go from 101 to 509. And that is the strange thing because the error message mentions some bus numbers that are not even in my base case. I don't know if you can open files sent by users, but in case you can I have attached my base case. Thank you. ________________________________ From: bounce-123506270-83594...@list.cornell.edu <bounce-123506270-83594...@list.cornell.edu> on behalf of Ray Zimmerman <r...@cornell.edu> Sent: Tuesday, April 9, 2019 5:32 PM To: MATPOWER discussion forum Subject: Re: Problems creating a new case in Matpower - 14 Generator Model And can you confirm that bus_nbr is always a number between 1 and 59? Ray On Apr 9, 2019, at 3:36 PM, Tayeb Meridji <tayeb.meri...@mail.mcgill.ca<mailto:tayeb.meri...@mail.mcgill.ca>> wrote: Hi, Thanks for the reply. To change the load i used: mpc.bus(bus_nbr, PD) = MW. Thanks for the tips. Regards. Get Outlook for Android<https://aka.ms/ghei36> ________________________________ From: bounce-123504563-83594...@list.cornell.edu<mailto:bounce-123504563-83594...@list.cornell.edu> <bounce-123504563-83594...@list.cornell.edu<mailto:bounce-123504563-83594...@list.cornell.edu>> on behalf of Ray Zimmerman <r...@cornell.edu<mailto:r...@cornell.edu>> Sent: Tuesday, April 9, 2019 11:27:35 AM To: MATPOWER discussion forum Subject: Re: Problems creating a new case in Matpower - 14 Generator Model You haven’t shown how you “attempt to change the load at a given bus” but I suspect there is a problem there. Your original case has 59 rows in the bus matrix. Make sure that mpc.bus still has only 59 rows after updating the load. Two other tips … * If you aren’t already, you might consider using the scale_load()<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/scale_load.html> or apply_changes()<http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/apply_changes.html> functions to update your loads. * I notice that your line ratings are all set to 9999. If you actually mean “unconstrained” you should simply use 0, rather than a large number. This will eliminate line constraints completely and speed up the optimization. Ray On Apr 2, 2019, at 10:35 AM, Tayeb Meridji <tayeb.meri...@mail.mcgill.ca<mailto:tayeb.meri...@mail.mcgill.ca>> wrote: Hi, I am trying to create a new MAtpower case for the 14 generator system representing an equivalent of the Australian network. When I run an OPF of the created network, it seems to work fine. However, when I attempt to change the load at a given bus, I get an error message saying the following (see below). The problem is that the buses mentioned in the error message are not even part of the newly created network. So I am really puzzled. Please note that at the end of this message I have included the input of my new network. Your guidance would be greatly appreciated. Many thanks. Tayeb. Error messages: error: ext2int: bus 60 has an invalid BUS_TYPEext2int: bus 61 has an invalid BUS_TYPEext2int: bus 62 has an invalid BUS_TYPEext2int: bus 63 has an invalid BUS_TYPEext2int: bus 64 has an invalid BUS_TYPEext2int: bus 65 has an invalid BUS_TYPEext2int: bus 66 has an invalid BUS_TYPEext2int: bus 67 has an invalid BUS_TYPEext2int: bus 68 has an invalid BUS_TYPEext2int: bus 69 has an invalid BUS_TYPEext2int: bus 70 has an invalid BUS_TYPEext2int: bus 71 has an invali d BUS_TYPEext2int: bus 72 has an invalid BUS_TYPEext2int: bus 73 has an invalid BUS_TYPEext2int: bus 74 has an invalid BUS_TYPEext2int: bus 75 has an inval id BUS_TYPEext2int: bus 76 has an invalid BUS_TYPEext2int: bus 77 has an invalid BUS_TYPEext2int: bus 78 has an invalid BUS_TYPEext2int: bus 79 has an inva lid BUS_TYPEext2int: bus 80 has an invalid BUS_TYPEext2int: bus 81 has an invalid BUS_TYPEext2int: bus 82 has an invalid BUS_TYPEext2int: bus 83 has an inv alid BUS_TYPEext2int: bus 84 has an invalid BUS_TYPEext2int: bus 85 has an invalid BUS_TYPEext2int: bus 86 has an invalid BUS_TYPEext2int: bus 87 has an in valid BUS_TYPEext2int: bus 88 has an invalid BUS_TYPEext2int: bus 89 has an invalid BUS_TYPEext2int: bus 90 has an invalid BUS_TYPEext2int: bus 91 has an i nvalid BUS_TYPEext2int: bus 92 has an invalid BUS_TYPEext2int: bus 93 has an invalid BUS_TYPEext2int: bus 94 has an invalid BUS_TYPEext2int: bus 95 has an invalid BUS_TYPEext2int: bus 96 has an invalid BUS_TYPEext2int: bus 97 has an invalid BUS_TYPEext2int: bus 98 has an invalid BUS_TYPEext2int: bus 99 has an invalid BUS_TYPEext2int: bus 100 has an invalid BUS_TYPEext2int: bus 101 has an invalid BUS_TYPEext2int: bus 102 has an invalid BUS_TYPEext2int: bus 103 h as an invalid BUS_TYPEext2int: bus 104 has an invalid BUS_TYPEext2int: bus 105 has an invalid BUS_TYPEext2int: bus 106 has an invalid BUS_TYPEext2int: bus 107 has an invalid BUS_TYPEext2int: bus 108 has an invalid BUS_TYPEext2int: bus 109 has an invalid BUS_TYPEext2int: bus 110 has an invalid BUS_TYPEext2int: bus 111 has an invalid BUS_TYPEext2int: bus 112 has an invalid BUS_TYPEext2int: bus 113 has an invalid BUS_TYPEext2int: bus 114 has an invalid BUS_TYPEext 2int: bus 115 has an invalid BUS_TYPEext2int: bus 116 has an invalid BUS_TYPEext2int: bus 117 has an invalid BUS_TYPEext2int: bus 118 has an invalid BUS_TY PEext2int: bus 119 has an invalid BUS_TYPEext2int: bus 120 has an invalid BUS_TYPEext2int: bus 121 has an invalid BUS_TYPEext2int: bus 122 has an invalid B US_TYPEext2int: bus 123 has an invalid BUS_TYPEext2int: bus 124 has an invalid BUS_TYPEext2int: bus 125 has an invalid BUS_TYPEext2int: bus 126 has an inva lid BUS_TYPEext2int: bus 127 has an invalid BUS_TYPEext2int: bus 128 has an invalid BUS_TYPEext2int: bus 129 has an invalid BUS_TYPEext2int: bus 130 has an invalid BUS_TYPEext2int: bus 131 has an invalid BUS_TYPEext2int: bus 132 has an invalid BUS_TYPEext2int: bus 133 has an invalid BUS_TYPEext2int: bus 134 h as an invalid BUS_TYPEext2int: bus 135 has an invalid BUS_TYPEext2int: bus 136 has an invalid BUS_TYPEext2int: bus 137 has an invalid BUS_TYPEext2int: bus 138 has an invalid BUS_TYPEext2int: bus 139 has an invalid BUS_TYPEext2int: bus 140 has an invalid BUS_TYPEext2int: bus 141 has an invalid BUS_TYPEext2int: bus 142 has an invalid BUS_TYPEext2int: bus 143 has an invalid BUS_TYPEext2int: bus 144 has an invalid BUS_TYPEext2int: bus 145 has an invalid BUS_TYPEext 2int: bus 146 has an invalid BUS_TYPEext2int: bus 147 has an invalid BUS_TYPEext2int: bus 148 has an invalid BUS_TYPEext2int: bus 149 has an invalid BUS_TY PEext2int: bus 150 has an invalid BUS_TYPEext2int: bus 151 has an invalid BUS_TYPEext2int: bus 152 has an invalid BUS_TYPEext2int: bus 153 has an invalid B US_TYPEext2int: bus 154 has an invalid BUS_TYPEext2int: bus 155 has an invalid BUS_TYPEext2int: bus 156 has an invalid BUS_TYPEext2int: bus 157 has an inva lid BUS_TYPEext2int: bus 158 has an invalid BUS_TYPEext2int: bus 159 has an invalid BUS_TYPEext2int: bus 160 has an invalid BUS_TYPEext2int: bus 161 has an invalid BUS_TYPEext2int: bus 162 has an invalid BUS_TYPEext2int: bus 163 has an invalid BUS_TYPEext2int: bus 164 has an invalid BUS_TYPEext2int: bus 165 h as an invalid BUS_TYPEext2int: bus 166 has an invalid BUS_TYPEext2int: bus 167 has an invalid BUS_TYPEext2int: bus 168 has an invalid BUS_TYPEext2int: bus 169 has an invalid BUS_TYPEext2int: bus 170 has an invalid BUS_TYPEext2int: bus 171 has an invalid BUS_TYPEext2int: bus 172 has an invalid BUS_TYPEext2int: bus 173 has an invalid BUS_TYPEext2int: bus 174 has an invalid BUS_TYPEext2int: bus 175 has an invalid BUS_TYPEext2int: bus 176 has an invalid BUS_TYPEext 2int: bus 177 has an invalid BUS_TYPEext2int: bus 178 has an invalid BUS_TYPEext2int: bus 179 has an invalid BUS_TYPEext2int: bus 180 has an invalid BUS_TY PEext2int: bus 181 has an invalid BUS_TYPEext2int: bus 182 has an invalid BUS_TYPEext2int: bus 183 has an invalid BUS_TYPEext2int: bus 184 has an invalid B US_TYPEext2int: bus 185 has an invalid BUS_TYPEext2int: bus 186 has an invalid BUS_TYPEext2int: bus 187 has an invalid BUS_TYPEext2int: bus 188 has an inva lid BUS_TYPEext2int: bus 189 has an invalid BUS_TYPEext2int: bus 190 has an invalid BUS_TYPEext2int: bus 191 has an invalid BUS_TYPEext2int: bus 192 has an invalid BUS_TYPEext2int: bus 193 has an invalid BUS_TYPEext2int: bus 194 has an invalid BUS_TYPEext2int: bus 195 has an invalid BUS_TYPEext2int: bus 196 h as an invalid BUS_TYPEext2int: bus 197 has an invalid BUS_TYPEext2int: bus 198 has an invalid BUS_TYPEext2int: bus 199 has an invalid BUS_TYPEext2int: bus 200 has an invalid BUS_TYPEext2int: bus 201 has an invalid BUS_TYPEext2int: bus 202 has an invalid BUS_TYPEext2int: bus 203 has an invalid BUS_TYPEext2int: bus 204 has an invalid BUS_TYPEext2int: bus 205 has an invalid BUS_TYPEext2int: bus 206 has an invalid BUS_TYPEext2int: bus 207 has an invalid BUS_TYPEext 2int: bus 208 has an invalid BUS_TYPE error: called from ext2int at line 135 column 17 14 Gen Matpower Model function mpc = case14gen %CASE39 Power flow data for 14 Gen Australian equivalent system. % Please see CASEFORMAT for details on the case file format. % MATPOWER %% MATPOWER Case Format : Version 2 mpc.version = '2'; %%----- Power Flow Data -----%% %% system MVA base mpc.baseMVA = 100; %% bus data % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin mpc.bus = [ 101 3 0 0 0 0 1 1 0 15 1 1.06 0.94; 102 1 450 45 0 0 1 1.039 -1.4 330 1 1.06 0.94; 201 2 0 0 0 0 2 1 48.9 20 1 1.06 0.94; 202 2 0 0 0 0 2 1 37.9 20 1 1.06 0.94; 203 2 0 0 0 0 2 1 32.24 20 1 1.06 0.94; 204 2 0 0 0 0 2 1 39.21 20 1 1.06 0.94; 205 2 390 39 0 0 2 1.055 44.02 330 1 1.06 0.94; 206 1 130 13 0 0 2 1.0472 41.45 330 1 1.06 0.94; 207 1 1880 188 0 0 2 1.0221 27.69 330 1 1.06 0.94; 208 1 210 21 0 0 2 1.0324 26.25 330 1 1.06 0.94; 209 1 0 0 0 0 2 1.0379 30.41 330 1 1.06 0.94; 210 1 0 0 0 0 2 1.0594 26.68 500 1 1.06 0.94; 211 1 1700 170 0 0 2 1.0081 18.87 330 1 1.06 0.94; 212 1 1660 166 0 0 2 1.0084 19.13 330 1 1.06 0.94; 213 1 0 0 0 0 2 1.0426 22.67 500 1 1.06 0.94; 214 1 0 0 0 0 2 1.0257 18.75 330 1 1.06 0.94; 215 1 480 48 0 0 2 1.0433 33.11 330 1 1.06 0.94; 216 1 1840 184 0 0 2 1.013 16.1 330 1 1.06 0.94; 217 1 1260 126 0 0 2 1.0021 9.39 330 1 1.06 0.94; 301 2 0 0 0 0 3 1 -3.39 20 1 1.06 0.94; 302 2 0 0 0 0 3 1 -18.02 20 1 1.06 0.94; 303 1 0 0 0 0 3 1.0417 -10.85 500 1 1.06 0.94; 304 1 0 0 0 0 3 1.0093 -21.75 500 1 1.06 0.94; 305 1 0 0 0 0 3 1.0164 -22.76 500 1 1.06 0.94; 306 1 1230 123 0 0 3 1.016 -24.92 500 1 1.06 0.94; 307 1 650 65 0 0 3 1.0182 -24.97 500 1 1.06 0.94; 308 1 655 66 0 0 3 1.0359 -35.41 500 1 1.06 0.94; 309 1 195 20 0 0 3 1.0291 -9.34 330 1 1.06 0.94; 310 1 0 0 0 0 3 1.0198 -24.81 330 1 1.06 0.94; 311 1 0 0 0 0 3 1.0136 -17.66 330 1 1.06 0.94; 312 1 115 12 0 0 3 1.0383 -23.59 220 1 1.06 0.94; 313 2 2405 240 0 0 3 1.015 -30.23 220 1 1.06 0.94; 314 1 250 25 0 0 3 1.0158 -28.52 220 1 1.06 0.94; 315 1 0 0 0 0 3 1.0376 -39.01 275 1 1.06 0.94; 401 2 0 0 0 0 4 1 74.34 20 1 1.06 0.94; 402 2 0 0 0 0 4 1 107.94 20 1 1.06 0.94; 403 2 0 0 0 0 4 1 113.34 20 1 1.06 0.94; 404 2 0 0 0 0 4 1 106.53 20 1 1.06 0.94; 405 1 989.9999 99 0 0 4 1.0312 99.54 275 1 1.06 0.94; 406 1 740 74 0 0 4 1.0342 104.79 275 1 1.06 0.94; 407 1 0 0 0 0 4 1.0393 107.11 275 1 1.06 0.94; 408 1 150 15 0 0 4 1.0301 100.38 275 1 1.06 0.94; 409 1 260 26 0 0 4 0.9855 76.61 275 1 1.06 0.94; 410 1 530 53 0 0 4 1.0301 68.14 275 1 1.06 0.94; 411 1 575 58 0 0 4 0.993 62.25 275 1 1.06 0.94; 412 2 1255 126 0 0 4 1 60.97 275 1 1.06 0.94; 413 1 0 0 0 0 4 1.0398 60.95 275 1 1.06 0.94; 414 1 0 0 0 0 4 1.0418 60.2 330 1 1.06 0.94; 415 1 0 0 0 0 4 1.0474 56.66 330 1 1.06 0.94; 416 1 0 0 0 0 4 1.0552 50.27 330 1 1.06 0.94; 501 2 0 0 0 0 5 1 -55.32 20 1 1.06 0.94; 502 2 0 0 0 0 5 1 -53.95 15 1 1.06 0.94; 503 2 0 0 0 0 5 1 -56.28 15 1 1.06 0.94; 504 1 300 60 0 0 5 1.0478 -63.31 275 1 1.06 0.94; 505 1 0 0 0 0 5 1.0221 -60.87 275 1 1.06 0.94; 506 1 0 0 0 0 5 1.0204 -62.18 275 1 1.06 0.94; 507 2 1000 200 0 0 5 1.015 -63.45 275 1 1.06 0.94; 508 1 800 160 0 0 5 1.0104 -64.66 275 1 1.06 0.94; 509 2 200 40 0 0 5 1.03 -45.74 275 1 1.06 0.94; ]; %% generator data % bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10 ramp_30 ramp_q apf mpc.gen = [ 101 300.8091 311.4755 581.128 -581.128 1 100 1 1200 0 0 0 0 0 0 0 0 0 0 0 0; 201 3600 573.6404 1743.648 -1743.648 1 100 1 3600 0 0 0 0 0 0 0 0 0 0 0 0; 202 2500 663.3815 1210.9 -1210.9 1 100 1 2500 0 0 0 0 0 0 0 0 0 0 0 0; 203 1500 531.1833 968.72 -968.72 1 100 1 2000 0 0 0 0 0 0 0 0 0 0 0 0; 204 2950.2 734.1567 1743.648 -1743.648 1 100 1 3600 0 0 0 0 0 0 0 0 0 0 0 0; 301 4200 996.5973 2034.256 -2034.256 1 100 1 4200 0 0 0 0 0 0 0 0 0 0 0 0; 302 939.9 154.5849 581.127 -581.127 1 100 1 1200 0 0 0 0 0 0 0 0 0 0 0 0; 401 1400 514.8802 774.836 -774.836 1 100 1 1600 0 0 0 0 0 0 0 0 0 0 0 0; 402 837 177.8523 435.846 -435.846 1 100 1 900 0 0 0 0 0 0 0 0 0 0 0 0; 403 1400 209.1907 774.836 -774.836 1 100 1 1600 0 0 0 0 0 0 0 0 0 0 0 0; 404 1549.8 326.8041 871.692 -871.692 1 100 1 1800 0 0 0 0 0 0 0 0 0 0 0 0; 501 600 50.694 290.564 -290.564 1 100 1 600 0 0 0 0 0 0 0 0 0 0 0 0; 502 800 160.2667 600 -600 1 100 1 800 0 0 0 0 0 0 0 0 0 0 0 0; 503 436 100.8554 290.652 -290.652 1 100 1 600 0 0 0 0 0 0 0 0 0 0 0 0; 205 0 -68.2624 430 -220 1.055 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; 313 0 71.5048 600 -200 1.015 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; 412 0 58.1613 1100 -330 1 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; 507 0 22.6476 320 -180 1.015 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; 509 0 10.569 400 -150 1.03 100 1 0 0 0 0 0 0 0 0 0 0 0 0 0; ]; %% branch data % fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax mpc.branch = [ 102 217 0.0084 0.0667 0.817 9999 9999 9999 0 0 1 -360 360; 102 217 0.0084 0.0667 0.817 9999 9999 9999 0 0 1 -360 360; 102 217 0.0078 0.062 0.76 9999 9999 9999 0 0 1 -360 360; 102 217 0.0078 0.062 0.76 9999 9999 9999 0 0 1 -360 360; 102 309 0.0045 0.0356 0.437 9999 9999 9999 0 0 1 -360 360; 102 309 0.0045 0.0356 0.437 9999 9999 9999 0 0 1 -360 360; 102 309 0.0109 0.0868 0.76 9999 9999 9999 0 0 1 -360 360; 205 206 0.0096 0.076 0.931 9999 9999 9999 0 0 1 -360 360; 205 206 0.0096 0.076 0.931 9999 9999 9999 0 0 1 -360 360; 205 416 0.0037 0.046 0.73 9999 9999 9999 0 0 1 -360 360; 205 416 0.0037 0.046 0.73 9999 9999 9999 0 0 1 -360 360; 206 207 0.0045 0.0356 0.437 9999 9999 9999 0 0 1 -360 360; 206 207 0.0045 0.0356 0.437 9999 9999 9999 0 0 1 -360 360; 206 212 0.0066 0.0527 0.646 9999 9999 9999 0 0 1 -360 360; 206 212 0.0066 0.0527 0.646 9999 9999 9999 0 0 1 -360 360; 206 215 0.0066 0.0527 0.646 9999 9999 9999 0 0 1 -360 360; 206 215 0.0066 0.0527 0.646 9999 9999 9999 0 0 1 -360 360; 207 208 0.0018 0.014 0.171 9999 9999 9999 0 0 1 -360 360; 207 208 0.0018 0.014 0.171 9999 9999 9999 0 0 1 -360 360; 207 209 0.0008 0.0062 0.076 9999 9999 9999 0 0 1 -360 360; 208 211 0.0031 0.0248 0.304 9999 9999 9999 0 0 1 -360 360; 208 211 0.0031 0.0248 0.304 9999 9999 9999 0 0 1 -360 360; 208 211 0.0031 0.0248 0.304 9999 9999 9999 0 0 1 -360 360; 209 212 0.0045 0.0356 0.437 9999 9999 9999 0 0 1 -360 360; 210 213 0.001 0.0145 1.54 9999 9999 9999 0 0 1 -360 360; 210 213 0.001 0.0145 1.54 9999 9999 9999 0 0 1 -360 360; 211 212 0.0014 0.0108 0.133 9999 9999 9999 0 0 1 -360 360; 211 212 0.0014 0.0108 0.133 9999 9999 9999 0 0 1 -360 360; 211 214 0.0019 0.0155 0.19 9999 9999 9999 0 0 1 -360 360; 212 217 0.007 0.0558 0.684 9999 9999 9999 0 0 1 -360 360; 214 216 0.001 0.0077 0.095 9999 9999 9999 0 0 1 -360 360; 214 217 0.0049 0.0388 0.475 9999 9999 9999 0 0 1 -360 360; 215 216 0.0051 0.0403 0.494 9999 9999 9999 0 0 1 -360 360; 215 216 0.0051 0.0403 0.494 9999 9999 9999 0 0 1 -360 360; 215 217 0.0072 0.0574 0.703 9999 9999 9999 0 0 1 -360 360; 215 217 0.0072 0.0574 0.703 9999 9999 9999 0 0 1 -360 360; 216 217 0.0051 0.0403 0.494 9999 9999 9999 0 0 1 -360 360; 303 304 0.002 0.028 0.74 9999 9999 9999 0 0 1 -360 360; 303 304 0.002 0.028 0.74 9999 9999 9999 0 0 1 -360 360; 303 305 0.0011 0.016 1.7 9999 9999 9999 0 0 1 -360 360; 303 305 0.0011 0.016 1.7 9999 9999 9999 0 0 1 -360 360; 304 305 0.0003 0.004 0.424 9999 9999 9999 0 0 1 -360 360; 305 306 0.0002 0.003 0.32 9999 9999 9999 0 0 1 -360 360; 305 307 0.0003 0.0045 0.447 9999 9999 9999 0 0 1 -360 360; 305 307 0.0003 0.0045 0.447 9999 9999 9999 0 0 1 -360 360; 306 307 0.0001 0.0012 0.127 9999 9999 9999 0 0 1 -360 360; 307 308 0.0023 0.0325 3.445 9999 9999 9999 0 0 1 -360 360; 307 308 0.0023 0.0325 3.445 9999 9999 9999 0 0 1 -360 360; 309 310 0.0135 0.10695 0.58267 9999 9999 9999 0 0 1 -360 360; 309 310 0.0135 0.10695 0.58267 9999 9999 9999 0 0 1 -360 360; 309 310 0.0135 0.10695 0.58267 9999 9999 9999 0 0 1 -360 360; 310 311 0 -0.0337 0 9999 9999 9999 0 0 1 -360 360; 310 311 0 -0.0337 0 9999 9999 9999 0 0 1 -360 360; 312 313 0.006 0.045 0.3 9999 9999 9999 0 0 1 -360 360; 312 313 0.006 0.045 0.3 9999 9999 9999 0 0 1 -360 360; 312 313 0.006 0.045 0.3 9999 9999 9999 0 0 1 -360 360; 313 314 0.001 0.01 0.26 9999 9999 9999 0 0 1 -360 360; 313 314 0.001 0.01 0.26 9999 9999 9999 0 0 1 -360 360; 315 509 0.007 0.05 0.19 9999 9999 9999 0 0 1 -360 360; 315 509 0.007 0.05 0.19 9999 9999 9999 0 0 1 -360 360; 405 406 0.0039 0.0475 0.381 9999 9999 9999 0 0 1 -360 360; 405 406 0.0039 0.0475 0.381 9999 9999 9999 0 0 1 -360 360; 405 408 0.0054 0.05 0.189 9999 9999 9999 0 0 1 -360 360; 405 409 0.018 0.122 0.79 9999 9999 9999 0 0 1 -360 360; 405 409 0.018 0.122 0.79 9999 9999 9999 0 0 1 -360 360; 405 409 0.018 0.122 0.79 9999 9999 9999 0 0 1 -360 360; 406 407 0.0006 0.0076 0.062 9999 9999 9999 0 0 1 -360 360; 406 407 0.0006 0.0076 0.062 9999 9999 9999 0 0 1 -360 360; 407 408 0.0042 0.0513 0.412 9999 9999 9999 0 0 1 -360 360; 408 410 0.0165 0.192 0.67333 9999 9999 9999 0 0 1 -360 360; 408 410 0.0165 0.192 0.67333 9999 9999 9999 0 0 1 -360 360; 408 410 0.0165 0.192 0.67333 9999 9999 9999 0 0 1 -360 360; 409 411 0.0103 0.0709 0.46 9999 9999 9999 0 0 1 -360 360; 409 411 0.0103 0.0709 0.46 9999 9999 9999 0 0 1 -360 360; 410 411 0.0043 0.0532 0.427 9999 9999 9999 0 0 1 -360 360; 410 412 0.0043 0.0532 0.427 9999 9999 9999 0 0 1 -360 360; 410 412 0.0043 0.0532 0.427 9999 9999 9999 0 0 1 -360 360; 410 412 0.0043 0.0532 0.427 9999 9999 9999 0 0 1 -360 360; 410 412 0.0043 0.0532 0.427 9999 9999 9999 0 0 1 -360 360; 410 413 0.004 0.0494 0.4 9999 9999 9999 0 0 1 -360 360; 410 413 0.004 0.0494 0.4 9999 9999 9999 0 0 1 -360 360; 411 412 0.0012 0.0152 0.122 9999 9999 9999 0 0 1 -360 360; 411 412 0.0012 0.0152 0.122 9999 9999 9999 0 0 1 -360 360; 414 415 0.002 0.025 0.39 9999 9999 9999 0 0 1 -360 360; 414 415 0.002 0.025 0.39 9999 9999 9999 0 0 1 -360 360; 415 416 0.0037 0.046 0.73 9999 9999 9999 0 0 1 -360 360; 415 416 0.0037 0.046 0.73 9999 9999 9999 0 0 1 -360 360; 504 507 0.023 0.15 0.56 9999 9999 9999 0 0 1 -360 360; 504 507 0.023 0.15 0.56 9999 9999 9999 0 0 1 -360 360; 504 508 0.026 0.019 0.87 9999 9999 9999 0 0 1 -360 360; 504 508 0.026 0.019 0.87 9999 9999 9999 0 0 1 -360 360; 505 507 0.0016 0.017 0.03 9999 9999 9999 0 0 1 -360 360; 505 507 0.0016 0.017 0.03 9999 9999 9999 0 0 1 -360 360; 505 508 0.0025 0.028 0.17 9999 9999 9999 0 0 1 -360 360; 506 507 0.0016 0.017 0.03 9999 9999 9999 0 0 1 -360 360; 506 507 0.0016 0.017 0.03 9999 9999 9999 0 0 1 -360 360; 506 508 0.003 0.028 0.14 9999 9999 9999 0 0 1 -360 360; 507 508 0.002 0.019 0.09 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 507 509 0.09 0.66 0.3 9999 9999 9999 0 0 1 -360 360; 101 102 0 0.009 0 9999 9999 9999 0 0 1 -360 360; 201 206 0 0.004 0 9999 9999 9999 0 0 1 -360 360; 202 209 0 0.00576 0 9999 9999 9999 0 0 1 -360 360; 203 208 0 0.00765 0 9999 9999 9999 0 0 1 -360 360; 204 215 0 0.004 0 9999 9999 9999 0 0 1 -360 360; 209 210 0 0.00000503 0 9999 9999 9999 0 0 1 -360 360; 213 214 0 0.00000503 0 9999 9999 9999 0 0 1 -360 360; 301 303 0 0.00343 0 9999 9999 9999 0 0 1 -360 360; 302 312 0 0.01127 0 9999 9999 9999 0 0 1 -360 360; 304 313 0 0.016 0 9999 9999 9999 0 0 1 -360 360; 305 311 0 0.00048600 0 9999 9999 9999 0 0 1 -360 360; 305 314 0 0.012 0 9999 9999 9999 0 0 1 -360 360; 308 315 0 0.0135 0 9999 9999 9999 0 0 1 -360 360; 401 410 0 0.00845 0 9999 9999 9999 0 0 1 -360 360; 402 408 0 0.017 0 9999 9999 9999 0 0 1 -360 360; 403 407 0 0.00845 0 9999 9999 9999 0 0 1 -360 360; 404 405 0 0.0085 0 9999 9999 9999 0 0 1 -360 360; 413 414 0 0.00002133 0 9999 9999 9999 0 0 1 -360 360; 501 504 0 0.0255 0 9999 9999 9999 0 0 1 -360 360; 502 505 0 0.016 0 9999 9999 9999 0 0 1 -360 360; 503 506 0 0.025 0 9999 9999 9999 0 0 1 -360 360; ]; %%----- OPF Data -----%% %% generator cost data % 1 startup shutdown n x1 y1 ... xn yn % 2 startup shutdown n c(n-1) ... c0 mpc.gencost = [ 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; 2 0 0 3 0.01 0.3 0.2; ];