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
|
procedure mkbbcal (band, w1, w2, dw, weff, mab)
string band {prompt="Bandpass"}
real w1 {prompt="Starting wavelength (um)"}
real w2 {prompt="Ending wavelength (um)"}
real dw {prompt="Wavelength step (um)"}
real weff {prompt="Effective wavelength (um)"}
real mab {prompt="M_AB of Vega at eff. wavelength"}
begin
string b
real x1, x2, dx, x, m
b = band
x1 = w1
x2 = w2
dx = dw
x = weff
m = mab
printf ("# type blackbody\n")
printf ("# units microns\n")
printf ("# band %s\n", b)
printf ("# weff %g\n\n", x)
for (x=x1; x<=x2; x=x+dx)
printf ("%.5g\t%.6g\t%g\n", x, m, dx)
end
|
