Python for IPSA (Power flow analysis)

[Dave Angel]

On 05/28/2013 06:00 AM, Debbie wrote:
> Hi there,
> I am new to Python,

Welcome.  Could you tell us a little about yourself, such as whether 
you've experience in a few other languages, or if Python is your first 
programming experience?  Also, what version of Python (presumably 2.7 or 
2.6) and what OS ?

> and wondering if you could help me with python based coding for the IPSA (Power system analysis software). I have a electrical distribution network with generators, buses and loads, on which I am performing the load flow analysis every 1 hour for a period of 1 year.
>
> The code to perform instantaneous load/power flow analysis is given below. I need to modify it such that I can perform this load flow analysis every 1 hour for a period of 1 year. Please help.
>
> from ipsa import *

I have no idea what functionality is in ipsa, so this whole message is a 
guess.  First question is whether the ReadFile() below and the 
DoLoadFlow() read the whole year's data, or data for a particular hour. 
  And if the latter, how do you get the next set of data?


>
> ipsasys = IscInterface()
> net = ipsasys.ReadFile("refinery.iif")
> bok = net.DoLoadFlow();
> if bok:
>      busbars = net.GetBusbars()
>      print "Load Flow results:"
>      print ""
>      print "BusName     Vmag(kV)"
>      print "===================="
>      for bus in busbars.itervalues():
>          name = bus.GetName()
>          vm = bus.GetVoltageMagnitudekV()
>          res = "%-8s  %10.5f" % (name, vm)
>          print res
> else:
>      print "Load Flow failed!"
>
> Regards,
> Debbie
>

First you want to make a function to print out a particular hour's data. 
  That might turn out to be something like:

def one_hour(net):
     busbars = net.GetBusbars()
     print "Load Flow results:"
     print ""
     print "BusName     Vmag(kV)"
     print "===================="
     for bus in busbars.itervalues():
         name = bus.GetName()
         vm = bus.GetVoltageMagnitudekV()
         res = "%-8s  %10.5f" % (name, vm)
         print res

(Just pasted from your code, I added in the probable parameter to the 
definition)

Now your main function might be something like:


def main(filename):
     ipsasys = IscInterface()
     net = ipsasys.ReadFile(filename)

     #if net gets you data for one hour, figure out how
     #to get the data for the whole year, in the form
     #of a list or an iterator called nets

     bok = net.DoLoadFlow();
     if bok:
         for net in nets:
             one_hour(net)

     else:
         print "Load Flow failed!"


#and your top-level code is:
main("refinery.iif")


As to ordering in the source file, put the import first:
    from ipsa import *

then your function definitions one_hour() and main(), then
your top-level code.

-- 
DaveA