This set of subroutines should be arranged into into one folder.
MatLab File ModelDataPrepare.m produces the output to txt-file
     X        Y        Z       I       D        F    VGPLat  VGPLong
  +20.394   +4.906  +46.538  +65.74  +13.52   51.047  +69.73  288.13
  +20.854   +4.121  +46.656  +65.51  +11.18   51.270  +69.81  292.76
  +21.285   +3.278  +46.853  +65.31   +8.75   51.566  +69.89  297.57
..........

All necessary settings are at the header section of ModelDataPrepare.m
%%%%%%%%%%%%%%%%%%%%%%% Data Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MTimesteps=2000;%
Lat=67; % degrees
Long=135; % degrees
MaxDeg=7; %% maximal degree of spherical harmonic expansion
modeltype='QC';
TimeStep = 10; %% years
cor_type = 'gaussian';
%%%%%%%%%%%%%%%%%%%%%%% OutPut File Parameters %%%%%%%%%%%%%%%%%%%%%%%%%%%%

Other subroutines are auxiliary.
The corresponding theory can be found in Variat.pdf

