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)