Many thanks, prof. Zimmerman. I will give it a try. Best regards. Get Outlook for Android<https://aka.ms/ghei36>
________________________________ From: [email protected] <[email protected]> on behalf of Ray Zimmerman <[email protected]> Sent: Tuesday, April 16, 2019 11:28:15 AM To: MATPOWER discussion forum Subject: Re: Problems creating a new case in Matpower - 14 Generator Model I think there is some confusion about bus numbering. There are external bus numbers and there are row indices into the bus matrix. In the original case file, these need not be the same. We call the bus numbers the external bus numbers. They need not be consecutive, i.e. numbering from 101 to 509 is fine. However, these are not row indices into the bus matrix (which has only 59 rows). So assigning a value to bus(101, PD) simply creates a bunch of new blank rows in the bus matrix and assigns a value to the PD column of row 101. But now you have an invalid bus matrix. Specifically, lots of rows with the BUS_I column equal to zero. Internally, MATPOWER calls ext2int() to convert from external bus numbers to consecutive integers (1–59 in your case). After this conversion, bus numbers and row indices into the bus matrix are identical. So the bottom line is, when changing the load, you need to find the row index that corresponds to the bus number you are trying to change and use that to index the bus matrix. I often do the following to create and external to internal bus mapping. nb = size(mpc.bus, 1); max_bus_num = max(mpc.bus(:, BUS_I)); e2i = sparse(mpc.bus(:, BUS_I), 1, 1:nb, max_bus_num, 1); % ext to int bus number map Then you could use something like the following to update your loads … mpc.bus(e2i(bus_nbr), PD) = MW Hope this helps, Ray On Apr 10, 2019, at 7:49 AM, Tayeb Meridji <[email protected]<mailto:[email protected]>> wrote: Hi Ray, One last note. So below i took a snapshot of the bug. - when i run the OPF on the original base case, it works fine; - when I attempt to change the load at bus, it allows me to do that (at least i think so since no error is shown). But when right after i attempt to run the OPF again, it gives me the following error message (see below). Thanks. <pastedImage.png> ________________________________ From: Tayeb Meridji Sent: Wednesday, April 10, 2019 8:43 AM To: MATPOWER discussion forum Subject: Re: Problems creating a new case in Matpower - 14 Generator Model 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: [email protected]<mailto:[email protected]> <[email protected]<mailto:[email protected]>> on behalf of Ray Zimmerman <[email protected]<mailto:[email protected]>> 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 <[email protected]<mailto:[email protected]>> 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: [email protected]<mailto:[email protected]> <[email protected]<mailto:[email protected]>> on behalf of Ray Zimmerman <[email protected]<mailto:[email protected]>> 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 <[email protected]<mailto:[email protected]>> 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; ];
