📄 psse2psat
字号:
}
$i = $nbus;
$nbus >= 0 && printf OUT "$format];\n\n", $busidx[$i], $kvb[$i],
$ein[$i],$ang[$i],$kae[$i],$kzo[$i];
# -----------------------------------------------------------------------
# write SW.con
# -----------------------------------------------------------------------
if ($nsw >= 0) {
print OUT "SW.con = [ ...\n";
$format = "%4d %7.2f %6.2f %8.5f 0.00 %8.5f %8.5f 1.1 0.9 %8.5f 1";
printf OUT "$format];\n\n", $swbus,$swsb,$kvbase{$swbus},
$swvs,$swqt,$swqb,$swpg;
}
# -----------------------------------------------------------------------
# write PV.con
# -----------------------------------------------------------------------
$npv >= 0 && printf OUT "PV.con = [ ...\n";
$format = "%4d %7.2f %6.2f " . "%8.5f " x 4 . "1.1 0.9 1";
for ($i = 0; $i <= $npv-1; $i++) {
#$h = $idxpv[$i];
printf OUT "$format;\n", $pvbus[$i],$pvsb[$i],$kvbase{$pvbus[$i]},
$pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];
}
#$h = $idxpv[$npv];
$i = $npv;
$npv >= 0 && printf OUT "$format];\n\n",$pvbus[$i],$pvsb[$i],
$kvbase{$pvbus[$i]},$pvpg[$i],$pvvs[$i],$pvqt[$i],$pvqb[$i];
# -----------------------------------------------------------------------
# write PQ.con
# -----------------------------------------------------------------------
$npq >= 0 && printf OUT "PQ.con = [ ...\n";
$format = "%4d $pbas %8.2f %8.5f %8.5f 1.1 0.9 1";
for ($i = 0; $i <= $npq-1; $i++) {
$h = $pqbus[$i];
printf OUT "$format;\n", $h,$kvbase{$h},$pqpl[$i],$pqql[$i];
}
$h = $pqbus[$npq];
$i = $npq;
$npq >= 0 && printf OUT "$format];\n\n",$h,$kvbase{$h},$pqpl[$i],
$pqql[$i];
# -----------------------------------------------------------------------
# write Pl.con
# -----------------------------------------------------------------------
$npl >= 0 && print OUT "Pl.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 %8.5f %8.5f 0 %8.5f %8.5f 0 0";
for ($i = 0; $i <= $npl-1; $i++) {
printf OUT "$format;\n", $plbus[$i],$kvbase{$plbus[$i]},$plyp[$i],
$plip[$i],$plyq[$i],$pliq[$i];
}
$h = $plbus[$npl];
$i = $npl;
$npl >= 0 && printf OUT "$format];\n\n",$plbus[$i],$kvbase{$plbus[$i]},
$plyp[$i],$plip[$i],$plyq[$i],$pliq[$i];
# -----------------------------------------------------------------------
# write Shunt.con
# -----------------------------------------------------------------------
$nsh >= 0 && print OUT "Shunt.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 %8.5f %8.5f";
for ($i = 0; $i <= $nsh-1; $i++) {
if ($psh[$i] != 0 or $qsh[$i] != 0) {
$h = abs($idxsh[$i]);
printf OUT "$format;\n", $h,$kvbase{$h},$psh[$i],$qsh[$i];
}
}
$h = $idxsh[$nsh];
$i = $nsh;
$nsh >= 0 && printf OUT "$format];\n\n",$h,$kvbase{$h},$psh[$i],$qsh[$i];
# -----------------------------------------------------------------------
# write Line.con
# -----------------------------------------------------------------------
$nline >= 0 && print OUT "Line.con = [ ...\n";
#$format = "%4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 3 .
# " 0 %8.5f %8.3f %8.3f %8.3f";
$format = "%4d %4d %8.2f %8.2f 60 0 %5.2f " . "%8.5f " x 5 .
" %8.3f %8.3f %8.3f";
for ($i = 0; $i <= $nline-1; $i++) {
if ($busidw{$bustt[$i]}) {
if ($kvbase{$bustt[$i]} != 0) {
$k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};
} else { # print warnings
$k = 0;
print "WARNING: Check voltage rating of bus ".$busidw{$bustt[$i]}."\n";
}
if ($k == 1) {$k = 0;}
printf OUT "$format;\n", $busfr[$i],$bustt[$i],$lbas[$i],
$kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],$ratio[$i],$phas[$i],
$ratea[$i],$rateb[$i],$ratec[$i];
} else {
print "WARNING: Bus ".$bustt[$i]." does not exists -> Line ".$i." ignored.\n";
}
}
$i = $nline;
if ($nline >= 0) {
if ($busidw{$bustt[$i]}) {
if ($kvbase{$bustt[$i]} != 0) {
$k = $kvbase{$busfr[$i]}/$kvbase{$bustt[$i]};
} else { # print warnings
$k = 0;
print "WARNING: Check voltage rating of bus ".$busidw{$bustt[$i]}."\n";
}
if ($k == 1) {$k = 0;}
printf OUT "$format];\n\n",$busfr[$i],$bustt[$i],$lbas[$i],
$kvbase{$busfr[$i]},$k,$rest[$i],$reat[$i],$susc[$i],$ratio[$i],$phas[$i],
$ratea[$i],$rateb[$i],$ratec[$i];
} else {
print "WARNING: Bus ".$bustt[$i]." does not exists -> Line ".$i." ignored\n";
}
}
# -----------------------------------------------------------------------
# write Ltc.con
# -----------------------------------------------------------------------
$nltc >= 0 && print OUT "Ltc.con = [ ...\n";
$format = "%4d %4d %8.2f %8.2f 60 1 0.001 0.1 " .
"%8.5f " x 6 . " %4d %4d";
for ($i = 0; $i <= $nltc-1; $i++) {
printf OUT "$format;\n", $ltc_busfr[$i],$ltc_busto[$i],$ltc_bas[$i],
$kvbase{$ltc_busfr[$i]},$ltc_vma[$i],$ltc_vmi[$i],$ltc_dm[$i],
$ltc_ref[$i],$ltc_x[$i],$ltc_r[$i],$ltc_busct[$i],$ltc_ctr[$i];
}
$i = $nltc;
if ($nltc >= 0) {
printf OUT "$format];\n\n", $ltc_busfr[$i],$ltc_busto[$i],$ltc_bas[$i],
$kvbase{$ltc_busfr[$i]},$ltc_vma[$i],$ltc_vmi[$i],$ltc_dm[$i],
$ltc_ref[$i],$ltc_x[$i],$ltc_r[$i],$ltc_busct[$i],$ltc_ctr[$i];
}
# -----------------------------------------------------------------------
# write Phs.con
# -----------------------------------------------------------------------
$nphs >= 0 && print OUT "Phs.con = [ ...\n";
$format = "%4d %4d %8.2f %8.2f %8.2f 60 0.001 10 50 " . "%8.5f " x 5;
for ($i = 0; $i <= $nphs-1; $i++) {
printf OUT "$format;\n",$phs_busfr[$i],$phs_busto[$i],$phs_bas[$i],
$kvbase{$ltc_busfr[$i]},$kvbase{$ltc_busto[$i]},$phs_ref[$i],
$phs_r[$i],$phs_x[$i],$phs_vma[$i],$phs_vmi[$i];
}
$i = $nphs;
if ($nphs >= 0) {
printf OUT "$format];\n\n",$phs_busfr[$i],$phs_busto[$i],
$phs_bas[$i],$kvbase{$ltc_busfr[$i]},$kvbase{$ltc_busto[$i]},
$phs_ref[$i],$phs_r[$i],$phs_x[$i],$phs_vma[$i],$phs_vmi[$i];
}
# -----------------------------------------------------------------------
# write Twt.con
# -----------------------------------------------------------------------
if ($flag_version == V29) {
$ntw >= 0 && print OUT "Twt.con = [ ...\n";
$format = "%4d %4d %4d $pbas 60 %5.2f %5.2f %5.2f " . "%8.5f " x 6 .
" 1 " . "%7.3f " x 9;
for ($i = 0; $i <= $ntw-1; $i++) {
printf OUT "$format;\n", $bustw1[$i],$bustw2[$i],$bustw3[$i],
$kvbase{$bustw1[$i]},$kvbase{$bustw2[$i]},$kvbase{$bustw3[$i]},
$r12[$i],$r13[$i],$r23[$i],$x12[$i],$x13[$i],$x23[$i],
$rt1a[$i],$rt1b[$i],$rt1c[$i],$rt2a[$i],$rt2b[$i],$rt2c[$i],
$rt3a[$i],$rt3b[$i],$rt3c[$i];
}
$i = $ntw;
if ($ntw >= 0) {
printf OUT "$format];\n\n", $bustw1[$i],$bustw2[$i],$bustw3[$i],
$kvbase{$bustw1[$i]},$kvbase{$bustw2[$i]},$kvbase{$bustw3[$i]},
$r12[$i],$r13[$i],$r23[$i],$x12[$i],$x13[$i],$x23[$i],
$rt1a[$i],$rt1b[$i],$rt1c[$i],$rt2a[$i],$rt2b[$i],$rt2c[$i],
$rt3a[$i],$rt3b[$i],$rt3c[$i];
}
}
#my ($nsvc,@svcbus,@svcbma,@svcbmi,@svcref);
# -----------------------------------------------------------------------
# write Svc.con
# -----------------------------------------------------------------------
$nsvc >= 0 && print OUT "Svc.con = [ ...\n";
$format = "%4d $pbas %8.2f 60 1 1 10 " . "%8.5f " x 3;
for ($i = 0; $i <= $nsvc-1; $i++) {
printf OUT "$format;\n",$svcbus[$i],$kvbase{$svcbus[$i]},
$svcref[$i],$svcbma[$i],$svcbmi[$i];
}
$i = $nsvc;
if ($nsvc >= 0) {
printf OUT "$format];\n\n",$svcbus[$i],$kvbase{$svcbus[$i]},
$svcref[$i],$svcbma[$i],$svcbmi[$i];
}
# -----------------------------------------------------------------------
# write Varname.bus
# -----------------------------------------------------------------------
for ($i = 0; $i <= $nbus; $i++) {
$busname[$i] =~ s/\'//g;
$busname[$i] = deblank ($busname[$i]);
}
$nbus >= 0 && print OUT "Varname.bus = { ...\n";
$h = ($nbus+1) % 5;
if ($h == 0) {$h = 5;}
if (($nbus+1) > 5) {
for ($i = 0; $i <= $nbus-$h; $i+=5) {
print OUT " '$busname[$i]'; '$busname[$i+1]'; " .
"'$busname[$i+2]'; '$busname[$i+3]'; '$busname[$i+4]';\n";
}
}
print OUT " ";
for ($i = $nbus-$h+1; $i <= $nbus-1; $i++) {
print OUT "'$busname[$i]'; ";
}
print OUT "'$busname[$nbus]'};\n\n";
# -----------------------------------------------------------------------
# close output data file
# -----------------------------------------------------------------------
close(OUT) || die "cannot close $ARGV[1]: $!\n";
if ($pti > 20) {$flag_version = V23;}
print "Coversion completed assuming PSS/E V". $flag_version ." data format.\n";
#########################################################################
#########################################################################
# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub assign {
my $param;
$param = substr($_[0],$_[1],$_[2]);
$param =~ s/\s*//g;
unless ($param) {$param = "0";}
#$param =~ s/^0*\.0*//;
unless ($param =~ /[1-9]/) {$param = $_[3];}
return $param;
}
# -----------------------------------------------------------------------
# function for formatting data and/or assigning default values
# -----------------------------------------------------------------------
sub deblank {
my $name = $_[0];
$name =~ s/( *)$//;
return $name;
}
# -----------------------------------------------------------------------
# function for writing a separator
# -----------------------------------------------------------------------
sub carret {
$verbose && print "-" x 50 . "\n";
}
# -----------------------------------------------------------------------
# function for reading comma separated data
# -----------------------------------------------------------------------
sub mysplit {
if ($_ =~ /\s*,\s*/) {
my $string = $_[0];
if (/,,,/) {
# this is to parse branch data only in PSS/E v26 data files
$string =~ s/\/(.*)//;
$flag_version = V26;
$string =~ s/,,,//;
my @mydata = split /\s+/, $string;
if ($mydata[0] eq "") { shift(@mydata); }
#print "PSS/E V26 data parsed correctly...\n";
return @mydata;
} else { }
$string =~ s/\/(.*)//;
my @mydata = split /\s*,\s*/, $string;
if ($mydata[0] eq "") { shift(@mydata); }
#print "Data parsed correctly...\n";
return @mydata;
} elsif ($_ =~ /\s+/) {
my $string = $_[0];
$string =~ s/\/(.*)//;
# if ($_ =~ /\'/) {
# my $diff = substr"$_",0,18;
# if ($diff =~ /\'/) {
# $flag_version = V26;
# print "$diff\n";
# } else {
# $flag_version = 3;
# print "$diff V23\n";
# }
# }
$string =~ s/\'//g;
my @mydata = split /\s+/, $string;
if ($mydata[0] eq "") { shift(@mydata); }
#print "Data parsed correctly...\n";
return @mydata;
}
}
# -----------------------------------------------------------------------
# function for finding multiple PQ loads at the same bus
# -----------------------------------------------------------------------
sub pqidx {
my $mypqbus = $_[0];
my $check = -1;
for ($i = 0; $i <= $npq; $i++) {
if ($pqbus[$i] == $mypqbus) {
$check = $i;
last;
}
}
return $check;
}
# -----------------------------------------------------------------------
# function for finding multiple PV generators at the same bus
# -----------------------------------------------------------------------
sub pvidx {
my $mypvbus = $_[0];
my $check = -1;
for ($i = 0; $i <= $npv; $i++) {
if ($pvbus[$i] == $mypvbus) {
$check = $i;
last;
}
}
return $check;
}
# -----------------------------------------------------------------------
# function for determining if the bus is isolated
# -----------------------------------------------------------------------
sub ididx {
my $myisbus = $_[0];
my $check = 0;
for ($i = 0; $i <= $nisl; $i++) {
if ($islbus[$i] == $myisbus) {
$check = 1;
last;
}
}
return $check;
}
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -