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
|
! $Id: x13f.f90,v 1.1 2006/05/16 20:24:13 airwin Exp $
! Pie chart demo.
!
! Copyright (C) 2004 Alan W. Irwin
!
! This file is part of PLplot.
!
! PLplot is free software; you can redistribute it and/or modify
! it under the terms of the GNU General Library Public License as
! published by the Free Software Foundation; either version 2 of the
! License, or (at your option) any later version.
!
! PLplot is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU Library General Public License for more details.
!
! You should have received a copy of the GNU Library General Public
! License along with PLplot; if not, write to the Free Software
! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
use plplot, PI => PL_PI
implicit none
real(kind=plflt) just, dx, dy, x(500), y(500), per(5)
integer i, j, dthet, theta0, theta1, theta
character*20 text(5)
data text / 'Maurice', 'Geoffrey', 'Alan', 'Rafael', 'Vince'/
! Process command-line arguments
call plparseopts(PL_PARSE_FULL)
per(1) = 10._plflt
per(2) = 32._plflt
per(3) = 12._plflt
per(4) = 30._plflt
per(5) = 16._plflt
! Ask user to specify the output device.
call plinit()
call plenv( 0._plflt, 10._plflt, 0._plflt, 10._plflt, 1, -2 )
call plcol0(2)
! n.b. all theta quantities scaled by 2*pi/500 to be integers to avoid
! floating point logic problems.
theta0 = 0
dthet = 1
do i = 0, 4
j = 0
x(j+1) = 5._plflt
y(j+1) = 5._plflt
j = j + 1
! n.b. the theta quantities multiplied by 2*pi/500 afterward so
! in fact per is interpreted as a percentage.
theta1 = theta0 + 5*per(i+1)
if (i .eq. 4) theta1 = 500
do theta = theta0, theta1, dthet
x(j+1) = 5 + 3*cos((2._plflt*pi/500._plflt)*theta)
y(j+1) = 5 + 3*sin((2._plflt*pi/500._plflt)*theta)
j = j + 1
enddo
call plcol0(i+1)
call plpsty( mod(i+3, 8) + 1 )
call plfill(x(:j), y(:j))
call plcol0(1)
call plline(x(:j), y(:j))
just = (2._plflt*pi/500._plflt)*(theta0 + theta1) / 2._plflt
dx = 0.25_plflt * cos(just)
dy = 0.25_plflt * sin(just)
if ((theta0 + theta1) .lt. 250 .or. &
(theta0 + theta1) .gt. 750) then
just = 0._plflt
else
just = 1._plflt
endif
call plptex( x(j/2+1)+dx, y(j/2+1)+dy, 1._plflt, 0._plflt, &
just, text(i+1) )
theta0 = theta - dthet
enddo
call plfont(2)
call plschr( 0._plflt, 1.3_plflt)
call plptex( 5._plflt, 9._plflt, 1._plflt, 0._plflt, &
0.5_plflt, 'Percentage of Sales' )
call plend
end
|
