aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/ptools/pselect.cl
blob: 86ebb9bf49ec9cb2b0ebc3183fec0ae67d116576 (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
# PSELECT - Select records from an ST table of an APPHOT/DAOPHOT text file
# based on the value of a boolean expression.

procedure pselect (infiles, outfiles, expr)

string	infiles		{prompt="Input apphot/daophot database(s)"}
string	outfiles	{prompt="Output apphot/daophot database(s)"}
string	expr		{prompt="Boolean expression for record selection"}

struct	*inlist, *outlist

begin
	# Local variable declarations.
	file	tmpin, tmpout
	int	nin, nout
	string	in, out, ex, inname, outname

	# Cache the istable parameters.
	cache ("istable")

	# Get the positional parameters.
	in = infiles
	out = outfiles
	ex = expr

	# Make temporary names.
	tmpin = mktemp ("tmp$")
	tmpout = mktemp ("tmp$")

	# Expand the file list names.
	files (in, sort=no, > tmpin)
	files (out, sort=no, > tmpout)

	# Compute the lengths of the input and output lists.
	inlist = tmpin
	for (nin = 0; fscan (inlist, inname) != EOF; nin = nin + 1)
	    ;
	inlist = ""
	outlist = tmpout
	for (nout = 0; fscan (outlist, outname) != EOF; nout = nout + 1)
	    ;
	outlist = ""

	# Delete the temporary files.
	delete (tmpin, ver-, >& "dev$null")
	delete (tmpout, ver-, >& "dev$null")

	# Quit if the number of input and output files is different.
	if (nin != nout) {
	    print ("ERROR: Input and output file lists are different lengths")
	    return
	}

	# Expand the file list names.
	files (in, sort=no, > tmpin)
	files (out, sort=no, > tmpout)

	# Loop over each file in the input and output lists selecting records.
	inlist = tmpin
	outlist = tmpout
	while (fscan (inlist, inname) != EOF && fscan (outlist, outname) !=
	    EOF) {
	    istable (inname)
	    if (istable.table) {
		tselect (inname, outname, ex)
	    } else if (istable.text) {
		txselect (inname, outname, ex)
	    } else {
		print ("ERROR: Cannot run PSELECT on file: " // inname)
	    }
	}

	delete (tmpin, ver-, >& "dev$null")
	delete (tmpout, ver-, >& "dev$null")
	inlist = ""
	outlist = ""
end