This is an R script written to calculate the relative terrain aspect parameter, which is utilized in the preferential snowfall downscaling scheme for mountainous terrain. The code reads a terrain aspect angle grid from file 'terrain_aspect_grid.txt' and writes the computed relative terrain aspect angle grid to file 'delta_psi_grid.txt'.