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
|
C
C gfortran griddata.f -L/usr/local/gr/lib -lGR
C
program griddata
integer GKS_K_MARKERTYPE_SOLID_CIRCLE
parameter (GKS_K_MARKERTYPE_SOLID_CIRCLE = -1)
double precision xd(100), yd(100), zd(100)
double precision x(200), y(200), z(200, 200)
double precision h(20)
integer i
call srand(0)
do 1 i = 1, 100
xd(i) = -2 + 4 * rand()
yd(i) = -2 + 4 * rand()
zd(i) = xd(i) * exp(-xd(i)*xd(i) - yd(i)*yd(i))
1 continue
call gr_setviewport(0.1D0, 0.95D0, 0.1D0, 0.95D0)
call gr_setwindow(-2.0D0, 2.0D0, -2.0D0, 2.0D0)
call gr_setspace(-0.5D0, 0.5D0, 0, 90)
call gr_setmarkersize(1.0D0)
call gr_setmarkertype(GKS_K_MARKERTYPE_SOLID_CIRCLE)
call gr_setcharheight(0.024D0)
call gr_settextalign(2, 0)
call gr_settextfontprec(3, 0)
call gr_gridit(100, xd, yd, zd, 200, 200, x, y, z)
do 2 i = 1, 20
h(i) = -0.5 + i / 19.0
2 continue
call gr_surface(200, 200, x, y, z, 5)
call gr_contour(200, 200, 20, x, y, h, z, 0)
call gr_polymarker(100, xd, yd)
call gr_axes(0.25D0, 0.25D0, -2D0, -2D0, 2, 2, 0.01D0)
call gr_updatews()
end
|
