This code supports the paper: Calafat, F. M., Frederikse, T., and Horsburgh, K. (2022). The Sources of Sea-Level Changes in the Mediterranean Sea since 1960, Journal of Geophysical Research Oceans, under review. Please cite this paper when using this code. Overview: The code implements spatiotemporal Bayesian hierarchical modeling of sea-level observations from tide gauges and satellite altimetry in combination with sea-level fingerprints of contemporary land-mass changes (GRD effects) and glacial isostatic adjustment (GIA). It makes estimates of sea-level changes as well as rates of change and quantifies the individual contribution from sterodynamic changes, GRD effects, and GIA to sea-level changes at any arbitrary location (observed or unobserved). Inference in the model is performed using Hamiltonian Monte Carlo as implemented by the Stan probabilistic programming language. Input data to run the Bayesian hierarchical model and needed to replicate the results of the paper are available at https://doi.org/10.5281/zenodo.6797928.