diff options
| author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
|---|---|---|
| committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
| commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
| tree | 4464880c571602d54f6ae114729bf62a89518057 /noao/obsutil/src/sptime/blazeang.cl | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'noao/obsutil/src/sptime/blazeang.cl')
| -rw-r--r-- | noao/obsutil/src/sptime/blazeang.cl | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/noao/obsutil/src/sptime/blazeang.cl b/noao/obsutil/src/sptime/blazeang.cl new file mode 100644 index 00000000..2c55b3a2 --- /dev/null +++ b/noao/obsutil/src/sptime/blazeang.cl @@ -0,0 +1,24 @@ +procedure blazeang (g, w) + +real g = 316 {prompt="l/mm"} +real w = 7500 {prompt="Blaze wavelength (A)"} +real phi = 46. {prompt="Camera-collimator angle (deg)"} +real m = 1 {prompt="Order"} +real n = 1. {prompt="Index of refraction"} +real prism = 22 {prompt="Prism angle (deg)"} + +begin + real dtor, val + + dtor = 3.14159 / 180. + + if (n <= 1.) { + val = g * w * m / cos (phi/2*dtor) / 2e7 + val = atan2 (val, sqrt (1 - val**2)) / dtor + } else { +# val = g * w * m / 1e7 / (n - 1.) +# val = atan2 (val, sqrt (1 - val**2)) / dtor + val = g * w * m / 1e7 / sin (dtor * prism) + 1 + } + printf ("%.4g\n", val) +end |