diff --git a/src/CalSurfG.f90 b/src/CalSurfG.f90
index 788a0a2..de7f308 100644
--- a/src/CalSurfG.f90
+++ b/src/CalSurfG.f90
@@ -1108,7 +1108,7 @@ subroutine CalSurfG(nx,ny,nz,nparpi,vels,iw,rw,col,dsurf, &
     igr=0
 !    print*,kmax
     call caldespersion(nx,ny,nz,vels,pvRc, &
-        iwave,igr,kmax,tRg,depz,minthk)
+        iwave,igr,kmaxRg,tRg,depz,minthk)
     igr=1
     print*,'Rayleigh wave group velocity depth kernel'
     call depthkernel(nx,ny,nz,vels,pvRg,sen_vsRg,sen_vpRg, &
@@ -1126,7 +1126,7 @@ subroutine CalSurfG(nx,ny,nz,nparpi,vels,iw,rw,col,dsurf, &
     iwave=1
     igr=0
     call caldespersion(nx,ny,nz,vels,pvLc, &
-        iwave,igr,kmax,tLg,depz,minthk)
+        iwave,igr,kmaxLg,tLg,depz,minthk)
     igr=1
     call depthkernel(nx,ny,nz,vels,pvLg,sen_vsLg,sen_vpLg, &
       sen_rhoLg,iwave,igr,kmaxLg,tLg,depz,minthk)
diff --git a/src/main.f90 b/src/main.f90
index e9a4334..c3d6bfb 100644
--- a/src/main.f90
+++ b/src/main.f90
@@ -174,9 +174,9 @@ program SurfTomo
                 if (checkstat > 0) stop 'error allocating RP'
                 read(10,*)(tRc(i),i=1,kmaxRc)
                 write(*,*)'Rayleigh wave phase velocity used,periods:(s)'
-                write(*,'(50f7.1)')(tRc(i),i=1,kmaxRc)
+                write(*,'(50f7.2)')(tRc(i),i=1,kmaxRc)
                 write(66,*)'Rayleigh wave phase velocity used,periods:(s)'
-                write(66,'(50f7.1)')(tRc(i),i=1,kmaxRc)
+                write(66,'(50f7.2)')(tRc(i),i=1,kmaxRc)
         endif
         read(10,*)kmaxRg
         if(kmaxRg.gt.0)then
@@ -184,9 +184,9 @@ program SurfTomo
                 if (checkstat > 0) stop 'error allocating RP'
                 read(10,*)(tRg(i),i=1,kmaxRg)
                 write(*,*)'Rayleigh wave group velocity used,periods:(s)'
-                write(*,'(50f7.1)')(tRg(i),i=1,kmaxRg)
+                write(*,'(50f7.2)')(tRg(i),i=1,kmaxRg)
                 write(66,*)'Rayleigh wave group velocity used,periods:(s)'
-                write(66,'(50f7.1)')(tRg(i),i=1,kmaxRg)
+                write(66,'(50f7.2)')(tRg(i),i=1,kmaxRg)
         endif
         read(10,*)kmaxLc
         if(kmaxLc.gt.0)then
@@ -194,9 +194,9 @@ program SurfTomo
                 if (checkstat > 0) stop 'error allocating RP'
                 read(10,*)(tLc(i),i=1,kmaxLc)
                 write(*,*)'Love wave phase velocity used,periods:(s)'
-                write(*,'(50f7.1)')(tLc(i),i=1,kmaxLc)
+                write(*,'(50f7.2)')(tLc(i),i=1,kmaxLc)
                 write(66,*)'Love wave phase velocity used,periods:(s)'
-                write(66,'(50f7.1)')(tLc(i),i=1,kmaxLc)
+                write(66,'(50f7.2)')(tLc(i),i=1,kmaxLc)
         endif
         read(10,*)kmaxLg
         if(kmaxLg.gt.0)then
@@ -204,9 +204,9 @@ program SurfTomo
                 if (checkstat > 0) stop 'error allocating RP'
                 read(10,*)(tLg(i),i=1,kmaxLg)
                 write(*,*)'Love wave group velocity used,periods:(s)'
-                write(*,'(50f7.1)')(tLg(i),i=1,kmaxLg)
+                write(*,'(50f7.2)')(tLg(i),i=1,kmaxLg)
                 write(66,*)'Love wave group velocity used,periods:(s)'
-                write(66,'(50f7.1)')(tLg(i),i=1,kmaxLg)
+                write(66,'(50f7.2)')(tLg(i),i=1,kmaxLg)
         endif
         read(10,*)ifsyn
         read(10,*)noiselevel
diff --git a/src/surfdisp96.f b/src/surfdisp96.f
index 78cc0b9..711c898 100644
--- a/src/surfdisp96.f
+++ b/src/surfdisp96.f
@@ -308,7 +308,7 @@ c -----         print *, itst,iq,t(k),t1a,t1b,cc0,cc1,gvel
         if(iverb(ifunc).eq.0)then
             iverb(ifunc) = 1
           write(LOT,*)'improper initial value in disper - no zero found'
-          write(*,*)'WARNING:improper initial value in disper - no zero found'
+          !write(*,*)'WARNING:improper initial value in disper - no zero found'
         write(LOT,*)'in fundamental mode '
         write(LOT,*)'This may be due to low velocity zone '
         write(LOT,*)'causing reverse phase velocity dispersion, '