file='./fort.16'
readcol, file, skip=1, w, nsi1, ksi, nc1, kc, neis1, keis

; descending order in wavelength:

; ------------
file='./Const/IceWarren2008.tab'
readcol, file, skip=3, wwaro, nwaro,  kwaro

wwar = reverse(wwaro)
nwar = reverse(nwaro)
kwar = reverse(kwaro)

nwar1 = nwar -1.


; ------------
file='./Const/IceHudgins1993.tab'
readcol, file, skip=3, T, wh1cm, nh,  kh
nh1 = nh - 1.
wh  = 1./wh1cm * 1e4
ii10  = where(T eq 10.)


; ------------
file ='nk_Ice.tab'
readcol, file, skip=3, wice, nice, kice




plot_oi, wwar, nwar1, xr=[0.09, 2000.], /xsty,ytitle = 'n(Ice) -1', /ysty, $
        xthick=4, ythick=4,  xtitle = '!9l !3(!9m!3m)',  /nodata
oplot, wh(ii10), nh1(ii10), thick=1 
oplot, wwar, nwar1, color=fsc_color('green'), linesty=0

oplot, wice, nice-1., color=fsc_color('orange')
oplot, w, neis1
oplot, w, neis1, psym=4

plot_oi, wwar, kwar, xr=[5.5e-2, 700.], /xsty,ytitle = 'k(Ice)',  $
        xthick=4, ythick=4,  xtitle = '!9l !3(mu!3m)'   ;,  /nodata,
	
oplot, wh(ii10), kh(ii10), thick=2
oplot, wwar, kwar, color=fsc_color('green'), linesty=0
oplot, wice, kice, color=fsc_color('brown')
oplot, weis, keis
oplot, weis, keis, psym=4
end