cosmo vol coordinates hdr

function [transformed, outside] = cosmo_vol_coordinates(ds, coords)