from numpy import *
from zPage import *
from zPlot import *
erre=0.25*(69.7+73)
erre2=erre*erre
elle=145.0
N=90.0/(elle*1.0e-3)
Sref=1743.0e-6
mu0=4.0e-7*pi
zbpol=45
elle1=zbpol+0.5*elle
elle2=0.5*elle-zbpol
cs1=elle1/sqrt(erre2+elle1**2)
cs2=elle2/sqrt(erre2+elle2**2)


freq=array([1001.,4967.,10088.,20068.,49917.,102493.])
cur= array([994.7,278.1,139.9,71.0,28.0,11.2])*1.0e-3
Vbpol=array([38.45,53.13,54.20,54.65,54.40,52.19])*1.0e-3
Vref=array([6.98,9.30,9.50,9.50,9.49,8.86])*1.0e-3

B=mu0*N*cur*(cs1+cs2)/2.0
Bp=2.0*pi*freq*B
Bpref=Vref/Sref
Z=zPage()
P=zPlot(Z)
P.plot2d(freq,[Bp,Bpref])