Rbulsara, thank you for mentioning the IEEE book, I will try and find it. It looks like a useful book. I know the PU system very well and I also teach it.
Z (pu) =Z(ohm)*(Vbase (kv))/(Sbase (MVA)^2). Clearly, if you make Vbase =1kV and Sbase=1MVA, then Z(pu)=Z(Ohm). Actually may be I exaggerated the extent of the descrepancies in the results, the Real Power flows balance exactly, the reactive powers by about up to 1 to 5%. The voltage solution is correct. I am beginning to think now that in the example which I thought had descrepancies, my (NEPLAN) results are correct and I will continue with my method for more examples and compare. I will continue like this until I am 100% sure since it will help greatly if I do not have to convert literally hundreds of input data by hand and in my opinion for no apparent reason.
mm