function [transformed, outside] = cosmo_vol_coordinates(ds, coords)
cosmo unflatten hdr
cosmo vol grid convert hdr