Initial commit

28 files changed, 1316 insertions, 0 deletions

diff --git a/noao/artdata/mkexamples/archdr.dat b/noao/artdata/mkexamples/archdr.dat
new file mode 100644
index 00000000..b52dcc32
--- /dev/null
+++ b/noao/artdata/mkexamples/archdr.dat
@@ -0,0 +1,19 @@
+OBSERVER= 'IRAF              '  /  observers
+OBSERVAT= 'KPNO              '  /  observatory
+INSTRUME= 'IRAF ARTDATA      '  /  instrument
+TELESCOP= 'kpcdf             '  /  telescope name
+DETECTOR= 'te1k              '  /  detector
+EXPTIME =                  60.  /  actual integration time
+DARKTIME=                  60.  /  total elapsed time
+IMAGETYP= 'comp              '  /  object, dark, bias, etc.
+DATE-OBS= '1991-11-26T12:11:30.00' /  date (dd/mm/yy) of obs.
+UT      = '12:11:30.00       '  /  universal time
+ST      = '09:04:54.00       '  /  sidereal time
+RA      = '06:37:02.00       '  /  right ascension
+DEC     = '06:09:03.00       '  /  declination
+EPOCH   =               1991.9  /  epoch of ra and dec
+ZD      = '48.760            '  /  zenith distance
+AIRMASS =                   0.  /  airmass
+APERTURE= '250micron slit    '  /  aperture
+GRATPOS =               4624.3  /  grating position
+CCDPROC = 'Nov 26  5:44 CCD processing done'
diff --git a/noao/artdata/mkexamples/ecarc.cl b/noao/artdata/mkexamples/ecarc.cl
new file mode 100644
index 00000000..174eba2f
--- /dev/null
+++ b/noao/artdata/mkexamples/ecarc.cl
@@ -0,0 +1,20 @@
+# ecarc - Echelle extracted thorium arc (calibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$archdr.dat"
+i = 10
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header=hdr,
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="extracted", width=20., scattered=0., f=590., gmm=31.6, blaze=63.,
+    theta=69., order=112, wavelength=5007.49, dispersion=2.61, cf=590.,
+    cgmm=226., cblaze=4.53, ctheta=-11.97, corder=1, cwavelength=6700.,
+    cdispersion=70., rv=0., z=no, continuum=20., temperature=0.,
+    lines="mkexamples$ecthorium.dat", sigma=0.05, >& "dev$null")
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/ecarc2d.cl b/noao/artdata/mkexamples/ecarc2d.cl
new file mode 100644
index 00000000..8261e84a
--- /dev/null
+++ b/noao/artdata/mkexamples/ecarc2d.cl
@@ -0,0 +1,20 @@
+# ecarc2d - Echelle thorium-argon slit spectrum
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$archdr.dat"
+i = 10
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header=hdr,
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="slit", width=20., scattered=10., f=590., gmm=31.6, blaze=63.,
+    theta=69., order=112, wavelength=5007.49, dispersion=2.61, cf=590.,
+    cgmm=226., cblaze=4.53, ctheta=-11.97, corder=1, cwavelength=6700.,
+    cdispersion=70., rv=0., z=no, continuum=20., temperature=0.,
+    lines="mkexamples$ecthorium.dat", sigma=0.05, >& "dev$null")
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/ecarcdc.cl b/noao/artdata/mkexamples/ecarcdc.cl
new file mode 100644
index 00000000..d188d52f
--- /dev/null
+++ b/noao/artdata/mkexamples/ecarcdc.cl
@@ -0,0 +1,21 @@
+# ecarcdc - Echelle thorium-argon arc (calibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$archdr.dat"
+i = 10
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header=hdr,
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="extracted", width=20., scattered=0., f=590., gmm=INDEF,
+    blaze=INDEF, theta=INDEF, order=112, wavelength=5007.49,
+    dispersion=2.61, cf=590., cgmm=226., cblaze=4.53, ctheta=-11.97,
+    corder=1, cwavelength=6700., cdispersion=70., rv=0., z=no,
+    continuum=20., temperature=0., lines="mkexamples$ecthorium.dat",
+    sigma=0.05, >& "dev$null")
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/echelle.cl b/noao/artdata/mkexamples/echelle.cl
new file mode 100644
index 00000000..4e6dc197
--- /dev/null
+++ b/noao/artdata/mkexamples/echelle.cl
@@ -0,0 +1,22 @@
+# echelle - Echelle absorption spectrum (calibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$objhdr.dat"
+i = 10
+k = 10 * i
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header=hdr,
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="extracted", width=20., scattered=0., f=590., gmm=INDEF,
+    blaze=INDEF, theta=INDEF, order=112, wavelength=5007.49,
+    dispersion=2.61, cf=590., cgmm=226., cblaze=4.53, ctheta=-11.97,
+    corder=1, cwavelength=6700., cdispersion=70., rv=0., z=no,
+    continuum=1000., temperature=5700., lines="", nrandom=k, peak=-0.5,
+    sigma=0.5, seed=i, >& "dev$null")
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/ecobj2d.cl b/noao/artdata/mkexamples/ecobj2d.cl
new file mode 100644
index 00000000..303dbb69
--- /dev/null
+++ b/noao/artdata/mkexamples/ecobj2d.cl
@@ -0,0 +1,29 @@
+# ecobj2d - Echelle object slit spectrum
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$objhdr.dat"
+i = 10
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header=hdr,
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="gaussian", width=4., scattered=25., f=590., gmm=31.6, blaze=63.,
+    theta=69., order=112, wavelength=5007.49, dispersion=2.61, cf=590.,
+    cgmm=226., cblaze=4.53, ctheta=-11.97, corder=1, cwavelength=6700.,
+    cdispersion=70., rv=0., z=no, continuum=500., temperature=7700.,
+    lines="", nrandom=500, peak=-0.5, sigma=0.5, seed=i, >& "dev$null")
+
+mkechelle (out, yes, ncols=512, nlines=512, norders=i,
+    title="Artificial Echelle Spectrum", header="",
+    list=no, make=yes, comments=b1, xc=235.5, yc=INDEF, pixsize=0.027,
+    profile="slit", width=20., scattered=10., f=590., gmm=31.6, blaze=63.,
+    theta=69., order=112, wavelength=5007.49, dispersion=2.61, cf=590.,
+    cgmm=226., cblaze=4.53, ctheta=-11.97, corder=1, cwavelength=6700.,
+    cdispersion=70., rv=0., z=no, continuum=200., temperature=5700.,
+    lines="", nrandom=20, peak=5.0, sigma=0.1, seed=i+1, >& "dev$null")
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/ecthorium.dat b/noao/artdata/mkexamples/ecthorium.dat
new file mode 100644
index 00000000..2624bca7
--- /dev/null
+++ b/noao/artdata/mkexamples/ecthorium.dat
@@ -0,0 +1,655 @@
+4588.4273 65.5462
+4589.8978 784.211
+4592.666 56.8004
+4593.6437 4.85515
+4595.4206 222.271
+4596.0967 64.2832
+4596.3078 7.56086
+4598.7621 119.88
+4603.1444 15.7948
+4607.9347 11.4279
+4608.6218 11.6267
+4609.372 7.02327
+4609.5673 1221.28
+4611.2433 5.10087
+4612.5435 4.90891
+4613.6044 17.7617
+4615.024 10.0449
+4615.3351 17.7089
+4619.4794 6.41386
+4621.1629 18.3723
+4624.314 6.47323
+4627.2973 9.49829
+4628.1991 16.5995
+4628.4409 35.8668
+4631.6382 7.29758
+4633.7628 18.198
+4637.233 85.9085
+4637.629 3.94716
+4638.6849 26.9969
+4640.0462 6.61454
+4641.2418 8.27625
+4644.7072 4.26447
+4646.6861 10.9471
+4647.2546 26.3496
+4650.2343 17.455
+4650.9195 4.13841
+4651.9895 6.17074
+4655.2088 5.39692
+4657.9041 836.603
+4659.5698 5.48127
+4663.2026 61.9759
+4666.0027 16.1667
+4666.5155 40.8868
+4666.7985 47.1163
+4668.172 81.6867
+4669.9853 117.237
+4673.0535 5.4883
+4673.6621 545.703
+4676.0555 82.2877
+4676.133 4.29782
+4680.2377 5.75225
+4680.646 5.52363
+4683.353 30.5339
+4686.1946 153.856
+4689.2523 41.9807
+4690.6219 49.3453
+4691.3452 8.5765
+4691.6354 46.0446
+4694.0914 54.4282
+4694.9323 12.2334
+4695.0381 168.923
+4695.4542 13.1503
+4700.7837 4.93863
+4702.3161 93.1766
+4703.3581 10.2278
+4703.9898 817.83
+4705.7606 21.5926
+4706.2383 3.97915
+4708.294 31.0609
+4708.5366 9.58955
+4712.4814 53.4481
+4712.8408 15.1739
+4715.4308 4.88338
+4719.9888 4.53178
+4720.4586 13.3998
+4721.2765 33.1405
+4721.591 22.633
+4722.0886 44.4139
+4723.4398 734.51
+4723.784 64.7118
+4726.8705 1064.39
+4728.133 6.98479
+4729.1293 92.2413
+4729.8795 7.27099
+4730.6614 5.23277
+4732.0532 180.301
+4735.9079 506.782
+4739.6764 45.7522
+4740.5292 63.9671
+4740.9585 48.6428
+4741.2967 14.9477
+4742.1174 17.1227
+4742.5535 5.73107
+4743.6923 23.7768
+4745.3322 16.1008
+4746.8239 3.38238
+4748.04 4.84055
+4749.2002 121.145
+4749.7935 9.05726
+4749.9713 47.3657
+4752.4141 99.979
+4757.2196 10.6699
+4761.1101 8.57456
+4764.3463 78.5935
+4764.8646 2661.86
+4765.5959 28.4925
+4766.6006 135.507
+4768.6814 12.8589
+4773.241 13.2449
+4774.2607 35.9558
+4775.313 14.2325
+4775.794 17.3488
+4777.1933 14.2463
+4778.297 563.167
+4779.7306 15.7418
+4780.422 10.1095
+4780.7505 8.35713
+4781.526 6.25546
+4782.7608 11.4082
+4783.8617 23.7332
+4784.0396 9.64609
+4786.1434 3.9864
+4786.531 31.3448
+4787.1476 53.7384
+4789.3868 271.327
+4792.0819 9.75035
+4793.2446 20.983
+4795.9131 29.3319
+4800.1741 10.6862
+4805.6064 6.48416
+4806.0205 1231.64
+4808.1351 498.92
+4809.6161 150.796
+4812.3755 16.1879
+4813.0075 18.7566
+4813.7204 51.2083
+4813.8963 88.0909
+4817.0206 7.78436
+4818.6458 7.23017
+4819.193 19.4464
+4819.5674 3.51181
+4820.4649 13.8117
+4820.8847 93.4916
+4821.2631 4.19723
+4821.385 7.22092
+4821.5878 111.817
+4821.8592 4.93727
+4822.8548 294.935
+4823.1692 8.09662
+4823.6058 110.695
+4823.9967 56.2398
+4826.7004 149.96
+4829.7973 18.3661
+4831.1225 423.617
+4831.5975 126.894
+4832.8029 71.37
+4833.8317 4.28069
+4840.4744 26.8947
+4840.8492 521.946
+4843.9395 65.1344
+4844.1653 12.8414
+4845.1626 32.4845
+4847.3208 7.04908
+4847.8095 577.432
+4847.9985 17.7145
+4848.3646 300.478
+4849.1429 27.1511
+4849.8617 28.7808
+4850.4408 55.1639
+4852.8696 40.0447
+4858.3344 15.3578
+4861.2167 44.0273
+4861.7173 34.9526
+4863.1724 227.657
+4865.4775 441.647
+4865.91 91.5246
+4867.5558 12.5734
+4868.5254 4.98434
+4868.8814 47.874
+4871.289 48.0372
+4872.0309 6.36617
+4872.9169 292.235
+4874.3662 116.992
+4876.2611 46.5136
+4876.4923 19.5098
+4877.0035 8.01475
+4877.8121 13.1595
+4878.0106 127.864
+4878.7344 565.658
+4879.3501 9.17893
+4879.8635 2699.79
+4880.8509 7.95572
+4881.2046 54.3629
+4881.8531 10.9968
+4882.2455 36.8137
+4887.9516 42.1399
+4889.0422 412.134
+4889.4903 19.0889
+4890.4582 20.1678
+4892.7592 24.1195
+4893.4436 15.9623
+4894.6906 16.6734
+4894.957 1132.26
+4898.4563 4.74024
+4898.8044 4.04661
+4899.2441 20.2471
+4902.0556 68.2129
+4902.7943 19.836
+4904.7533 127.776
+4907.2093 6.81641
+4909.8328 4.27015
+4910.1576 51.7022
+4910.7929 26.4061
+4911.3787 23.3431
+4912.5293 30.4426
+4919.8157 381.301
+4920.0226 45.5284
+4921.6134 25.5825
+4922.938 7.03937
+4924.4223 10.6983
+4925.9501 11.3703
+4927.298 9.44349
+4927.7803 74.3743
+4929.086 12.2534
+4929.9754 5.72947
+4932.0473 11.9265
+4933.2136 214.685
+4933.8521 14.704
+4936.7746 57.1255
+4937.8295 95.4003
+4938.5071 19.9503
+4939.2703 6.0869
+4939.6422 301.318
+4941.416 5.25372
+4943.0642 168.915
+4945.4614 456.276
+4946.6637 44.9907
+4947.5751 22.4592
+4949.4037 4.26514
+4950.2513 13.8656
+4950.6264 33.1245
+4952.6953 11.0807
+4954.6546 11.0296
+4955.1204 7.72677
+4956.7647 15.5955
+4957.2963 8.85339
+4958.0988 9.11654
+4960.4279 14.3728
+4961.726 24.3599
+4963.1881 7.4831
+4963.7664 6.96135
+4964.1157 8.74783
+4964.2533 4.31933
+4965.0812 1356.08
+4965.7315 71.1095
+4968.7556 21.4262
+4970.0766 30.2119
+4972.1597 95.446
+4973.398 8.93496
+4975.9439 3.991
+4980.1859 57.2975
+4980.9513 11.0146
+4982.4875 99.3208
+4985.3725 251.563
+4985.9469 6.5111
+4987.1505 11.3022
+4987.5582 9.65165
+4989.3086 78.4742
+4989.9557 6.99063
+4992.1257 6.07045
+4992.6372 5.15659
+4993.7488 38.5607
+4994.1058 11.8688
+4997.7993 3.54205
+4999.9398 6.91979
+5000.2516 5.5321
+5002.0972 624.716
+5003.6785 99.5414
+5009.3344 351.003
+5009.9424 7.14161
+5010.9532 6.36927
+5011.4774 4.75069
+5015.8892 19.6963
+5016.5356 11.4415
+5017.1673 698.923
+5017.5084 36.3313
+5017.6235 10.7912
+5018.0568 8.5097
+5019.3336 8.20173
+5019.8062 151.717
+5020.5335 4.04989
+5021.2531 4.14241
+5022.0051 33.1907
+5022.1698 11.5854
+5023.7059 12.9307
+5026.9532 3.96646
+5028.6556 325.768
+5029.0102 9.39148
+5029.6295 5.87679
+5029.8916 84.2817
+5035.3516 2.98555
+5039.2303 360.229
+5039.5265 19.1678
+5040.5393 6.15313
+5040.6805 10.461
+5041.1224 18.8219
+5041.5971 6.75403
+5043.5131 9.80939
+5044.7195 678.079
+5045.2481 89.4546
+5046.3518 13.9731
+5046.6372 8.14425
+5047.0434 191.836
+5048.8139 30.8297
+5048.9354 35.1897
+5049.5543 8.69854
+5049.796 240.343
+5050.7842 203.828
+5051.8887 76.8459
+5054.1776 19.4488
+5055.3473 69.0541
+5056.5358 3.91915
+5057.9866 43.4318
+5059.8611 118.26
+5060.0871 36.3014
+5061.2199 9.60687
+5061.6562 41.2795
+5062.0395 408.861
+5062.9325 18.4902
+5063.5157 40.647
+5063.9988 33.0461
+5064.602 265.528
+5064.9463 369.346
+5066.1355 57.7505
+5066.7773 29.6794
+5067.1371 24.0714
+5067.9737 968.248
+5069.3384 10.6411
+5073.0746 9.28652
+5074.6637 4.94064
+5077.7903 221.971
+5081.4462 36.1646
+5084.6987 5.1124
+5084.9935 20.9643
+5087.1131 7.51694
+5090.0513 10.254
+5090.5148 191.383
+5095.0639 20.4182
+5096.4848 193.562
+5098.0432 32.5449
+5098.932 11.7805
+5100.6211 223.131
+5101.1299 43.7184
+5101.3395 2.88295
+5109.7207 5.92458
+5110.5493 7.00854
+5110.8709 7.06356
+5111.2781 56.3139
+5115.0477 467.92
+5116.3016 4.97391
+5117.2923 10.4614
+5118.205 18.4299
+5122.4995 18.2438
+5125.4895 15.4357
+5125.9502 109.096
+5126.279 6.87144
+5128.4906 54.4392
+5130.2338 4.79469
+5131.0699 7.36472
+5133.1051 6.47732
+5134.746 110.644
+5136.121 4.46243
+5137.4733 29.6407
+5140.7736 45.7047
+5141.7827 335.653
+5143.2673 55.4042
+5143.9165 128.674
+5145.3083 191.136
+5146.0544 8.24467
+5148.2087 18.9162
+5149.207 9.03956
+5150.7064 4.93751
+5151.612 422.456
+5151.858 9.02414
+5154.243 378.949
+5158.6059 1084.03
+5159.4574 5.80872
+5159.6246 14.2469
+5160.7163 83.4544
+5161.5396 95.0381
+5162.2892 207.227
+5162.7493 10.0091
+5163.4604 233.684
+5165.066 4.84967
+5165.7728 37.3169
+5166.6553 7.52561
+5168.9225 21.3244
+5170.2227 3.38792
+5173.0794 17.3811
+5173.6757 25.3887
+5174.2217 4.48851
+5174.7988 12.5656
+5175.3248 46.0539
+5175.9115 14.3605
+5176.2292 30.8526
+5176.4037 24.8192
+5176.961 228.93
+5177.6228 30.8468
+5178.4726 6.05432
+5180.7209 3.79134
+5182.5269 9.67873
+5183.0965 6.48001
+5183.9896 13.7039
+5184.4541 16.2505
+5186.4132 20.5776
+5187.1693 9.6076
+5187.3374 80.8798
+5187.4583 11.2173
+5187.7491 277.819
+5190.872 19.9127
+5191.3463 3.51862
+5193.8256 39.5266
+5193.9316 6.88467
+5194.4576 80.7526
+5197.2361 7.67376
+5198.7999 130.392
+5199.1637 605.975
+5201.9993 8.19701
+5203.8474 10.9719
+5205.1522 29.1488
+5206.6722 4.99346
+5207.8015 5.52073
+5209.7246 76.0411
+5210.492 6.73147
+5211.2343 435.898
+5213.3492 103.41
+5214.7867 14.0768
+5215.6977 18.0725
+5216.6014 42.831
+5218.5271 11.4526
+5219.1137 181.613
+5220.7068 11.6886
+5220.9301 24.6622
+5221.271 89.5286
+5221.8618 5.7256
+5228.2246 28.1778
+5228.9951 6.50876
+5231.1597 1030.95
+5233.2254 25.7881
+5234.1087 13.1727
+5238.8137 112.491
+5239.552 108.328
+5240.1968 32.4576
+5245.3797 5.56783
+5247.1977 20.5721
+5247.6548 133.966
+5250.8771 17.9253
+5254.2588 12.0692
+5254.4699 12.9062
+5258.3602 414.673
+5260.1041 61.1699
+5261.472 18.3194
+5264.7865 9.64064
+5266.7145 215.526
+5267.1597 5.01696
+5272.927 29.2593
+5273.1317 9.19802
+5274.1188 69.5206
+5277.1433 14.8433
+5277.5002 22.9387
+5278.1166 40.8362
+5280.3546 7.23929
+5281.069 53.1053
+5281.6257 8.22585
+5286.887 11.7575
+5289.8929 4.02498
+5294.3997 36.0776
+5295.7567 40.0391
+5296.2811 50.8073
+5297.7431 67.7718
+5298.2824 62.5283
+5300.5234 34.4504
+5301.3968 4.29366
+5303.4845 9.75274
+5305.699 10.4784
+5306.9963 9.24972
+5307.4632 10.0681
+5309.5983 4.18063
+5310.2665 23.0678
+5312.0063 273.93
+5312.5288 106.566
+5312.9045 58.8661
+5315.2278 5.62052
+5317.4983 42.7434
+5320.7757 8.43051
+5325.1438 10.8274
+5326.2774 26.3838
+5326.9756 139.316
+5329.3745 10.7514
+5329.7088 6.99339
+5330.0804 20.0113
+5340.4982 7.20135
+5343.5812 484.175
+5345.3079 4.87436
+5347.9714 27.1696
+5351.1294 27.9601
+5351.8297 7.16019
+5354.6016 10.7326
+5355.6465 8.3333
+5358.3826 13.1422
+5359.8269 8.53931
+5360.1501 53.6849
+5361.1593 12.9852
+5362.5751 24.8142
+5369.2819 22.9514
+5369.447 24.2969
+5370.7096 21.5016
+5372.7067 15.7941
+5373.4943 8.9765
+5374.8219 65.5878
+5375.3587 20.0024
+5375.7688 5.72737
+5376.1304 6.71117
+5376.7809 19.7164
+5378.8349 44.3788
+5379.1105 92.3125
+5382.9276 12.0274
+5384.0356 10.4107
+5384.3682 4.78892
+5386.6148 178.392
+5387.8266 3.94538
+5388.0507 10.9853
+5390.4375 108.219
+5392.5726 25.4667
+5393.6047 3.81093
+5393.9744 31.1239
+5394.7608 173.502
+5397.5144 13.2397
+5398.2047 5.77056
+5398.7015 46.0962
+5398.9219 28.3115
+5399.1745 46.9412
+5399.6218 7.48724
+5400.1455 6.34557
+5402.6048 15.4506
+5403.1993 8.38744
+5407.6535 119.875
+5410.4885 8.56801
+5410.7687 119.322
+5415.5226 45.915
+5417.4926 350.868
+5419.1234 4.06181
+5421.3517 35.6623
+5421.8359 8.0436
+5424.0079 60.1299
+5425.6783 33.8477
+5426.4141 11.7312
+5431.112 141.196
+5435.1229 4.50729
+5435.892 13.6085
+5437.3877 10.8337
+5438.1488 7.10681
+5439.9891 29.8005
+5440.6015 21.9521
+5442.2477 6.35966
+5443.2326 17.9646
+5443.6938 7.1111
+5447.1536 13.7033
+5449.4733 14.3251
+5451.652 105.803
+5452.2187 132.988
+5454.3015 13.2046
+5457.2985 12.3395
+5457.7973 5.71497
+5461.7402 6.77458
+5462.6057 6.94785
+5464.2123 83.3744
+5467.1465 7.21199
+5469.0844 4.43358
+5470.7591 40.2228
+5473.4541 14.1218
+5479.0751 15.313
+5484.1367 8.1378
+5488.6363 14.6626
+5488.9151 8.21965
+5489.0833 5.74745
+5490.1267 4.08873
+5492.3362 11.3461
+5493.2043 35.802
+5494.3327 41.7922
+5495.8738 274.697
+5496.1382 96.6002
+5498.175 9.52024
+5499.2567 644.75
+5499.6461 26.91
+5500.3125 4.43996
+5501.2784 15.8737
+5501.9367 7.14131
+5504.3018 114.524
+5506.1152 55.5132
+5506.9034 5.49701
+5507.5376 26.0302
+5507.7676 4.82248
+5508.5564 12.7312
+5509.9938 158.421
+5514.8721 102.678
+5518.9841 15.0484
+5519.2615 5.45086
+5524.5824 27.0831
+5524.957 24.5668
+5527.3022 13.4549
+5528.2376 12.4776
+5533.2452 8.99002
+5534.4989 7.95551
+5535.9572 5.60565
+5537.1194 5.84807
+5537.5563 15.7154
+5538.608 23.1786
+5539.2614 388.735
+5539.9104 65.4704
+5539.933 78.3845
+5541.1446 17.2482
+5541.5833 12.5174
+5541.94 3.13988
+5542.8911 66.4701
+5546.122 4.80963
+5548.1774 189.504
+5551.3701 12.0377
+5552.6221 14.6156
+5554.0219 29.5953
+5554.1319 6.26146
+5555.5278 8.09759
+5557.0459 149.272
+5558.3407 206.791
+5558.702 263.229
+5559.6677 11.9343
+5559.8912 48.3765
+5567.9986 17.0399
+5570.9235 11.7199
+5571.1937 136.006
+5572.4866 145.08
+5573.3535 302.095
+5576.2053 63.5346
+5577.689 25.2927
+5578.5506 4.47167
+5579.3594 168.85
+5580.0795 14.7368
+5580.7436 5.80429
+5581.8644 7.58005
+5582.8565 6.70252
+5583.7635 13.0022
+5587.026 752.61
+5588.7244 16.4156
+5593.2794 2.48508
+5593.6193 28.1029
diff --git a/noao/artdata/mkexamples/galcluster.cl b/noao/artdata/mkexamples/galcluster.cl
new file mode 100644
index 00000000..8a6e1e7f
--- /dev/null
+++ b/noao/artdata/mkexamples/galcluster.cl
@@ -0,0 +1,39 @@
+# GALCLUSTER - Galaxy cluster
+
+file	image, dat
+
+image = s1
+dat = mktemp ("art")
+
+gallist (dat, 100, interactive=no, spatial="hubble", xmin=1., xmax=512.,
+    ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF, core_radius=50.,
+    base=0., sseed=i+1, luminosity="schecter", minmag=-7., maxmag=0.,
+    mzero=15., power=0.6, alpha=-1.24, mstar=-21.41, lseed=i+1, egalmix=0.8,
+    ar=0.7, eradius=20., sradius=1., absorption=1.2, z=0.05, sfile="",
+    nssample=100, sorder=10, lfile="", nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, dbinsize=0.5, ebinsize=0.1, pbinsize=20.,
+    graphics="stdgraph", cursor="")
+
+gallist (dat, 500, interactive=no, spatial="uniform", xmin=1.,
+    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
+    core_radius=50., base=0., sseed=i+1, luminosity="powlaw", minmag=-7.,
+    maxmag=0., mzero=15., power=0.6, alpha=-1.24, mstar=-21.41, lseed=i+1,
+    egalmix=0.4, ar=0.7, eradius=20., sradius=1., absorption=1.2, z=0.05,
+    sfile="", nssample=100, sorder=10, lfile="", nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, dbinsize=0.5, ebinsize=0.1, pbinsize=20.,
+    graphics="stdgraph", cursor="")
+
+starlist (dat, 100, "", "", interactive=no, spatial="uniform", xmin=1.,
+    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
+    core_radius=30., base=0., sseed=i, luminosity="powlaw", minmag=-7.,
+    maxmag=0., mzero=-4., power=0.6, alpha=0.74, beta=0.04, delta=0.294,
+    mstar=1.28, lseed=i, nssample=100, sorder=10, nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, graphics="stdgraph", cursor="")
+
+mkobjects (image, output="", ncols=512, nlines=512,
+    title="Example artificial galaxy cluster", header="artdata$stdheader.dat",
+    background=1000., objects=dat, xoffset=0., yoffset=0., star="moffat",
+    radius=1.0, beta=2.5, ar=1., pa=0., distance=1., exptime=1., magzero=7.,
+    gain=3., rdnoise=10., poisson=yes, seed=j, comments=b1)
+
+delete (dat, verify=no)
diff --git a/noao/artdata/mkexamples/galfield.cl b/noao/artdata/mkexamples/galfield.cl
new file mode 100644
index 00000000..168972fe
--- /dev/null
+++ b/noao/artdata/mkexamples/galfield.cl
@@ -0,0 +1,30 @@
+# GALFIELD - Uniform galaxy field
+
+file	image, dat
+
+image = s1
+dat = mktemp ("art")
+
+gallist (dat, 1000, interactive=no, spatial="uniform", xmin=1.,
+    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
+    core_radius=50., base=0., sseed=i+1, luminosity="powlaw", minmag=-7.,
+    maxmag=0., mzero=15., power=0.45, alpha=-1.24, mstar=-21.41, lseed=i+1,
+    egalmix=0.4, ar=0.7, eradius=10., sradius=1., absorption=1.2, z=0.05,
+    sfile="", nssample=100, sorder=10, lfile="", nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, dbinsize=0.5, ebinsize=0.1, pbinsize=20.,
+    graphics="stdgraph", cursor="")
+
+starlist (dat, 100, "", "", interactive=no, spatial="uniform", xmin=1.,
+    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
+    core_radius=30., base=0., sseed=i, luminosity="powlaw", minmag=-7.,
+    maxmag=0., mzero=-4., power=0.6, alpha=0.74, beta=0.04, delta=0.294,
+    mstar=1.28, lseed=i, nssample=100, sorder=10, nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, graphics="stdgraph", cursor="")
+
+mkobjects (image, output="", ncols=512, nlines=512,
+    title="Example artificial galaxy field", header="artdata$stdheader.dat",
+    background=400., objects=dat, xoffset=0., yoffset=0., star="moffat",
+    radius=1.0, beta=2.5, ar=1., pa=0., distance=1., exptime=1., magzero=5.5,
+    gain=5., rdnoise=10., poisson=yes, seed=j, comments=b1)
+
+delete (dat, verify=no)
diff --git a/noao/artdata/mkexamples/globular.cl b/noao/artdata/mkexamples/globular.cl
new file mode 100644
index 00000000..f70f3f1e
--- /dev/null
+++ b/noao/artdata/mkexamples/globular.cl
@@ -0,0 +1,23 @@
+# GLOBULAR - Globular cluster
+
+file	image, dat
+
+image = s1
+dat = mktemp ("art")
+
+starlist (dat, 5000, "", "", interactive=no, spatial="hubble",
+    xmin=1., xmax=512., ymin=1., ymax=512., xcenter=INDEF,
+    ycenter=INDEF, core_radius=30., base=0., sseed=i,
+    luminosity="bands", minmag=-7., maxmag=0., mzero=-4., power=0.6,
+    alpha=0.74, beta=0.04, delta=0.294, mstar=1.28, lseed=i,
+    nssample=100, sorder=10, nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, graphics="stdgraph", cursor="")
+
+mkobjects (image, output="", ncols=512, nlines=512,
+    title="Example artificial globular cluster",
+    header="artdata$stdheader.dat", background=1000., objects=dat,
+    xoffset=0., yoffset=0., star="moffat", radius=1.0, beta=2.5,
+    ar=1., pa=0., distance=1., exptime=1., magzero=7.,
+    gain=3., rdnoise=10., poisson=yes, seed=j, comments=b1)
+
+delete (dat, verify=no)
diff --git a/noao/artdata/mkexamples/henear.cl b/noao/artdata/mkexamples/henear.cl
new file mode 100644
index 00000000..cec48b98
--- /dev/null
+++ b/noao/artdata/mkexamples/henear.cl
@@ -0,0 +1,20 @@
+# henear - Helium-Neon-Argon spectrum (uncalibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$archdr.dat"
+
+for (k=1; k<=i; k+=1) {
+    mk1dspec (out, output="", ap=k, rv=0., z=no, ncols=512, naps=i,
+	wstart=4209.0+k, wend=7361.7+k, title="Helium-Neon-Argon Arc Example",
+	header="", continuum=0.5, slope=0., temperature=0.,
+	lines="mkexamples$henear2.dat", profile="gaussian", gfwhm=14,
+	comments=b1)
+}
+
+mkheader (out, hdr, append=no, verbose=no)
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/henear1.dat b/noao/artdata/mkexamples/henear1.dat
new file mode 100644
index 00000000..6dfd4bdf
--- /dev/null
+++ b/noao/artdata/mkexamples/henear1.dat
@@ -0,0 +1,56 @@
+# HENEAR TRUE WAVELENGTHS
+4259.361 456.992 
+4300.4 151.516 
+4426.01 304.683 
+4471.477 808.484 
+4510.733 197.966 
+4545.08 371.968 
+4579.39 334.567 
+4657.94 269.339 
+4713.143 365.173 
+4764.89 628.846 
+4806.07 308.678 
+4879.9 418.101 
+4921.929 644.809 
+4965.12 234.898 
+5015.675 2573.67 
+5187.746 292.801 
+5221.27 193.722 
+5400.562 326.122 
+5495.872 454.027 
+5572.548 335.12 
+5606.732 326.878 
+5650.703 283.252 
+5748.299 224.797 
+5764.418 269.445 
+5852.4878 20334.6 
+5875.618 15259.1 
+5944.8342 2891.43 
+6029.9971 1970.72 
+6074.3377 2573.29 
+6096.163 4871.86 
+6143.0623 7662.59 
+6163.5939 2296.45 
+6217.2813 1579.43 
+6266.495 4593.41 
+6304.7892 1604.55 
+6334.4279 3425.67 
+6382.9914 6427.61 
+6402.246 10708.7 
+6506.5279 6053.5 
+6532.8824 2271.5 
+6598.9529 3574.04 
+6678.2 30911.1 
+6717.0428 4317.68 
+6752.832 3451.33 
+6871.29 2877.43 
+6929.468 5863.05 
+6965.43 35198.7 
+7032.4127 14128.7 
+7065.188 32334.6 
+7107.496 1037.61 
+7147.041 3887.49 
+7173.939 1277.6 
+7206.986 1433.43 
+7245.167 4155.93 
+7281.349 11982.2 
diff --git a/noao/artdata/mkexamples/henear1d.cl b/noao/artdata/mkexamples/henear1d.cl
new file mode 100644
index 00000000..2595ad6a
--- /dev/null
+++ b/noao/artdata/mkexamples/henear1d.cl
@@ -0,0 +1,4 @@
+# HENEAR1D - 1D Helium-Neon-Argon Arc Spectrum
+
+i = 1
+cl < "mkexamples$henear.cl"
diff --git a/noao/artdata/mkexamples/henear2.dat b/noao/artdata/mkexamples/henear2.dat
new file mode 100644
index 00000000..a8fb3236
--- /dev/null
+++ b/noao/artdata/mkexamples/henear2.dat
@@ -0,0 +1,56 @@
+# HENEAR WAVELENGTHS WITH 10A QUADRATIC DISPERSION FUNCTION
+4268.73 456.992 
+4309.26 151.516 
+4433.41 304.683 
+4478.38 808.484 
+4517.22 197.966 
+4551.22 371.968 
+4585.18 334.567 
+4662.99 269.339 
+4717.69 365.173 
+4769.00 628.846 
+4809.85 308.678 
+4883.11 418.101 
+4924.84 644.809 
+4967.74 234.898 
+5017.97 2573.67 
+5189.10 292.801 
+5222.47 193.722 
+5401.09 326.122 
+5496.16 454.027 
+5572.69 335.12 
+5606.83 326.878 
+5650.75 283.252 
+5748.30 224.797 
+5764.42 269.445 
+5852.53 20334.6 
+5875.68 15259.1 
+5944.98 2891.43 
+6030.31 1970.72 
+6074.76 2573.29 
+6096.65 4871.86 
+6143.69 7662.59 
+6164.29 2296.45 
+6218.17 1579.43 
+6267.58 4593.41 
+6306.05 1604.55 
+6335.83 3425.67 
+6384.63 6427.61 
+6403.99 10708.7 
+6508.88 6053.5 
+6535.40 2271.5 
+6601.92 3574.04 
+6681.75 30911.1 
+6720.90 4317.68 
+6756.98 3451.33 
+6876.49 2877.43 
+6935.22 5863.05 
+6971.54 35198.7 
+7039.22 14128.7 
+7072.35 32334.6 
+7115.13 1037.61 
+7155.13 3887.49 
+7182.34 1277.6 
+7215.79 1433.43 
+7254.44 4155.93 
+7291.07 11982.2 
diff --git a/noao/artdata/mkexamples/heneardc.cl b/noao/artdata/mkexamples/heneardc.cl
new file mode 100644
index 00000000..8e1ab75f
--- /dev/null
+++ b/noao/artdata/mkexamples/heneardc.cl
@@ -0,0 +1,18 @@
+# heneardc - Helium-Neon-Argon spectrum (calibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$archdr.dat"
+
+for (k=1; k<=i; k+=1) {
+    mk1dspec (out, output="", ap=k, rv=0., z=no, ncols=512, naps=i,
+	wstart=4210.0, wend=7362.7, title="Helium-Neon-Argon Arc Example",
+	header=hdr, continuum=0.5, slope=0., temperature=0.,
+	lines="mkexamples$henear1.dat", profile="gaussian",
+	gfwhm=14, comments=b1)
+}
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=10., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/longslit.cl b/noao/artdata/mkexamples/longslit.cl
new file mode 100644
index 00000000..5be43bcd
--- /dev/null
+++ b/noao/artdata/mkexamples/longslit.cl
@@ -0,0 +1,77 @@
+# LONGSLIT - Long slit example
+
+real	w1, w2
+file	out, obj, bkg, dat, hdr
+
+out = s1
+hdr = s2
+obj = mktemp ("art")
+bkg = mktemp ("art")
+dat = mktemp ("art")
+
+w1 = 4209 + i
+w2 = 7361 + i
+
+if (i == 2) {		# Featureless hot continuum + sky
+    mk1dspec (obj, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=10000., lines="", nlines=0, peak=-0.5, profile="gaussian",
+	gfwhm=23.5, seed=1, comments=b1)
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=5800., lines="", nlines=20, peak=1., profile="gaussian",
+	gfwhm=12, seed=1, comments=b1)
+    print (obj, " 1 gauss 3 0 50 .002", > dat)
+    print (bkg, " 10 slit 90 0 50 0", >> dat)
+} else if (i == 3) {	# Sky only
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=5800., lines="", nlines=20, peak=1., profile="gaussian",
+	gfwhm=12, seed=1, comments=b1)
+    print (bkg, " 10 slit 90 0 50 0", > dat)
+} else if (i == 4) {	# Flat field
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=8000., lines="", nlines=0, peak=-0.5, profile="gaussian",
+	gfwhm=23.5, seed=1, comments=b1)
+    print (bkg, " 10 slit 90 0 50 0", > dat)
+} else if (i == 5) {	# HE-NE-AR
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=0.5, slope=0.,
+	temperature=0., lines="mkexamples$henear2.dat", profile="gaussian",
+	gfwhm=14, comments=b1)
+    print (bkg, " 100 slit 90 0 50 0", > dat)
+} else if (i == 6) {	# Galaxy absorption line spectrum + sky
+    mk1dspec (obj, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=7000., lines="", nlines=50, peak=-0.5, profile="gaussian",
+	gfwhm=23.5, seed=i+1, comments=b1)
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=5800., lines="", nlines=20, peak=1., profile="gaussian",
+	gfwhm=12, seed=1, comments=b1)
+    print (obj, " 10 gauss 30 0 50 .002", > dat)
+    print (bkg, " 10 slit 200 0 50 0", >> dat)
+} else {		# Star absorption line spectrum + sky
+    mk1dspec (obj, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=7000., lines="", nlines=50, peak=-0.5, profile="gaussian",
+	gfwhm=23.5, seed=i+1, comments=b1)
+    mk1dspec (bkg, output="", ap=1, rv=0., z=no, title="", header="",
+	ncols=512, naps=1, wstart=w1, wend=w2, continuum=1000., slope=0.,
+	temperature=5800., lines="", nlines=20, peak=1., profile="gaussian",
+	gfwhm=12, seed=1, comments=b1)
+    print (obj, " 1 gauss 3 0 50 .002", > dat)
+    print (bkg, " 10 slit 90 0 50 0", >> dat)
+}
+
+mk2dspec (out, output="", model=dat, comments=b1, ncols=100, nlines=512,
+    title="Example artificial long slit image",  header=hdr)
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=3., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
+
+imdelete (obj, verify=no, >& "dev$null")
+imdelete (bkg, verify=no, >& "dev$null")
+delete (dat, verify=no, >& "dev$null")
diff --git a/noao/artdata/mkexamples/lsarc.cl b/noao/artdata/mkexamples/lsarc.cl
new file mode 100644
index 00000000..2ca4afdc
--- /dev/null
+++ b/noao/artdata/mkexamples/lsarc.cl
@@ -0,0 +1,5 @@
+# LSARC - Longslit Helium-Neon-Argon Arc Spectrum
+
+i = 5
+s2 = "mkexamples$archdr.dat"
+cl < "mkexamples$longslit.cl"
diff --git a/noao/artdata/mkexamples/lsgal.cl b/noao/artdata/mkexamples/lsgal.cl
new file mode 100644
index 00000000..0c9aaad7
--- /dev/null
+++ b/noao/artdata/mkexamples/lsgal.cl
@@ -0,0 +1,5 @@
+# LSGAL - Long slit galaxy spectrum (uncalibrated)
+
+i = 6
+s2 = "mkexamples$objhdr.dat"
+cl < "mkexamples$longslit.cl"
diff --git a/noao/artdata/mkexamples/lsobj.cl b/noao/artdata/mkexamples/lsobj.cl
new file mode 100644
index 00000000..fe00bf9b
--- /dev/null
+++ b/noao/artdata/mkexamples/lsobj.cl
@@ -0,0 +1,5 @@
+# LSOBJ - Long slit object spectrum (uncalibrated)
+
+i = 1
+s2 = "mkexamples$objhdr.dat"
+cl < "mkexamples$longslit.cl"
diff --git a/noao/artdata/mkexamples/mkexamples.cl b/noao/artdata/mkexamples/mkexamples.cl
new file mode 100644
index 00000000..ed9150e4
--- /dev/null
+++ b/noao/artdata/mkexamples/mkexamples.cl
@@ -0,0 +1,68 @@
+# MKEXAMPLES -- Make ARTDATA examples.
+# The example script files use variable s1 to pass the image name,
+# variable s2 to pass the image header file, variable i to pass an object
+# seed, variable j to pass a noise seed, and variable b1 to pass the comment
+# flag.
+
+procedure mkexamples (name, image)
+
+string	name		{prompt="Example name"}
+string	image		{prompt="Image name"}
+
+int	oseed=1		{prompt="Object seed"}
+int	nseed=1		{prompt="Noise seed"}
+bool	comments=no	{prompt="Add comments to image?"}
+bool	verbose=yes	{prompt="Print operation?"}
+bool	errors=yes	{prompt="Report errors?"}
+bool	list=no		{prompt="List example script file only?"}
+
+begin
+	string	name1, name2, example
+
+
+	# Get and check parameters.
+	if ($nargs < 1 && mode != "h") {
+	    name1 = "mkexamples"
+	    name2 = ""
+	    while (name1 != name2) {
+		example = "mkexamples$" // name1 // ".men"
+		if (!access (example))
+		    break
+		type (example)
+		name2 = name1
+		name1 = name
+	    }
+	    if (name1 == name2)
+		return
+	} else
+	    name1 = name
+
+	example = "mkexamples$" // name1 // ".cl"
+	if (!access (example)) {
+	    if (errors)
+	        error (2, "Unknown example " // name1) 
+	    return
+	}
+
+	# Make or list the example.
+	if (list)
+	    page (example)
+	else {
+	    s1 = image 
+	    s2 = "artdata$stdheader.dat"
+	    i = oseed
+	    j = nseed
+	    b1 = comments
+	    if (s1 == "" || name1 == "")
+		return
+	    if ((access (s1) || access (s1//"."//envget("imtype")))) {
+		if (errors)
+		    error (1, "Image " // s1 // " already exists")
+		return
+	    }
+
+	    if (verbose)
+                print ("Creating example ", name1, " in image ", s1, " ...")
+	    cl (< example)
+	}
+end
diff --git a/noao/artdata/mkexamples/mkexamples.men b/noao/artdata/mkexamples/mkexamples.men
new file mode 100644
index 00000000..feb89380
--- /dev/null
+++ b/noao/artdata/mkexamples/mkexamples.men
@@ -0,0 +1,10 @@
+	    		MKEXAMPLE Menu
+
+     galcluster - Galaxy cluster
+       globular - Globular cluster
+       galfield - Galaxy field
+      starfield - Starfield
+
+       onedspec - Menu of one dimensional spectra
+       twodspec - Menu of two dimensional spectra
+     threedspec - Menu of three dimensional spectra
diff --git a/noao/artdata/mkexamples/multifiber.cl b/noao/artdata/mkexamples/multifiber.cl
new file mode 100644
index 00000000..af48c4e7
--- /dev/null
+++ b/noao/artdata/mkexamples/multifiber.cl
@@ -0,0 +1,31 @@
+# MULTIFIBER - Multifiber example
+
+file	out, obj, dat
+
+out = s1
+obj = mktemp ("art")
+dat = mktemp ("art")
+
+mk1dspec (obj, output="", ap=1, rv=0., z=no, title="", header="",
+    ncols=512, naps=1, wstart=4000., wend=8000., continuum=1000., slope=0.,
+    temperature=5700., lines="", nlines=50, peak=-0.5, profile="gaussian",
+    gfwhm=24, seed=i, comments=b1)
+
+print (obj, " .1 gauss 3 0 10.1 .002", > dat)
+print (obj, " .2 gauss 3 0 20.2 .002", >> dat)
+print (obj, " .3 gauss 3 0 30.3 .002", >> dat)
+print (obj, " .4 gauss 3 0 40.4 .002", >> dat)
+print (obj, " .5 gauss 3 0 50.5 .002", >> dat)
+print (obj, " .6 gauss 3 0 60.6 .002", >> dat)
+print (obj, " .7 gauss 3 0 70.7 .002", >> dat)
+print (obj, " .8 gauss 3 0 80.8 .002", >> dat)
+print (obj, " .9 gauss 3 0 90.9 .002", >> dat)
+mk2dspec (out, output="", model=dat, comments=b1, ncols=100, nlines=512,
+    title="Example artificial multifiber image", header="artdata$stdheader.dat")
+
+mknoise (out, output="", title="", header="", ncols=512, nlines=512,
+    background=0., gain=1., rdnoise=3., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
+
+imdelete (obj, verify=no)
+delete (dat, verify=no)
diff --git a/noao/artdata/mkexamples/objhdr.dat b/noao/artdata/mkexamples/objhdr.dat
new file mode 100644
index 00000000..eabebc5f
--- /dev/null
+++ b/noao/artdata/mkexamples/objhdr.dat
@@ -0,0 +1,19 @@
+OBSERVER= 'IRAF              '  /  observers
+OBSERVAT= 'KPNO              '  /  observatory
+INSTRUME= 'IRAF ARTDATA      '  /  instrument
+TELESCOP= 'kpcdf             '  /  telescope name
+DETECTOR= 'te1k              '  /  detector
+EXPTIME =                1200.  /  actual integration time
+DARKTIME=                1200.  /  total elapsed time
+IMAGETYP= 'object            '  /  object, dark, bias, etc.
+DATE-OBS= '1991-11-26T12:19:55.00' /  date (dd/mm/yy) of obs.
+UT      = '12:19:55.00       '  /  universal time
+ST      = '09:13:15.00       '  /  sidereal time
+RA      = '06:37:02.00       '  /  right ascension
+DEC     = '06:08:52.00       '  /  declination
+EPOCH   =               1991.9  /  epoch of ra and dec
+ZD      = '44.580            '  /  zenith distance
+AIRMASS =                   0.  /  airmass
+APERTURE= '250micron slit    '  /  aperture
+GRATPOS =               4624.3  /  grating position
+CCDPROC = 'Nov 26  5:44 CCD processing done'
diff --git a/noao/artdata/mkexamples/onedspec.men b/noao/artdata/mkexamples/onedspec.men
new file mode 100644
index 00000000..903eb45a
--- /dev/null
+++ b/noao/artdata/mkexamples/onedspec.men
@@ -0,0 +1,10 @@
+	    		ONEDSPEC Menu
+
+         henear - Helium-neon-argon spectrum (uncalibrated)
+       heneardc - Helium-neon-argon spectrum (calibrated)
+
+          ecarc - Echelle thorium-argon spectrum (uncalibrated)
+        ecarcdc - Echelle thorium-argon spectrum (calibrated)
+
+       spectrum - Absorption spectrum (calibrated)
+        echelle - Echelle absorption spectrum (calibrated)
diff --git a/noao/artdata/mkexamples/spectrum.cl b/noao/artdata/mkexamples/spectrum.cl
new file mode 100644
index 00000000..a5045b8f
--- /dev/null
+++ b/noao/artdata/mkexamples/spectrum.cl
@@ -0,0 +1,19 @@
+# spectrum - Object spectrum (calibrated)
+
+file	out, hdr
+
+out = s1
+hdr = "mkexamples$objhdr.dat"
+
+for (k=1; k<=i; k+=1) {
+    x = k * 100.
+    mk1dspec (out, output="", ap=k, rv=0., z=no, ncols=512, naps=i,
+	wstart=4210.0, wend=7362.7, title="Artificial Spectrum",
+	header=hdr, continuum=x, slope=0., temperature=5700.,
+        lines="", nlines=50, peak=-0.5, profile="gaussian", gfwhm=24,
+	seed=i, comments=b1)
+}
+
+mknoise (out, output="", ncols=512, nlines=512, title="", header="",
+    background=0., gain=1., rdnoise=3., poisson=yes, seed=j, cosrays="",
+    ncosrays=0, energy=30000., radius=0.5, ar=1., pa=0., comments=b1)
diff --git a/noao/artdata/mkexamples/starfield.cl b/noao/artdata/mkexamples/starfield.cl
new file mode 100644
index 00000000..3e1c19b9
--- /dev/null
+++ b/noao/artdata/mkexamples/starfield.cl
@@ -0,0 +1,22 @@
+# STARFIELD - Uniform star field
+
+file	image, dat
+
+image = s1
+dat = mktemp ("art")
+
+starlist (dat, 1000, "", "", interactive=no, spatial="uniform", xmin=1.,
+    xmax=512., ymin=1., ymax=512., xcenter=INDEF, ycenter=INDEF,
+    core_radius=30., base=0., sseed=i, luminosity="powlaw", minmag=-7.,
+    maxmag=0., mzero=-4., power=0.6, alpha=0.74, beta=0.04, delta=0.294,
+    mstar=1.28, lseed=i, nssample=100, sorder=10, nlsample=100, lorder=10,
+    rbinsize=10., mbinsize=0.5, graphics="stdgraph", cursor="")
+
+mkobjects (image, output="", ncols=512, nlines=512,
+    title="Example artificial uniform star field",
+    header="artdata$stdheader.dat",
+    background=1000., objects=dat, xoffset=0., yoffset=0., star="moffat",
+    radius=1.0, beta=2.5, ar=1., pa=0., distance=1., exptime=1., magzero=7.,
+    gain=3., rdnoise=10., poisson=yes, seed=j, comments=b1)
+
+delete (dat, verify=no)
diff --git a/noao/artdata/mkexamples/threedspec.men b/noao/artdata/mkexamples/threedspec.men
new file mode 100644
index 00000000..8aba4892
--- /dev/null
+++ b/noao/artdata/mkexamples/threedspec.men
@@ -0,0 +1,3 @@
+	    		THREEDSPEC Menu
+
+	None available
diff --git a/noao/artdata/mkexamples/twodspec.men b/noao/artdata/mkexamples/twodspec.men
new file mode 100644
index 00000000..08ebff20
--- /dev/null
+++ b/noao/artdata/mkexamples/twodspec.men
@@ -0,0 +1,10 @@
+	    		TWODSPEC Menu
+
+        ecarc2d - Echelle thorium-argon slit spectrum
+        ecobj2d - Echelle object slit spectrum
+
+          lsarc - Long slit helium-neon-argon spectrum
+	  lsobj - Long slit object spectrum
+	  lsgal - Long slit galaxy spectrum
+
+     multifiber - Multifiber spectrum