Add VDI profiles based on lpagg

src/demandlib/bdew/heat_building.py

@@ -241,7 +241,7 @@ def get_sf_values(self, filename="shlp_hour_factors.csv"):
             + (["weekday"] if not residential else [])
         )
 
-        sf_mat = sf_mat.drop(drop_cols, 1)
+        sf_mat = sf_mat.drop(labels=drop_cols, axis=1)
 
         # Determine the h values
         length = len(self.temperature)

src/demandlib/vdi/__init__.py

@@ -1,2 +1,3 @@
 from .dwd_try import find_try_region  # noqa: F401
+from .dwd_try import read_dwd_weather_file  # noqa: F401
 from .regions import Region  # noqa: F401

src/demandlib/vdi/dwd_try.py

@@ -46,8 +46,14 @@ def find_try_region(longitude, latitude):
     return int(try_map.loc[try_map.contains(my_point), "TRY_code"])
 
 
-def read_dwd_weather_file(weather_file_path):
-    """Read and interpolate "DWD Testreferenzjahr" files."""
+def read_dwd_weather_file(weather_file_path=None, try_region=None):
+    """Read and interpolate 'DWD Testreferenzjahr' files."""
+    if weather_file_path is None:
+        weather_file_path = os.path.join(
+            os.path.dirname(__file__),
+            "resources_weather",
+            "TRY2010_{:02d}_Jahr.dat".format(try_region),
+        )
     # The comments in DWD files before the header are not commented out.
     # Thus we have to search for the line with the header information:
     header_row = None