Add constants.f90

.gitignore

@@ -1,2 +1,5 @@
 *.o
 *.mod
+*.f90~
+#Makefile#
+test

caotic_pendulum.f90

@@ -1,12 +1,10 @@
 module caotic_pendulum
   use precision
+  use constants
   use functions
   use solvers
   implicit none
 
-  real(dp), parameter :: pi = 3.14159265358979323846264338327_dp 
-  real(dp), parameter :: eps = 1.d-10
-
 contains
 
   subroutine solve(Nstep, Nperiods, x0, v0)

constants.f90

@@ -0,0 +1,7 @@
+module constants
+  use precision
+
+  real(dp), parameter :: pi = acos(-1.0_dp) 
+  real(dp), parameter :: eps = 1.d-10
+
+end module constants

make.dep

@@ -1,6 +1,8 @@
 caotic_pendulum.o: precision.o
+caotic_pendulum.o: constants.o
 caotic_pendulum.o: functions.o
 caotic_pendulum.o: solvers.o
+constants.o: precision.o
 functions.o: precision.o
 functions.o: precision.o
 solvers.o: precision.o

test.f90

@@ -1,10 +1,10 @@
 program test
     use precision
+    use constants
     use functions
     use solvers
     implicit none
 
-    real(dp), parameter        :: pi = acos(-1.0_dp)
     real(dp)                   :: err, t, dt, x0, v0, Qmin, Qmax, dQ, AA, BB, xex
     real(dp), allocatable      :: u(:), u0(:), param(:)
     integer                    :: i, j, k, NQ, Nstep, Ncicli, num_args, funit            ! funit = file unit