diff --git a/src/radial_derivatives_even.f90 b/src/radial_derivatives_even.f90
index fc9e3780..2ba9be0a 100644
--- a/src/radial_derivatives_even.f90
+++ b/src/radial_derivatives_even.f90
@@ -17,10 +17,10 @@ subroutine get_ddr_even(f,df,ddf,n_f_max,n_f_start,n_f_stop,   &
               &            n_r_max,n_cheb_max,dr_fac,work1,work2, &
               &            chebt_odd, chebt_even)
       !
-      !  Returns first rarial derivative df and second radial
+      !  Returns first radial derivative df and second radial
       !  derivative ddf of the input function f.
       !  Array f(n_f_max,*) may contain several functions numbered by
-      !  the first index. The subroutine calculates the derivaties of
+      !  the first index. The subroutine calculates the derivatives of
       !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
       !  to a Chebychev representation using n_r_max radial grid points.
       !  The cheb transforms have to be initialized by calling
@@ -73,10 +73,9 @@ subroutine get_drNS_even(f,df,n_f_max,n_f_start,n_f_stop, &
               &             n_r_max,n_cheb_max,dr_fac,work1, &
               &             chebt_odd, chebt_even)
       !
-      !  Returns first rarial derivative df and second radial
-      !  derivative ddf of the input function f.
+      !  Returns first rarial derivative df.
       !  Array f(n_f_max,*) may contain several functions numbered by
-      !  the first index. The subroutine calculates the derivaties of
+      !  the first index. The subroutine calculates the derivative of
       !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
       !  to a Chebychev representation using n_r_max radial grid points.
       !  The cheb transforms have to be initialized by calling
@@ -116,10 +115,10 @@ subroutine get_ddrNS_even(f,df,ddf,n_f_max,n_f_start,n_f_stop, &
               &              n_r_max,n_cheb_max,dr_fac,work1,     &
               &              chebt_odd, chebt_even)
       !
-      !  Returns first rarial derivative df and second radial
+      !  Returns first radial derivative df and second radial
       !  derivative ddf of the input function f.
       !  Array f(n_f_max,*) may contain several functions numbered by
-      !  the first index. The subroutine calculates the derivaties of
+      !  the first index. The subroutine calculates the derivatives of
       !  the functions f(n_f_start,*) to f(n_f_stop) by transforming
       !  to a Chebychev representation using n_r_max radial grid points.
       !  The cheb transforms have to be initialized by calling