summaryrefslogtreecommitdiffstats
path: root/benchmarks/dc39.dat
blob: 66f2d337b975ff359675e322a56bc1f2cf2089b7 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
BEGIN NEW DATA CASE
C     BENCHMARK DC-39
C      Test of supporting routine "ARRDAT" which reads ZnO data, and punches the
C      associated Type-92 EMTP branch cards.  See Rule Book, end of Sect. XIX-I.
C      1st of 3 subcases has manual selection of a single exponential: NEXP = 1.
C      As for the (i, v) data points,  these come from the evaluation of the ZnO
C      function   i = 2500 * ( v / 778000 ) ** 26   as used by 1st case of DC-37
C      Note that paralleling (D4 = number of columns) is used, although the more
C      modern and flexible approach would be to delay this until usage (variable
C      COL  of columns 76-80  of the first Type-92 ZnO miscellaneous data card).
C         Sept, 1997,  double the size by adding a second subcase that processes
C         the same 3 data subsets only using kilo scaling (KV and KA input data)
ZNO FITTER    { Special-request to transfer to  SUBROUTINE ARRDAT  of overlay 40
$ERASE
BRANCH  RECA        RECB        RECC    { Name the 3 terminal pairs for punching
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
           1           3           2           0     412500.     380000.
C         A1          A2          A3          A4          A5        AMIN
     707.107     192000.        .962         3.0     .892857        .001
         1.0     1164.8            { First point of (I,V) characteristic for ZnO
         2.0     1181.6            { Except for rounding,  these define a  nice,
         5.0     1198.4            { smooth,  exponential function with exponent
         10.     1209.6            { equal to  26.  Scaling is set to draw  2500
         20.     1232.0            { amperes at 778 kV. How?  Consult Vladimir's
         50.     1260.0            { explanation, Rule Book p. 103k, Section 3.0
        100.     1288.0
        200.     1323.84
        500.     1388.8
       1000.     1442.56
       2000.     1512.0
       3000.     1556.8             { last point of (I,V) characteristic for ZnO
BLANK card bounds points of the ZnO characteristic
$PUNCH
C C  Rating = 192000.0   V-mult =  9.62000E-01   I-mult =  3.00000E+00   Gapless
C 92RECA                                 5555.
C C             V-reference              V-flashover
C    4.1250000000000000E+05   9.2121212121212130E-01
C C              Multiplier                 Exponent                    V-min
C    2.9479544296115170E+04   2.6530262418533470E+01   5.4505063612285050E-01
C                      9999
C    5.9605957177725360E+05   2.6530262418533470E+01   4.8665227581674000E-01
C                      9999
C 92RECB        RECA                     5555.
C 92RECC        RECA                     5555.
C       2nd of 3 subsets has manual selection of the 2 exponentials (NEXP = 2).
C       Otherwise, the configuration of interest is the same as for 1st subset.
C       Specifically, the i-v characteristic is the same (segmented i-v of 1st).
BRANCH  RECA        RECB        RECC        { 2nd subcase of 3 begins with names
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
           2           3                       1     412500.     380000.
C         A1          A2          A3          A4          A5        AMIN
     707.107     192000.        .962         3.0     .892857        .001
         1.0     1164.8
         2.0     1181.6
         5.0     1198.4
         10.     1209.6
         20.     1232.0
         50.     1260.0
        100.     1288.0
BLANK card bounds points of the ZnO characteristic - end of 1st segment, begin 2
        200.     1323.84
        500.     1388.8
       1000.     1442.56
       2000.     1512.0
       3000.     1556.8
BLANK card bounds points of the ZnO characteristic --- end 2nd and final segment
$PUNCH
C C  Rating = 192000.0   V-mult =  9.62000E-01   I-mult =  3.00000E+00   Gapless
C 92RECA                                 5555.
C C             V-reference              V-flashover
C    4.1250000000000000E+05   9.2121212121212130E-01
C C              Multiplier                 Exponent                    V-min
C    5.0558478867726780E+06   4.6419997332462390E+01   6.3275408479727320E-01
C    1.2276715303900460E+04   1.6677590344522540E+01   8.1674801890784120E-01
C                      9999
C    9.7386264053136350E+08   4.6419997332462390E+01   5.6495891388983900E-01
C    8.1268335203236510E+04   1.6677590344522540E+01   7.2923918591799820E-01
C                      9999
C 92RECB        RECA                     5555.
C 92RECC        RECA                     5555.
C                  3rd of 3 subsets illustrate automatic ZnO segment selection.
C                  For the tolerance chosen (ERRLIM=.05), 5 exponentials result.
BRANCH  RECA        RECB        RECC        { 3rd subcase of 3 beings with names
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
          -1           3         .05           1     412500.     380000.
C         A1          A2          A3          A4          A5        AMIN
     707.107     192000.        .962         3.0     .892857        .001
         1.0     1164.8
         2.0     1181.6
         5.0     1198.4
         10.     1209.6
         20.     1232.0
         50.     1260.0
        100.     1288.0
        200.     1323.84
        500.     1388.8
       1000.     1442.56
       2000.     1512.0
       3000.     1556.8
BLANK card bounds points of the ZnO characteristic
$PUNCH
C C  Rating = 192000.0   V-mult =  9.62000E-01   I-mult =  3.00000E+00   Gapless
C 92RECA                                 5555.
C C             V-reference              V-flashover
C    4.1250000000000000E+05   9.2121212121212130E-01
C C              Multiplier                 Exponent                    V-min
C    7.5028204133564620E+06   4.8403950953194400E+01   6.3949942977557550E-01
C    2.4489276643997050E+09   6.8408843374214610E+01   7.4876091788815410E-01
C    1.1046477599319610E+06   3.9475540404878120E+01   7.6623729035951440E-01
C    8.4663437752101640E+04   2.7934580347853540E+01   8.0046407236184350E-01
C    1.6489996944083740E+04   1.8756547209277270E+01   8.3673818423150000E-01
C    1.1130526844643490E+04   1.4446009330513530E+01   9.1284753595740940E-01
C                      9999
C    1.8095670102550080E+09   4.8403950953194400E+01   5.7098154237113100E-01
C    5.7007069310942390E+12   6.8408843374214610E+01   6.6853642686286360E-01
C    9.6857800318865540E+07   3.9475540404878120E+01   6.8414032835852490E-01
C    2.0071630292080190E+06   2.7934580347853540E+01   7.1469995025677850E-01
C    1.3815991082985590E+05   1.8756547209277270E+01   7.4708754495838430E-01
C    5.7216427799055640E+04   1.4446009330513530E+01   8.1504231241232490E-01
C                      9999
C 92RECB        RECA                     5555.
C 92RECC        RECA                     5555.
BLANK card ends data subcases within "ARRDAT"
BEGIN NEW DATA CASE
C      The preceding subase of this disk file involved 3 data sets using units
C      of volts and amps.  Now we repeat these three,  only using kilovolts and
C      kiloamps.  Output (including punched cards) is changed only by roundoff.
ZNO FITTER    { Special-request to transfer to  SUBROUTINE ARRDAT  of overlay 40
$ERASE
INPUT KV AND KA  { Input is to be in kilovolts and kiloamps (not volts and amps)
BRANCH  RECA        RECB        RECC    { Name the 3 terminal pairs for punching
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
           1           3           2           0       412.5        380.
C         A1          A2          A3          A4          A5        AMIN
     707.107        192.        .962         3.0     .892857        .001
        .001     1.1648
        .002     1.1816
        .005     1.1984
         .01     1.2096
         .02     1.2320
         .05     1.2600
         .10     1.2880
         .20     1.32384
         .50     1.3888
        1.00     1.44256
        2.00     1.5120
        3.00     1.5568
BLANK card bounds points of the ZnO characteristic
$PUNCH
C       2nd of 3 data subsets.  We do not need another  INPUT KV AND KA  declaration
C       as long as control remains within  ZNO FITTER.
BRANCH  RECA        RECB        RECC        { 2nd subcase of 3 begins with names
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
           2           3                       1       412.5        380.
C         A1          A2          A3          A4          A5        AMIN
     707.107        192.        .962         3.0     .892857        .001
        .001     1.1648
        .002     1.1816 
        .005     1.1984
         .01     1.2096
         .02     1.2320
         .05     1.2600
         .10     1.2880
BLANK card bounds points of the ZnO characteristic - end of 1st segment, begin 2
         .20     1.32384
         .50     1.3888
        1.00     1.44256
        2.00     1.5120
        3.00     1.5568
BLANK card bounds points of the ZnO characteristic --- end 2nd and final segment
$PUNCH
C       3rd of 3 data subsets of 2nd subcase.  The  INPUT KV AND KA  of the
C       first of the 3 data subsets remains in effect.  Answers should be
C       identical to the corresponding ones of the preceding subcase (which
C       involves no such kilo scaling).
BRANCH  RECA        RECB        RECC        { 3rd subcase of 3 beings with names
C       NEXP      IPHASE      ERRLIM      IPRZNO        VREF      VFLASH
          -1           3         .05           1       412.5        380.
C         A1          A2          A3          A4          A5        AMIN
     707.107        192.        .962         3.0     .892857        .001
        .001     1.1648
        .002     1.1816 
        .005     1.1984
         .01     1.2096
         .02     1.2320
         .05     1.2600
         .10     1.2880
         .20     1.32384
         .50     1.3888
        1.00     1.44256
        2.00     1.5120
        3.00     1.5568
BLANK card bounds points of the ZnO characteristic
$PUNCH
BLANK card ends data subcases within "ARRDAT"
BEGIN NEW DATA CASE
BLANK