summaryrefslogtreecommitdiffstats
path: root/benchmarks/dc57.dat
blob: d962ea9895cbe0d9288627a528d27c909afe6de2 (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
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
BEGIN NEW DATA CASE
C     BENCHMARK DC-57   (former DC-71 of "M39." vintage)
C     The basic  DC-30  problem has been augmented by batch-mode SPY commands as
C     an illustration of  $SPY  usage.  See  DC-56  (this is 2nd of 2 examples).
C     Unlike  DC-56,  this present case uses two sets of in-line  $SPY  data, so
C     it  is more general.    For SPY use,  the  user may desire conversion from
C     universal character plotting to computer-dependent vector-graphic plotting
C     (usually in a separate, third window).   Applicable  to  Apollo for such a
C     change  are the comment cards beginning with "C ---." (see two locations).
C       1st of 2 subcases is the original for Salford EMTP.  It involves one
C       vector plot,  and no  ROLL.  For ROLL, see the 2nd subcase.  Execution
C       with  DISK  will result in screen output that is limited to the vector
C       plots (1 from 1st subcase and 4 from 2nd).  Execution using  BOTH  will
C       differ in that text will be outputted to the screen until the 1st vector
C       plot.  Then, in effect, ATP switches from  BOTH  to  DISK.  March, 2002.
$DEPOSIT, D4FACT=2.0  { Use SPY DEPOSIT to change CALCOMP PLOT time from STARTUP
$DEPOSIT, NOCALC=-1    { We want screen graphics during DISK in spite of STARTUP
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 9 output columns)
  1.0E-6  1.2E-3
       1       1       1       0       1      -1
       5       5      20      20     100     100
TACS HYBRID
 1DUMMY   +UNITY
       1.0
       1.0    0.5E-3
90BUS2
90BUS3
99VSW                     BUS2  -     BUS3
99DRIVE              ABS  VSW
99BREAK       1.5E+8*     TIMEX             +     1.0E+5
98GRID  51+UNITY                                                    BREAK DRIVE
33DUMMY TIMEX UNITY BUS2  BUS3  VSW   DRIVE BREAK GRID
BLANK card ending the last TACS data
  GEN   BUS1                       15.
  BUS1                                   2.9
  BUS1  BUS2                       0.1
  BUS2                                   0.1
  BUS3                                  .017
  BUS3                            490.
  BUS2  BUS2R              24.34
  BUS3  BUS3R BUS2  BUS2R
BLANK card ending electric network branch cards
  BUS2  BUS3         -1.                1.E9                                   2
11BUS2R BUS3R                  20.                              GRID          13
$SPY         { Request for batch-mode SPY connection: in-line SPY commands begin
BRANCH          { SPY request for a display of the program branch table (List 2)
ALL                                       { Display all rows of the branch table
SPY        { Exit the loop of the branch table display;  back to  "SPY:"  prompt
SOURCE
TACS
ALL
SPY   { Exit the loop of the TACS source table display;  back to  "SPY:"  prompt
time   { Extra, do-nothing command makes the last something other than spy
$SPYEND        { Bound on in-line SPY commands;  back to batch-mode program data
BLANK card ending switch cards
14GEN         66500.       50.   -2.0508                           -1.
$SPY         { Request for batch-mode SPY connection: in-line SPY commands begin
SOURCE
ELEC
ALL
SPY  { Exit the loop of the electric source table display; back to "SPY:" prompt
SWITCH     { Show ATP storage related to switches
ALL        { Display all rows of switch table
SPY        { Exit the loop of the switch table display;  back to  "SPY:"  prompt
BREAK      { Initiate a request for a later SPY break
16, -750    { Next SPY break will be at .75E-3 seconds, which is step number 750
GO             { Begin EMTP execution;  no "@" read until break point is reached
PLOT                        { Having just broken at step 750, enter SPY plotting
C ---.    For 3-window SPY usage and vector plotting,  comment out the MODE card
C ---.    that follows (assuming LTEK=1 for vector plots is set by STARTUP file)
C MODE         { Toggle from vector (Apollo default) to character mode of plotting
SET COLUMN                       { Request a change in the output plotting width
131   { The default width of 79 is now being changed to 131 for maximum accuracy
CHOICE          { Request a list of the output variables, available for plotting
NAME                        { Plot command to begin inputting the variable names
END               { Step over node voltage input, and begin branch voltage input
BUS2  BUS3 { The 1st and only variable to be plotted will be this branch voltage
LAST      { Exit the loop over plot input variables;  back to basic plot command
LABEL         { Define labeling for plots,  beginning with super title (next):
Super Title Line
Vertical axis
1st of two lines of case title
2nd and final line of case title
end
TIMEUNITS  { Default time units of seconds are not convenient, so change to:
4          { 4 ==> msec (scale by 1000 in labeling the time axis)
ALL TIME                      { Perform the plot using all available time points
SPY                                   { Exit plotting;  back to  "SPY:"  command
BREAK     { Request a pause in the simulation that will soon resume (after "GO")
16, -1001           { Next SPY activity will recommence at time step number 1001
GO             { Begin EMTP execution;  no "@" read until break point is reached
TIME       { Trivial SPY command before program stop, just to prove SPY activity
$SPYEND        { Bound on in-line SPY commands;  back to batch-mode program data
BLANK card ending all electric sources
C       Total network loss  P-loss  by summing injections =  -1.885346136987E-07
C      Step      Time       BUS2        BUS2R       BUS3R       BUS2R       BUS3
C                           BUS3        BUS3R
C
C                           TACS        TACS        TACS        TACS        TACS
C                           DUMMY       TIMEX       UNITY       BUS2        BUS3
C ***  Phasor I(0) = -1.5049840E+01       Switch  "BUS2  "  to  "BUS3  "  closed
C       0       0.0          0.0         0.0  64751.3498  64751.3498  64751.3498
C                            0.0         0.0         1.0         0.0         0.0
C ***          Open  switch  "BUS2  "  to  "BUS3  "  after  1.00000000E-06  sec.
C       1     .1E-5          0.0         0.0  64752.0751  64752.0751  64752.0751
C                     .001998002       .1E-5         1.0  64752.0751  64752.0751
C $OPEN,  UNIT=42  STATUS=UNKNOWN  FILE=DEBUG.LIS
 1                         { Request for the output of all network node voltages
C ---.        All following cards will only be read for single-window execution:
C Last step: 1200  .0012  105797.943 105797.943 -37370.689 68427.2546 -37370.689
C Last step:              .909282119      .0012       1.0  68427.2546 -37370.689
C maxima : 177392.698 177392.698  102764.956  84422.9848  102764.956  84422.9848
C maxima : .909282119      .0012         1.0  84422.9848  102764.956  177392.698
C times max : .497E-3    .497E-3      .83E-4     .478E-3      .83E-4     .478E-3
C times max :   .0012      .0012         0.0     .478E-3      .83E-4     .497E-3
  PRINTER PLOT
 184 .1 0.0 1.2         BUS2  BUS3               { Axis limits: (-0.386,  1.774)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C        2nd of 2 subcases is appended  18 March 2002.  This is modification of
C        the first.  A big change was made for Salford, as Salford VECPLT in
C        the UTPF was reformed.  Previously,  SPY PLOT  worked,  but it was not
C        a ROLL-ing vector plot as Mingw32 ATP demonstrated in its separate plot
C        window for DISLIN graphics.  Now, Salford EMTP, too, can have a ROLLing
C        vector  SPY PLOT  just as Mingw32 ATP did.  There is a price, however:
C        Execution must be using  DISK  rather than  BOTH  (text output to the
C        screen is not possible).  Also,  NOCALC = -1  is required to enable
C        CALCOMP PLOT  graphics  (screen graphics)  during  DISK.  Then, during
C        all plotting,  the screen is held in graphic mode.   Previously, TSTALL
C        was made positive to slow text output.  Since text and graphics are
C        mixed for the ROLL-ing SPY PLOT  in the dT loop, this still works.  The
C        graphics continue to be slowed,  even though text can not be seen on
C        the screen.  Using default 50% overlap of plots, the first begins at
C        1/4 msec,  and extends through 1.0 msec.  The 2nd then is (0.5, 1.5),
C        and the 3rd covers (1.0, 2.0),  of which the simulation covers 80% when
C        T-max is reached.  Finally, real  CALCOMP PLOT  produces a 4th plot.
C        All 4 should be in the  .PS  (PostScript) output,  and will be found in
C        the  ATPHPGL.*  for  HP-GL  (one disk file per plot).
$DEPOSIT, KSLOWR=5    { Use SPY DEPOSIT to change ROLL buffer from STARTUP value
PRINTED NUMBER WIDTH, 12, 2,  { Request maximum precision (for 9 output columns)
  1.0E-6  1.8E-3
       1       1       1       0       1      -1
       5       5      20      20     100     100
TACS HYBRID
 1DUMMY   +UNITY
       1.0
       1.0    0.5E-3
90BUS2
90BUS3
99VSW                     BUS2  -     BUS3
99DRIVE              ABS  VSW
99BREAK       1.5E+8*     TIMEX             +     1.0E+5
98GRID  51+UNITY                                                    BREAK DRIVE
33DUMMY TIMEX UNITY BUS2  BUS3  VSW   DRIVE BREAK GRID
BLANK card ending the last TACS data
  GEN   BUS1                       15.
  BUS1                                   2.9
  BUS1  BUS2                       0.1
  BUS2                                   0.1
  BUS3                                  .017
  BUS3                            490.
  BUS2  BUS2R              24.34
  BUS3  BUS3R BUS2  BUS2R
BLANK card ending electric network branch cards
  BUS2  BUS3         -1.                1.E9                                   2
11BUS2R BUS3R                  20.                              GRID          13
$SPY         { Request for batch-mode SPY connection: in-line SPY commands begin
BRANCH          { SPY request for a display of the program branch table (List 2)
ALL                                       { Display all rows of the branch table
SPY        { Exit the loop of the branch table display;  back to  "SPY:"  prompt
SOURCE
TACS
ALL
SPY   { Exit the loop of the TACS source table display;  back to  "SPY:"  prompt
time   { Extra, do-nothing command makes the last something other than spy
$SPYEND        { Bound on in-line SPY commands;  back to batch-mode program data
BLANK card ending switch cards
14GEN         66500.       50.   -2.0508                           -1.
$SPY         { Request for batch-mode SPY connection: in-line SPY commands begin
SOURCE
ELEC
ALL
SPY  { Exit the loop of the electric source table display; back to "SPY:" prompt
SWITCH
ALL
SPY        { Exit the loop of the switch table display;  back to  "SPY:"  prompt
BREAK      { Initiate a request for a later SPY break
16, -250    { Next SPY break will be at .25E-3 seconds, which is step number 250
GO             { Begin EMTP execution;  no "@" read until break point is reached
PLOT                        { Having just broken at step 750, enter SPY plotting
C ---.    For 3-window SPY usage and vector plotting,  comment out the MODE card
C ---.    that follows (assuming LTEK=1 for vector plots is set by STARTUP file)
C MODE         { Toggle from vector (Apollo default) to character mode of plotting
SET COLUMN                       { Request a change in the output plotting width
131   { The default width of 79 is now being changed to 131 for maximum accuracy
CHOICE          { Request a list of the output variables, available for plotting
NAME                        { Plot command to begin inputting the variable names
END               { Step over node voltage input, and begin branch voltage input
BUS2  BUS3 { The 1st and only variable to be plotted will be this branch voltage
LAST      { Exit the loop over plot input variables;  back to basic plot command
LABEL         { Define labeling for plots,  beginning with super title (next):
Super Title Line
Vertical label
FLUSH    { Erase any multi-line case title now being stored (preceding subcase)
1st of two lines of case title
2nd and final line of case title
END           { Terminate this multi-line case title
TIMEUNITS  { Default time units of seconds are not convenient, so change to:
4          { 4 ==> msec (scale by 1000 in labeling the time axis)
TIME 0.0 1.0    { Plot 4 times as far as points now exist (0, 1/4 msec)
ROLLV           { Switch to ROLL-ing plot, Vector mode (not C ==> character)
SPY                                   { Exit plotting;  back to  "SPY:"  command
BREAK     { Request a pause in the simulation that will soon resume (after "GO")
16, -1001           { Next SPY activity will recommence at time step number 1001
GO             { Begin EMTP execution;  no "@" read until break point is reached
TIME       { Trivial SPY command before program stop, just to prove SPY activity
$SPYEND        { Bound on in-line SPY commands;  back to batch-mode program data
BLANK card ending all electric sources
C    The following slows output just prior to entry into dT loop (and plotting).
C    This serves to hold the 1st plot on the screen for several seconds, & also
C    to slow following ROLL-ing SPY PLOTs of this subcase enough so the user is
C    able to study structure.  Without TSTALL, the average PC is too fast, so:
$DEPOSIT, TSTALL=0.1      { 1/10 second delay will occur after each printed line
 1                         { Request for the output of all network node voltages
  CALCOMP PLOT  { Switch to vector plotting (preceding subcase had PRINTER PLOT)
 2Arbitrary 78-character case title text of which this is an example, I hope.
  First of two lines of 78-byte graph subheading text.
  Second and final such line of graph subheading text.
 184 .2 0.0 2.0         BUS2  BUS3              Graph heading---Vertical axis la
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK