Receiver function calculation in seismology