Revision history [back]

Picking up on the thread in your question, here is how I'd add load to a bunch of buses and measure the response at the slack bus:

electrical_distance = {}

buses = subsystem_info('bus', ["NUMBER", "TYPE", "AREA"])
buses_in_qld = filter(qld_only, buses)

for busno, type, area in buses_in_qld:

  # rolls back to original saved case on each for-loop iteration.
  with transaction():

    # add one MW load to this bus
    add_load_to(busno, {"pl": 1}) 

    # solve probably with fixed taps, shunts, fixed everything
    psspy.fnsl(**options)

    # difference at slack bus
    slack_diff = original_slack - get_p(slack_busno)

    electrical_distance[busno] = slack_diff

I've tied some code up into helpful functions. You'd have to write add_load_to, define some options for fnsl and write a function get_p.

The transaction function has already been written: http://www.whit.com.au/blog/2011/09/turn-back-time-using-python-and-psse/

Picking up on the thread in your question, here is how I'd add load to a bunch of buses and measure the response at the slack bus:

electrical_distance = {}

buses = subsystem_info('bus', ["NUMBER", "TYPE", "AREA"])
buses_in_qld = filter(qld_only, buses)

for busno, type, area in buses_in_qld:

  # rolls back to original saved case on each for-loop iteration.
  with transaction():

    # add one MW load to this bus
    add_load_to(busno, {"pl": 1}) 

    # solve probably with fixed taps, shunts, fixed everything
    psspy.fnsl(**options)

    # difference at slack bus
    slack_diff = get_p(slack_busno) - original_slack - get_p(slack_busno)
 
    electrical_distance[busno] = slack_diff

I've tied some code up into helpful functions. You'd have to write add_load_to, define some options for fnsl and write a function get_p.

The transaction function has already been written: http://www.whit.com.au/blog/2011/09/turn-back-time-using-python-and-psse/