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]> 
> 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] 
> <[email protected]> on behalf of Ray Zimmerman 
> <[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;
>> 
>> ];
> 

Reply via email to