For headobs.py, the coverage testing showed that the read_metadata() method was started but not really implemented yet. The code also has the aquifer names hard-coded into it. I propose:
complete the read_metadata() method and move the code for that out of get_data(), get_data() can call read_metadata() as needed
get rid of the get_layer() method since this preprocessor doesn't have a modflow object passed to it
replace the code in get_data() that does the aquifer names by importing a dictionary
using from mapgwm.lookups import aquifer_cd_name
in lookups.py add the code to read in the aquifer names yaml file into a dictionary that can be imported.
For headobs.py, the coverage testing showed that the read_metadata() method was started but not really implemented yet. The code also has the aquifer names hard-coded into it. I propose: