def test_read_ccmp(self):
wind = filament.Wind()
wind.read_from_ccmp(self.windfile)
len(wind.lon) == 1440
len(wind.lat) == 628
wind.lon[0] == 0.125
wind.lon[1] == 0.375
wind.lat[0] == -78.375
wind.lat[1] == -78.125
wind.u.shape == (1, 628, 1440)
wind.u.mean() == -0.11586497
def test_read_ascat_domain(self):
wind = filament.Wind()
r = wind.read_ascat(self.windfile, self.domain)
self.assertTrue(r)
self.assertEqual(len(wind.lon), 6134)
self.assertEqual(len(wind.lat), 6134)
self.assertEqual(wind.lon[-1], -47.16028)
self.assertEqual(wind.lat[600], 4.96835)
self.assertEqual(wind.speed.min(), 7.01)
self.assertEqual(len(wind.speed.compressed()), 6134)
self.assertFalse(np.ma.is_masked(wind.u[999]))
self.assertFalse(np.ma.is_masked(wind.v[1999]))
self.assertEqual(wind.u[124], -7.4646841483642)
def test_read_ascat(self):
wind = filament.Wind()
r = wind.read_ascat(self.windfile)
self.assertTrue(r)
self.assertEqual(len(wind.lon), 136513)
self.assertEqual(len(wind.lat), 136513)
self.assertAlmostEqual(wind.lon[20], -52.24672999)
self.assertEqual(wind.lat[200], 1.31254)
self.assertEqual(wind.speed.mean(), 8.85660589984728)
self.assertEqual(len(wind.speed.compressed()), 62205)
self.assertTrue(np.ma.is_masked(wind.u[999]))
self.assertTrue(np.ma.is_masked(wind.v[1999]))
self.assertEqual(wind.u[22222], 3.603894954004558)
def test_read_quikcat_domain(self):
wind = filament.Wind()
r = wind.read_from_quikscat(self.windfile, self.domain)
self.assertTrue(r)
self.assertTrue(wind.lon.shape == (29198, ))
self.assertTrue(wind.lat.shape == (29198, ))
self.assertAlmostEqual(wind.lon.min(), -124.99933, places=5)
self.assertAlmostEqual(wind.lon.mean(), -115.39173, places=5)
self.assertAlmostEqual(wind.lat[0], 49.949474, places=6)
self.assertAlmostEqual(wind.lat.max(), 49.999313, places=6)
self.assertAlmostEqual(wind.u.max(), 13.594069, places=6)
self.assertTrue(np.ma.is_masked(wind.u[0]))
self.assertTrue(wind.u.mean(), 2.214242427089597)
self.assertTrue(wind.v.mean(), -5.515122331629585)
self.assertTrue(np.ma.is_masked(wind.angle[1]))
def test_read_quikcat(self):
wind = filament.Wind()
r = wind.read_from_quikscat(self.windfile)
self.assertTrue(r)
self.assertTrue(wind.lon.shape == (242935, ))
self.assertTrue(wind.lat.shape == (242935, ))
self.assertAlmostEqual(wind.lon.min(), -179.99792, places=5)
self.assertAlmostEqual(wind.lon.mean(), -120.590515, places=6)
self.assertAlmostEqual(wind.lat[0], 70.397896, places=6)
self.assertAlmostEqual(wind.lat.max(), 89.866196, places=6)
self.assertAlmostEqual(wind.u.max(), 29.148708, places=6)
self.assertTrue(np.ma.is_masked(wind.u[0]))
self.assertTrue(wind.u.mean(), 1.5385874209615187)
self.assertTrue(wind.v.mean(), 1.2612295160130862)
self.assertTrue(np.ma.is_masked(wind.angle[1]))
def test_read_ascat_domain(self):
wind = filament.Wind()
r = wind.read_ascat(self.windfile, self.domain2)
self.assertFalse(r)