def test_all_data_single_pressure_iastate(): """Test Iowa State request for all data at a single pressure level.""" df = IAStateUpperAir.request_all_data(datetime(1999, 5, 4, 0), pressure=500) idx = df.loc[df.station == 'KDDC'].index[0] assert (df['time'][idx] == datetime(1999, 5, 4, 0)) assert (df['station'][idx] == 'KDDC') assert_almost_equal(df['pressure'][idx], 500.0, 2) assert_almost_equal(df['height'][idx], 5606.0, 2) assert_almost_equal(df['temperature'][idx], -17.2, 2) assert_almost_equal(df['dewpoint'][idx], -22.9, 2) assert_almost_equal(df['u_wind'][idx], 30.834774, 2) assert_almost_equal(df['v_wind'][idx], 35.47135, 2) assert_almost_equal(df['speed'][idx], 47.0, 1) assert_almost_equal(df['direction'][idx], 221.0, 1) assert (df.units['pressure'] == 'hPa') assert (df.units['height'] == 'meter') assert (df.units['temperature'] == 'degC') assert (df.units['dewpoint'] == 'degC') assert (df.units['u_wind'] == 'knot') assert (df.units['v_wind'] == 'knot') assert (df.units['speed'] == 'knot') assert (df.units['direction'] == 'degrees') assert (df.units['station'] is None) assert (df.units['time'] is None)
def test_iastate(): """Test that we are properly parsing data from the Iowa State archive.""" df = IAStateUpperAir.request_data(datetime(1999, 5, 4, 0), 'OUN') assert (df['time'][0] == datetime(1999, 5, 4, 0)) assert (df['station'][0] == 'KOUN') assert_almost_equal(df['pressure'][6], 872.7, 2) assert_almost_equal(df['height'][6], 1172.0, 2) assert_almost_equal(df['temperature'][6], 18.2, 2) assert_almost_equal(df['dewpoint'][6], 15.1, 2) assert_almost_equal(df['u_wind'][6], 4.631, 2) assert_almost_equal(df['v_wind'][6], 37.716, 2) assert_almost_equal(df['speed'][6], 38.0, 1) assert_almost_equal(df['direction'][6], 187.0, 1) assert (df.units['pressure'] == 'hPa') assert (df.units['height'] == 'meter') assert (df.units['temperature'] == 'degC') assert (df.units['dewpoint'] == 'degC') assert (df.units['u_wind'] == 'knot') assert (df.units['v_wind'] == 'knot') assert (df.units['speed'] == 'knot') assert (df.units['direction'] == 'degrees') assert (df.units['station'] is None) assert (df.units['time'] is None)
def test_high_alt_iastate(): """Test Iowa State data that starts at pressure less than 925 hPa.""" df = IAStateUpperAir.request_data(datetime(2010, 12, 9, 12), 'BOI') assert (df['time'][0] == datetime(2010, 12, 9, 12)) assert (df['station'][0] == 'KBOI') assert_almost_equal(df['pressure'][0], 919.0, 2) assert_almost_equal(df['height'][0], 871.0, 2) assert_almost_equal(df['temperature'][0], -0.1, 2) assert_almost_equal(df['dewpoint'][0], -0.2, 2) assert_almost_equal(df['u_wind'][0], 2.598, 2) assert_almost_equal(df['v_wind'][0], 1.500, 2) assert_almost_equal(df['speed'][0], 3.0, 1) assert_almost_equal(df['direction'][0], 240.0, 1) assert (df.units['pressure'] == 'hPa') assert (df.units['height'] == 'meter') assert (df.units['temperature'] == 'degC') assert (df.units['dewpoint'] == 'degC') assert (df.units['u_wind'] == 'knot') assert (df.units['v_wind'] == 'knot') assert (df.units['speed'] == 'knot') assert (df.units['direction'] == 'degrees') assert (df.units['station'] is None) assert (df.units['time'] is None)
def test_interp(): """Test the interpolation flag for the IA State retrieval.""" df = IAStateUpperAir.request_data(datetime(2011, 4, 14, 18), 'OUN', interp_nans=True) assert not df['temperature'].isnull().any() assert not df['dewpoint'].isnull().any()
def test_high_alt_iastate(): """Test Iowa State data that starts at pressure less than 925 hPa.""" df = IAStateUpperAir.request_data(datetime(2010, 12, 9, 12), 'BOI') assert(df['time'][0] == datetime(2010, 12, 9, 12)) assert(df['station'][0] == 'KBOI') assert_almost_equal(df['pressure'][0], 919.0, 2) assert_almost_equal(df['height'][0], 871.0, 2) assert_almost_equal(df['temperature'][0], -0.1, 2) assert_almost_equal(df['dewpoint'][0], -0.2, 2) assert_almost_equal(df['u_wind'][0], 2.598, 2) assert_almost_equal(df['v_wind'][0], 1.500, 2) assert_almost_equal(df['speed'][0], 3.0, 1) assert_almost_equal(df['direction'][0], 240.0, 1) assert(df.units['pressure'] == 'hPa') assert(df.units['height'] == 'meter') assert(df.units['temperature'] == 'degC') assert(df.units['dewpoint'] == 'degC') assert(df.units['u_wind'] == 'knot') assert(df.units['v_wind'] == 'knot') assert(df.units['speed'] == 'knot') assert(df.units['direction'] == 'degrees') assert(df.units['station'] is None) assert(df.units['time'] is None)
def test_all_data_single_pressure_iastate(): """Test Iowa State request for all data at a single pressure level.""" df = IAStateUpperAir.request_all_data(datetime(1999, 5, 4, 0), pressure=500) idx = df.loc[df.station == 'KDDC'].index[0] assert(df['time'][idx] == datetime(1999, 5, 4, 0)) assert(df['station'][idx] == 'KDDC') assert_almost_equal(df['pressure'][idx], 500.0, 2) assert_almost_equal(df['height'][idx], 5606.0, 2) assert_almost_equal(df['temperature'][idx], -17.2, 2) assert_almost_equal(df['dewpoint'][idx], -22.9, 2) assert_almost_equal(df['u_wind'][idx], 30.834774, 2) assert_almost_equal(df['v_wind'][idx], 35.47135, 2) assert_almost_equal(df['speed'][idx], 47.0, 1) assert_almost_equal(df['direction'][idx], 221.0, 1) assert(df.units['pressure'] == 'hPa') assert(df.units['height'] == 'meter') assert(df.units['temperature'] == 'degC') assert(df.units['dewpoint'] == 'degC') assert(df.units['u_wind'] == 'knot') assert(df.units['v_wind'] == 'knot') assert(df.units['speed'] == 'knot') assert(df.units['direction'] == 'degrees') assert(df.units['station'] is None) assert(df.units['time'] is None)
def test_iastate(): """Test that we are properly parsing data from the Iowa State archive.""" df = IAStateUpperAir.request_data(datetime(1999, 5, 4, 0), 'OUN') assert(df['time'][0] == datetime(1999, 5, 4, 0)) assert(df['station'][0] == 'KOUN') assert_almost_equal(df['pressure'][6], 872.7, 2) assert_almost_equal(df['height'][6], 1172.0, 2) assert_almost_equal(df['temperature'][6], 18.2, 2) assert_almost_equal(df['dewpoint'][6], 15.1, 2) assert_almost_equal(df['u_wind'][6], 4.631, 2) assert_almost_equal(df['v_wind'][6], 37.716, 2) assert_almost_equal(df['speed'][6], 38.0, 1) assert_almost_equal(df['direction'][6], 187.0, 1) assert(df.units['pressure'] == 'hPa') assert(df.units['height'] == 'meter') assert(df.units['temperature'] == 'degC') assert(df.units['dewpoint'] == 'degC') assert(df.units['u_wind'] == 'knot') assert(df.units['v_wind'] == 'knot') assert(df.units['speed'] == 'knot') assert(df.units['direction'] == 'degrees') assert(df.units['station'] is None) assert(df.units['time'] is None)
def test_all_data_all_levels_iastate(): """Test Iowa State request for all data at all levels.""" df = IAStateUpperAir.request_all_data(datetime(1999, 5, 4, 0)) idx = df.loc[df.station == 'KDDC'].index[12] assert(df['time'][idx] == datetime(1999, 5, 4, 0)) assert(df['station'][idx] == 'KDDC') assert_almost_equal(df['pressure'][idx], 700.0, 2) assert_almost_equal(df['height'][idx], 2969.0, 2) assert_almost_equal(df['temperature'][idx], 5.8, 2) assert_almost_equal(df['dewpoint'][idx], -11.7, 2) assert_almost_equal(df['u_wind'][idx], 33.234, 2) assert_almost_equal(df['v_wind'][idx], 33.234, 2) assert_almost_equal(df['speed'][idx], 47.0, 1) assert_almost_equal(df['direction'][idx], 225.0, 1) assert(df.units['pressure'] == 'hPa') assert(df.units['height'] == 'meter') assert(df.units['temperature'] == 'degC') assert(df.units['dewpoint'] == 'degC') assert(df.units['u_wind'] == 'knot') assert(df.units['v_wind'] == 'knot') assert(df.units['speed'] == 'knot') assert(df.units['direction'] == 'degrees') assert(df.units['station'] is None) assert(df.units['time'] is None)
def test_all_data_all_levels_iastate(): """Test Iowa State request for all data at all levels.""" df = IAStateUpperAir.request_all_data(datetime(1999, 5, 4, 0)) idx = df.loc[df.station == 'KDDC'].index[12] assert (df['time'][idx] == datetime(1999, 5, 4, 0)) assert (df['station'][idx] == 'KDDC') assert_almost_equal(df['pressure'][idx], 700.0, 2) assert_almost_equal(df['height'][idx], 2969.0, 2) assert_almost_equal(df['temperature'][idx], 5.8, 2) assert_almost_equal(df['dewpoint'][idx], -11.7, 2) assert_almost_equal(df['u_wind'][idx], 33.234, 2) assert_almost_equal(df['v_wind'][idx], 33.234, 2) assert_almost_equal(df['speed'][idx], 47.0, 1) assert_almost_equal(df['direction'][idx], 225.0, 1) assert (df.units['pressure'] == 'hPa') assert (df.units['height'] == 'meter') assert (df.units['temperature'] == 'degC') assert (df.units['dewpoint'] == 'degC') assert (df.units['u_wind'] == 'knot') assert (df.units['v_wind'] == 'knot') assert (df.units['speed'] == 'knot') assert (df.units['direction'] == 'degrees') assert (df.units['station'] is None) assert (df.units['time'] is None)
def read_upper_air_iastate(date): try: data = IAStateUpperAir().request_all_data(date) except ValueError: #print("ERROR: no data available for ",date) #print(" Trying previous synoptic time") return None #date = date - timedelta(hours = 12) #data = IAStateUpperAir().request_all_data(date) data = add_station_lat_lon(data) return data
def test_no_future_data_iastate(): """Test Iowa State data when future data are requested.""" with pytest.raises(ValueError): IAStateUpperAir.request_data(datetime(2999, 12, 9, 12), 'BOI')
def test_no_data_iastate(): """Test Iowa State data when no data are available.""" with pytest.raises(ValueError): IAStateUpperAir.request_data(datetime(2010, 12, 9, 1), 'BOI')
def test_no_future_data_with_pressure_iastate(): """Test Iowa State data when future data are requested with the pressure argument.""" with pytest.raises(ValueError): IAStateUpperAir.request_all_data(datetime(2999, 12, 9, 12), pressure=850.)
# Set a date and time for upper-air observations (should only be 00 or 12 # UTC for the hour). # # Request all data from Iowa State using the Siphon package. The result is # a pandas DataFrame containing all of the sounding data from all # available stations. # # Set date for desired UPA data today = datetime.utcnow() # Go back one day to ensure data availability date = datetime(today.year, today.month, today.day, 0) - timedelta(days=1) # Request data using Siphon request for data from Iowa State Archive data = IAStateUpperAir.request_all_data(date) ###################################################################### # Subset Observational Data # ------------------------- # # From the request above will give all levels from all radisonde sites # available through the service. For plotting a pressure surface map there # is only need to have the data from that level. Below the data is subset # and a few parameters set based on the level chosen. Additionally, the # station information is obtained and latitude and longitude data is added # to the DataFrame. # level = 500
# Import the Iowa state simple web service upper air object from siphon.simplewebservice.iastate import IAStateUpperAir # Create the datetime and station variables you'll need request_time = datetime(2011, 4, 14, 18) station = 'OUN' # Make the request for the data df = IAStateUpperAir.request_data(request_time, station) # Attach units to the data sounding = pandas_dataframe_to_unit_arrays(df)