thanks!!!!

Is there a way to check if the case solve or not and output that status to a text file?

agreed -the keyword function is just my personal preference - I tend to forget what those values represent after a period of time.

If I understands your question correctly, assuming you want to use the PSSE API python method to take a branch out of service in the powerflow model.

Here's the standard python syntax per PSSE_API document - i = frombus, j = tobus, ckt = circuit id, intgar is an array of 6 elements, so using the keyword function, intgar1 is where the branch status can be set from 1 to 0 where 1 = in service.

psspy.branchdata3( i=12345,j=34567,ckt = '1',intgar1 = 0)

This is how I run my dynamic simulations with multiple fault files and converted models.

Here's a simple way to do it - Others may have better codes than this, but this works well for me. Of course you could add additional details in between as you wish but it should get what you are asking for. workdir and flt are folders but you could specifically assign what workdir path is or completely forget it and have everything input/output from your main folder.

Reading and storing all fault files in flt folder

    mypath_ = work_dir + "/flt"
    flt_file = [a for a in listdir(mypath_) if isfile(join(mypath_, a))]

    #splitting up file extensions
    flt_infile = []
    for r in flt_file:
            flt_infile.append(os.path.splitext(r)[0])

    print 'FAULT_DEF = ' , flt_infile
    print (" \n")

Reading and storing all converted cases convlf folder

    mypath = work_dir + "/convlf"
    clf_file = [g for g in listdir(mypath) if isfile(join(mypath, g))]

    #splitting up the .sav extension
    clf_infile = []
    for t in clf_file:
            clf_infile.append(os.path.splitext(t)[0])

    print 'CONVLF_DEF = ' , clf_infile
    print (" \n")

Looping through all cases and flt files

    flt_flag = 0     #1 = idv files, 0 = python files
    dist_run = 20.0     #20 seconds run

    for x in clf_infile:
        for y in flt_infile:

            psspy.case(work_dir + "/convlf/" + x + ".sav")
            psspy.rstr(work_dir + "/snap/" + snapfile + ".snp")

            psspy.strt(option = 0,
                       outfile = work_dir + "/bin/" + x + "-" + y + ".out"
                       )

            # run disturbance in steady state for 0.1 sec
            psspy.run(option = 0,
                      tpause = flat_start,
                      nprt = 0,
                      nplt = 1,
                      crtplt = 0
                      )

            # run disturbance and clear fault
            if flt_flag == 0:
                    execfile(work_dir + "/flt/" + y + ".py")

            else:
                    psspy.runrspnsfile(work_dir + "/flt/" + y + ".idv")

            # finish 20 second run
            psspy.run(option = 0, 
                      tpause = dist_run,
                      nprt = 0,
                      nplt = 1,
                      crtplt = 0
                      )

I'm trying to figure out how I can clear only the island that the tripping element is causing during the tripping of the fault element right after distclearfault().

When I use psspy.tree(1,0) and psspy.tree(2,1), it clears ALL islands within the case which I do not want.

That works great if there are island buses, but it is throwing me an ierr = 1 which indicates that there are no island buses. So how would I check if there are island buses or not and if so then invoke the psspy.island()? Ok nvm I figured it out. Thanks for your help!

What's the best way to check for island buses? I found the psspy.tree(1,0) but unfortunately, if there are no island buses, it throws an error when I change the option to 2 to disconnect all island buses. What I would like to do is check if there are island buses not equal to zero, then execute the psspy.tree(2,1).

I'm trying to verify the line constants for my short circuit model - I have 8 different substations and each T line has different length - When I model the structures in PSSE Line Prop, do I need to model all structures and set the circuit location in ft based on my reference substation or do I need to build the structure individually and get the R, X values for just that circuit.

The impedance of each circuit is different with all structures modeled vs individually. What's the correct way?

This code prints out all channels available - where is the actual scan?

Thank you.

I need to scan for low voltages and other identifiers - How do I automate this in python to scan the voltages and if the voltages are below say .70 then flag it as unstable?

The stability script does change the work directory under variable work_dir. So whatever the user has for the work_dir variable as is the working directory.

Thank You!!!

How would I iterate this code? I have like over 400 stability postscript files in a folder. What I would like to do is tell python to go into that directory and start the conversion for all *.ps files 2 pdf in that directory. Thank you!

The small script above worked great alone but when I embed it within my stability run script at the end, it is complaining that it can't locate the idv file. Any suggestions?