;
; Script for visualizating magnetic field (in the x-y plane)
; vectors show x and y field, while color coding correspond to z field
;
if !d.name eq 'PS' then begin
  device,xsize=14,ysize=14,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
if !d.name eq 'X' then begin
  window,0,xsize=640*9/10,ysize=480*9/10
end
;
siz=1.4
!p.charsize=siz
!p.multi=[0,1,1]
loadct,3
;
; mv idl.ps ~/rmhd/other/anvar/cluster/visualize_Pm=30.ps
;
t0=2517.55
urms0=0.123
Brms0=0.063
to=1/urms0
;
; Restore data
;
file='store_slice_t= 806.29_zpt=99.sav'
restore,file
;
; Set some auxillary variables
;
nx=256
xx=(indgen(nx)-(nx-1)/2.)*!pi/((nx-1)/2.)
;
; Reduce number of mesh points for clarity
;
red_fac=8
new_nx=nx/red_fac
bb=fltarr(new_nx,new_nx,3)
x=fltarr(new_nx)
ii=0
for i=0,nx-1,red_fac do begin
    jj=0
    for j=0,nx-1,red_fac do begin
        bb(ii,jj,*)=b_save(i,j,*)
        jj=jj+1
    end
    x(ii)=xx(i)
    ii=ii+1
end
;
; Plot reults
;
;contour,bb(*,*,2),x,x,nlevels=60,/fill
contour,b_save(*,*,2),xx,xx,nlevels=60,/fill,$
  xtitle='!8x!6',ytitle='!8y!6',ticklen=-0.02
velovect,bb(*,*,0),bb(*,*,1),x,x,/overplot,length=3
;
loadct,5
;
END