imod.wq.BasicFlow creates a new dimension if presented with array of floats
in this statement:
imod.wq.BasicFlow(ibound=xr.ones_like(like, dtype=int), top=0., bottom=0.-thick.cumsum(), starting_head=0.)
where thick as a np.array, results in:
<xarray.BasicFlow>
Dimensions: (bottom: 35, layer: 35, x: 499, y: 1)
Coordinates:
* layer (layer) int32 1 2 3 4 5 6 7 8 9 ... 28 29 30 31 32 33 34 35
* y (y) float64 0.5
* x (x) float64 5.0 15.0 25.0 ... 4.965e+03 4.975e+03 4.985e+03
dx int32 10
dy float64 -1.0
* bottom (bottom) float64 -1.0 -2.0 -3.0 -4.0 ... -240.0 -265.0 -290.0
Data variables:
ibound (layer, y, x) int32 1 1 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1 1
top float64 0.0
starting_head float64 0.0
inactive_head float64 1e+30
with a new dimension created for bottom, instead of it being automatically assigned to the layer coord.
Edited by Huite Bootsma