change reprojection to Dem instead of output floodmap

flood_adapt/adapter/sfincs_adapter.py

@@ -8,7 +8,6 @@
 import hydromt_sfincs.utils as utils
 import numpy as np
 import pandas as pd
-import rioxarray as rxr
 import xarray as xr
 from cht_tide.read_bca import SfincsBoundary
 from cht_tide.tide_predict import predict
@@ -523,6 +522,7 @@ def write_geotiff(self, zsmax, demfile: Path, floodmap_fn: Path):
             UnitTypesLength("meters")
         )
         dem = dem_conversion * self.sf_model.data_catalog.get_rasterdataset(demfile)
+        dem = dem.rio.reproject(self.sf_model.crs)
 
         # determine conversion factor for output floodmap
         floodmap_units = self.site.attrs.sfincs.floodmap_units
@@ -536,11 +536,6 @@ def write_geotiff(self, zsmax, demfile: Path, floodmap_fn: Path):
             hmin=0.01,
             floodmap_fn=str(floodmap_fn),
         )
-        # reproject floodmap to sfincs projection
-        floodmap = rxr.open_rasterio(floodmap_fn, masked=True).squeeze()
-        crs = self.sf_model.crs
-        floodmap = floodmap.rio.reproject(crs)
-        floodmap.rio.to_raster(Path(floodmap_fn))
 
     def downscale_hmax(self, zsmax, demfile: Path):
         # read DEM and convert units to metric units used by SFINCS