Add -Wall compiler optnio

Author: Ben Wells

CMakeLists.txt

@@ -22,8 +22,8 @@ get_filename_component (Fortran_COMPILER_NAME ${CMAKE_Fortran_COMPILER} NAME)
 
 if (Fortran_COMPILER_NAME MATCHES "gfortran.*")
   # gfortran
-  set (CMAKE_Fortran_FLAGS_RELEASE "-funroll-all-loops -fno-f2c -O3")
-  set (CMAKE_Fortran_FLAGS_DEBUG   "-fno-f2c -O0 -g")
+  set (CMAKE_Fortran_FLAGS_RELEASE "-funroll-all-loops -fno-f2c -O3 -Wall")
+  set (CMAKE_Fortran_FLAGS_DEBUG   "-fno-f2c -O0 -g -Wall")
 elseif (Fortran_COMPILER_NAME MATCHES "ifort.*")
   # ifort (untested)
   set (CMAKE_Fortran_FLAGS_RELEASE "-f77rtl -O3")

src/main.f90

@@ -23,7 +23,7 @@ program pme
    !*********************************************************************************
    implicit none
    !---------------------------------------------------------------------------------
-   integer :: nodes, i, reports, reportid
+   integer :: nodes, reports, reportid
    double precision, dimension(:), allocatable :: u, x, theta, x_dot
    double precision :: output_t, delta_t, total_t, t_init, rzero, tzero
    double precision :: report_step, report_time
@@ -177,8 +177,8 @@ subroutine initial_conditions(u, x, nodes, m, Q, t_init, rzero, tzero)
    one_over_m = 1.0d0/m
 
    ! Calculate quantities for self-similar solution.
-   gamman = gamma(one_over_m + (3.0d0/2.0d0))
-   gammad = gamma(one_over_m + 1.0d0)
+   gamman = GammaFunction(one_over_m + (3.0d0/2.0d0))
+   gammad = GammaFunction(one_over_m + 1.0d0)
 
    rzero = Q*gamman/(SQRT(pi)*gammad)
    tzero = ((rzero**2)*m)/(2.0d0*(m + 2.0d0))

src/special_functions.f90

@@ -2,7 +2,7 @@ module special_functions
 
 contains
 
-   double precision function gamma(xx)
+   double precision function GammaFunction(xx)
       !*********************************************************************************
       !**                                                                             **
       !**  This function calculates the gamma funtion for the value xx.               **
@@ -36,11 +36,11 @@ double precision function gamma(xx)
          ser = ser + cof(j)/y
       end do
 
-      gamma = tmp + dlog(stp*ser/x)
-      gamma = dexp(gamma)
+      GammaFunction = tmp + dlog(stp*ser/x)
+      GammaFunction = dexp(GammaFunction)
 
       return
 
-   end function gamma
+   end function GammaFunction
 
 end module special_functions