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
|
include <ctype.h>
define VLIGHT 2.997925e18 # Speed of light in Angstroms/sec
# Options
define OPTION "|ext|mags|fnu|flam|bands|stars|"
define EXT 1 # Extinction
define MAGS 2 # Standard star magnitudes
define FNU 3 # Standard star fluxes
define FLAM 4 # Standard star fluxes
define BANDS 5 # Standard star band passes
define STARS 6 # Standard stars
# T_LCALIB -- List information in calibration file:
# 1) Extinction vs wavelength
# 2) Magnitude vs wavelength
# 3) F-nu vs wavelength
# 4) F-lambda vs wavelength
# 5) Bandpass vs wavelength
# 6) Standard stars
procedure t_lcalib ()
int i, nwaves, fd
real fnu, flam, fnuzero
pointer sp, str, file, waves, bands, mags, extns
int getline(), open(), clgwrd()
real clgetr()
errchk ext_load, getcalib
begin
call smark (sp)
call salloc (str, SZ_LINE, TY_CHAR)
call salloc (file, SZ_LINE, TY_CHAR)
#Switch on the option.
switch (clgwrd ("option", Memc[str], SZ_LINE, OPTION)) {
case EXT:
call ext_load (waves, extns, nwaves)
do i = 1, nwaves {
call printf ("%6f %12.5g\n")
call pargr (Memr[waves+i-1])
call pargr (Memr[extns+i-1])
}
call ext_free (waves, extns)
case MAGS:
call getcalib (waves, bands, mags, nwaves)
do i = 1, nwaves {
call printf ("%6f %12.5g\n")
call pargr (Memr[waves+i-1])
call pargr (Memr[mags+i-1])
}
call freecalib (waves, bands, mags)
case FNU:
fnuzero = clgetr ("fnuzero")
call getcalib (waves, bands, mags, nwaves)
do i = 1, nwaves {
fnu = fnuzero * 10. ** (-0.4 * Memr[mags+i-1])
call printf ("%6f %12.5g\n")
call pargr (Memr[waves+i-1])
call pargr (fnu)
}
call freecalib (waves, bands, mags)
case FLAM:
fnuzero = clgetr ("fnuzero")
call getcalib (waves, bands, mags, nwaves)
do i = 1, nwaves {
fnu = fnuzero * 10. ** (-0.4 * Memr[mags+i-1])
flam = fnu * VLIGHT / Memr[waves+i-1] ** 2
call printf ("%6f %12.5g\n")
call pargr (Memr[waves+i-1])
call pargr (flam)
}
call freecalib (waves, bands, mags)
case BANDS:
call getcalib (waves, bands, mags, nwaves)
do i = 1, nwaves {
call printf ("%6f %12.5g\n")
call pargr (Memr[waves+i-1])
call pargr (Memr[bands+i-1])
}
call freecalib (waves, bands, mags)
case STARS:
call clgstr ("caldir", Memc[str], SZ_LINE)
call sprintf (Memc[file], SZ_LINE, "%sstandards.men")
call pargstr (Memc[str])
fd = open (Memc[file], READ_ONLY, TEXT_FILE)
while (getline (fd, Memc[file]) != EOF)
call putline (STDERR, Memc[file])
call close (fd)
default:
call eprintf ("Unknown option: %s\n")
call pargstr (Memc[str])
}
call sfree (sp)
end
|
