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
|
SUBROUTINE TCNSMT (IERROR)
C
C LATEST REVISION JUNE 1984
C
C PURPOSE TO PROVIDE A SIMPLE DEMONSTRATION OF
C CONRECSMTH AND TO TEST CONRECSMTH ON A SINGLE
C PROBLEM
C
C USAGE CALL TCNSMT (IERROR)
C
C ARGUMENTS
C
C ON OUTPUT IERROR
C AN INTEGER VARIABLE
C = 0, IF THE TEST WAS SUCCESSFUL,
C = 1, OTHERWISE
C
C I/O IF THE TEST IS SUCCESSFUL, THE MESSAGE
C
C CONRECSMTH TEST SUCCESSFUL . . . SEE PLOT TO
C VERIFY PERFORMANCE
C
C IS PRINTED ON UNIT 6.
C IN ADDITION, TWO FRAMES CONTAINING THE CONTOUR
C PLOT ARE PRODUCED ON THE MACHINE GRAPHICS
C DEVICE. IN ORDER TO DETERMINE IF THE TEST
C WAS SUCCESSFUL, IT IS NECESSARY TO EXAMINE
C THESE PLOTS.
C
C PRECISION SINGLE
C
C
C LANGUAGE FORTRAN
C
C ALGORITHM THE FUNCTION
C Z(X,Y) = X + Y + 1./((X-.1)**2+Y**2+.09)
C -1./((X+.1)**2+Y**2+.09)
C FOR X = -1. TO +1. IN INCREMENTS OF .1 AND
C Y = -1.2 TO +1.2 IN INCREMENTS OF .1
C IS COMPUTED.
C TCNSMT CALLS SUBROUTINES EZCNTR, CONREC, AND
C WTSTR TO DRAW TWO LABELLED CONTOUR PLOTS OF THE
C ARRAY Z.
C
C PORTABILITY ANSI FORTRAN77 STANDARD
C
C
C Z CONTAINS THE VALUES TO BE PLOTTED.
C
REAL Z(21,25)
C
C SPECIFY COORDINATES FOR PLOT TITLES. ON AN ABSTRACT GRID WHERE
C THE INTEGER COORDINATES RANGE FROM 0.0 TO 1.0, THE VALUES TX AND TY
C DEFINE THE CENTER OF THE TITLE STRING.
C
c DATA TX/0.42676/, TY/0.97656/
TX = 0.42676
TY = 0.97656
C
C
C INITIALIZE ERROR PARAMETER
C
IERROR = 0
C
C FILL TWO DIMENSIONAL ARRAY TO BE PLOTTED
C
DO 20 I=1,21
X = .1*FLOAT(I-11)
DO 10 J=1,25
Y = .1*FLOAT(J-13)
Z(I,J) = X+Y+1./((X-.10)**2+Y**2+.09)-
1 1./((X+.10)**2+Y**2+.09)
10 CONTINUE
20 CONTINUE
C
C SELECT NORMAIZATION TRANS NUMBER TO WRITE TITLES
C
CALL GSELNT (0)
C
C ENTRY EZCNTR REQUIRES ONLY THE ARRAY NAME AND ITS DIMENSIONS
C
C THE TITLE FOR THIS PLOT IS
C
C DEMONSTRATION PLOT FOR EZCNTR ENTRY OF CONRECSMTH
C
CALL WTSTR (TX,TY,
1 'DEMONSTRATION PLOT FOR EZCNTR ENTRY OF CONRECSMTH',
2 2,0,0)
CALL EZCNTR (Z,21,25)
C
C
C ENTRY CONREC ALLOWS USER SPECIFICATION OF PLOT PARAMETERS, IF DESIRED
C
C IN THIS EXAMPLE, THE LOWEST CONTOUR LEVEL (-4.5), THE HIGHEST CONTOUR
C LEVEL (4.5), AND THE INCREMENT BETWEEN CONTOUR LEVELS (0.3) ARE
C SPECIFIED.
C
C THE TITLE FOR THIS PLOT IS
C
C DEMONSTRATION PLOT FOR CONREC ENTRY OF CONRECSMTH
C
CALL WTSTR (TX,TY,
1 'DEMONSTRATION PLOT FOR CONREC ENTRY OF CONRECSMTH',
2 2,0,0)
CALL CONREC (Z,21,21,25,-4.5,4.5,.3,0,0,0)
c CALL NEWFM
C
c WRITE (6,1001)
RETURN
C
c 1001 FORMAT (' CONRECSMTH TEST SUCCESSFUL',24X,
c 1 'SEE PLOT TO VERIFY PERFORMANCE')
C
C
C---------------------------------------------------------------------
C
C REVISION HISTORY
C
C JUNE 1984 CONVERTED TO FORTRAN 77 AND GKS
C
C---------------------------------------------------------------------
END
|
