;;; restore,'/wkspace2/brandenb/f77/ndynd/pencil_test/idlsave.dat'
;;; print,'time=',time
xxx=x(l1:l2)
yyy=y(m1:m2)
ccc=cos(yyy)
sss=sin(yyy)
aaa=pc_noghost(aa)
m=51 ;;; & plot,a(*,m),yst=0,yr=[0,.005] & oplot,aaa(*,m,2),col=122,li=2
;
!p.multi=0
bbb=pc_noghost(bb) & fac=1.
;;; m=41 & plot,b(*,m),yst=0,yr=[0,.00015]
oplot,fac*bbb(*,m,2),col=122
;
;  try and reproduce Bphi by -(1/r)*(d/dtheta)*Ar
;
dt=8e-6 & alp=1e3 & fac=1.
Bphi_theo=ampl*(1.-xxx)*(xxx-.2)*sss(m)+fac*dt*alp*ampl*sss(m)*2./xxx^3
oplot,Bphi_theo,col=55,li=2
print,Bphi_theo
;
;  result: bbb is correct, and b is not!
;
;;; stop
;;; m=51 & plot,aaa(*,m,2)-a(*,m)
END