diff --git a/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py index 03937ee..9310615 100644 --- a/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py +++ b/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py @@ -13,7 +13,7 @@ # Load the shapefile area driver = ogr.GetDriverByName('ESRI Shapefile') -shpfile = driver.Open('shpfile/france.shp') +shpfile = driver.Open('shpfile/italy.shp') AREA = shpfile.GetLayer() # Check if point is inside of shapefile area @@ -55,8 +55,8 @@ def parse_data(filepath): # Retrieve flight level values flightlevels = df.index.get_level_values('lv_AMSL0') # Filter to a specific flight level: - # FL100 = 10000 feet = 3048 meters - df = df.loc[(flightlevels == 3048)] + # FL450 = 45000 feet = 13716 meters + df = df.loc[(flightlevels == 13716)] # Get longitude values from index lons = df.index.get_level_values('lon_0') # Map longitude range from (0 to 360) into (-180 to 180) @@ -98,10 +98,10 @@ def plot_data(data): c=color, s=10, cmap='Spectral_r', - edgecolors='gray', + # edgecolors='gray', linewidths=0.1 ) - plt.title('Clear-Air Turbulence % at FL100') + plt.title('Clear-Air Turbulence % at FL450') plt.colorbar() plt.show() diff --git a/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py index e651dde..95536d1 100644 --- a/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py +++ b/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py @@ -13,7 +13,7 @@ # Load the shapefile area driver = ogr.GetDriverByName('ESRI Shapefile') -shpfile = driver.Open('shpfile/france.shp') +shpfile = driver.Open('shpfile/italy.shp') AREA = shpfile.GetLayer() # Check if point is inside of shapefile area diff --git a/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py index 22c56a9..4d7ed0b 100644 --- a/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py +++ b/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py @@ -13,7 +13,7 @@ # Load the shapefile area driver = ogr.GetDriverByName('ESRI Shapefile') -shpfile = driver.Open('shpfile/france.shp') +shpfile = driver.Open('shpfile/italy.shp') AREA = shpfile.GetLayer() # Check if point is inside of shapefile area @@ -100,7 +100,7 @@ def plot_data(data): c=color, s=10, cmap='coolwarm', - edgecolors='gray', + # edgecolors='gray', linewidths=0.1 ) plt.title('Temperature (C) at 20 hPa') diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.prj b/examples/working_with_grib_data/cropped_area/shpfile/france.prj deleted file mode 100644 index f45cbad..0000000 --- a/examples/working_with_grib_data/cropped_area/shpfile/france.prj +++ /dev/null @@ -1 +0,0 @@ -GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] \ No newline at end of file diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.shp b/examples/working_with_grib_data/cropped_area/shpfile/france.shp deleted file mode 100644 index 166a4d0..0000000 Binary files a/examples/working_with_grib_data/cropped_area/shpfile/france.shp and /dev/null differ diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.shx b/examples/working_with_grib_data/cropped_area/shpfile/france.shx deleted file mode 100644 index a16c3c8..0000000 Binary files a/examples/working_with_grib_data/cropped_area/shpfile/france.shx and /dev/null differ diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.dbf b/examples/working_with_grib_data/cropped_area/shpfile/italy.dbf similarity index 62% rename from examples/working_with_grib_data/cropped_area/shpfile/france.dbf rename to examples/working_with_grib_data/cropped_area/shpfile/italy.dbf index 2cec854..708d516 100644 Binary files a/examples/working_with_grib_data/cropped_area/shpfile/france.dbf and b/examples/working_with_grib_data/cropped_area/shpfile/italy.dbf differ diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.prj b/examples/working_with_grib_data/cropped_area/shpfile/italy.prj new file mode 100644 index 0000000..a30c00a --- /dev/null +++ b/examples/working_with_grib_data/cropped_area/shpfile/italy.prj @@ -0,0 +1 @@ +GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]] \ No newline at end of file diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.shp b/examples/working_with_grib_data/cropped_area/shpfile/italy.shp new file mode 100644 index 0000000..1a5ceff Binary files /dev/null and b/examples/working_with_grib_data/cropped_area/shpfile/italy.shp differ diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.shx b/examples/working_with_grib_data/cropped_area/shpfile/italy.shx new file mode 100644 index 0000000..9208c59 Binary files /dev/null and b/examples/working_with_grib_data/cropped_area/shpfile/italy.shx differ