vectorpotentialdipole_verify2.py File Reference

Go to the source code of this file.

Namespaces
namespace	vectorpotentialdipole_verify2

Variables
	vectorpotentialdipole_verify2.testset = int(sys.argv[1])
bool	vectorpotentialdipole_verify2.calcderivatives = True
bool	vectorpotentialdipole_verify2.plotmagnitude = False
str	vectorpotentialdipole_verify2.outfilename = "./vecpotdip_verify2_"+str(testset)+".png"
int	vectorpotentialdipole_verify2.RE = 6371000.
int	vectorpotentialdipole_verify2.epsilon = 1.e-15
list	vectorpotentialdipole_verify2.BGB = [0.,0.,0.]
int	vectorpotentialdipole_verify2.tilt_angle_phi = 0.
int	vectorpotentialdipole_verify2.tilt_angle_theta = 0.
int	vectorpotentialdipole_verify2.line_theta = 0.
	vectorpotentialdipole_verify2.line_start = np.array([0,0,0])
int	vectorpotentialdipole_verify2.line_phi = np.array([0,45,80,90,110,135])*math.pi/180.
float	vectorpotentialdipole_verify2.step = 0.1
int	vectorpotentialdipole_verify2.linewidth = 2
int	vectorpotentialdipole_verify2.linthresh = 1.e-10
int	vectorpotentialdipole_verify2.fontsize = 20
	vectorpotentialdipole_verify2.dip = fieldmodels.dipole(0,0,0,tilt_angle_phi,tilt_angle_theta)
	vectorpotentialdipole_verify2.mdip = fieldmodels.dipole(80*RE,0,0,tilt_angle_phi,180.-tilt_angle_theta)
	vectorpotentialdipole_verify2.imfpot = fieldmodels.IMFpotential(radius_z=10, radius_f=40, IMF=BGB)
	vectorpotentialdipole_verify2.fig = plt.figure()
	vectorpotentialdipole_verify2.nsubplots = len(line_theta)
	vectorpotentialdipole_verify2.axes = fig.get_axes()
int	vectorpotentialdipole_verify2.radii = np.arange(0.1,100,step)*RE
	vectorpotentialdipole_verify2.nr = len(radii)
int	vectorpotentialdipole_verify2.radiiRE = radii/RE
	vectorpotentialdipole_verify2.ax = axes[i]
	vectorpotentialdipole_verify2.transform
	vectorpotentialdipole_verify2.transAxes
	vectorpotentialdipole_verify2.bbox
int	vectorpotentialdipole_verify2.xv = line_start[0]*RE + radii*np.sin(line_phi[i])*np.cos(line_theta[i])
int	vectorpotentialdipole_verify2.yv = line_start[1]*RE + radii*np.sin(line_phi[i])*np.sin(line_theta[i])
int	vectorpotentialdipole_verify2.zv = line_start[2]*RE + radii*np.cos(line_phi[i])
	vectorpotentialdipole_verify2.B1 = np.zeros([nr,4])
	vectorpotentialdipole_verify2.B2 = np.zeros([nr,4])
	vectorpotentialdipole_verify2.B3 = np.zeros([nr,4])
	vectorpotentialdipole_verify2.B4 = np.zeros([nr,4])
list	vectorpotentialdipole_verify2.colors = ['r','k','b','magenta']
list	vectorpotentialdipole_verify2.coords = ['x','y','z','mag']
	vectorpotentialdipole_verify2.plotrange = range(3)
	vectorpotentialdipole_verify2.c
	vectorpotentialdipole_verify2.linestyle
	vectorpotentialdipole_verify2.label
	vectorpotentialdipole_verify2.zorder
	vectorpotentialdipole_verify2.linthreshy
	vectorpotentialdipole_verify2.ylims = np.array(ax.get_ylim())
	vectorpotentialdipole_verify2.handles
	vectorpotentialdipole_verify2.labels
float	vectorpotentialdipole_verify2.step2 = 0.00001
	vectorpotentialdipole_verify2.dB1 = np.zeros([nr,3,3])
	vectorpotentialdipole_verify2.dB2 = np.zeros([nr,3,3])
	vectorpotentialdipole_verify2.dB3 = np.zeros([nr,3,3])
	vectorpotentialdipole_verify2.dB4 = np.zeros([nr,3,3])
	vectorpotentialdipole_verify2.kk = kkk
tuple	vectorpotentialdipole_verify2.cdbx = (dip.getX(xv[j]+step2*RE,yv[j],zv[j],0,k,0) - dip.getX(xv[j]-step2*RE,yv[j],zv[j],0,k,0) + imfpot.get(xv[j]+step2*RE,yv[j],zv[j],0,k,0) - imfpot.get(xv[j]-step2*RE,yv[j],zv[j],0,k,0))/(2*step2*RE)
tuple	vectorpotentialdipole_verify2.cdby = (dip.getX(xv[j],yv[j]+step2*RE,zv[j],0,k,0) - dip.getX(xv[j],yv[j]-step2*RE,zv[j],0,k,0))/(2*step2*RE)
tuple	vectorpotentialdipole_verify2.cdbz = (dip.getX(xv[j],yv[j],zv[j]+step2*RE,0,k,0) - dip.getX(xv[j],yv[j],zv[j]-step2*RE,0,k,0))/(2*step2*RE)