In this paper, we introduce, model, and solve a clustered resource allocation and routing problem for humanitarian aid distribution in the event of an earthquake and subsequent tsunami. First, for the preparedness stage, we build a set of clusters to identify, classify, sort, focus, and prioritize the aid distribution. The clusters are built with k-means method and a modified version of the capacitated p-median model. Each cluster has a set of beneficiaries and candidate delivery aid points. Second, vehicle routes are strategically determined to visit the clusters for the response stage. A mixed integer linear programming model is presented to determine efficient vehicle routes, minimizing the aid distribution times. A vulnerability index is added to our model to prioritize aid distribution. A case study is solved for the city of Iquique, Chile.