### By Xiaohui Yuan (Yuan et al., 1997)

- psmove/psmove2/psmove3 are for Ps/M1/M2, respectively

- To run: psmove \n input.rf \n moved.rf

- To make: use f77 or gfortran (ifort may have problem)
f77 -o psmove psmove.f subs.f saclib.f
f77 -o psmove2 psmove2.f subs.f saclib.f
f77 -o psmove3 psmove3.f subs.f saclib.f
  OR:
gfortran -o psmove psmove.f subs.f saclib.f
gfortran -o psmove2 psmove2.f subs.f saclib.f
gfortran -o psmove3 psmove3.f subs.f saclib.f

- Line 45 in psmove.f is for moveout phase:
  1 for Ps, 2 for M1 (PpPs), 3 for M2 (PpSs+PsPs)

- May need to modify Line 380 for ray parameter,
  which is currently set as 0.06 s/km
