Vlasiator ebf0dd394 on dev (v5.4.0 + 1054 commits)
Loading...
Searching...
No Matches
vectorpotentialdipole_verify.py File Reference

Go to the source code of this file.

Namespaces

namespace  vectorpotentialdipole_verify

Variables

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