From b18347ffc9db9641e215995edea1c04c363b2bdf Mon Sep 17 00:00:00 2001 From: Angelo Rossi Date: Wed, 21 Jun 2023 12:04:16 +0000 Subject: Initial commit. --- benchmarks/dc18.dat | 384 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 384 insertions(+) create mode 100644 benchmarks/dc18.dat (limited to 'benchmarks/dc18.dat') 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 -- cgit v1.2.3