;$Id: pvectors.pro,v 1.1 2006/10/07 12:25:53 brandenb Exp $
;
;  plots magnetic field in the y-z plane, for an extended z-range
;
pc_potentialfield_exterior_z,aa,aaa=aaa,zzz=zzz
bbb=curl(aaa)
;
l=11
contour,reform(bbb(l,m1:m2,nn1:nn2,0)),nlev=22,y(m1:m2),zzz(nn1:nn2),/fil
vel_a,reform(bbb(l,m1:m2,nn1:nn2,1)),reform(bbb(l,m1:m2,nn1:nn2,2)),$
  nvec=800,y(m1:m2),zzz(nn1:nn2),/over
oplot,x,x-x+.5,col=255
oplot,x,x-x-.5,col=255
;
;  interfaces
;
n=n_elements(zzz)
nnn1=(n-nz)/2
nnn2=(n+nz)/2-1
;
tvscl,bbb(l1:l2,m1:m2,nnn1),0
tvscl,bbb(l1:l2,m1:m2,nnn2),1
;
END