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
|
SUBROUTINE TCONAN (IERROR)
C
C LATEST REVISION JULY 1984
C
C PURPOSE TO PROVIDE A SIMPLE DEMONSTRATION OF
C CONRAN, THE STANDARD ENTRY POINT OF THE
C CONRAN PACKAGE.
C
C THIS SAME SUBROUTINE CAN BE USED TO PRODUCE
C DEMO PLOTS OF THE SMOOTH VERSION OF CONRAN
C BY LOADING DASHSMTH INSTEAD OF DASHCHAR.
C
C USAGE CALL TCONAN (IERROR)
C
C ARGUMENTS
C
C ON OUTPUT IERROR
C AN INTEGER VARIABLE
C .EQ. 0, IF THE TEST WAS SUCCESSFUL,
C .NE. 0, OTHERWISE
C IF NOT ZERO THE NUMBER PRODUCED WILL
C CORRESPOND TO THE ERROR NUMBERS IN
C THE CONRAN LISTING.
C
C I/O IF THE TEST IS SUCCESSFUL, THE MESSAGE
C
C CONRAN 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 AND TRIANGULATION OF THE DATA ARE PRODUCED
C ON THE DEFAULT GRAPHICS DEVICE UNLESS THE USER
C SPECIFIES OTHERWISE VIA JCL.
C
C PRECISION SINGLE
C
C REQUIRED LIBRARY CONRAN
C FILES CONTERP
C CONCOM
C
C LANGUAGE FORTRAN77
C
C ALGORITHM A SPARSE DATA SET IS DEFINED VIA DATA
C STATEMENTS. OPTIONS ARE SET TO PRODUCE A
C TITLE AND DISPLAY THE TRIANGULATION GENERATED
C BY THE INTERPOLATING ROUTINES. BY DEFAULT
C A MESSAGE AT THE BOTTEM OF THE PLOT AND A
C PERIMETER ARE ALSO PRODUCED. THIS ROUTINE
C TAKES ADVANTAGE OF THE PORT ERROR HANDLING
C ROUTINES TO DETERMINE IF CONRAN TERMINATED
C NORMALLY.
C
C PORTABILITY ANSI FORTRAN77 STANDARD
C
C COMMON /RANINT/ IRANMJ, IRANMN, IRANTX
C SET UP THE SCRATCH SPACES REQUIRED BY CONRAN
C
DIMENSION WK(221),IWK(744),SCR(1600)
C
C SET UP THE ARRAYS TO DEFINE THE DATA SET
C
DIMENSION XD(17),YD(17),ZD(17)
C
C DEFINE THE DATA SET
C
DATA XD(1),XD(2),XD(3),XD(4),XD(5),XD(6),XD(7),XD(8),
1 XD(9),XD(10),XD(11),XD(12),XD(13),XD(14),XD(15),
2 XD(16),XD(17)
3 /3.,3.,10.,18.,18.,10.,10.,5.,1.,15.,20.,
4 5.,15.,10.,7.,13.,16./
C
DATA YD(1),YD(2),YD(3),YD(4),YD(5),YD(6),YD(7),YD(8),
1 YD(9),YD(10),YD(11),YD(12),YD(13),YD(14),YD(15),
2 YD(16),YD(17)
3 /3.,18.,18.,3.,18.,10.,1.,5.,10.,5.,10.,
4 15.,15.,15.,20.,20.,8./
C
DATA ZD(1),ZD(2),ZD(3),ZD(4),ZD(5),ZD(6),ZD(7),ZD(8),
1 ZD(9),ZD(10),ZD(11),ZD(12),ZD(13),ZD(14),ZD(15),
2 ZD(16),ZD(17)
3 /25.,25.,25.,25.,25.,-5.,1.,1.,1.,1.,1.,
4 1.,1.,1.,1.,1.,25./
C
C SET UP PARAMETER FOR NUMBER OF INPUT POINTS
C
DATA NDP/17/
call conbdn
C
C SET UP TITLE FOR PLOT
C
CALL CONOP4('TLE=ON','DEMONSTRATION PLOT FOR CONRAN',29, 0)
C
CALL CONRAN(XD,YD,ZD,NDP,WK,IWK,SCR)
C
RETURN
END
|
