fill masked arrays with NaNs at initialization
minor bug fix io.nc
axes now also recognize "indexing" and convert to "position_index" as
appropriate