Roger Dugan

For power flow, three different ways to represent power vs. Voltage: 1 - Constant kVA 2 - Constant impedance 3 - Constant P. |V| objective (set PU = property)

The single-phase loads are all defined with the incorrect voltage rating. Instead of 12.66, it should be 12.66/sqrt(3) = 7.3 kV. The way it is currently defined the loads declared with Phases=1 are all low. The 33-bus system was set up as a positive-sequence model, so it is a little tricky to make it into an unbalanced model with 1-phase loads. You have to pay careful attention to the definition of single-phase elements. To do this right, you would need R0 and X0 defined for the line segments. Since...

If you don't change modes from Daily mode, the monitors will continue to collect for each Solve command. Solving multiple days for 86400 points is going to create a lot of points. Are you sure you want to do that?

The command Show Voltage LL Nodes gives one report of LL voltages. You can use Python to move from Line to Line and get the voltages for the active CktElement.

When you set a Bus active, you can use the puVmagAngle property to get the pu voltage and angle (degrees) for each node at the bus. Try to find the OpenDSS_COM.chm file. It has this information in it.

We don't maintain the 33-bus system. I wouldn't be surprised if someone somewhere has it running in OpenDSS, but we don't supply it with OpenDSS.

I would guess a syntax error such as a missing space in the InvControl definition or the PVSystem definition.

I haven't traced out your circuit completely, but you should be able to make it work. I noticed one thing when I executed your code: the loads are floating. Nothing appears to be grounded, which may make interpreting the power flows a little difficult. The voltages a node to the system Y reference and the phase angle of the currents can be different. The sum of the powers usually works but what is reported for the individual terminals may seem weird. This is the Blondel theorem in action. If you...