The program read.mk12wm13.f is a fortran program that asks you for the coordinates and depth of a point in the mantle and gives you as output percent perturbation in mu over rho, kappa over rho, bulk sound velocity, S velocity, P velocity, as predicted by MK12WM13 (Su and Dziewonski, PEPI 1997). You can compile it by simply doing f77 read.mk12wm13.f mk12wm13_kappa.dat and mk12wm13_mu.dat are just the MK12WM13 model. foanis05.mod is PREM and velpert.h is a header file that accompanies read.mk12wm13.f