diff --git a/seispy/psrayp.py b/seispy/psrayp.py
index 1d7e6177..db9d243d 100644
--- a/seispy/psrayp.py
+++ b/seispy/psrayp.py
@@ -50,7 +50,6 @@ def taup_rayp(self, this_dis=50, this_dep=10, taup='taup_time'):
 
     def get_rayp(self):
         for i in range(self.dis.shape[0]):
-            print('{}'.format(self.dis[i]))
             for j in range(self.dep.shape[0]):
                 self.rayp[i, j, :] = self.taup_rayp(this_dis=self.dis[i], this_dep=self.dep[j])[1]
 
diff --git a/seispy/rfcorrect.py b/seispy/rfcorrect.py
index afe89389..fdf80de0 100644
--- a/seispy/rfcorrect.py
+++ b/seispy/rfcorrect.py
@@ -716,7 +716,6 @@ def psrf_1D_raytracing(stadatar, YAxisRange, velmod='iasp91', srayp=None, sphere
         for i in range(stadatar.ev_num):
             tps[i], x_s, x_p, raylength_s[i], raylength_p[i] = xps_tps_map(
                 dep_mod, stadatar.rayp[i], stadatar.rayp[i], is_raylen=True, sphere=sphere, phase=phase)
-            print(type(stadatar.bazi[i]), type(rad2deg(x_s)))
             pplat_s[i], pplon_s[i] = latlon_from(stadatar.stla, stadatar.stlo, stadatar.bazi[i], rad2deg(x_s))
             pplat_p[i], pplon_p[i] = latlon_from(stadatar.stla, stadatar.stlo, stadatar.bazi[i], rad2deg(x_p))
     elif isinstance(srayp, str) or isinstance(srayp, np.lib.npyio.NpzFile):
diff --git a/setup.py b/setup.py
index 1aaa05b1..dafe2770 100644
--- a/setup.py
+++ b/setup.py
@@ -6,7 +6,7 @@
     long_description = fh.read()
 
 
-VERSION = "1.3.7"
+VERSION = "1.3.8"
 setup(name='python-seispy',
       version=VERSION,
       author='Mijian Xu',