Check out the 'help get_profile'. Much like you pass the previous profiles object on the load profile line, you can pass the previous profiles object on the other lines:
profiles = getprofiles('wind_profile1',iwind); profiles = getprofiles('wind_profile2',profiles,iwind); profiles = getprofiles('load_profile_30',profiles) See if that works! On Thu, Sep 21, 2017 at 7:24 PM, Joshua Sebben <joshkeep...@gmail.com> wrote: > Stephen, > > Thank you for your reply. > > casefile = 'case30_modified_MedPen'; > mpc = loadcase(casefile); > xgd = loadxgendata('xgd_uc_30_MedPen', mpc); > [iwind, mpc, xgd] = addwind('wind_uc_30_MedPen', mpc, xgd); > profiles = getprofiles('wind_profile1', iwind); > profiles = getprofiles('wind_profile2', iwind); > profiles = getprofiles('wind_profile3', iwind); > profiles = getprofiles('wind_profile4', iwind); > profiles = getprofiles('wind_profile5', iwind); > profiles = getprofiles('wind_profile6', iwind); > profiles = getprofiles('wind_profile7', iwind); > profiles = getprofiles('load_profile_30', profiles); > nt = size(profiles(1).values, 1); % number of periods > > I created the 7 wind profiles with each of the row counts changed to 1 to > 7 for their respective generators. > > This would then change my main code assuming to the snippet above? > This again also has the problem of only changing the last generator in the > list as each of the previous wind profile changes are overwritten by the > new getprofile? > > Regards, > Josh > > On 22 September 2017 at 04:02, Stephen Suffian <ssuff...@villanova.edu> > wrote: > >> I also couldn't figure it out, but I got around it by adding a separate >> profile for each generator and changing the rows value. So if a wind >> generator is in row 5, you would have the struct look like this below. >> >> >> windprofile = struct( ... >> 'type', 'mpcData', ... >> 'table', CT_TGEN, ... >> 'rows', 5, ... >> 'col', PMAX, ... >> 'chgtype', CT_REL, ... >> 'values', [] ); >> >> And then for row 6 you would add another profile with the struct looking >> tlike this: >> >> >> windprofile = struct( ... >> 'type', 'mpcData', ... >> 'table', CT_TGEN, ... >> 'rows', 6, ... >> 'col', PMAX, ... >> 'chgtype', CT_REL, ... >> 'values', [] ); >> >> It is a bit of a work around, but I believe it should work (I did >> something similar to set the Pmax for each conventional generator in a >> profile, so I imagine it will work the same). >> >> On Thu, Sep 21, 2017 at 1:52 PM, Stephen Suffian < >> stephen.suff...@gmail.com> wrote: >> >>> I also couldn't figure it out, but I got around it by adding a separate >>> profile for each generator and changing the rows value. So if a wind >>> generator is in row 5, you would have the struct look like this below. >>> >>> >>> windprofile = struct( ... >>> 'type', 'mpcData', ... >>> 'table', CT_TGEN, ... >>> 'rows', 5, ... >>> 'col', PMAX, ... >>> 'chgtype', CT_REL, ... >>> 'values', [] ); >>> >>> And then for row 6 you would add another profile with the struct looking >>> tlike this: >>> >>> >>> windprofile = struct( ... >>> 'type', 'mpcData', ... >>> 'table', CT_TGEN, ... >>> 'rows', 6, ... >>> 'col', PMAX, ... >>> 'chgtype', CT_REL, ... >>> 'values', [] ); >>> >>> It is a bit of a work around, but I believe it should work (I did >>> something similar to set the Pmax for each conventional generator in a >>> profile, so I imagine it will work the same). >>> >>> On Thu, Sep 21, 2017 at 5:39 AM, Joshua Sebben <joshkeep...@gmail.com> >>> wrote: >>> >>>> By the way, >>>> >>>> >>>> >>>> I have also tried setting the row count to [1 2 3 4 5 6 7] for my 7 >>>> extra generators that I want to add, however when I run my code I get an >>>> error: >>>> >>>> >>>> >>>> Error using apply_profile (line 148) >>>> >>>> apply_profile: third dimension of profile.values should match length of >>>> pro= >>>> >>>> file.rows >>>> >>>> >>>> >>>> Error in loadmd (line 508) >>>> >>>> optab =3D apply_profile(profiles(p), optab); >>>> >>>> >>>> >>>> Error in Test (line 33) >>>> >>>> mdi =3D loadmd(mpc, transmat, xgd, [], [], profiles); >>>> >>>> >>>> Regards, >>>> >>>> Josh >>>> >>>> On 20 September 2017 at 22:17, Joshua Sebben <joshkeep...@gmail.com> >>>> wrote: >>>> >>>>> Currently working on adding wind generator units to my model. I am >>>>> working off the 30 bus example case in MOST. I am trying to run the >>>>> following code >>>>> >>>>> casefile = 'case30'; >>>>> mpc = loadcase(casefile); >>>>> xgd = loadxgendata('xgd_uc', mpc); >>>>> [iwind, mpc, xgd] = addwind('wind_uc_30', mpc, xgd); >>>>> profiles = getprofiles('wind_profile', iwind); >>>>> profiles = getprofiles('load_profile', profiles); >>>>> nt = size(profiles(1).values, 1); % number of periods >>>>> >>>>> %%----- Full Transition Probabilities ----- >>>>> transmat = ex_transmat(nt); >>>>> mdi = loadmd(mpc, transmat, xgd, [], [], profiles); >>>>> mdo = most(mdi, mpopt); >>>>> if verbose >>>>> ms = most_summary(mdo); >>>>> end >>>>> >>>>> However the added generator units of which there are 7 of them in >>>>> wind_uc_30 don't seem to get the profile set to them in the output. >>>>> Instead >>>>> only the first Generator unit gets the profile while the rest of them are >>>>> set at PMAX. >>>>> >>>>> windprofile = struct( ... >>>>> 'type', 'mpcData', ... >>>>> 'table', CT_TGEN, ... >>>>> 'rows', 1, ... >>>>> 'col', PMAX, ... >>>>> 'chgtype', CT_REL, ... >>>>> 'values', [] ); >>>>> >>>>> windprofile.values(:, :, 1) = [ >>>>> 0.80; >>>>> 0.65; >>>>> 0.60; >>>>> 0.82; >>>>> 1.00; >>>>> 0.70; >>>>> 0.50; >>>>> 0.85; >>>>> 1.00; >>>>> 1.10; >>>>> 1.06; >>>>> 0.95; >>>>> >>>>> Above is a snippet from the wind_profile. I am assuming it has >>>>> something to do with the row count, however, I haven't been able to find a >>>>> solution. >>>>> >>>>> Could I please get some help. to ensure all the generator units >>>>> receive the load profile. >>>>> >>>>> Thankyou, >>>>> >>>>> -- >>>>> J.Sebben >>>>> >>>> >>>> >>>> >>>> -- >>>> J.Sebben >>>> >>> >>> >> > > > -- > J.Sebben >