function index_ice_266K,wave_in ;----------------------------------------------------------------- ; ; purpose: ; Get the refractive index for ice at -7 C (266K) (Warren, Appl. ; Opt., V23, p1206, 1984) at the specified wavelength. ; The data are stored in the routine for simplicity and ease. ; ; input: ; wave_in..........wavelength, microns ; ; output: ; index_ice_266K...complex refractive index ; ; source: Mark Hervig ; ;----------------------------------------------------------------- numw = 232 ; # wavelengths wave = fltarr(numw) ; wavelength, microns dxr = fltarr(numw) ; real index dxi = fltarr(numw) ; imaginary index waves = $ '4.43000E-02 4.86000E-02 5.39000E-02 6.05000E-02 6.89000E-02 '+$ '7.75000E-02 9.18000E-02 0.107800 0.121000 0.131900 0.139300 '+$ '0.145000 0.148500 0.155000 0.162500 0.173700 0.18210 0.193000 '+$ '0.350000 0.460000 0.470000 0.480000 0.490000 0.50000 0.510000 '+$ '0.520000 0.530000 0.540000 0.550000 0.560000 0.57000 0.580000 '+$ '0.590000 0.600000 0.610000 0.620000 0.630000 0.64000 0.650000 '+$ '0.660000 0.670000 0.680000 0.690000 0.700000 0.71000 0.720000 '+$ '0.730000 0.740000 0.750000 0.760000 0.770000 0.78000 0.790000 '+$ '0.800000 0.810000 0.820000 0.830000 0.840000 0.85000 0.860000 '+$ '0.870000 0.880000 0.890000 0.900000 0.910000 0.92000 0.930000 '+$ '0.940000 0.950000 0.960000 0.970000 0.980000 0.99000 1.00000 '+$ '1.01000 1.02000 1.03000 1.04000 1.05000 1.06000 1.07000 '+$ '1.08000 1.09000 1.10000 1.11000 1.12000 1.13000 1.14000 '+$ '1.15000 1.16000 1.17000 1.18000 1.19000 1.20000 1.21000 '+$ '1.22000 2.37000 2.39000 2.41000 2.43000 2.46000 2.50000 '+$ '2.52000 2.55000 2.56500 2.58000 2.60000 2.62000 2.67500 '+$ '2.72500 2.77800 2.81700 2.83300 2.84900 2.86500 2.88200 '+$ '2.89900 2.91500 2.93300 2.95000 2.96700 2.98500 3.00300 '+$ '3.02100 3.04000 3.05800 3.07700 3.09600 3.11500 3.13500 '+$ '3.15500 3.17500 3.19500 3.21500 3.23600 3.25700 3.27900 '+$ '3.30000 3.32200 3.34500 3.36700 3.39000 3.41300 3.43600 '+$ '3.46000 3.48400 3.50900 3.53400 3.55900 3.62400 3.73200 '+$ '3.77500 3.84700 3.96900 4.09900 4.23900 4.34800 4.38700 '+$ '4.44400 4.50500 4.54700 4.56000 4.58000 4.71900 4.90400 '+$ '5.00000 5.10000 5.20000 5.26300 5.40000 5.55600 5.71400 '+$ '5.74700 5.78000 5.81400 5.84800 5.88200 6.06100 6.13500 '+$ '6.25000 6.28900 6.32900 6.36900 6.41000 6.45200 6.49400 '+$ '6.57900 6.66700 6.75700 6.89700 7.04200 7.14300 7.24600 '+$ '7.35300 7.46300 7.56700 7.69200 7.81200 7.93700 8.06500 '+$ '8.19700 8.33300 8.47500 8.69600 8.92900 9.09100 9.25900 '+$ '9.52400 9.80400 10.00000 10.20000 10.31000 10.42000 10.53000 '+$ '10.6400 10.7500 10.8700 11.0000 11.1100 11.3600 11.6300 '+$ '11.9000 12.2000 12.5000 12.8200 12.9900 13.1600 13.3300 '+$ '13.5100 13.7000 13.8900 14.0800' reads,waves,wave dxrs = $ '0.83440 0.840400 0.861600 0.893000 1.02330 1.20040 1.41640 '+$ '1.40220 1.39630 1.34140 1.32030 1.52130 1.61170 1.63830 '+$ '1.55640 1.45990 1.42920 1.40420 1.32490 1.31510 1.31450 '+$ '1.31400 1.31350 1.31300 1.31260 1.31220 1.31180 1.31140 '+$ '1.31100 1.31060 1.31030 1.31000 1.30970 1.30940 1.30910 '+$ '1.30880 1.30850 1.30830 1.30800 1.30780 1.30760 1.30730 '+$ '1.30710 1.30690 1.30670 1.30650 1.30620 1.30600 1.30580 '+$ '1.30570 1.30550 1.30530 1.30510 1.30490 1.30470 1.30450 '+$ '1.30440 1.30420 1.30400 1.30380 1.30370 1.30350 1.30330 '+$ '1.30320 1.30300 1.30280 1.30270 1.30250 1.30230 1.30220 '+$ '1.30200 1.30180 1.30170 1.30150 1.30130 1.30120 1.30100 '+$ '1.30080 1.30060 1.30050 1.30030 1.30010 1.30000 1.29980 '+$ '1.29960 1.29950 1.29930 1.29910 1.29890 1.29870 1.29850 '+$ '1.29840 1.29820 1.29800 1.29780 1.29760 1.24510 1.24270 '+$ '1.24000 1.23730 1.23270 1.22580 1.22200 1.21550 1.21180 '+$ '1.20790 1.20210 1.19570 1.17410 1.14730 1.10770 1.06740 '+$ '1.04760 1.02650 1.00360 0.982000 0.965000 0.959600 0.972700 '+$ '0.99170 1.00670 1.02190 1.04270 1.07600 1.12950 1.21270 '+$ '1.32510 1.42600 1.49660 1.55100 1.59990 1.63630 1.65020 '+$ '1.64280 1.62690 1.61280 1.59240 1.57330 1.55770 1.54130 '+$ '1.52650 1.51140 1.49730 1.48450 1.47210 1.46120 1.45130 '+$ '1.44210 1.43370 1.41550 1.39420 1.38730 1.37730 1.36450 '+$ '1.35410 1.34460 1.33880 1.33810 1.33850 1.34050 1.34290 '+$ '1.34420 1.34630 1.34420 1.33450 1.32900 1.32330 1.31800 '+$ '1.31430 1.30620 1.29720 1.28900 1.28730 1.28600 1.28510 '+$ '1.28540 1.28810 1.30160 1.30900 1.31720 1.31890 1.32040 '+$ '1.32200 1.32240 1.32150 1.32040 1.31810 1.31710 1.31810 '+$ '1.31950 1.31930 1.31900 1.31910 1.31800 1.31630 1.31540 '+$ '1.31540 1.31550 1.31450 1.31190 1.30680 1.29930 1.29250 '+$ '1.28390 1.27400 1.26720 1.25990 1.24510 1.22240 1.19910 '+$ '1.17150 1.15530 1.13700 1.11810 1.10130 1.09080 1.08730 '+$ '1.09250 1.10650 1.14780 1.20200 1.25820 1.32310 1.38570 '+$ '1.44480 1.47170 1.49620 1.51650 1.53330 1.54900 1.56280 '+$ '1.57320 ' reads,dxrs,dxr dxis = $ '0.164000 0.223000 0.279000 0.392000 0.517000 0.524000 0.380000 '+$ '0.251000 0.247000 0.180000 0.344000 0.389000 0.298000 0.156000 '+$ '1.00000E-01 3.00000E-05 4.00000E-07 1.64000E-08 3.76500E-09 '+$ '1.53000E-09 1.55000E-09 1.64000E-09 1.78000E-09 1.91000E-09 '+$ '2.14000E-09 2.26000E-09 2.54000E-09 2.93000E-09 3.11000E-09 '+$ '3.29000E-09 3.52000E-09 4.04000E-09 4.88000E-09 5.73000E-09 '+$ '6.89000E-09 8.58000E-09 1.04000E-08 1.22000E-08 1.43000E-08 '+$ '1.66000E-08 1.89000E-08 2.09000E-08 2.40000E-08 2.90000E-08 '+$ '3.44000E-08 4.03000E-08 4.30000E-08 4.92000E-08 5.87000E-08 '+$ '7.08000E-08 8.58000E-08 1.02000E-07 1.18000E-07 1.34000E-07 '+$ '1.40000E-07 1.43000E-07 1.45000E-07 1.51000E-07 1.83000E-07 '+$ '2.15000E-07 2.65000E-07 3.35000E-07 3.92000E-07 4.20000E-07 '+$ '4.44000E-07 4.74000E-07 5.11000E-07 5.53000E-07 6.02000E-07 '+$ '7.55000E-07 9.26000E-07 1.12000E-06 1.33000E-06 1.62000E-06 '+$ '2.00000E-06 2.25000E-06 2.33000E-06 2.33000E-06 2.17000E-06 '+$ '1.96000E-06 1.81000E-06 1.74000E-06 1.73000E-06 1.70000E-06 '+$ '1.76000E-06 1.82000E-06 2.04000E-06 2.25000E-06 2.29000E-06 '+$ '3.04000E-06 3.84000E-06 4.77000E-06 5.76000E-06 6.71000E-06 '+$ '8.66000E-06 1.02000E-05 6.18000E-04 6.80000E-04 7.30000E-04 '+$ '7.82000E-04 7.82000E-04 9.25000E-04 9.20000E-04 8.92000E-04 '+$ '8.70000E-04 8.90000E-04 1.01000E-03 1.35000E-03 3.42000E-03 '+$ '7.92000E-03 2.00000E-02 3.80000E-02 5.20000E-02 6.80000E-02 '+$ '9.23000E-02 0.127000 0.169000 0.221000 0.276000 0.312000 0.347000 '+$ '0.388000 0.438000 0.493000 0.554000 0.612000 0.625000 0.593000 '+$ '0.539000 0.491000 0.438000 0.372000 0.300000 0.238000 0.193000 '+$ '0.158000 0.121000 1.03000E-01 8.36000E-02 6.68000E-02 5.40000E-02 '+$ '4.22000E-02 3.42000E-02 2.74000E-02 2.20000E-02 1.86000E-02 '+$ '1.52000E-02 1.26000E-02 1.06000E-02 8.02000E-03 6.85000E-03 '+$ '6.60000E-03 6.96000E-03 9.16000E-03 1.11000E-02 1.45000E-02 '+$ '2.00000E-02 2.30000E-02 2.60000E-02 2.90000E-02 2.93000E-02 '+$ '3.00000E-02 2.85000E-02 1.73000E-02 1.29000E-02 1.20000E-02 '+$ '1.25000E-02 1.34000E-02 1.40000E-02 1.75000E-02 2.40000E-02 '+$ '3.50000E-02 3.80000E-02 4.20000E-02 4.60000E-02 5.20000E-02 '+$ '5.70000E-02 6.90000E-02 7.00000E-02 6.70000E-02 6.50000E-02 '+$ '6.40000E-02 6.20000E-02 5.90000E-02 5.70000E-02 5.60000E-02 '+$ '5.50000E-02 5.70000E-02 5.80000E-02 5.70000E-02 5.50000E-02 '+$ '5.50000E-02 5.40000E-02 5.20000E-02 5.20000E-02 5.20000E-02 '+$ '5.20000E-02 5.00000E-02 4.70000E-02 4.30000E-02 3.90000E-02 '+$ '3.70000E-02 3.90000E-02 4.00000E-02 4.20000E-02 4.40000E-02 '+$ '4.50000E-02 4.60000E-02 4.70000E-02 5.10000E-02 6.50000E-02 '+$ '7.50000E-02 8.80000E-02 0.108000 0.134000 0.168000 0.204000 '+$ '0.248000 0.280000 0.341000 0.379000 0.409000 0.422000 0.422000 '+$ '0.403000 0.389000 0.374000 0.354000 0.335000 0.315000 0.294000 '+$ '0.271000 ' reads,dxis,dxi ;---- interpolate index to wave_in if (wave_in ge wave(0) and wave_in le wave(numw-1) ) then begin for i = 0,(numw-2) do begin if (wave(i) le wave_in and wave(i+1) gt wave_in) then begin term = (wave_in-wave(i)) /(wave(i+1)-wave(i)) dr = term * (dxr(i+1)-dxr(i)) + dxr(i) di = term * (dxi(i+1)-dxi(i)) + dxi(i) result = complex(dr,di) goto,jump1 endif endfor print ,wave_in,' not in range in index_ice_266K.pro' result = complex(0.0,0.0) endif jump1: return,result end