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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
|
include <fset.h>
# T_GALACTIC -- convert between equatorial and galactic coordinates.
procedure t_galactic ()
char fname[SZ_FNAME]
int filelist, prt_coords, in_coords
bool streq(), clgetb()
int clpopni(), clgfil(), btoi(), clgwrd()
begin
# Input can come from the standard input, a file, or a list of files.
# The following procedure makes both cases look like a list of files.
filelist = clpopni ("input")
# Get output option
in_coords = clgwrd ("in_coords", fname, SZ_FNAME,
"|equatorial|galactic|")
prt_coords = btoi (clgetb ("print_coords"))
# Process each coordinate list. If reading from the standard input,
# set up the standard output to flush after every output line, so that
# converted coords come back immediately when working interactively.
while (clgfil (filelist, fname, SZ_FNAME) != EOF) {
if (streq (fname, "STDIN"))
call fseti (STDOUT, F_FLUSHNL, YES)
else
call fseti (STDOUT, F_FLUSHNL, NO)
switch (in_coords) {
case 1:
call in_equatorial (STDOUT, fname, prt_coords)
case 2:
call in_galactic (STDOUT, fname, prt_coords)
}
}
call clpcls (filelist)
end
# IN_EQUITORIAL -- convert a list of equatorial coordinates read from the
# named file to galactic coordinates, writing the results on the output file.
procedure in_equatorial (out, listfile, prt_coords)
int out # output stream
char listfile[SZ_FNAME] # input file
int prt_coords # print coordinates in output file?
int in
double ra, dec, year, lii, bii
int fscan(), nscan(), open()
errchk open, fscan, printf
begin
in = open (listfile, READ_ONLY, TEXT_FILE)
# Read successive RA,DEC coordinate pairs from the standard input,
# converting and printing the result on the standard output.
while (fscan (in) != EOF) {
call gargd (ra)
call gargd (dec)
call gargd (year)
# If there is something wrong with the input coords, print warning
# and skip the conversion. If year is not given with entry
# assume 1950.0.
switch (nscan()) {
case 2:
year = 1950.0
case 3:
default:
call eprintf ("Bad entry in coordinate list\n")
next
}
# Call routine to perform conversion and write lII, bII to
# the standard output.
call ast_galactic (ra, dec, year, lii, bii)
if (prt_coords == YES ) {
call fprintf (out, "%13.2h %12.1h %8.2f")
call pargd (ra)
call pargd (dec)
call pargd (year)
}
call fprintf (out, "%13.4f %9.4f" )
call pargd (lii)
call pargd (bii)
call fprintf (out, "\n")
}
call close (in)
end
# IN_GALACTIC -- convert a list of galactic coordinates read from the
# named file to equatorial coordinates, writing the results on the output
# file.
procedure in_galactic (out, listfile, prt_coords)
int out # output stream
char listfile[SZ_FNAME] # input file
int prt_coords # print coordinates in output file?
int in
double ra, dec, year, lii, bii
int fscan(), nscan(), open()
errchk open, fscan, printf
begin
in = open (listfile, READ_ONLY, TEXT_FILE)
# Read successive RA,DEC coordinate pairs from the standard input,
# converting and printing the result on the standard output.
while (fscan (in) != EOF) {
call gargd (lii)
call gargd (bii)
call gargd (year)
# If there is something wrong with the input coords, print warning
# and skip the conversion. If year is not given with entry
# assume 1950.0.
switch (nscan()) {
case 2:
year = 1950.0
case 3:
default:
call eprintf ("Bad entry in coordinate list\n")
next
}
# Call routine to perform conversion and write lII, bII to
# the standard output.
call ast_galtoeq (lii, bii, ra, dec, year)
if (prt_coords == YES ) {
call fprintf (out, "%13.4f %9.4f" )
call pargd (lii)
call pargd (bii)
}
call fprintf (out, "%13.2h %12.1h %8.2f")
call pargd (ra)
call pargd (dec)
call pargd (year)
call fprintf (out, "\n")
}
call close (in)
end
|
