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
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
|
BEGIN NEW DATA CASE
C BENCHMARK DC-18
C Artificially concocted sample test of various logical operators of TACS.
C There is no electric network. Note the unusual use of the TSTART field
C (columns 61-70) and TSTOP (columns 71-80) of the TACS source cards,
C which turn the constant source amplitudes on and off during the several
C steps of this test. A study of the operations involved will show that
C the variables RESL1 and RESL2 are residuals which should always
C be zero by logical identity. See DCPR24.DAT for DIAGNOSTIC version.
C For a DIAGNOSTIC version of the very same data case, refer to DCPR-24.
0.1 2.0
1 1 0 0 1
TACS STAND ALONE
DUMMY +UNITY
11LGCL1 1. 0.25
11FST 1. 0.05 0.15
11SCND 1. 0.25 0.35
99LGCL4 = FST + SCND
C 99LGCL4 FST + SCND
99LGCL2 .NOT. LGCL4
98NAND LGCL1 .NAND.LGCL2
98RESL1 .NOT. LGCL1 .OR. LGCL4 - NAND
98NOR LGCL1 .NOR. LGCL2
98RESL2 .NOT. LGCL1 .AND. LGCL4 - NOR
33LGCL1 LGCL2 LGCL4 NAND NOR RESL1 RESL2 FST SCND
C Step Time TACS TACS TACS TACS TACS TACS TACS TACS TACS
C LGCL1 LGCL2 LGCL4 NAND NOR RESL1 RESL2 FST SCND
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 1 0.1 0.0 0.0 1.0 1.0 1.0 0.0 0.0 1.0 0.0
C 2 0.2 0.0 1.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0
C 3 0.3 1.0 0.0 1.0 1.0 0.0 0.0 0.0 0.0 1.0
C 4 0.4 1.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
BLANK card ending all TACS data cards
C 20 2.0 1.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C Variable maxima : 1.0 1.0 1.0 1.0 1.0 0.0 0.0 1.0 1.0
C Times of maxima : 0.3 0.2 0.1 0.1 0.1 0.0 0.0 0.1 0.3
C Variable minima : 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C Times of minima : 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
PRINTER PLOT
143 .2 0.0 1.0 LGCL1 NOR { Axis limits: (0.000, 1.000)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 2nd subcase of DC-18 tests various free-format supplemental variables.
C There is no connection between this and the original test case data.
C Algebraic and trigonometric computation is done 2 ways: once in high-
C level expressions (TEST1, TEST2, TEST3), and a second time using a
C reordered low-level equivalents (X1, .... X11). The two are compared
C (residual RESID should in theory be zero). As for the following
C allocation of TACS tables, this is just a modification of the default
C sizes that were overflowed: 20 90 100 20 30 250 350 60 4
C The complaint was about Table 6, so this was expanded at expense of 1.
C 11 March 2001, TACS table number 9 was added for Type-53 devices.
C This data involves none, so value is immaterial. Default is shown.
TACS WARN LIMIT, 5, 1.0, { No use other than to illustrate the declaration
ABSOLUTE TACS DIMENSIONS { Table #: 5 6 7 8 9
C 10 90 100 20 30 400 350 60 4
C Expand TACS Table 1 from 57 to 130 on 1 April 2007. Copy use from DC-2
C without worrying about probably waste that might be involved:
C 57 256 285 36 85 713 998 171 --- default
130 65 80 20 70 325 120 115 { From DC2.DAT
.02 2.0 { Even though no dynamics, vary time to produce nice plots!
1 1 0 0 1 -1
5 5
TACS STAND ALONE
99TEST1 = 10.0 * ( UNITY + TIMEX ) ** 2 + 50. { High-level, single-card result
99X1 = 1.0 + TIMEX { Expression within parentheses on preceding card
99X2 = X1 * X1 { Mimic exponentiation for special case of square
99X3 = 10 * X2 + 50 { 3rd of low-level, 3-card equivalent to TEST1
99TEST2 = 1.E2 * COS ( 2.0 * PI * TEST1 / 100. ) { High-level d1-card result
99X4 = PI * TEST1 { Beginning terms of trigonometric argument
99X5 = X4 * .02 { Complete trig argument using modified form of 2/100
99X6 = COS ( X5 )
99X7 = X6 * 100. { 4th and final line of low-level equivalent gives TEST2
99TEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( TEST2 ) )
99X8 = ABS ( X7 )
99X9 = SQRT ( X8 )
99X10 = X9 * 5.5 * TIMEX
99X11 = X10 + 10 { 4th and final line of low-level equivalent gives TEST3
99RESID = ABS ( TEST1 - X3 ) + ABS ( TEST2 - X7 ) + ABS ( TEST3 - X11 )
99NORMAL = GAUSS ( 0.0 ) { Illustrate Gaussian random # generator. October 2010
1TEST4 +TEST3
1.0
0.0 1.0
99TEST5 = 2.0 * PI * TEST1 / 100.
99TEST6 = SIN ( TEST5 )
99UNITY? = 1.E-4 * ( TEST2 * TEST2 + 1.E4 * TEST6 * TEST6 )
33TEST1 X3 TEST2 X7 TEST3 X11 TEST4 UNITY?RESID NORMAL
C Column headings for the 10 EMTP output variables follow. These are divided among the 5 possible classes as follows ....
C Next 10 output variables belong to TACS (with "TACS" an internally-added upper name of pair).
C Step Time TACS TACS TACS TACS TACS TACS TACS TACS TACS TACS
C TEST1 X3 TEST2 X7 TEST3 X11 TEST4 UNITY? RESID NORMAL
C 0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
C 1 .02 60.404 60.404 -79.383758 -79.383758 10.9800732 10.9800732 .109800732 1.0 0.0 .890624441
C 2 .04 60.816 60.816 -77.783088 -77.783088 11.9402839 11.9402839 .339004303 1.0 .30198E-13 -.26382016
BLANK card ending all TACS data cards
C 95 1.9 134.1 134.1 -54.112125 -54.112125 86.8711868 86.8711868 88.6795836 1.0 0.0 -1.4747572
C 100 2.0 140. 140. -80.901699 -80.901699 108.939909 108.939909 98.5281276 1.0 0.0 -.47311878
C Variable maxima : 140. 140. 99.9938856 99.9938856 108.939909 108.939909 98.5281276 1.0 .33396E-12 2.73027422
C Times of maxima : 2.0 2.0 1.24 1.24 2.0 2.0 2.0 .08 1.82 0.3
C Variable minima : 0.0 0.0 -80.901699 -80.901699 0.0 0.0 0.0 0.0 0.0 -2.7809904
C Times of minima : 0.0 0.0 2.0 2.0 0.0 0.0 0.0 0.0 0.0 1.7
CALCOMP PLOT { Needed for 2nd stacked case; 1st ended with PRINTER PLOT use
143 .2 0.0 2.0 TEST1 TEST2 TEST3 TEST4
PRINTER PLOT
143 .4 0.0 2.0 TEST1 TEST2 TEST4 { Axis limits: (-0.809, 1.400)
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 3rd of 5 subcases is identical to 2nd except that here the pocket
C calculator is used in place of Dube's logic for supplemental variables
C of TACS. Answers should be identical except for roundoff error.
C The plots will be omitted for simplicity.
C TACS ASSEMBLY LANGUAGE { Temporary request for use of pocket calculator
TACS POCKET CALCULATOR { 12 January 2001, this new line replaces preceding
TACS POCKET CALCULATOR ON { 1 February 2002, add this equivalent declaration
C The two preceding declarations have exactly the same effect. Either or both
C can be repeated any number of times. Either is simply a request to use the
C pocket calculator rather than Dube's logic. Making the request more than
C once has no effect. The ON alternative became effective 1 February 2002.
C 19 March 2009, add request to prohibit optimization by the pocket calculator.
C Addition here is purely illustrative. The answer changes only slightly, with
C near-zero RESID changing to exact zeros on two time steps (for Salford ATP).
C This is variable (TACS, RESID) which had nonzero values .14211E-13 and
C .14211E-13 for steps 60 and 70 when there was optimization. WSM.
$PARAMETER NOOPT=1 { Set optimization of pocket calculator compilation to none
ABSOLUTE TACS DIMENSIONS
C 10 90 100 20 30 400 350 60
C Expand TACS Table 1 from 57 to 130 on 1 April 2007. Copy use from DC-2
C without worrying about probably waste that might be involved:
C 57 256 285 36 85 713 998 171 --- default
130 65 80 20 70 325 120 115 { From DC2.DAT
.02 2.0 { Even though no dynamics, vary time to produce nice plots!
1 1 0 0 1 -1
5 5
TACS STAND ALONE
99TEST1 = 10.0 * ( UNITY + TIMEX ) ** 2 + 50. { High-level, single-card result
99X1 = 1.0 + TIMEX { Expression within parentheses on preceding card
99X2 = X1 * X1 { Mimic exponentiation for special case of square
99X3 = 10 * X2 + 50 { 3rd of low-level, 3-card equivalent to TEST1
99TEST2 = 1.E2 * COS ( 2.0 * PI * TEST1 / 100. ) { High-level d1-card result
99X4 = PI * TEST1 { Beginning terms of trigonometric argument
99X5 = X4 * .02 { Complete trig argument using modified form of 2/100
99X6 = COS ( X5 )
99X7 = X6 * 100. { 4th and final line of low-level equivalent gives TEST2
99TEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( TEST2 ) )
99X8 = ABS ( X7 )
99X9 = SQRT ( X8 )
99X10 = X9 * 5.5 * TIMEX
99X11 = X10 + 10 { 4th and final line of low-level equivalent gives TEST3
99RESID = ABS ( TEST1 - X3 ) + ABS ( TEST2 - X7 ) + ABS ( TEST3 - X11 )
1TEST4 +TEST3
1.0
0.0 1.0
99TEST5 = 2.0 * PI * TEST1 / 100.
99TEST6 = SIN ( TEST5 )
99UNITY? = 1.E-4 * ( TEST2 * TEST2 + 1.E4 * TEST6 * TEST6 )
33TEST1 X3 TEST2 X7 TEST3 X11 TEST4 UNITY?RESID
BLANK card ending all TACS data cards
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 4th of 5 subcases documents the data of MATHCOMP as first mentioned
C in the April, 1997 newsletter. The data shown here demonstrates use
C of the pocket calculator rather than compiled TACS, however. This is
C the data that was used to estimate improved speed compared with Dube's
C code for handling supplemental variables (see July, 2001, newsletter).
C TACS POCKET CALCULATOR { Use pocket calculator rather than Dube's logic (default)
C 2 February 2002, remove the preceding line by commenting. This has no
C effect on the answer. In fact, the TPC declaration of the 3rd subcase will
C remain in effect. The TPC declaration has an effect here if and only if no
C preceding subcase had one. Any TPC declaration remains in effect for all
C later data unless cancelled by a TPC OFF declaration (see DCN21.DAT).
C $DEPOSIT, KOMPAR=1 { Use SPY DEPOSIT to change STARTUP value so elapsed times
C For benchmark purposes, enable the preceding statement. Regardless
C of the value of KOMPAR in STARTUP, this will gives elapsed times at
C the end of execution.
C Add NLS (following 7 noncomment cards) on 8 March 2003 to prove it works.
C In fact, dimensions are unchanged from LISTSIZE.333 as used for RUN.BAT
NEW LIST SIZES { Request for a change to program table sizes (VARDIM input)
DEFAULT 3.0 { Any blank field will be converted to 3 times default value
BLANK card for lists 1-10.
BLANK card for lists 11-20.
BLANK card for lists 21-30.
200 300 { Optional card for Lists 31, etc.
240000 742 { Offsets for supporting programs (non-simulation overlays)
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
ABSOLUTE TACS DIMENSIONS
10 90 100 20 60 1600 350 120
C .00002 2.0 { Original line gives 100K steps of July, 2001, newsletter
.002 2.0 { For standard test cases, speed execution by factor of 100
1 -11 0 0 1 -1
C 5 5 20 20 100 100 1000 1000 10000 10000
5 5 20 20 100 100 { Preceding was for 100K steps
TACS STAND ALONE
99TEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99TEST2 = 1.E2 * COS ( 2.0 * 3.14159 * TEST1 / 100. )
99TEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( TEST2 ) )
99TEST5 = 2.0 * PI * TEST1 / 100.
99TEST6 = 2.0 * SIN ( TEST5 ) { Show that in-line comments are allowed
99TWOA = 2.E-4 * ( TEST2 * TEST2 + 2500. * TEST6 * TEST6 )
99BEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99BEST2 = 1.E2 * COS ( 2.0 * 3.14159 * BEST1 / 100. )
99BEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( BEST2 ) )
99BEST5 = 2.0 * PI * BEST1 / 100.
99BEST6 = 2.0 * SIN ( BEST5 )
99TWOB = 2.E-4 * ( BEST2 * BEST2 + 2500. * BEST6 * BEST6 )
99ZEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99ZEST2 = 1.E2 * COS ( 2.0 * 3.14159 * ZEST1 / 100. )
99ZEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( ZEST2 ) )
99ZEST5 = 2.0 * PI * ZEST1 / 100.
99ZEST6 = 2.0 * SIN ( ZEST5 )
99TWOC = 2.E-4 * ( ZEST2 * ZEST2 + 2500. * ZEST6 * ZEST6 )
99REST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99REST2 = 1.E2 * COS ( 2.0 * 3.14159 * REST1 / 100. )
99REST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( REST2 ) )
99REST5 = 2.0 * PI * REST1 / 100.
99REST6 = 2.0 * SIN ( REST5 )
99TWOD = 2.E-4 * ( REST2 * REST2 + 2500. * REST6 * REST6 )
99DEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99DEST2 = 1.E2 * COS ( 2.0 * 3.14159 * DEST1 / 100. )
99DEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( DEST2 ) )
99DEST5 = 2.0 * PI * DEST1 / 100.
99DEST6 = 2.0 * SIN ( DEST5 )
99TWOE = 2.E-4 * ( DEST2 * DEST2 + 2500. * DEST6 * DEST6 )
99FEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99FEST2 = 1.E2 * COS ( 2.0 * 3.14159 * FEST1 / 100. )
99FEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( FEST2 ) )
99FEST5 = 2.0 * PI * FEST1 / 100.
99FEST6 = 2.0 * SIN ( FEST5 )
99TWOF = 2.E-4 * ( FEST2 * FEST2 + 2500. * FEST6 * FEST6 )
99GEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99GEST2 = 1.E2 * COS ( 2.0 * 3.14159 * GEST1 / 100. )
99GEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( GEST2 ) )
99GEST5 = 2.0 * PI * GEST1 / 100.
99GEST6 = 2.0 * SIN ( GEST5 )
99TWOG = 2.E-4 * ( GEST2 * GEST2 + 2500. * GEST6 * GEST6 )
99HEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99HEST2 = 1.E2 * COS ( 2.0 * 3.14159 * HEST1 / 100. )
99HEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( HEST2 ) )
99HEST5 = 2.0 * PI * HEST1 / 100.
99HEST6 = 2.0 * SIN ( HEST5 )
99TWOH = 2.E-4 * ( HEST2 * HEST2 + 2500. * HEST6 * HEST6 )
99LEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99LEST2 = 1.E2 * COS ( 2.0 * 3.14159 * LEST1 / 100. )
99LEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( LEST2 ) )
99LEST5 = 2.0 * PI * LEST1 / 100.
99LEST6 = 2.0 * SIN ( LEST5 )
99TWOI = 2.E-4 * ( LEST2 * LEST2 + 2500. * LEST6 * LEST6 )
99MEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99MEST2 = 1.E2 * COS ( 2.0 * 3.14159 * MEST1 / 100. )
99MEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( MEST2 ) )
99MEST5 = 2.0 * PI * MEST1 / 100.
99MEST6 = 2.0 * SIN ( MEST5 )
99TWOJ = 2.E-4 * ( MEST2 * MEST2 + 2500. * MEST6 * MEST6 )
33MEST1 MEST2 MEST3 MEST5 MEST6 TWOJ { Outputs are for 10th of 10 data sets
BLANK card ending all TACS data cards
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 5th of 5 subcases documents the use of IF blocks within TACS. This
C began 14 March 2001 when most-primitive IF-THEN-ELSE-ENDIF first
C produced correct answers. The pocket calculator is required (do not
C try to execute after removing the TPC declaration below).
C DIAGNOSTIC 0 9 0 0 0 0 0 0 0 0 0 0 0 0 0 9 { Make comment 22 Apr 07
TACS POCKET CALCULATOR { Use pocket calculator for TACS supplemental variables
.02 .14
1 -1
TACS STAND ALONE
C Begin with simplest of block IF statements: IF-THEN-ELSE-ENDIF:
IF( TIMEX .LE. 0.08 ) THEN { If simulation time T does not exceed 4 * dT:
99X1 = 1.0 + TIMEX { Formula for steps 1 thru 4 will produce a ramp up
ELSE { Alternatively (if simulation time T does exceed 4 * dT):
99X1 = 1.16 - TIMEX { Formula for steps 5, 6, and 7 will produce ramp down
ENDIF { Terminate 5-line block that chooses among 2 formulas for supplemental X1
C As 2nd illustration, add ELSEIF to preceding illustration, and
C make the resulting signal more interesting. First, ramp it up. Then,
C in the middle, will be a flat top. Finally, ramp it back down to
C where it began. This enhancement became operational March 16th:
IF( TIMEX .LE. 0.04 ) THEN { If simulation time T does not exceed 2 * dT:
99X2 = 1.0 + TIMEX { For steps 1 and 2, assign values 1.02 and 1.04
ELSEIF( TIMEX .LE. 0.08 ) THEN { If simulation time T is either 3*dT or 4*dT:
99X2 = 1.06 { For steps 3 and 4, assign flat top value of 1.06
ELSE { Alternatively (if simulation time T does exceed 4 * dT):
99X2 = 1.06 - ( TIMEX - 0.08 ) { For steps 5 thru 7: 1.04, 1.02, and 1.0
ENDIF { Terminate 7-line block that chooses among 3 formulas for supplemental X2
28SOUR28 DEV69D { Cols. 41-46 carries user name
C Add the preceding Type-28 TACS source, which is defined in user-supplied TACS
C source code of SUBROUTINE DEVT69. This is disconnected from rest of data.
C Addition here on 30 January 2003 simply confirms function DEV69D of the UTPF.
C Since the Type-28 source is only defined within the dT loop, step 0 will
C show value zero. If the user is unhappy with this, the following initial
C condition (now commented out) will provide continuity:
C 77SOUR28 .75 { Initial condition on the Type-28 function avoids zero
1TEST4 +X2
1.0
0.0 1.0
33X1 X2 TEST4 SOUR28
BLANK card ending all TACS data cards
BLANK card ending plot cards
BEGIN NEW DATA CASE
C 6th of 6 subcases documents is the same as the 4th except that List 32
C is made too small for the data. Added 9 March 2003, this demonstrates
C rejection if the new variable limit LSIZ32 is inadequate -- but only
C for F95 (Lahey). For F77, the 20 below is ignored. For F77, the size
C of tables of the pocket calculator remain fixed at 200. For F77, this
C data demonstrates that any attempt to restrict List 31 is ignored. To
C simplify, drop extrema and limit output to T-max = 10 * dT.
NEW LIST SIZES { Request for a change to program table sizes (VARDIM input)
DEFAULT 3.0 { Any blank field will be converted to 3 times default value
BLANK card for lists 1-10.
BLANK card for lists 11-20.
BLANK card for lists 21-30.
200 20 { Optional card for Lists 31, etc.
240000 742 { Offsets for supporting programs (non-simulation overlays)
PRINTED NUMBER WIDTH, 10, 2, { Request maximum precision (for 8 output columns)
ABSOLUTE TACS DIMENSIONS
10 90 100 20 60 1600 350 120
C .00002 2.0 { Original line gives 100K steps of July, 2001, newsletter
.002 .02 { 10 steps is more than enough to demonstrate LSIZ32 limit
1 -11 0 0 0 -1
5 5
TACS STAND ALONE
99TEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99TEST2 = 1.E2 * COS ( 2.0 * 3.14159 * TEST1 / 100. )
99TEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( TEST2 ) )
99TEST5 = 2.0 * PI * TEST1 / 100.
99TEST6 = 2.0 * SIN ( TEST5 ) { Show that in-line comments are allowed
99TWOA = 2.E-4 * ( TEST2 * TEST2 + 2500. * TEST6 * TEST6 )
99BEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99BEST2 = 1.E2 * COS ( 2.0 * 3.14159 * BEST1 / 100. )
99BEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( BEST2 ) )
99BEST5 = 2.0 * PI * BEST1 / 100.
99BEST6 = 2.0 * SIN ( BEST5 )
99TWOB = 2.E-4 * ( BEST2 * BEST2 + 2500. * BEST6 * BEST6 )
99ZEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99ZEST2 = 1.E2 * COS ( 2.0 * 3.14159 * ZEST1 / 100. )
99ZEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( ZEST2 ) )
99ZEST5 = 2.0 * PI * ZEST1 / 100.
99ZEST6 = 2.0 * SIN ( ZEST5 )
99TWOC = 2.E-4 * ( ZEST2 * ZEST2 + 2500. * ZEST6 * ZEST6 )
99REST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99REST2 = 1.E2 * COS ( 2.0 * 3.14159 * REST1 / 100. )
99REST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( REST2 ) )
99REST5 = 2.0 * PI * REST1 / 100.
99REST6 = 2.0 * SIN ( REST5 )
99TWOD = 2.E-4 * ( REST2 * REST2 + 2500. * REST6 * REST6 )
99DEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99DEST2 = 1.E2 * COS ( 2.0 * 3.14159 * DEST1 / 100. )
99DEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( DEST2 ) )
99DEST5 = 2.0 * PI * DEST1 / 100.
99DEST6 = 2.0 * SIN ( DEST5 )
99TWOE = 2.E-4 * ( DEST2 * DEST2 + 2500. * DEST6 * DEST6 )
99FEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99FEST2 = 1.E2 * COS ( 2.0 * 3.14159 * FEST1 / 100. )
99FEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( FEST2 ) )
99FEST5 = 2.0 * PI * FEST1 / 100.
99FEST6 = 2.0 * SIN ( FEST5 )
99TWOF = 2.E-4 * ( FEST2 * FEST2 + 2500. * FEST6 * FEST6 )
99GEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99GEST2 = 1.E2 * COS ( 2.0 * 3.14159 * GEST1 / 100. )
99GEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( GEST2 ) )
99GEST5 = 2.0 * PI * GEST1 / 100.
99GEST6 = 2.0 * SIN ( GEST5 )
99TWOG = 2.E-4 * ( GEST2 * GEST2 + 2500. * GEST6 * GEST6 )
99HEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99HEST2 = 1.E2 * COS ( 2.0 * 3.14159 * HEST1 / 100. )
99HEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( HEST2 ) )
99HEST5 = 2.0 * PI * HEST1 / 100.
99HEST6 = 2.0 * SIN ( HEST5 )
99TWOH = 2.E-4 * ( HEST2 * HEST2 + 2500. * HEST6 * HEST6 )
99LEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99LEST2 = 1.E2 * COS ( 2.0 * 3.14159 * LEST1 / 100. )
99LEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( LEST2 ) )
99LEST5 = 2.0 * PI * LEST1 / 100.
99LEST6 = 2.0 * SIN ( LEST5 )
99TWOI = 2.E-4 * ( LEST2 * LEST2 + 2500. * LEST6 * LEST6 )
99MEST1 = 10.0 * ( 1.0 + TIMEX ) ** 2 + 50.
99MEST2 = 1.E2 * COS ( 2.0 * 3.14159 * MEST1 / 100. )
99MEST3 = 10.0 + 5.5 * TIMEX * SQRT ( ABS ( MEST2 ) )
99MEST5 = 2.0 * PI * MEST1 / 100.
99MEST6 = 2.0 * SIN ( MEST5 )
99TWOJ = 2.E-4 * ( MEST2 * MEST2 + 2500. * MEST6 * MEST6 )
33MEST1 MEST2 MEST3 MEST5 MEST6 TWOJ { Outputs are for 10th of 10 data sets
BLANK card ending all TACS data cards
BLANK card ending plot cards
BEGIN NEW DATA CASE
BLANK
|
