def test_find_magequator(self): expected = {'Bmin': array([ 1030.456337, 3444.077016 ])} Bmin_loci = array([[ 2.99935449, 0.005511 , -0.032353 ], [ 2.00289871, -0.00734881, 0.045382]]) actual = ib.find_magequator(self.ticks, self.loci, omnivals=self.omnivals) numpy.testing.assert_almost_equal(expected['Bmin'], actual['Bmin'], decimal=6) numpy.testing.assert_almost_equal(Bmin_loci, actual['loci'].data, decimal=6)
def _find_loss_cone(self, position, time): foot = irb.find_footpoint(time, position, extMag='T89')['Bfoot'] eq = irb.find_magequator(time, position, extMag='T89')['Bmin'] pitch = 90 #for particles mirroring at 100km return np.rad2deg( np.arcsin(np.sqrt(np.sin(np.deg2rad(pitch))**2 * eq / foot)))
def test_find_magequator(self): expected = {'Bmin': array([ 1030.50271503, 3444.0670669 ])} Bmin_loci = array([[ 2.99935449, 0.0055078 , -0.032328 ], [ 2.00289871, -0.00734881, 0.04538073]]) actual = ib.find_magequator(self.ticks, self.loci, omnivals=self.omnivals) numpy.testing.assert_almost_equal(expected['Bmin'], actual['Bmin'], decimal=6) numpy.testing.assert_almost_equal(Bmin_loci, actual['loci'].data, decimal=6)
def to_equatorial(position, time, pitch): """ take in spacepy coord class and ticktock class """ blocal = irb.get_Bfield(time, position, extMag='T89')['Blocal'] beq = irb.find_magequator(time, position, extMag='T89')['Bmin'] eq_pitch = np.arcsin(np.sqrt(np.sin(np.deg2rad(pitch))**2 * beq / blocal)) return np.rad2deg(eq_pitch)
def test_find_magequator(self): expected = {'Bmin': array([1030.50271503, 3444.0670669])} Bmin_loci = array([[2.99935449, 0.0055078, -0.032328], [2.00289871, -0.00734881, 0.04538073]]) actual = ib.find_magequator(self.ticks, self.loci, omnivals=self.omnivals) numpy.testing.assert_almost_equal(expected['Bmin'], actual['Bmin'], decimal=6) numpy.testing.assert_almost_equal(Bmin_loci, actual['loci'].data, decimal=6)
def test_find_magequator(self): expected = {'Bmin': array([1030.4563374, 3444.07701499])} Bmin_loci = [[2.999354, 0.00551125, -0.03235315], [2.00289878, -0.0073488, 0.04538181]] actual = ib.find_magequator(self.ticks, self.loci, omnivals=self.omnivals) numpy.testing.assert_almost_equal(expected['Bmin'], actual['Bmin'], decimal=6) numpy.testing.assert_almost_equal(Bmin_loci, actual['loci'].data, decimal=6)