|
Vlasiator ebf0dd394 on dev (v5.4.0 + 1054 commits)
|
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 | |
Definition at line 112 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.axes = fig.get_axes() |
Definition at line 83 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.B1 = np.zeros([nr,3]) |
Definition at line 95 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.B2 = np.zeros([nr,3]) |
Definition at line 96 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.B3 = np.zeros([nr,3]) |
Definition at line 97 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.B4 = np.zeros([nr,3]) |
Definition at line 98 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.c |
Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.
| list vectorpotentialdipole_compare_with_data.colors = ['r','g','b','magenta','k'] |
Definition at line 51 of file vectorpotentialdipole_compare_with_data.py.
| list vectorpotentialdipole_compare_with_data.coords = ['x','y','z'] |
Definition at line 109 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.cut = res[0].data |
Definition at line 124 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.fig = plt.figure() |
Definition at line 78 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.fontsize = 20 |
Definition at line 70 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.handles |
Definition at line 148 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.label |
Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.labels |
Definition at line 148 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.line_phi = float(sys.argv[1]) |
Definition at line 41 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.line_start = np.array([0,0,0]) |
Definition at line 65 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.line_theta = 0. * math.pi/180. |
Definition at line 63 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.linestyle |
Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.linewidth = 2 |
Definition at line 68 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.linthresh = 1.e-10 |
Definition at line 69 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.linthreshy |
Definition at line 135 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.nr = len(radii) |
Definition at line 86 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.nsubplots = 3 |
Definition at line 80 of file vectorpotentialdipole_compare_with_data.py.
| str vectorpotentialdipole_compare_with_data.outfilename = "./vecpotdip_compare_"+str(int(line_phi))+".png" |
Definition at line 47 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.pr_coords = res[2].data |
Definition at line 126 of file vectorpotentialdipole_compare_with_data.py.
| vectorpotentialdipole_compare_with_data.pr_dist = res[1].data |
Definition at line 125 of file vectorpotentialdipole_compare_with_data.py.
Definition at line 127 of file vectorpotentialdipole_compare_with_data.py.
Definition at line 85 of file vectorpotentialdipole_compare_with_data.py.
Definition at line 87 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.RE = 6371000. |
Definition at line 59 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| float vectorpotentialdipole_compare_with_data.step = 0.1 |
Definition at line 67 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.tilt_angle_phi = 0. |
Definition at line 61 of file vectorpotentialdipole_compare_with_data.py.
| int vectorpotentialdipole_compare_with_data.tilt_angle_theta = 0. |
Definition at line 62 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| list vectorpotentialdipole_compare_with_data.times = [0,10,50,100,200,500] |
Definition at line 50 of file vectorpotentialdipole_compare_with_data.py.
Definition at line 121 of file vectorpotentialdipole_compare_with_data.py.
| list vectorpotentialdipole_compare_with_data.vlsvobj = [] |
Definition at line 55 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.ylims = np.array(ax.get_ylim()) |
Definition at line 141 of file vectorpotentialdipole_compare_with_data.py.
| 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.
| vectorpotentialdipole_compare_with_data.zorder |
Definition at line 114 of file vectorpotentialdipole_compare_with_data.py.
| 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.