diff --git a/code/ysolve.f90 b/code/ysolve.f90
index 0f3d650..aa221e5 100644
--- a/code/ysolve.f90
+++ b/code/ysolve.f90
@@ -132,7 +132,7 @@
        END SUBROUTINE calc_diff
       END INTERFACE
 
-      INTEGER :: i,j,k,savenum
+      INTEGER :: i
       REAL :: pflame,pflold,delf=0.
       REAL :: residue = 0.
 
@@ -280,7 +280,7 @@
      END SUBROUTINE write_sd
 
      SUBROUTINE save_final_field
-      INTEGER :: i,j,k
+      INTEGER :: i
 
       OPEN (305,FILE='sfield.bin',form='unformatted',status='unknown')
       DO i=1,nx
@@ -363,9 +363,7 @@
 
      SUBROUTINE BASE_FLAME_PROFILE(x)
       INTEGER :: i, ifl, si
-      REAL :: xi
       REAL :: x(nx)
-      REAL :: max_ysum
 
       ! initial flame thickness (0.2 pi n grids)
       ifl=INT(2.0*pi*0.1/hx)
@@ -519,7 +517,6 @@
       REAL, INTENT(INOUT),DIMENSION(:,:) :: ri,r1,r2
       REAL, INTENT(OUT),DIMENSION(:,:) :: f
 
-      INTEGER :: i,j,k
       REAL :: at,bt
 
        CALL calc_diff(ri)
@@ -562,9 +559,8 @@
       REAL, INTENT(OUT),DIMENSION(:,:) :: f
       REAL, DIMENSION(nsp,nx) :: ux, dux, d2ux
       REAL, DIMENSION(nx) :: dxdm
-      INTEGER :: i,j,k
+      INTEGER :: i
       REAL :: wrate,Ly,Dy,Lt,Dt
-      REAL :: wrate1, wrate2
 
 ! x-direction
       DO i=1,nx
@@ -604,8 +600,8 @@
       REAL, INTENT(IN),DIMENSION(:,:) :: r1
       REAL, INTENT(OUT),DIMENSION(:,:) :: f
       REAL, DIMENSION(3,nx) :: ux, dux, d2ux
-      INTEGER :: i,j,k
-      REAL :: wrate,Ly,Dy,Lt,Dt
+      INTEGER :: i
+      REAL :: Ly,Dy,Lt,Dt
       REAL :: wrate1, wrate2
 
 ! x-direction
@@ -620,7 +616,6 @@
       CALL d2fnonp(nx,hx,ux(:,:),d2ux(:,:),3,1)
 
       DO i=1,nx
-        wrate=rate_1step(ux(1,i), ux(3,i))
         wrate1=rate1_2step(ux(1,i), ux(2,i), ux(3,i))
         wrate2=rate2_2step(ux(2,i), ux(3,i))