1 files changed, 84 insertions, 0 deletions
diff --git a/noao/imred/dtoi/hdicfit/hdicgsample.x b/noao/imred/dtoi/hdicfit/hdicgsample.x
new file mode 100644
index 00000000..12eef158
--- /dev/null
+++ b/noao/imred/dtoi/hdicfit/hdicgsample.x
@@ -0,0 +1,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
+