Thanks! I also saw it after I fixed the issue... On Mon, Mar 16, 2020 at 12:37 PM Ray Daniel Zimmerman <[email protected]> wrote:
> Right, I had forgotten about that, but I see it does say in the help for > scenarios_ACTIVSg200 … > > This change table provides zonal load changes for 8760 hours > for case_ACTIVSg200 for the year 2017. Requires swapping the BUS_AREA > and ZONE columns in the BUS matrix of the MATPOWER case struct before > applying via APPLY_CHANGES. > > -- Ray > > > On Mar 12, 2020, at 1:56 PM, Steven G <[email protected]> wrote: > > Thanks for the suggestion! The problem was caused by the settings in > case_ACTIVSg200, there is only 1 area but 6 zones (indexed from 2 to 7). > After adding area numbers `mpc.bus(:, BUS_AREA) = mpc.bus(:, ZONE) - 1; `, > there is no issue creating the mdi structure, but a new problem showed up > that UC has no solution.. > > On Thu, Mar 12, 2020 at 12:29 PM Ray Daniel Zimmerman <[email protected]> > wrote: > >> From the error message, it would seem that your base case does not have >> any loads defined. You can use case_info() on it to see. >> >> Ray >> >> >> On Mar 11, 2020, at 3:03 PM, Steven G <[email protected]> wrote: >> >> Thanks Dr. Zimmerman! >> >> I tried your suggestions, but got the following error message. >> >> Error using scale_load (line 248) >> scale_load: impossible to make zone 1 load equal 144.2 by >> scaling non-existent loads >> >> Error in apply_changes (line 352) >> [mpc.bus, mpc.gen] = scale_load(dmd, mpc.bus, >> mpc.gen, load_zone, opt); >> >> Error in most (line 411) >> mpc = apply_changes(label, mpc, >> mdi.tstep(t).OpCondSched(j).tab); >> >> Error in run_SCUC (line 44) >> mdo = most(mdi, mpopt); >> >> >> On Tue, Mar 10, 2020 at 4:13 PM Ray Daniel Zimmerman <[email protected]> >> wrote: >> >>> I think you want to define your load profile something like this … >>> >>> profiles = struct( ... >>> 'type', 'mpcData', ... >>> 'table', CT_TAREALOAD, ... >>> 'rows', [1:na], ... >>> 'col', CT_LOAD_ALL_P, ... >>> 'chgtype', CT_REP, ... >>> 'values', [] ); >>> profiles.values(:, 1, :) = area_load; >>> >>> Hope this helps, >>> >>> Ray >>> >>> >>> >>> On Mar 9, 2020, at 8:00 PM, Steven G <[email protected]> wrote: >>> >>> Dear Matpower Users, >>> This post is a follow-up of my previous post: >>> https://www.mail-archive.com/[email protected]/msg07377.html >>> >>> I'm trying to solve SCUC for the 200-bus system (`case_ACTIVSg200`),and >>> created the following example (code attached as .m file, and also pasted at >>> the end). However, Matlab returns error message "Unrecognized function or >>> variable 'PR_TAREALOAD'.", since the variable `PR_TAREALOAD` is not defined >>> in the function `loadmd.m`. >>> >>> In addition, I'm aware of this way of defining profile.load is a bad >>> idea, but I just haven't not figured out the "right way" of doing it using >>> `apply_changes.m` or `apply_profile.m`. >>> >>> Thanks, >>> Steven >>> >>> clear; clc; >>> >>> define_constants; >>> casename = 'case_ACTIVSg200'; >>> >>> " OutmailID: 124442730, List: 'matpower-l', MemberID: 78095970 SCRIPT: >>> "Raw Settings >>> mpc = loadcase(casename); >>> nb = size(mpc.bus, 1); nl = size(mpc.branch, 1); ng = size(mpc.gen, 1); >>> nt = 24; % 24 hours >>> na = 6; % 6 areas >>> scenarios = scenarios_ACTIVSg200; % get a change table >>> >>> % PROFILES can take the form of >>> % a struct with the following fields: >>> % .wind: 3-dim array (NT x NJ_MAX x num wind sites) >>> % .load: 3-dim array (NT x NJ_MAX x num load zones) >>> >>> load_data = scenarios(1:nt*na, end); % last column: loads of 6 areas >>> area_load = reshape(load_data, na, nt)'; % nt-by-na matrix >>> profiles.load = zeros(nt,1,na); >>> profiles.load(:,1,:) = area_load; >>> profiles.wind = []; >>> >>> >>> " TCL MERGE ERROR ( 03/10/2020 16:13:04 ): "invalid command name "Raw" >>> Construct MOST struct >>> mdi = loadmd(mpc, nt, [], [], [], profiles); >>> >>> " OutmailID: 124442730, List: 'matpower-l', MemberID: 78095970 SCRIPT: >>> "Solve SCUC >>> mpopt = mpoption(mpopt, 'most.dc_model', 0); >>> >>> " TCL MERGE ERROR ( 03/10/2020 16:13:04 ): "invalid command name "Solve" >>> Issues: >>> % 1. Missing xgd data for the 200-bus synthetic system >>> % 2. According to the comments in `loadmd.m`, profiles.load is from an >>> % older version >>> <Undefined_PR_TAREALOAD.m> >>> >>> >>> <run_SCUC.m> >> >> >> >
