in well_data.py, remove has-site-cod sites;
can plot confined and unconfined sites separately