create definition assign_to_top_active_layer()
Often it is usefull to assign a 2d raster to the top_active cells in the 3D model. In my script i usually do something like this:
raster_2d_bot = imod.idf.open(path_raster_2d) raster_test = ~np.isnan(bnd.where((top > raster_2d_bot) & (bot <= raster_2d_bot))) raster_3d_bot = xr.full_like(bnd, raster_2d_bot).where(raster_test)
Maybe something like that can be standardized. or is it allready?