diff --git a/src/driver.F90 b/src/driver.F90
index 97edc5e..3a8f13e 100644
--- a/src/driver.F90
+++ b/src/driver.F90
@@ -36,6 +36,7 @@ program amun
   use config   , only : fdt
 #endif /* FORCE */
   use coordinates, only : initialize_coordinates, finalize_coordinates
+  use equations     , only : initialize_equations
   use evolution, only : evolve, find_new_timestep, n, t, dt, dtn
 #ifdef FORCE
   use forcing  , only : init_forcing, clear_forcing
@@ -221,6 +222,10 @@ program amun
 !
     call initialize_mesh(.true.)
 
+! initialize module EQUATIONS
+!
+    call initialize_equations()
+
 ! initialize the integrals module
 !
     call init_integrals(.true.)
diff --git a/src/makefile b/src/makefile
index 22c54bd..121d1fc 100644
--- a/src/makefile
+++ b/src/makefile
@@ -208,9 +208,9 @@ boundaries.o    : boundaries.F90 blocks.o config.o error.o interpolation.o     \
 config.o        : config.F90 error.o
 constants.o     : constants.F90
 coordinates.o   : coordinates.F90 parameters.o
-driver.o        : driver.F90 blocks.o config.o coordinates.o evolution.o       \
-                  forcing.o integrals.o io.o mesh.o mpitools.o parameters.o    \
-                  random.o
+driver.o        : driver.F90 blocks.o config.o coordinates.o equations.o       \
+                  evolution.o forcing.o integrals.o io.o mesh.o mpitools.o     \
+                  parameters.o random.o
 equations.o     : equations.F90 parameters.o variables.o
 error.o         : error.F90
 evolution.o     : evolution.F90 blocks.o boundaries.o config.o coordinates.o   \