ds= requester(
"SW_EXPT_EFIB_TCT02", #Mag B, high resolution, 50Hz B (Magnetic field)
measurements_E, #Magnetic field in NEC coordinates
True,
asynchronous=False,
show_progress=False)
velocity=np.array([ds["VsatN"], ds["VsatE"], ds["VsatC"]]) #We need the velocities of the satellite in NEC so we can get the unit vector and then get our Electric field
velocity_unit=unit_array(velocity) #function that gives the unitary vector
ENEC=np.multiply(np.array([ds["Evx"], ds["Evy"], ds["Evz"]]), velocity_unit)
Gets my E, however, when I plot this it does show the velocities when plotted in lompe.plot by itself but not orange quivers are produced.
model = lompe.Emodel(grid, (cmod.hall, cmod.pedersen)) model.add_data(Eswarm) model.run_inversion(l1 = 1, l2 = 10)
CassandraAuri
commented
3 months ago
Hello, I am using swarm to calculate E
Eused=np.array([ENEC[1],ENEC[0]]) #East North coords=np.array([ ds['Longitude'].to_numpy(),ds['Latitude'].to_numpy()]) Eswarm = lompe.Data(Eused, coordinates = coords, datatype = 'Efield', error = 5e-3, iweight=1.0 )
Gets my E, however, when I plot this it does show the velocities when plotted in lompe.plot by itself but not orange quivers are produced.
model = lompe.Emodel(grid, (cmod.hall, cmod.pedersen)) model.add_data(Eswarm) model.run_inversion(l1 = 1, l2 = 10)https://github.com/CassandraAuri/Physics_Work/blob/Deployment-version/Aurora_Work/Data_handling_with_Lompe.ipynb Here is the jupyter notebook