diff options
Diffstat (limited to 'noao/onedspec/sensfunc/sfmove.x')
| -rw-r--r-- | noao/onedspec/sensfunc/sfmove.x | 166 |
1 files changed, 166 insertions, 0 deletions
diff --git a/noao/onedspec/sensfunc/sfmove.x b/noao/onedspec/sensfunc/sfmove.x new file mode 100644 index 00000000..4451938e --- /dev/null +++ b/noao/onedspec/sensfunc/sfmove.x @@ -0,0 +1,166 @@ +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 |