include	<gset.h>
include "sensfunc.h"


# SF_MOVE -- Move point, star, or wavelength.

procedure sf_move (gp, stds, nstds, key, istd, ipt, shift)

pointer	gp			# GIO pointer
pointer	stds[nstds]		# Standard star data
int	nstds			# Number of standard stars
int	key			# Delete point, star, or wavelength
int	istd			# Index of standard star
int	ipt			# Index of point
real	shift

int	i, j, n, wcs, mark, mdel, cdel, color, stridx()
real	wave, szmark, szmdel
pointer	x, y, z, w, gio

begin
	gio = GP_GIO(gp)
	mdel = GP_MDEL(gp)
	cdel = GP_CDEL(gp)
	szmdel = GP_SZMDEL(gp)
	szmark = GP_SZMARK(gp)

	# Move points in each displayed graph.
	for (wcs = 1; GP_GRAPHS(gp,wcs) != EOS; wcs = wcs + 1) {
	    if (stridx (GP_GRAPHS(gp,wcs), "ars") == 0)
	        next

	    call gseti (gio, G_WCS, wcs)
	    call sf_data (stds, nstds, GP_GRAPHS(gp,wcs))
	    switch (key) {
	    case 'p':
		if (istd != nstds-1) {
		    mark = GP_MARK(gp)
		    color = GP_CMARK(gp)
		} else {
		    mark = GP_MADD(gp)
		    color = GP_CADD(gp)
		}
	        x = STD_X(stds[istd])+ipt-1
	        y = STD_Y(stds[istd],1)+ipt-1
	        w = STD_WTS(stds[istd])+ipt-1
		if (Memr[w] != 0.) {
		    call gseti (gio, G_PMLTYPE, 0)
	            call gmark (gio, Memr[x], Memr[y], mark, szmark, szmark)
		    call gseti (gio, G_PMLTYPE, 1)
		    call gseti (gio, G_PLCOLOR, color)
	            call gmark (gio, Memr[x], Memr[y]+shift, mark, szmark,
			szmark)
		} else {
		    call gseti (gio, G_PMLTYPE, 0)
	            call gmark (gio, Memr[x], Memr[y], mdel, szmdel, szmdel)
		    call gseti (gio, G_PMLTYPE, 1)
		    call gseti (gio, G_PLCOLOR, cdel)
	            call gmark (gio, Memr[x], Memr[y]+shift, mdel, szmdel,
			szmdel)
		}
	    case 's':
		if (istd != nstds-1) {
		    mark = GP_MARK(gp)
		    color = GP_CMARK(gp)
		} else {
		    mark = GP_MADD(gp)
		    color = GP_CADD(gp)
		}
	        n = STD_NWAVES(stds[istd])
	        x = STD_X(stds[istd])
	        y = STD_Y(stds[istd])
	        w = STD_WTS(stds[istd])
		do i = 1, n {
		    if (Memr[w] != 0.) {
		        call gseti (gio, G_PMLTYPE, 0)
	                call gmark (gio, Memr[x], Memr[y], mark, szmark, szmark)
		        call gseti (gio, G_PMLTYPE, 1)
		        call gseti (gio, G_PLCOLOR, color)
	                call gmark (gio, Memr[x], Memr[y]+shift, mark, szmark,
			    szmark)
		    } else {
		        call gseti (gio, G_PMLTYPE, 0)
	                call gmark (gio, Memr[x], Memr[y], mdel, szmdel, szmdel)
		        call gseti (gio, G_PMLTYPE, 1)
		        call gseti (gio, G_PLCOLOR, cdel)
	                call gmark (gio, Memr[x], Memr[y]+shift, mdel, szmdel,
			    szmdel)
		    }
		    x = x + 1
		    y = y + 1
		    w = w + 1
		}
	    case 'w':
	        wave = Memr[STD_WAVES(stds[istd])+ipt-1]
	        do i = 1, nstds {
	            if (STD_FLAG(stds[i]) != SF_INCLUDE)
		        next
		    if (i != nstds-1) {
			mark = GP_MARK(gp)
			color = GP_CMARK(gp)
		    } else {
			mark = GP_MADD(gp)
			color = GP_CADD(gp)
		    }
	            n = STD_NWAVES(stds[i])
	            x = STD_X(stds[i])
	            y = STD_Y(stds[i])
		    z = STD_WAVES(stds[i])
	            w = STD_WTS(stds[i])
		    do j = 1, n {
		        if (Memr[z] == wave) {
			    if (Memr[w] != 0.) {
		                call gseti (gio, G_PMLTYPE, 0)
	                        call gmark (gio, Memr[x], Memr[y], mark, szmark,
				    szmark)
		                call gseti (gio, G_PMLTYPE, 1)
		                call gseti (gio, G_PLCOLOR, color)
	                        call gmark (gio, Memr[x], Memr[y]+shift, mark,
				    szmark, szmark)
			    } else {
		                call gseti (gio, G_PMLTYPE, 0)
	                        call gmark (gio, Memr[x], Memr[y], mdel, szmdel,
				    szmdel)
		                call gseti (gio, G_PMLTYPE, 1)
		                call gseti (gio, G_PLCOLOR, cdel)
	                        call gmark (gio, Memr[x], Memr[y]+shift, mdel,
				    szmdel, szmdel)
			    }
		        }
			x = x + 1
			y = y + 1
			z = z + 1
			w = w + 1
		    }
	        }
	    }
	}

	# Now add the shift to the data.
	switch (key) {
	case 'p':
	    y = STD_SENS(stds[istd])+ipt-1
	    Memr[y] = Memr[y] + shift
	case 's':
	    n = STD_NWAVES(stds[istd])
	    y = STD_SENS(stds[istd])
	    call aaddkr (Memr[y], shift, Memr[y], n)
	    STD_SHIFT(stds[istd]) = STD_SHIFT(stds[istd]) + shift
	case 'w':
	    wave = Memr[STD_WAVES(stds[istd])+ipt-1]
	    do i = 1, nstds {
	        if (STD_FLAG(stds[i]) != SF_INCLUDE)
		    next
	        n = STD_NWAVES(stds[i])
		z = STD_WAVES(stds[i])
	        y = STD_SENS(stds[i])
		do j = 1, n {
		    if (Memr[z] == wave)
	                Memr[y] = Memr[y] + shift
		    w = w + 1
		    y = y + 1
		}
	    }
	}
end