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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
define MAX_FILES 50
# COLUMNS -- Convert a multicolumn file into a series of single column files.
# One output file `columns.n' is produced for each column in the input file.
#
# usage: COLUMNS input_file ncols
#
# This routine was originally written to allow SDAS to treat multicolumn
# tables as simple CL lists. Each column in the table is referenced in
# SDAS by a different parameter, pointing in the .par file to
# a different list.
procedure t_columns()
int numcols, infd, nchar, nfile, ip
pointer sp, fname, ofile, oroot, line, word, fnum, ofd
int clgeti(), open(), getline(), itoc(), ctowrd()
begin
call smark (sp)
call salloc (fname, SZ_FNAME, TY_CHAR)
call salloc (ofile, SZ_FNAME, TY_CHAR)
call salloc (oroot, SZ_FNAME, TY_CHAR)
call salloc (line, SZ_LINE, TY_CHAR)
call salloc (word, SZ_LINE, TY_CHAR)
call salloc (fnum, SZ_FNAME, TY_CHAR)
call salloc (ofd, MAX_FILES, TY_INT)
# Get the number of columns and the input file name
call clgstr ("filename", Memc[fname], SZ_FNAME)
numcols = clgeti ("numcols")
if (numcols > MAX_FILES)
call error (1, "too many files (columns)")
call clgstr ("outroot", Memc[oroot], SZ_FNAME)
# Open all the files
infd = open (Memc[fname], READ_ONLY, TEXT_FILE)
for (nfile=1; nfile <= numcols; nfile=nfile+1) {
nchar = itoc (nfile, Memc[fnum], 2)
call strcpy (Memc[oroot], Memc[ofile], SZ_FNAME)
call strcat (Memc[fnum], Memc[ofile], SZ_FNAME)
Memi[ofd+nfile-1] = open (Memc[ofile], NEW_FILE, TEXT_FILE)
}
# Separate each line of the input file
while (getline (infd, Memc[line]) != EOF) {
if ((Memc[line] != '#') && (Memc[line] != '\n')) {
ip = 1
for (nfile=1; nfile <= numcols; nfile=nfile+1) {
nchar = ctowrd (Memc[line], ip, Memc[word], SZ_LINE)
call strcat ("\n", Memc[word], SZ_LINE)
call putline (Memi[ofd+nfile-1], Memc[word])
}
}
}
# Close the files.
call close (infd)
for (nfile=1; nfile <= numcols; nfile=nfile+1)
call close (Memi[ofd+nfile-1])
call sfree (sp)
end
|
