Netcdf4 yields masked arrays. This gave sg problems: now fixed by converting them to ndarrays with Nan's. Added functionality where if coord name is inferred from ax name, the coord name is called ax_crd: e.g. X_crd, to distinguish (new Levitus data).
updated version in setup.py