;  $Id: pbb.pro,v 1.1 2004/06/20 14:09:35 brandenb Exp $
;
;  puts B-field (without ghost zones) into bb array
;  the name pbb is chosen, because cbb does already exist
;
bb=pc_noghost(curl(aa))
;
;  minmax diagnostics
;
for i=0,2 do print,i,minmax(bb(*,*,*,i))
;
;  for convenience, also calcuate coordinates without ghost zones
;
xxx=x(l1:l2)
yyy=y(m1:m2)
zzz=z(n1:n2)
;
;  plot example
;
plot,zzz,bb(0,0,*,1),yr=[-1,1]*.2
oplot,xxx,bb(*,0,0,1),col=111
tvscl,bb(*,0,*,1)
;wait,.1
;contour,reform(bb(*,10,*,1)),/fill,nlev=60
;
END