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

Go to the source code of this file.

Namespaces

namespace  vectorpotentialdipole_compare_with_data

Variables

 vectorpotentialdipole_compare_with_data.line_phi = float(sys.argv[1])
str vectorpotentialdipole_compare_with_data.outfilename = "./vecpotdip_compare_"+str(int(line_phi))+".png"
str vectorpotentialdipole_compare_with_data.inputLocation = "/turso/group/spacephysics/vlasiator/data/L0/2D/BCH/bulk/"
list vectorpotentialdipole_compare_with_data.times = [0,10,50,100,200,500]
list vectorpotentialdipole_compare_with_data.colors = ['r','g','b','magenta','k']
list vectorpotentialdipole_compare_with_data.timefulls = [str(time).rjust(7, '0') for time in times]
list vectorpotentialdipole_compare_with_data.file_names = [inputLocation+"bulk."+timefull+".vlsv" for timefull in timefulls]
list vectorpotentialdipole_compare_with_data.vlsvobj = []
int vectorpotentialdipole_compare_with_data.RE = 6371000.
int vectorpotentialdipole_compare_with_data.tilt_angle_phi = 0.
int vectorpotentialdipole_compare_with_data.tilt_angle_theta = 0.
int vectorpotentialdipole_compare_with_data.line_theta = 0. * math.pi/180.
 vectorpotentialdipole_compare_with_data.line_start = np.array([0,0,0])
float vectorpotentialdipole_compare_with_data.step = 0.1
int vectorpotentialdipole_compare_with_data.linewidth = 2
int vectorpotentialdipole_compare_with_data.linthresh = 1.e-10
int vectorpotentialdipole_compare_with_data.fontsize = 20
 vectorpotentialdipole_compare_with_data.dip = fieldmodels.dipole(0,0,0,tilt_angle_phi,tilt_angle_theta)
 vectorpotentialdipole_compare_with_data.mdip = fieldmodels.dipole(80*RE,0,0,tilt_angle_phi,180.-tilt_angle_theta)
 vectorpotentialdipole_compare_with_data.fig = plt.figure()
int vectorpotentialdipole_compare_with_data.nsubplots = 3
 vectorpotentialdipole_compare_with_data.axes = fig.get_axes()
int vectorpotentialdipole_compare_with_data.radii = np.arange(0.1,45,step)*RE
 vectorpotentialdipole_compare_with_data.nr = len(radii)
int vectorpotentialdipole_compare_with_data.radiiRE = radii/RE
int vectorpotentialdipole_compare_with_data.xv = line_start[0]*RE + radii*np.sin(line_phi)*np.cos(line_theta)
int vectorpotentialdipole_compare_with_data.yv = line_start[1]*RE + radii*np.sin(line_phi)*np.sin(line_theta)
int vectorpotentialdipole_compare_with_data.zv = line_start[2]*RE + radii*np.cos(line_phi)
 vectorpotentialdipole_compare_with_data.B1 = np.zeros([nr,3])
 vectorpotentialdipole_compare_with_data.B2 = np.zeros([nr,3])
 vectorpotentialdipole_compare_with_data.B3 = np.zeros([nr,3])
 vectorpotentialdipole_compare_with_data.B4 = np.zeros([nr,3])
list vectorpotentialdipole_compare_with_data.coords = ['x','y','z']
 vectorpotentialdipole_compare_with_data.ax = axes[k]
 vectorpotentialdipole_compare_with_data.c
 vectorpotentialdipole_compare_with_data.linestyle
 vectorpotentialdipole_compare_with_data.label
 vectorpotentialdipole_compare_with_data.zorder
list vectorpotentialdipole_compare_with_data.vf = vlsvobj[i]
 vectorpotentialdipole_compare_with_data.res = pt.calculations.cut_through_step(vf, [xv[0],0,zv[0]], [xv[-1],0,zv[-1]])
 vectorpotentialdipole_compare_with_data.cut = res[0].data
 vectorpotentialdipole_compare_with_data.pr_dist = res[1].data
 vectorpotentialdipole_compare_with_data.pr_coords = res[2].data
int vectorpotentialdipole_compare_with_data.pr_Re = np.array(pr_dist)/RE
list vectorpotentialdipole_compare_with_data.pr_B = vf.read_variable("B", operator=coords[k],cellids=cut)
 vectorpotentialdipole_compare_with_data.linthreshy
 vectorpotentialdipole_compare_with_data.ylims = np.array(ax.get_ylim())
 vectorpotentialdipole_compare_with_data.handles
 vectorpotentialdipole_compare_with_data.labels