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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <ctype.h>
include <plset.h>
# PL_ASCIIDUMP -- Dump a two dimensional region of a mask as a printable ASCII
# character array on the output stream. Intended as a simple debugging tool;
# see also PL_SAVEIM.
procedure pl_asciidump (pl, vs, ve, outfd)
pointer pl #I mask descriptor
long vs[ARB] #I ll vector (only first two elements used)
long ve[ARB] #I ur vector (only first two elements used)
int outfd #I output file
pointer sp, pv, cv
int npix, ch, i
long v[PL_MAXDIM]
errchk pl_valid
begin
call pl_valid (pl)
npix = ve[1] - vs[1] + 1
call smark (sp)
call salloc (pv, npix, TY_INT)
call salloc (cv, npix, TY_CHAR)
# Output mask.
call amovl (vs, v, PL_MAXDIM)
v[2] = ve[2]
while (v[2] >= vs[2]) {
call pl_glpi (pl, v, Memi[pv], 0, npix, PIX_SRC)
do i = 1, npix {
ch = mod (Memi[pv+i-1], 128)
if (ch < 32)
ch = ch + 32
if (ch <= 32 || ch == 127)
ch = '.'
Memc[cv+i-1] = ch
}
call fprintf (outfd, "%3d ")
call pargi (v[2])
call write (outfd, Memc[cv], npix)
call putci (outfd, '\n')
v[2] = v[2] - 1
}
# Label the columns.
call fprintf (outfd, "%5t")
do i = 1, npix
call putci (outfd, TO_DIGIT(mod (i,10)))
call fprintf (outfd, "\n")
call fprintf (outfd, "%5t")
do i = 1, npix
if (mod (i, 10) == 0)
call putci (outfd, TO_DIGIT(i / 10))
else
call putci (outfd, ' ')
call fprintf (outfd, "\n")
call sfree (sp)
end
|
