diff options
Diffstat (limited to 'benchmarks/dc18.dat')
| -rw-r--r-- | benchmarks/dc18.dat | 384 |
1 files changed, 384 insertions, 0 deletions
diff --git a/benchmarks/dc18.dat b/benchmarks/dc18.dat new file mode 100644 index 0000000..f833ce5 --- /dev/null +++ b/benchmarks/dc18.dat @@ -0,0 +1,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
|