How to find zero impedance branch lists

answered 2017-11-29 15:33:25 -0500

perolofl
2947 ●1 ●8 ●16

I didn't find any API returning ZILs so I wrote my own:

def get_zil(sid=-1):
    """ Returns ibus, jbus and id for zero impedance lines in three lists
    """
    ierr, thrshz = psspy.prmdat('THRSHZ')
    ierr, iarray = psspy.abrnint(sid, string=['FROMNUMBER','TONUMBER'])
    ierr, carray = psspy.abrnchar(sid,entry=1, string='ID')
    ierr, cmplxarray = psspy.abrncplx(sid, string=['RX'])
    ibl = []
    jbl = []
    idl = []
    for ibus,jbus,id,z in zip(iarray[0],iarray[1],carray[0],cmplxarray[0]):
        if z.real==0. and abs(z.imag)<=thrshz:
            ibl.append(ibus)
            jbl.append(jbus)
            idl.append(id)
    return ibl,jbl,idl

Create lists with:

ibusl,jbusl,idl = get_zil()