You’ve got the right idea. You can simply set the branch status to 0 and then use find_islands() and to find the isolated buses and set their BUS_TYPE to NONE (4). You may then be left with multiple islands, also found via find_islands() or extract_islands(). To solve a power flow, each island needs its own reference bus (which may have been eliminated even for the main system). The load / generation mismatch could also be potential very large in some islands, and that would show up in the slack generation (negative, or unreasonably large positive).
So I think the main challenge in automating everything is in selecting a new reference bus for islands without one. I’d probably go with the first PV bus, and if there are no PV buses, then that whole island blacks out and you can set all its buses to type NONE. Hope this helps, Ray On Dec 4, 2021, at 12:58 PM, prem panigrahi <prempanigr...@outlook.com<mailto:prempanigr...@outlook.com>> wrote: Dear Ray sir, is their any direct process to remove bus along with their branches . I am using mpc.branch method to remove a bus with its branches. But I am facing difficulties when I have to remove more than 100 buses. mpc.branch([6161;7479;7480;7488;7489;7490], BR_STATUS) = 0; results=rundcpf(mpc,opt) [groups, isolated] = find_islands(mpc) mpc_list = extract_islands(mpc) [GROUPS, ISOLATED] = case_info(mpc) mpc.bus([ISOLATED],BUS_TYPE)=4 results=rundcpf(mpc,opt) Sent from Mail<https://go.microsoft.com/fwlink/?LinkId=550986> for Windows