- use chunks and dask delayed to keep memory usage efficient
- get_geotop now gets raw data
- add deprecationwarning to get_geotop_raw_within_extent()

- now takes geotop_ds as input to convert to layers (not extent)

- do not use sortby but check diff of y dimension and flip if necessary

- use .values instead of .data to load data into memory if using dask
- undo onno changes for now w missing strat code

- 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

- add function name in all cases
- add line number when level is DEBUG

- use .values to deal with dask arrays
- use get_geotop
- update docstring

…top_to_ml_layers

Add

* fix for #199

* fix

* return da in some situations

---------

Co-authored-by: Ruben Caljé <r.calje@artesia-water.nl>

Small fix for Error: Assigning CRS to a GeoDataFrame without a geometry column is not supported

add fontsize to plot_layers in DatasetCrossSection

* fix for #198

* codacy fix

---------

Co-authored-by: Ruben Caljé <r.calje@artesia-water.nl>

… certain datasets (#202)

* add option in get_ds_empty to specify dimensions

* fix error for vertex grids

---------

Co-authored-by: Ruben Caljé <r.calje@artesia-water.nl>

- fix comments
- add docstrings

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

Improve logging and docstrings

- 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

- replace Exception with more specific types
- add LayerError and MissingValueError
- fix some typos

Improve geotop functionality

Improve reading model output

…a variable with dtype object present in the dataset (#209)

* fix for #208

* better log message

- if steady state select first timestep from evaporation data
- squeeze evaporation data when building stress_period_data

- do not nearest extrpolate if no nans in data
- do not raise exception, just log warning

…d.gwf.drn, but then with a default value of 0 for chd)

The defaults remain the same (although questionable), so this will not alter existing scripts

…ded from a file)

add contour-array flopy func

Updates based on model in Zeeland

Solves #223

…p-from-cache

…of-cachedir-in-write_and_run_model

- add custom crs definition string

And implement in transport notebook

- use nlmod.plot.get_map() where possible
- simplify plotting with nlmod.plot.map_array()
- other nb improvements (e.g. descriptions, text, etc).
- rename kwarg background to backgroundmap in get_map
- improve default axis ticks in get_map
- some typos
- rename layer coord to source_layer in get_concentration_at_gw_surface()

…workspace-instead-of-cachedir-in-write_and_run_model

227 write dataset to model workspace instead of cachedir in write and run model

…t-packages

210 nlmod adding rch and evt packages