pro show_kernels
  plot,[1],[1],xr=[9,500],yr=[-0.005,0.07],/xlog,xtit='N',ytit='log(rho_sim/rho_ini)',xst=1,yst=1,chars=2,charth=2,xth=2,yth=2
  oplot,[1,1000],[0,0],col=96

  name=['WendlandC2','Cubic Spline','WendlandC6']
  kernel=['wc2','cs','wc6']
  sym=[1,2,4,5]
  kernel_min=[13,12,28]
  kernel_max=[128,128,400]
  FOR k=0,N_ELEMENTS(kernel)-1 DO BEGIN
     oplot,[200],[0.065-k*0.003],psym=sym[k],symsize=2
     xyouts,220,0.065-k*0.003,name[k],charsize=2,charthick=2
  END

  FOR ingb=13,400 DO BEGIN
     FOR k=0,N_ELEMENTS(kernel)-1 DO BEGIN
        IF(ingb GE kernel_min[k] AND ingb LE kernel_max[k]) THEN BEGIN
           snr=kernel[k]+'_'+strtrim(string(ingb,form='(i)'),1)+'_000'
           readnew,snr,h,'HEAD'
           readnew,snr,mass,'MASS'
           rho_ini = TOTAL(mass,/double) / h.BOXSIZE^3 
           readnew,snr,rho,'RHO'
           rho_sim = TOTAL(rho) / N_ELEMENTS(rho)
           oplot,[ingb],[alog10(rho_sim/rho_ini)],psym=sym[k]
        END
     END
  END
  
end