Hello! I'm modifying the matpower case30.m in order to evaluate the load curtailment in each PQ bus. So, I added a "fictional generator" to those buses, with a real power output limited by 0 and the real power demand. For example, for bus number i:
0 <= Pgi <= Pdi I did the same for the reactive power. The problem is that, after running the opf, some of those fictional generatores have a positive reactive power output, even when the real power output is zero. My goal is that when there isn't load curtailment in one bus (Pgi = 0), the reactive power output would be also null, i.e., I don't want to manipulate the real and reactive power separately, but as the apparent power. Any suggestions ? Cheers, José Miguel C. Costa
