;set the ever-useful angstrom plotting symbol in IDL character specification terms
ang = '!3!sA!r!u!9 %!3!n'

;set continuum level
cont=9500.

;*CAK*****
READCOL,'spectre_3928_10.dat',cax,cay,numline=5000
cay=cay/cont


;***********HA************
READCOL,'spectre_6558_10.dat',hax,hay,numline=5000
hay=hay/cont


!P.thick=2
!p.charsize=2
!p.charthick=3

window,0,xsize=800,ysize=400

toggle,file='profiles_ha_ca.ps'
plot,cax,cay,pos=[0.2,0.7,0.45,0.9],/norm,$
    xrange=[3933.7-1.5,3933.7+1.5],xst=10,yst=11,$
    ytit='Solar Intensity',xtit='Wavelength ('+ang+')',$
    xtickv=[3933,3934,3935],xticks=2,$
    tit='(a) Ca !4II!3 K'
plots,[3933.45,3932.7],[720./cont,600./cont]
plots,[3933.85,3933.99],[620./cont,540./cont]

xyouts,3932.6,540./cont,/data,'K!D2V!N'
xyouts,3934.1,515./cont,/data,'K!D2R!N'
xyouts,3933.2,446./cont,/data,'K!D3!N'
xyouts,3932.6,1230./cont,/data,'K!D1!N'
xyouts,3934.4,1230./cont,/data,'K!D1!N'

plot,hax,hay,pos=[0.60,0.7,0.90,0.9],/norm,/noerase,$
    xrange=[6562.8-2.5,6562.8+2.5],xst=10,yst=11,$
    ytit='Solar Intensity',xtit='Wavelength ('+ang+')',$
    xtickv=[6561,6563,6565],xticks=2,$
    tit='(b) H!4a!3'
toggle
$gv profiles_ha_ca.ps



















print, caxmin & READ,'enter element number: ',caxminele
0
caxmin=caxmin(caxminele)

carange=900

size_gauss=(2*carange)-1
height=1

cags=gauss_put(size_gauss,0,caheight,carange-1,64)

cafwhm=where (cags GT 0.49*caheight AND cags LT 0.51*caheight)

print,cafwhm
cafwhm_x1=0
cafwhm_x2=0

READ,'input lower value: ',cafwhm_x1 & READ,'input higher value: ',cafwhm_x2
823
974


plot,cax(caxmin-carange:caxmin+carange-1),cay(caxmin-carange:caxmin+carange-1),$
  ystyle=8,ytit='Solar Intensity',xtit='Wavelength',pos=[0.2,0.2,0.9,0.95],/norm,$
  charsize=1.4,charthick=3.5,thick=3.5,$
  xtickv=[3932,3933,3934,3935],xticks=3,$
  ytickv=[400,600,800,1000,1200,1400],yticks=5,yrange=[0,1400]
AXIS,yrange=[0,3.5],/yaxis,/save,ytit='Normalised Filter Transmission',$
  charsize=1.4,charthick=2.0,ytickv=[0,0.5,1],yticks=2
oplot,cax(caxmin-carange:caxmin+carange-1),cags,thick=2.0

;horline,height/2.
verline,cax(cafwhm_x1+caxmin-carange),line=2
verline,cax(cafwhm_x2+caxmin-carange),line=2
print,cax(cafwhm_x2+caxmin-carange)-cax(cafwhm_x1+caxmin-carange)

;***********HA************
READCOL,'ha2.txt',hax,hay,nlines=2500

;plot,x,y

haxmin=where(hay EQ min(hay))

print, haxmin & READ,'enter element number: ',haxminele 
6
haxmin=haxmin(haxminele)
harange=913

size_gauss=(2*harange)-1
height=1

hags=gauss_put(size_gauss,0,height,harange-1,45)

hafwhm=where (hags GT 0.49*height AND hags LT 0.51*height)

print,hafwhm
hafwhm_x1=0
hafwhm_x2=0

READ,'input lower value: ',hafwhm_x1 & READ,'input higher value: ',hafwhm_x2
859
965


plot,hax(haxmin-harange:haxmin+harange-1),hay(haxmin-harange:haxmin+harange-1),$
  ystyle=8,ytit='Solar Intensity',xtit='Wavelength',pos=[0.2,0.2,0.9,0.95],/norm,$
  charsize=2,charthick=2,thick=2,xrange=[6561,6564.5],$
  xticks=3,xtickv=[6561,6562,6563,6564],$
  yrange=[-2000,8000],ytickv=[2000,4000,6000,8000],yticks=3
AXIS,yrange=[0,3.5],/yaxis,/save,ytit='Normalised Filter Transmission',$
  charsize=2,charthick=2,ytickv=[0,0.5,1],yticks=2
oplot,hax(haxmin-harange:haxmin+harange-1),hags,thick=2

;horline,height/2.
verline,hax(hafwhm_x1+haxmin-harange),line=2
verline,hax(hafwhm_x2+haxmin-harange),line=2
print,hax(hafwhm_x2+haxmin-harange)-hax(hafwhm_x1+haxmin-harange)

toggle
plot,cax(caxmin-carange:caxmin+carange-1),cay(caxmin-carange:caxmin+carange-1),$
  ystyle=8,ytit='Solar Intensity',pos=[0.06,0.60,0.42,0.90],/norm,$
  charsize=1.4,charthick=3.5,thick=3.5,$
  xtickv=[3932,3933,3934,3935],xticks=3,$
  ytickv=[400,600,800,1000,1200,1400],yticks=5,yrange=[0,1400],tit='(a) Ca II K3'
AXIS,yrange=[0,3.5],/yaxis,/save,$
  charsize=1.4,charthick=2.0,ytickv=[0,0.5,1],yticks=2
oplot,cax(caxmin-carange:caxmin+carange-1),cags,thick=2.0

plot,hax(haxmin-harange:haxmin+harange-1),hay(haxmin-harange:haxmin+harange-1),$
  ystyle=8,pos=[0.54,0.60,0.90,0.90],/norm,$
  charsize=1.4,charthick=3.5,thick=3.5,$
  xticks=3,xtickv=[6561,6562,6563,6564],$
  yrange=[-2000,8000],ytickv=[2000,4000,6000,8000],yticks=3,/noerase,tit='(b) H-alpha'
AXIS,yrange=[0,3.5],/yaxis,/save,ytit='Normalised Filter Transmission',$
  charsize=1.4,charthick=2.0,ytickv=[0,0.5,1],yticks=2
oplot,hax(haxmin-harange:haxmin+harange-1),hags,thick=2.0