readcol,"IRCALFILTER-J.txt",lambda_j,t_j
readcol,"IRCALFILTER-H.txt",lambda_h,t_h
readcol,"IRCALFILTER-K.txt",lambda_k,t_k
readcol,"IRCALFILTER-KsBarr.txt",lambda_ks,t_ks

readcol,"IRCALFILTER-BrG.txt",lambda_brg,t_brg
readcol,"IRCALFILTER-H2.txt",lambda_h2,t_h2
readcol,"IRCALFILTER-HCONT.txt",lambda_hcont,t_hcont

readcol,"IRCALFILTER-JCH4.txt",lambda_jch4,t_jch4
readcol,"IRCALFILTER-KCH4.txt",lambda_kch4,t_kch4

;!p.multi=[0,1,2]
plot,lambda_j,t_j,/nodata,xr=[1000,2500],yr=[0,1],xtitle="Wavelength (nm)",ytitle="Transmission Fraction"
labeloplot,lambda_j,t_j,"J"
labeloplot,lambda_h,t_h,"H"
labeloplot,lambda_k,t_k,"K",xoffset=250,yoffset=-0.05
labeloplot,lambda_ks,t_ks,"Ks",color=getcolor('gray')


labeloplot,lambda_brg,t_brg,"BrG",color=getcolor('red')
labeloplot,lambda_h2,t_h2,textoidl("H_2"),color=getcolor('orange')
labeloplot,lambda_hcont,t_hcont,textoidl("H_{cont}"),color=getcolor('cyan')
labeloplot,lambda_kch4,t_kch4,textoidl("CH_4 2.1"),color=getcolor('yellow')


  
  lambda = findgen(1500)+1000
  
  tr_j = interpol(t_j,lambda_j,lambda) >0
  tr_jch4 = interpol(t_jch4,lambda_jch4,lambda) >0
  labeloplot,lambda,tr_jch4*tr_j,color=getcolor('cyan'),textoidl("CH_4-1.2 \times J")



; OTHER FILTERS

names = ['[FeII]', textoidl('H_{cont}'), '2.2/0.04',textoidl('K_{cont}')]
lambda_c = [1644,1570,2192,2270]
fwhms = [16,20,47,20]
colors=[getcolor('blue'),getcolor('blue'),getcolor('magenta'),getcolor('purple')]

sigmas = fltarr(4)
peak=0.8
for i = 0,3 do begin
	sigmas[i] = fwhms[i] / ( 2 * SQRT(2* ALOG(2)) )

	fakedata = gaussian(lambda,[peak,lambda_c[i],sigmas[i]])
	labeloplot,lambda,fakedata,color=colors(i),names[i],lines=2
endfor

end