Initial commit

1 files changed, 28 insertions, 0 deletions

diff --git a/noao/lib/onedstds/blackbody/mkbbcal.cl b/noao/lib/onedstds/blackbody/mkbbcal.cl
new file mode 100644
index 00000000..42575276
--- /dev/null
+++ b/noao/lib/onedstds/blackbody/mkbbcal.cl
@@ -0,0 +1,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