Source Code for Module pypower.t.t_case_int

  1  # Copyright (C) 2004-2011 Power System Engineering Research Center (PSERC) 
  2  # Copyright (C) 2011 Richard Lincoln 
  3  # 
  4  # PYPOWER is free software: you can redistribute it and/or modify 
  5  # it under the terms of the GNU General Public License as published 
  6  # by the Free Software Foundation, either version 3 of the License, 
  7  # or (at your option) any later version. 
  8  # 
  9  # PYPOWER is distributed in the hope that it will be useful, 
 10  # but WITHOUT ANY WARRANTY], without even the implied warranty of 
 11  # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 
 12  # GNU General Public License for more details. 
 13  # 
 14  # You should have received a copy of the GNU General Public License 
 15  # along with PYPOWER. If not, see <http://www.gnu.org/licenses/>. 
 16   
 17  """Case data in internal format. 
 18  """ 
 19   
 20  from numpy import array 
 21   
 22   
 23 -def t_case_int(): 
 24      """Case data in internal format used to test C{ext2int} and C{int2ext}. 
 25      """ 
 26      ppc = {} 
 27   
 28      ## PYPOWER Case Format : Version 2 
 29      ppc['version'] = '2' 
 30   
 31      ##-----  Power Flow Data  -----## 
 32      ## system MVA base 
 33      ppc['baseMVA'] = 100.0 
 34   
 35      ## bus data 
 36      # bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 
 37      ppc['bus'] = array([ 
 38          [0, 3, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 39          [1, 2, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 40          [2, 2, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 41          [3, 1, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 42          [4, 1, 90,  30, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 43          [5, 1, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 44          [6, 1, 100, 35, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 45          [7, 1, 0,   0,  0, 0, 1, 1, 0, 345, 1, 1.1, 0.9], 
 46          [8, 1, 125, 50, 0, 0, 1, 1, 0, 345, 1, 1.1, 0.9] 
 47      ]) 
 48   
 49      ## generator data 
 50      # bus, Pg, Qg, Qmax, Qmin, Vg, mBase, status, Pmax, Pmin, Pc1, Pc2, 
 51      # Qc1min, Qc1max, Qc2min, Qc2max, ramp_agc, ramp_10, ramp_30, ramp_q, apf 
 52      ppc['gen'] = array([ 
 53          [0, 0,   0, 300, -300, 1, 100, 1, 250, 90, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 54          [1, 163, 0, 300, -300, 1, 100, 1, 300, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], 
 55          [2, 85,  0, 300, -300, 1, 100, 1, 270, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 
 56      ], float) 
 57   
 58      ## branch data 
 59      # fbus, tbus, r, x, b, rateA, rateB, rateC, ratio, angle, status, angmin, angmax 
 60      ppc['branch'] = array([ 
 61          [0, 3, 0,      0.0576, 0,     0,   250, 250, 0, 0, 1, -360, 360], 
 62          [3, 4, 0.017,  0.092,  0.158, 0,   250, 250, 0, 0, 1, -360, 360], 
 63          [4, 5, 0.039,  0.17,   0.358, 150, 150, 150, 0, 0, 1, -360, 360], 
 64          [2, 5, 0,      0.0586, 0,     0,   300, 300, 0, 0, 1, -360, 360], 
 65          [5, 6, 0.0119, 0.1008, 0.209, 40,  150, 150, 0, 0, 1, -360, 360], 
 66          [6, 7, 0.0085, 0.072,  0.149, 250, 250, 250, 0, 0, 1, -360, 360], 
 67          [7, 1, 0,      0.0625, 0,     250, 250, 250, 0, 0, 1, -360, 360], 
 68          [7, 8, 0.032,  0.161,  0.306, 250, 250, 250, 0, 0, 1, -360, 360], 
 69          [8, 3, 0.01,   0.085,  0.176, 250, 250, 250, 0, 0, 1, -360, 360] 
 70      ]) 
 71   
 72      ##-----  OPF Data  -----## 
 73      ## area data 
 74      # area refbus 
 75      ppc['areas'] = array([ 
 76          [1, 4] 
 77      ], float) 
 78   
 79      ## generator cost data 
 80      # 1 startup shutdown n x1 y1 ... xn yn 
 81      # 2 startup shutdown n c(n-1) ... c0 
 82      ppc['gencost'] = array([ 
 83          [1, 0, 0, 4,  0, 0, 100, 2000, 200, 4403.5, 270, 6363.5], 
 84          [1, 0, 0, 4,  0, 0, 100, 2500, 200, 5500,   250, 7250.0], 
 85          [2, 0, 0, 2, 15, 0,   0,    0,   0,    0,     0,    0.0] 
 86      ]) 
 87   
 88      ppc['A'] = array([ 
 89          [1, 2, 3, 4,  5,  7,  8,  9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 24, 22, 21, 28, 26, 25, 29, 30], 
 90          [2, 4, 6, 8, 10, 14, 16, 18, 20, 22, 24, 26, 28, 30, 34, 36, 38, 40, 48, 44, 42, 56, 52, 50, 58, 60] 
 91      ], float) 
 92   
 93      ppc['N'] = array([ 
 94          [30, 29, 28, 27, 26, 24, 23, 22, 21, 20, 19, 18, 17, 16, 14, 13, 12, 11, 7,  9,  10, 3, 5,  6,  2, 1], 
 95          [60, 58, 56, 54, 52, 48, 46, 44, 42, 40, 38, 36, 34, 32, 28, 26, 24, 22, 14, 18, 20, 6, 10, 12, 4, 2] 
 96      ], float) 
 97   
 98      return ppc 
 99