diff --git a/src/scheme.F90 b/src/scheme.F90
index 249913c..8f9eea1 100644
--- a/src/scheme.F90
+++ b/src/scheme.F90
@@ -30,12 +30,12 @@ module scheme
 
   contains
 !
-!======================================================================
+!==============================================================================
 !
-! update: subroutine sweeps over all directions and integrates
-!         the increment of the solution
+! update: subroutine sweeps over all directions and integrates the directional
+!         derivatives of the flux in order to get the increment of solution
 !
-!======================================================================
+!==============================================================================
 !
   subroutine update(u, du, dxi, dyi, dzi)
 
@@ -58,7 +58,7 @@ module scheme
 !
     real, dimension(nvars,ngrids) :: ul, fl
 !
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
     du(:,:,:,:) = 0.0
 
@@ -175,16 +175,16 @@ module scheme
     end do
 #endif /* NDIMS == 3 */
 
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
   end subroutine update
 #ifdef HLL
 !
-!======================================================================
+!===============================================================================
 !
-! hll: subroutine to compute flux approximated by HLL method
+! hll: subroutine computes the approximated flux using HLL method
 !
-!======================================================================
+!===============================================================================
 !
   subroutine hll(m, n, uc, f)
 
@@ -206,7 +206,7 @@ module scheme
     real, dimension(n)   :: cl, cr
     real                 :: al, ar, ap, div
 !
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
 #ifdef CONREC
 ! reconstruct left and right states of conserved variables
@@ -270,14 +270,14 @@ module scheme
 !
     f(:,2:n) = - fx(:,2:n) + fx(:,1:n-1)
 
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
   end subroutine hll
 #endif /* HLL */
 !
 !===============================================================================
 !
-! fluxspeed: subroutine to compute fluxes and speeds
+! fluxspeed: subroutine computes fluxes and speeds for a given set of equations
 !
 !===============================================================================
 !
@@ -299,7 +299,7 @@ module scheme
     integer :: i
     real    :: cs
 !
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
 ! sweep over all points
 !
@@ -330,13 +330,13 @@ module scheme
 #endif /* ISO */
     enddo
 
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
   end subroutine fluxspeed
 !
 !===============================================================================
 !
-! cons2prim: subroutine to calculate primitive variables from conservative ones
+! cons2prim: subroutine converts primitive variables to conservative
 !
 !===============================================================================
 !
@@ -357,7 +357,7 @@ module scheme
     integer :: i
     real    :: dni
 !
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
     do i = 1, n
       dni    = 1.0 / u(1,i)
@@ -373,13 +373,13 @@ module scheme
 #endif /* ADI */
     enddo
 
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
   end subroutine cons2prim
 !
 !===============================================================================
 !
-! prim2cons: subroutine to calculate conservative variables from primitive ones
+! prim2cons: subroutine converts primitive variables to conservative
 !
 !===============================================================================
 !
@@ -399,7 +399,7 @@ module scheme
 !
     integer :: i
 !
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
     do i = 1, n
       u(1,i) = q(1,i)
@@ -412,10 +412,10 @@ module scheme
 #endif /* ADI */
     enddo
 
-!----------------------------------------------------------------------
+!-------------------------------------------------------------------------------
 !
   end subroutine prim2cons
 
-!======================================================================
+!===============================================================================
 !
 end module