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+'/cg-coefs-'+xCase+'.dat')
x = v.freq/(v.lval+0.5)
IF !d.name EQ 'X' THEN window, xsize = 1600, ysize = 1200
IF !d.name EQ 'X' THEN !p.charsize = 2 ELSE !p.charsize = 1.5
!p.title = xDir
!p.multi = [0, 3, 2]
!x.title = 'log(freq/L)'
plot, x, v.cg[0], psym = 3, yr = [350, 370], /ystyle, ytitle = 'CG[1]'
plot, x, v.cg[1], psym = 3, yr = [-1, 1]*.1, /ystyle, ytitle = 'CG[2]'
plot, x, v.cg[2], psym = 3, yr = [-10, -4],  /ystyle, ytitle = 'CG[3]'
plot, x, v.cg[3], psym = 3, yr = [-1, 1]*.1, /ystyle, ytitle = 'CG[4]'
plot, x, v.cg[4], psym = 3, yr = [-3, 2],    /ystyle, ytitle = 'CG[6]'
plot, x, v.cg[5], psym = 3, yr = [-1, 1]*.1, /ystyle, ytitle = 'CG[5]'
pause
FOR k = 0, 5 DO plot, x, v.scg[k] > 1e-4, psym = 3, /ylog, ytitle = 'sigma-CG['+string(k+1, format = '(i0)')+']'

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