addwind

ADDWIND   Adds wind generators and corresponding xGenData to existing data.

   [IDX, NEW_MPC] = ADDWIND(WIND, MPC)
   [IDX, NEW_MPC, NEW_XGD] = ADDWIND(WIND, MPC)
   [IDX, NEW_MPC, NEW_XGD] = ADDWIND(WIND, MPC, XGD)

   Given a WindUnitData structure, or the name of a file containing such
   a structure, this function adds the specified wind generators to an
   existing MATPOWER case and xGenData struct.

   Inputs:
       WIND : a WindUnitData struct or the name of an M-file
                   or MAT-file that returns one, with the following fields
           .gen     : rows to be appended to the GEN matrix from MPC
           .gencost : (optional) rows to be added to the GENCOST matrix
                      from MPC, default is zero cost
           .xgd_table : xGenData table struct or filename providing data for
                        the wind units being added. See LOADXGENDATA for
                        more information on the xGenData table format.
       MPC : MATPOWER case struct to which wind generators will be added
       XGD : (optional) xGenData struct corresponding to the generators
             already in MPC, to which the new xGenData for the wind units
             will be added.

   Output:
       IDX : Generator indices of newly added wind units.
       NEW_MPC : MPC with wind units appended to MPC.GEN and MPC.GENCOST
                 MPC.GENFUEL (= 'wind').
       NEW_XGD : XGD with xGenData for new wind units appended.

   See also LOADXGENDATA.

0001 function [idx, new_mpc, new_xgd] = addwind(wind, mpc, xgd)
0002 %ADDWIND   Adds wind generators and corresponding xGenData to existing data.
0003 %
0004 %   [IDX, NEW_MPC] = ADDWIND(WIND, MPC)
0005 %   [IDX, NEW_MPC, NEW_XGD] = ADDWIND(WIND, MPC)
0006 %   [IDX, NEW_MPC, NEW_XGD] = ADDWIND(WIND, MPC, XGD)
0007 %
0008 %   Given a WindUnitData structure, or the name of a file containing such
0009 %   a structure, this function adds the specified wind generators to an
0010 %   existing MATPOWER case and xGenData struct.
0011 %
0012 %   Inputs:
0013 %       WIND : a WindUnitData struct or the name of an M-file
0014 %                   or MAT-file that returns one, with the following fields
0015 %           .gen     : rows to be appended to the GEN matrix from MPC
0016 %           .gencost : (optional) rows to be added to the GENCOST matrix
0017 %                      from MPC, default is zero cost
0018 %           .xgd_table : xGenData table struct or filename providing data for
0019 %                        the wind units being added. See LOADXGENDATA for
0020 %                        more information on the xGenData table format.
0021 %       MPC : MATPOWER case struct to which wind generators will be added
0022 %       XGD : (optional) xGenData struct corresponding to the generators
0023 %             already in MPC, to which the new xGenData for the wind units
0024 %             will be added.
0025 %
0026 %   Output:
0027 %       IDX : Generator indices of newly added wind units.
0028 %       NEW_MPC : MPC with wind units appended to MPC.GEN and MPC.GENCOST
0029 %                 MPC.GENFUEL (= 'wind').
0030 %       NEW_XGD : XGD with xGenData for new wind units appended.
0031 %
0032 %   See also LOADXGENDATA.
0033 
0034 %   MOST
0035 %   Copyright (c) 2013-2016, Power Systems Engineering Research Center (PSERC)
0036 %   by Ray Zimmerman, PSERC Cornell
0037 %
0038 %   This file is part of MOST.
0039 %   Covered by the 3-clause BSD License (see LICENSE file for details).
0040 %   See http://www.pserc.cornell.edu/matpower/ for more info.
0041 
0042 %% define named indices into data matrices
0043 [PW_LINEAR, POLYNOMIAL, MODEL, STARTUP, SHUTDOWN, NCOST, COST] = idx_cost;
0044 
0045 %% define fuel type for wind
0046 WIND_FUEL = 'wind';
0047 
0048 %% input arg handling
0049 if nargin < 3
0050     xgd = [];
0051 end
0052 if ischar(wind)
0053     infile = sprintf(' in file: ''%s''', wind);
0054     wind = loadgenericdata(wind, 'struct', 'gen', 'wind');
0055 else
0056     infile = '';
0057 end
0058 
0059 %% add to MPC
0060 nw = size(wind.gen, 1);
0061 if isfield(wind, 'gencost')
0062     wind_gencost = wind.gencost;
0063 else        %% use zero cost by default
0064     wind_gencost = repmat([POLYNOMIAL 0 0 2 0 0], nw, 1);
0065 end
0066 [new_mpc, idx] = addgen2mpc(mpc, wind.gen, wind_gencost, WIND_FUEL);
0067 
0068 %% handle xGenData
0069 if nargout > 2      %% output NEW_XGD requested
0070     if isfield(wind, 'xgd_table')
0071         wind_xgd = loadxgendata(wind.xgd_table, wind.gen);
0072     else
0073         error('addwind: missing XGD_TABLE field in WIND');
0074     end
0075 
0076     if isempty(xgd)     %% no input XGD provided
0077         new_xgd = wind_xgd;
0078     else                %% input XGD provided
0079         new_xgd = xgd;
0080         fields = fieldnames(xgd);
0081         for f = 1:length(fields)    %% append rows of every field in xgd
0082             ff = fields{f};
0083             %% dims of wind_xgd fields already checked by loadxgendata
0084             if size(xgd.(ff), 1) ~= size(mpc.gen, 1)
0085                 error('addwind: # of rows in XGD.%s (%d) does not match MPC.GEN (%d)', ...
0086                     ff, size(xgd.(ff), 1), size(mpc.gen, 1));
0087             end
0088             new_xgd.(ff) = [xgd.(ff); wind_xgd.(ff)];
0089         end
0090     end
0091 end