Thank you Dr. Zimmerman. I am more inclined to the third answer. The challenge for me right now is to update the [*Ybus*] matrix whenever any generator hits its reactive limit, and that I can do by adding the generator synchronous reactance in series to the diagonal element in the [Ybus] matrix, corresponding to that generator.
Thanks again Aaqib
