function mpc = TwoBus_CSI
%% 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 = [
     1	     3      0	0	0	0	1       1	0	138     1       1.05	0.95; 
     2       1      0	0	0	0	1       1	0	138 	1       1.05	0.95; 
     ];

%% generator data
%	bus   Pg	 Qg	      Qmax	Qmin	Vg      mBase	status	Pmax	Pmin	Unit Code
mpc.gen = [
     1	   0 	  0       100	-25     1.05	100      1      200      0   %	U400
     1	   0 	  0       100	-25     1.05	100      1      200      0   %	U400
     1	   0 	  0       100	-25     1.05	100      1      200      0   %	U400
     1	   0 	  0       100	-25     1.05	100      1      200      0   %	U400
     1	   0 	  0       100	-25     1.05	100      1      200      0   %	U400
     2	   0	  0        0    -22     1.05	100      1       0     -528 %	Flexible load
    ];

%% branch data
%	fbus   tbus	   r	   x	    b	 Con  LTE  STE  ratio	angle  status  angmin	angmax
mpc.branch = [
    1       2	 0.0001	 0.2112	 0.0572  303   0   0     0       0       1     -360	 360; % 138KV,1
    ];

% R_PU and rate(A) were calculated based on Maximum Tc = 100, WindSpeed = 0.6, WindAngle = 90, Ta = 40, SolarRadiationAngle = 0, hour = 3864, Diameter = 28.1 

%% Branch reliability data
mpc.branch_FRRD = [
    0.24 16; % 1 - 2
    ];

%% generator cost data
% Nodel Startup cost  shutdowncost    No.Coeff
   mpc.gencost = [							
    2	1500            0               3           0           130     400; % U200		
    2	1500            0               3           0           130     400; % U200	
    2	1500            0               3           0           130     400; % U200	
    2	1500            0               3           0           130     400; % U200	
    2	1500            0               3           0           130     400; % U200	
    2	1500            0               2           1500    	665     0; % Dispatchable load  
    ];
