boundary box can be given to EURAD plotting method
At the moment, the full map for the EURAD-IM output is plotted.
The areas of interest are nevertheless North Rhine-Westphalia and Berlin-Brandenburg.
Therefore, now a hardwired solution cuts the German data within the "eurad_vis.py".
The data should be cut using the values of the CTM output files.
Furthermore, it should be possible to specify a bounding box and the cities to be plotted (via an input file) (so that you can also zoom to North Rhine-Westphalia and Berlin-Brandenburg within Germany).
Adjustments to the source code of eurad_vis.py to show only a map for Germany (hardwired and with a new file h09.csv, which includes coordinates for Düsseldorf and Berlin):
l.44
- if self.is_imported():
- self.grid = int(self.dataset.attrs["dx"])
+ # is_imported does not work as expected...
+ #if self.is_imported():
+ self.grid = int(self.dataset.attrs["dx"])
l.219
- if self.grid == 3 or self.grid == 9:
- # print("3")
+ if self.grid == 3:
self.add_cities_boolean = True
csv_file = "h03.csv"
+ if self.grid == 9:
+ self.add_cities_boolean = True
+ csv_file = "h09.csv"
l.288
- cenlat = self.dataset.cenlat
- cenlon = self.dataset.cenlon
+ cenlat = 51.0
+ cenlon = 10.5
- lat_values = self.dataset.lat.values
- lon_values = self.dataset.lon.values
+ lat_values = self.dataset.lat.values[150:250,230:305]
+ lon_values = self.dataset.lon.values[150:250,230:305]
l. 332
- area = self.dataset[s].isel(time=t, z=self.height)
+ area = self.dataset[s].isel(time=t, z=self.height)[150:250,230:305]
fig = plt.figure()
ax = plt.axes(projection=proj)
+ ax.set_extent([6.0, 15.0, 47.0, 55.0], crs=crs)