- add kwarg delayed, if False, load data into memory, else return data array with delayed dask arrays
- add kwarg chunked, if True, chunk data array with chunks="auto"
- add x,y data to vertex grid data array

- fix comments
- add docstrings

- fix where command for dry/noflow
- improve support for loading output without gwf or ds

- new folder mfoutput
- add new flopy binary read functions that support multithreading (binaryfile.py)
- separate reading budget and head output
- modify gwf.output and gwt.output to use new methods
- split logic into multiple reusable functions
- add support for grb files
- add/improve tests
- add method to obtain dims, coords from modelgrid object

- codacy
- fix comments @OnnoEbbens
- some additional fixes