diff --git a/bfps/base.py b/bfps/base.py
index 52451d94797e56bc2015a402e7b2ae13cfc63d36..559089d1d614c746132ed85884168900237085fc 100644
--- a/bfps/base.py
+++ b/bfps/base.py
@@ -71,7 +71,7 @@ class base(object):
             elif type(self.parameters[key[i]]) == str:
                 src_txt += 'MPI_Bcast((void*)(' + key[i] + '), {0}, MPI_CHAR, {1}, MPI_COMM_WORLD);\n'.format(self.string_length, self.iorank)
             else:
-                src_txt += 'MPI_Bcast((void*)(&' + key[i] + '), 1, MPI_REAL8, {0}, MPI_COMM_WORLD);\n'.format(self.iorank)
+                src_txt += 'MPI_Bcast((void*)(&' + key[i] + '), 1, MPI_DOUBLE, {0}, MPI_COMM_WORLD);\n'.format(self.iorank)
         src_txt += ('MPI_Allreduce((void*)(&err_while_reading), (void*)(&errr), 1, MPI_INTEGER, MPI_SUM, MPI_COMM_WORLD);\n'
                   + 'if (errr > 0)\n{\n'
                   + 'fprintf(stderr, "Error reading parameters.\\nAttempting to exit.\\n");\n'
diff --git a/src/fluid_solver.cpp b/src/fluid_solver.cpp
index 26ec1305d91b52506b81c9593dc05e03fec934c7..b2082d30f266845d6e75bd30a8dc8195a656f53a 100644
--- a/src/fluid_solver.cpp
+++ b/src/fluid_solver.cpp
@@ -420,14 +420,14 @@ FLUID_SOLVER_DEFINITIONS(
         FFTW_MANGLE_FLOAT,
         float,
         fftwf_complex,
-        MPI_REAL4,
-        MPI_COMPLEX8)
+        MPI_FLOAT,
+        MPI_C_FLOAT_COMPLEX)
 //FLUID_SOLVER_DEFINITIONS(
 //        FFTW_MANGLE_DOUBLE,
 //        double,
 //        fftw_complex,
-//        MPI_REAL8,
-//        MPI_COMPLEX16)
+//        MPI_DOUBLE,
+//        MPI_C_DOUBLE_COMPLEX)
 /*****************************************************************************/
 
 
diff --git a/src/fluid_solver_base.cpp b/src/fluid_solver_base.cpp
index 931fd69a0eed8f9d8001fd28abc2086b6dad70e4..86d566ca893742f01e017ec8a642d67bd79ab1ee 100644
--- a/src/fluid_solver_base.cpp
+++ b/src/fluid_solver_base.cpp
@@ -59,7 +59,7 @@ void fluid_solver_base<rnumber>::cospectrum(cnumber *a, cnumber *b, double *spec
             (void*)cospec_local,
             (void*)spec,
             this->nshells,
-            MPI_REAL8, MPI_SUM, this->cd->comm);
+            MPI_DOUBLE, MPI_SUM, this->cd->comm);
     for (int n=0; n<this->nshells; n++)
     {
         spec[n] *= 12.5663706144*pow(this->kshell[n], 2) / this->nshell[n];
@@ -211,7 +211,7 @@ fluid_solver_base<R>::fluid_solver_base( \
             (void*)(kshell_local), \
             (void*)(this->kshell), \
             this->nshells, \
-            MPI_REAL8, MPI_SUM, this->cd->comm); \
+            MPI_DOUBLE, MPI_SUM, this->cd->comm); \
     for (int n=0; n<this->nshells; n++) \
     { \
         this->kshell[n] /= this->nshell[n]; \
@@ -450,14 +450,14 @@ FLUID_SOLVER_BASE_DEFINITIONS(
         FFTW_MANGLE_FLOAT,
         float,
         fftwf_complex,
-        MPI_REAL4,
-        MPI_COMPLEX8)
+        MPI_FLOAT,
+        MPI_C_FLOAT_COMPLEX)
 //FLUID_SOLVER_BASE_DEFINITIONS(
 //        FFTW_MANGLE_DOUBLE,
 //        double,
 //        fftw_complex,
-//        MPI_REAL8,
-//        MPI_COMPLEX16)
+//        MPI_DOUBLE,
+//        MPI_C_DOUBLE_COMPLEX)
 /*****************************************************************************/
 
 
diff --git a/src/slab_field_particles.cpp b/src/slab_field_particles.cpp
index e21d8fb6687db50de84ed3efffe987df9e290575..0fd1bdfaa94b79968a88ac232eeb009ae61d20b8 100644
--- a/src/slab_field_particles.cpp
+++ b/src/slab_field_particles.cpp
@@ -77,7 +77,7 @@ slab_field_particles<rnumber>::slab_field_particles(
             tbound,
             this->lbound,
             nprocs,
-            MPI_REAL8,
+            MPI_DOUBLE,
             MPI_SUM,
             this->fs->rd->comm);
     std::fill_n(tbound, nprocs, 0.0);
@@ -86,7 +86,7 @@ slab_field_particles<rnumber>::slab_field_particles(
             tbound,
             this->ubound,
             nprocs,
-            MPI_REAL8,
+            MPI_DOUBLE,
             MPI_SUM,
             this->fs->rd->comm);
     delete[] tbound;
@@ -140,7 +140,7 @@ void slab_field_particles<rnumber>::synchronize_single_particle(int p)
                 MPI_Send(
                         this->state + p*this->ncomponents,
                         this->ncomponents,
-                        MPI_REAL8,
+                        MPI_DOUBLE,
                         r,
                         p*this->computing[p],
                         this->fs->rd->comm);
@@ -148,7 +148,7 @@ void slab_field_particles<rnumber>::synchronize_single_particle(int p)
                 MPI_Recv(
                         this->state + p*this->ncomponents,
                         this->ncomponents,
-                        MPI_REAL8,
+                        MPI_DOUBLE,
                         this->computing[p],
                         p*this->computing[p],
                         this->fs->rd->comm,
@@ -176,7 +176,7 @@ void slab_field_particles<rnumber>::synchronize()
             tstate,
             this->state,
             this->array_size,
-            MPI_REAL8,
+            MPI_DOUBLE,
             MPI_SUM,
             this->fs->rd->comm);
     fftw_free(tstate);
@@ -291,7 +291,7 @@ void slab_field_particles<rnumber>::read()
     MPI_Bcast(
             this->state,
             this->array_size,
-            MPI_REAL8,
+            MPI_DOUBLE,
             0,
             this->fs->rd->comm);
     // initial assignment of particles
@@ -325,7 +325,7 @@ ptrdiff_t slab_field_particles<rnumber>::buffered_local_size()
 template <class rnumber>
 void slab_field_particles<rnumber>::rFFTW_to_buffered(rnumber *src, rnumber *dst)
 {
-    const MPI_Datatype MPI_RNUM = (sizeof(rnumber) == 4) ? MPI_REAL4 : MPI_REAL8;
+    const MPI_Datatype MPI_RNUM = (sizeof(rnumber) == 4) ? MPI_FLOAT : MPI_DOUBLE;
     const ptrdiff_t bsize = this->buffer_size*this->fs->rd->slice_size;
     MPI_Request *mpirequest = new MPI_Request;
     /* do big copy of middle stuff */
diff --git a/src/tracers.cpp b/src/tracers.cpp
index bb82010a97dd0c2ae66ff21db4f8ff8b4b8e5007..e1496c85a42ade745732fcea7a1993805f4de9ce 100644
--- a/src/tracers.cpp
+++ b/src/tracers.cpp
@@ -54,7 +54,7 @@ void tracers<rnumber>::jump_estimate(double *jump)
             tjump,
             jump,
             this->nparticles,
-            MPI_REAL8,
+            MPI_DOUBLE,
             MPI_SUM,
             this->fs->rd->comm);
     delete[] tjump;
@@ -125,8 +125,8 @@ TRACERS_DEFINITIONS(
         FFTW_MANGLE_FLOAT,
         float,
         fftwf_complex,
-        MPI_REAL4,
-        MPI_COMPLEX8)
+        MPI_FLOAT,
+        MPI_C_FLOAT_COMPLEX)
 /*****************************************************************************/