restore, 'routines.sav'

 xDir = 'MDI/2001'  & xCase = 'mdi-2001'

;xDir = 'MDI/2002'  & xCase = 'mdi-2002'
;xDir = 'MDI/2010'  & xCase = 'mdi-2010'

;xDir = 'GONG/2002' & xCase = 'gong-2002'
;xDir = 'GONG/2010' & xCase = 'gong-2010'

;xDir = 'HMI/2010'  & xCase = 'hmi-2010'


v = read_struct(xDir+'/multiplets-'+xCase+'.dat')
IF !d.name EQ 'X' THEN window, xsize = 1600, ysize = 1200
!p.multi = [0, 3, 2]
IF !d.name EQ 'X' THEN !p.charsize = 2 ELSE !p.charsize = 1.5

!p.title = xDir
plot, v.lval, v.freq, psym = 3, /ynoz, xtitle = 'l', ytitle = 'freq'
!x.title = 'freq'
plot, v.freq, v.sfreq, psym = 3, /ylog, ytitle = 'sigma-freq'
plot, v.freq, v.dfreq, psym = 3,        ytitle = 'delta-freq'

plot, v.freq, v.xfwhm  > 1,      psym = 3, /ylog, ytitle = 'xFWHM'
plot, v.freq, v.dalpha > (-.01), psym = 3,        ytitle = 'dAlpha'
plot, v.freq, v.ramp,            psym = 3, /ylog, ytitle = 'rAmp'

pause

plot, v.freq, v.fwhm_r > 1,   psym = 3, /ylog, ytitle = 'FWHM_ridge'
plot, v.freq, v.alpha_r,      psym = 3,        ytitle = 'alpha_ridge'
plot, v.freq, v.amp_r > 1e-4, psym = 3, /ylog, ytitle = 'amp_ridge'

plot, v.freq, abs(v.fwhm_m), psym = 3, /ylog, ytitle = 'FWHM_mode, FWHM_e'
oplot, v.freq, v.fwhm_e, psym = 3
plot, v.freq, v.alpha_m, psym = 3,        ytitle = 'alpha_mode'
plot, v.freq, v.amp_m > 1e-4,   psym = 3, /ylog, ytitle = 'amp_mode'

pause

plot, v.freq, v.sfwhm_r >  1e-2, psym = 3, /ylog, ytitle = 'sigma-FWHM_ridge'
plot, v.freq, v.salpha_r > 1e-3, psym = 3, /ylog, ytitle = 'sigma-alpha_ridge'
plot, v.freq, v.samp_r,   psym = 3, /ylog, ytitle = 'sigma-amp_ridge'

plot, v.freq, v.sfwhm_m  > 1e-2, psym = 3, /ylog, ytitle = 'sigma-FWHM_mode'
plot, v.freq, v.salpha_m > 1e-3, psym = 3, /ylog, ytitle = 'sigma-alpha_mode'
plot, v.freq, v.samp_m   > 1e-5, psym = 3, /ylog, ytitle = 'sigma-amp_mode'

!p.title = ''
!x.title = ''
!p.multi = 0
!p.charsize = 0
END