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
|
include <gset.h>
include <pkg/rg.h>
include <pkg/gtools.h>
include "hdicfit.h"
# ICG_SAMPLE -- Mark sample.
procedure icg_sampled (ic, gp, gt, x, npts, pltype)
pointer ic # ICFIT pointer
pointer gp # GIO pointer
pointer gt # GTOOLS pointer
double x[npts] # Ordinates of graph
int npts # Number of data points
int pltype # Plot line type
pointer rg
int i, axis, pltype1
real xl, xr, yb, yt, dy
real x1, x2, y1, y2, y3
int gstati(), stridxs(), gt_geti()
pointer rg_xrangesd()
begin
if (stridxs ("*", Memc[IC_SAMPLE(ic)]) > 0)
return
# Find axis along which the independent data is plotted.
if (IC_AXES(ic,IC_GKEY(ic),1) == 'x')
axis = 1
else if (IC_AXES(ic,IC_GKEY(ic),2) == 'x')
axis = 2
else
return
if (gt_geti (gt, GTTRANSPOSE) == YES)
axis = mod (axis, 2) + 1
pltype1 = gstati (gp, G_PLTYPE)
call gseti (gp, G_PLTYPE, pltype)
rg = rg_xrangesd (Memc[IC_SAMPLE(ic)], x, npts)
switch (axis) {
case 1:
call ggwind (gp, xl, xr, yb, yt)
dy = yt - yb
y1 = yb + dy / 100
y2 = y1 + dy / 20
y3 = (y1 + y2) / 2
do i = 1, RG_NRGS(rg) {
x1 = x[RG_X1(rg, i)]
x2 = x[RG_X2(rg, i)]
if ((x1 > xl) && (x1 < xr))
call gline (gp, x1, y1, x1, y2)
if ((x2 > xl) && (x2 < xr))
call gline (gp, x2, y1, x2, y2)
call gline (gp, x1, y3, x2, y3)
}
case 2:
call ggwind (gp, yb, yt, xl, xr)
dy = yt - yb
y1 = yb + dy / 100
y2 = y1 + dy / 20
y3 = (y1 + y2) / 2
do i = 1, RG_NRGS(rg) {
x1 = x[RG_X1(rg, i)]
x2 = x[RG_X2(rg, i)]
if ((x1 > xl) && (x1 < xr))
call gline (gp, y1, x1, y2, x1)
if ((x2 > xl) && (x2 < xr))
call gline (gp, y1, x2, y2, x2)
call gline (gp, y3, x1, y3, x2)
}
}
call gseti (gp, G_PLTYPE, pltype1)
call rg_free (rg)
end
|
