k_factor31 = 100e6;   % Conversion factor
k_factor31r = 1000e6; % Conversion factor
 
R31 = 0.0112013;    % R [pu] - Winding kV, system MVA
X31 = 0.367610;     % X [pu] - Winding kV, system MVA
 
R31r = 113400.007813;   % R [W]  - Load loss 
X31r = 0.118200;        % X [pu] - Z pu
 
V1_31 = 132e3;       % Nominal V1
V2_31 = 11e3;        % Nominal V2
 
V1_31r = 125e3;      % Nominal V1
V2_31r = 11.5e3;     % Nominal V2
 
MVA = 31.5e6;     	% MVA-base
 
%% PSSE-modell fra Statnett
Z31_base = V1_31^2/k_factor31;    % I Ohm (kfaktor pga System MVA)
R31_ohm = R31*Z31_base;
X31_ohm = X31*Z31_base;
Z31_ohm = abs(R31_ohm + 1j*X31_ohm);
I31_nom = sqrt(MVA / Z31_ohm);
 
%% Konvertert modell
 
Z31r_base = V1_31r^2 / (3*MVA); % I Ohm (3*MVA pga Winding MVA)
R31r_pu = R31r / (MVA * 3);    % pu av MVA-rating  
R31r_ohm = R31r_pu * Z31r_base;
X31r_ohm = 3*X31r*Z31_base;
Z31r_ohm = abs(R31r_ohm + 1j*X31r_ohm);
I31r_nom = 3*MVA / V1_31r;
 
Z_compare = [Z31_base Z31r_base];
R_compare = [R31_ohm R31r_ohm];
X_compare = [X31_ohm X31r_ohm];
Z_compare = [Z31_ohm Z31r_ohm];
I_compare = [I31_nom I31r_nom];
 
fprintf('\n########### Comparison ##############\n')
fprintf('               PSSE 31     PSSE 31r\n')
 
fprintf('Z base  : %12f %12f \nR [ohm] : %12f %12f \nX [ohm] : %12f %12f \nI [Amps]: %12f %12f\n', ...
    [Z_compare; R_compare; X_compare; I_compare].')

a19mYWN0b3IzMSA9IDEwMGU2OyAgICUgQ29udmVyc2lvbiBmYWN0b3IKa19mYWN0b3IzMXIgPSAxMDAwZTY7ICUgQ29udmVyc2lvbiBmYWN0b3IKClIzMSA9IDAuMDExMjAxMzsgICAgJSBSIFtwdV0gLSBXaW5kaW5nIGtWLCBzeXN0ZW0gTVZBClgzMSA9IDAuMzY3NjEwOyAgICAgJSBYIFtwdV0gLSBXaW5kaW5nIGtWLCBzeXN0ZW0gTVZBCgpSMzFyID0gMTEzNDAwLjAwNzgxMzsgICAlIFIgW1ddICAtIExvYWQgbG9zcyAKWDMxciA9IDAuMTE4MjAwOyAgICAgICAgJSBYIFtwdV0gLSBaIHB1CgpWMV8zMSA9IDEzMmUzOyAgICAgICAlIE5vbWluYWwgVjEKVjJfMzEgPSAxMWUzOyAgICAgICAgJSBOb21pbmFsIFYyCgpWMV8zMXIgPSAxMjVlMzsgICAgICAlIE5vbWluYWwgVjEKVjJfMzFyID0gMTEuNWUzOyAgICAgJSBOb21pbmFsIFYyCgpNVkEgPSAzMS41ZTY7ICAgICAJJSBNVkEtYmFzZQoKJSUgUFNTRS1tb2RlbGwgZnJhIFN0YXRuZXR0ClozMV9iYXNlID0gVjFfMzFeMi9rX2ZhY3RvcjMxOyAgICAlIEkgT2htIChrZmFrdG9yIHBnYSBTeXN0ZW0gTVZBKQpSMzFfb2htID0gUjMxKlozMV9iYXNlOwpYMzFfb2htID0gWDMxKlozMV9iYXNlOwpaMzFfb2htID0gYWJzKFIzMV9vaG0gKyAxaipYMzFfb2htKTsKSTMxX25vbSA9IHNxcnQoTVZBIC8gWjMxX29obSk7CgolJSBLb252ZXJ0ZXJ0IG1vZGVsbAoKWjMxcl9iYXNlID0gVjFfMzFyXjIgLyAoMypNVkEpOyAlIEkgT2htICgzKk1WQSBwZ2EgV2luZGluZyBNVkEpClIzMXJfcHUgPSBSMzFyIC8gKE1WQSAqIDMpOyAgICAlIHB1IGF2IE1WQS1yYXRpbmcgIApSMzFyX29obSA9IFIzMXJfcHUgKiBaMzFyX2Jhc2U7ClgzMXJfb2htID0gMypYMzFyKlozMV9iYXNlOwpaMzFyX29obSA9IGFicyhSMzFyX29obSArIDFqKlgzMXJfb2htKTsKSTMxcl9ub20gPSAzKk1WQSAvIFYxXzMxcjsKClpfY29tcGFyZSA9IFtaMzFfYmFzZSBaMzFyX2Jhc2VdOwpSX2NvbXBhcmUgPSBbUjMxX29obSBSMzFyX29obV07ClhfY29tcGFyZSA9IFtYMzFfb2htIFgzMXJfb2htXTsKWl9jb21wYXJlID0gW1ozMV9vaG0gWjMxcl9vaG1dOwpJX2NvbXBhcmUgPSBbSTMxX25vbSBJMzFyX25vbV07CgpmcHJpbnRmKCdcbiMjIyMjIyMjIyMjIENvbXBhcmlzb24gIyMjIyMjIyMjIyMjIyNcbicpCmZwcmludGYoJyAgICAgICAgICAgICAgIFBTU0UgMzEgICAgIFBTU0UgMzFyXG4nKQogICAgICAgICAKZnByaW50ZignWiBiYXNlICA6ICUxMmYgJTEyZiBcblIgW29obV0gOiAlMTJmICUxMmYgXG5YIFtvaG1dIDogJTEyZiAlMTJmIFxuSSBbQW1wc106ICUxMmYgJTEyZlxuJywgLi4uCiAgICBbWl9jb21wYXJlOyBSX2NvbXBhcmU7IFhfY29tcGFyZTsgSV9jb21wYXJlXS4nKQ==