aboutsummaryrefslogtreecommitdiff
path: root/noao/onedspec/ecidentify/ecshow.x
blob: e8fb5acc65d51f59b7978b65e6d06164a7c8fad2 (plain) (blame)
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
include	<pkg/center1d.h>
include	"ecidentify.h"

# EC_SHOW -- Show parameter information.

procedure ec_show (ec, file)

pointer	ec			# ID pointer
char	file[ARB]		# File

char	line[SZ_LINE]
int	fd

int	open(), ecf_geti()
double	ecf_getd()
errchk	open()

begin
	fd = open (file, APPEND, TEXT_FILE)

	call sysid (line, SZ_LINE)
	call fprintf (fd, "%s\n")
	call pargstr (line)

	call fprintf (fd, "image %s\n")
	    call pargstr (Memc[EC_IMAGE(ec)])
	switch (EC_FTYPE(ec)) {
	case EMISSION:
	    call fprintf (fd, "ftype emission\n")
	case ABSORPTION:
	    call fprintf (fd, "ftype absorption\n")
	}
	switch (EC_LABELS(ec)) {
	case 2:
	    call fprintf (fd, "labels index\n")
	case 3:
	    call fprintf (fd, "labels pixel\n")
	case 4:
	    call fprintf (fd, "labels user\n")
	default:
	    call fprintf (fd, "labels none\n")
	}
	call fprintf (fd, "maxfeatures %d\n")
	    call pargi (EC_MAXFEATURES(ec))
	call fprintf (fd, "match %g\n")
	    call pargr (EC_MATCH(ec))
	call fprintf (fd, "zwidth %g\n")
	    call pargr (EC_ZWIDTH(ec))
	call fprintf (fd, "fwidth %g\n")
	    call pargr (EC_FWIDTH(ec))
	call fprintf (fd, "database %s\n")
	    call pargstr (Memc[EC_DATABASE(ec)])
	call fprintf (fd, "coordlist %s\n")
	    call pargstr (Memc[EC_COORDLIST(ec)])
	call fprintf (fd, "cradius %g\n")
	    call pargr (EC_CRADIUS(ec))
	call fprintf (fd, "threshold %g\n")
	    call pargr (EC_THRESHOLD(ec))
	call fprintf (fd, "minsep %g\n")
	    call pargr (EC_MINSEP(ec))
	if (EC_ECF(ec) != NULL) {
	    call fprintf (fd, "function = %s\n")
	        call ecf_gets ("function", line, SZ_LINE)
		call pargstr (line)
	    call fprintf (fd, "xorder = %d, yorder = %d\n")
		call pargi (ecf_geti ("xorder"))
		call pargi (ecf_geti ("yorder"))
	    call fprintf (fd,
		"niterate = %d, lowreject = %g, highreject = %g\n")
		call pargi (ecf_geti ("niterate"))
		call pargd (ecf_getd ("low"))
		call pargd (ecf_getd ("high"))
	    call fprintf (fd, "Fit at first pixel = %0.8g\n")
		call pargd (Memd[EC_FITDATA(ec)])
	}

	call close (fd)
end