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

Variables

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

Variable Documentation

◆ ax

vectorpotentialdipole_compare_with_data.ax = axes[k]

Definition at line 112 of file vectorpotentialdipole_compare_with_data.py.

◆ axes

vectorpotentialdipole_compare_with_data.axes = fig.get_axes()

Definition at line 83 of file vectorpotentialdipole_compare_with_data.py.

◆ B1

vectorpotentialdipole_compare_with_data.B1 = np.zeros([nr,3])

Definition at line 95 of file vectorpotentialdipole_compare_with_data.py.

◆ B2

vectorpotentialdipole_compare_with_data.B2 = np.zeros([nr,3])

Definition at line 96 of file vectorpotentialdipole_compare_with_data.py.

◆ B3

vectorpotentialdipole_compare_with_data.B3 = np.zeros([nr,3])

Definition at line 97 of file vectorpotentialdipole_compare_with_data.py.

◆ B4

vectorpotentialdipole_compare_with_data.B4 = np.zeros([nr,3])

Definition at line 98 of file vectorpotentialdipole_compare_with_data.py.

◆ c

vectorpotentialdipole_compare_with_data.c

Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.

◆ colors

list vectorpotentialdipole_compare_with_data.colors = ['r','g','b','magenta','k']

Definition at line 51 of file vectorpotentialdipole_compare_with_data.py.

◆ coords

list vectorpotentialdipole_compare_with_data.coords = ['x','y','z']

Definition at line 109 of file vectorpotentialdipole_compare_with_data.py.

◆ cut

vectorpotentialdipole_compare_with_data.cut = res[0].data

Definition at line 124 of file vectorpotentialdipole_compare_with_data.py.

◆ dip

vectorpotentialdipole_compare_with_data.dip = fieldmodels.dipole(0,0,0,tilt_angle_phi,tilt_angle_theta)

Definition at line 74 of file vectorpotentialdipole_compare_with_data.py.

◆ fig

vectorpotentialdipole_compare_with_data.fig = plt.figure()

Definition at line 78 of file vectorpotentialdipole_compare_with_data.py.

◆ file_names

list vectorpotentialdipole_compare_with_data.file_names = [inputLocation+"bulk."+timefull+".vlsv" for timefull in timefulls]

Definition at line 53 of file vectorpotentialdipole_compare_with_data.py.

◆ fontsize

vectorpotentialdipole_compare_with_data.fontsize = 20

Definition at line 70 of file vectorpotentialdipole_compare_with_data.py.

◆ handles

vectorpotentialdipole_compare_with_data.handles

Definition at line 148 of file vectorpotentialdipole_compare_with_data.py.

◆ inputLocation

str vectorpotentialdipole_compare_with_data.inputLocation = "/turso/group/spacephysics/vlasiator/data/L0/2D/BCH/bulk/"

Definition at line 49 of file vectorpotentialdipole_compare_with_data.py.

◆ label

vectorpotentialdipole_compare_with_data.label

Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.

◆ labels

vectorpotentialdipole_compare_with_data.labels

Definition at line 148 of file vectorpotentialdipole_compare_with_data.py.

◆ line_phi

int vectorpotentialdipole_compare_with_data.line_phi = float(sys.argv[1])

Definition at line 41 of file vectorpotentialdipole_compare_with_data.py.

◆ line_start

vectorpotentialdipole_compare_with_data.line_start = np.array([0,0,0])

Definition at line 65 of file vectorpotentialdipole_compare_with_data.py.

◆ line_theta

int vectorpotentialdipole_compare_with_data.line_theta = 0. * math.pi/180.

Definition at line 63 of file vectorpotentialdipole_compare_with_data.py.

◆ linestyle

vectorpotentialdipole_compare_with_data.linestyle

Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.

◆ linewidth

vectorpotentialdipole_compare_with_data.linewidth = 2

Definition at line 68 of file vectorpotentialdipole_compare_with_data.py.

◆ linthresh

int vectorpotentialdipole_compare_with_data.linthresh = 1.e-10

Definition at line 69 of file vectorpotentialdipole_compare_with_data.py.

◆ linthreshy

vectorpotentialdipole_compare_with_data.linthreshy

Definition at line 135 of file vectorpotentialdipole_compare_with_data.py.

◆ mdip

vectorpotentialdipole_compare_with_data.mdip = fieldmodels.dipole(80*RE,0,0,tilt_angle_phi,180.-tilt_angle_theta)

Definition at line 75 of file vectorpotentialdipole_compare_with_data.py.

◆ nr

vectorpotentialdipole_compare_with_data.nr = len(radii)

Definition at line 86 of file vectorpotentialdipole_compare_with_data.py.

◆ nsubplots

int vectorpotentialdipole_compare_with_data.nsubplots = 3

Definition at line 80 of file vectorpotentialdipole_compare_with_data.py.

◆ outfilename

str vectorpotentialdipole_compare_with_data.outfilename = "./vecpotdip_compare_"+str(int(line_phi))+".png"

Definition at line 47 of file vectorpotentialdipole_compare_with_data.py.

◆ pr_B

list vectorpotentialdipole_compare_with_data.pr_B = vf.read_variable("B", operator=coords[k],cellids=cut)

Definition at line 128 of file vectorpotentialdipole_compare_with_data.py.

◆ pr_coords

vectorpotentialdipole_compare_with_data.pr_coords = res[2].data

Definition at line 126 of file vectorpotentialdipole_compare_with_data.py.

◆ pr_dist

vectorpotentialdipole_compare_with_data.pr_dist = res[1].data

Definition at line 125 of file vectorpotentialdipole_compare_with_data.py.

◆ pr_Re

int vectorpotentialdipole_compare_with_data.pr_Re = np.array(pr_dist)/RE

Definition at line 127 of file vectorpotentialdipole_compare_with_data.py.

◆ radii

int vectorpotentialdipole_compare_with_data.radii = np.arange(0.1,45,step)*RE

Definition at line 85 of file vectorpotentialdipole_compare_with_data.py.

◆ radiiRE

int vectorpotentialdipole_compare_with_data.radiiRE = radii/RE

Definition at line 87 of file vectorpotentialdipole_compare_with_data.py.

◆ RE

int vectorpotentialdipole_compare_with_data.RE = 6371000.

Definition at line 59 of file vectorpotentialdipole_compare_with_data.py.

◆ res

vectorpotentialdipole_compare_with_data.res = pt.calculations.cut_through_step(vf, [xv[0],0,zv[0]], [xv[-1],0,zv[-1]])

Definition at line 123 of file vectorpotentialdipole_compare_with_data.py.

◆ step

float vectorpotentialdipole_compare_with_data.step = 0.1

Definition at line 67 of file vectorpotentialdipole_compare_with_data.py.

◆ tilt_angle_phi

int vectorpotentialdipole_compare_with_data.tilt_angle_phi = 0.

Definition at line 61 of file vectorpotentialdipole_compare_with_data.py.

◆ tilt_angle_theta

int vectorpotentialdipole_compare_with_data.tilt_angle_theta = 0.

Definition at line 62 of file vectorpotentialdipole_compare_with_data.py.

◆ timefulls

list vectorpotentialdipole_compare_with_data.timefulls = [str(time).rjust(7, '0') for time in times]

Definition at line 52 of file vectorpotentialdipole_compare_with_data.py.

◆ times

list vectorpotentialdipole_compare_with_data.times = [0,10,50,100,200,500]

Definition at line 50 of file vectorpotentialdipole_compare_with_data.py.

◆ vf

list vectorpotentialdipole_compare_with_data.vf = vlsvobj[i]

Definition at line 121 of file vectorpotentialdipole_compare_with_data.py.

◆ vlsvobj

list vectorpotentialdipole_compare_with_data.vlsvobj = []

Definition at line 55 of file vectorpotentialdipole_compare_with_data.py.

◆ xv

int vectorpotentialdipole_compare_with_data.xv = line_start[0]*RE + radii*np.sin(line_phi)*np.cos(line_theta)

Definition at line 91 of file vectorpotentialdipole_compare_with_data.py.

◆ ylims

vectorpotentialdipole_compare_with_data.ylims = np.array(ax.get_ylim())

Definition at line 141 of file vectorpotentialdipole_compare_with_data.py.

◆ yv

int vectorpotentialdipole_compare_with_data.yv = line_start[1]*RE + radii*np.sin(line_phi)*np.sin(line_theta)

Definition at line 92 of file vectorpotentialdipole_compare_with_data.py.

◆ zorder

vectorpotentialdipole_compare_with_data.zorder

Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.

◆ zv

int vectorpotentialdipole_compare_with_data.zv = line_start[2]*RE + radii*np.cos(line_phi)

Definition at line 93 of file vectorpotentialdipole_compare_with_data.py.