diff --git a/src/simsoptpp/surfacerzfourier.cpp b/src/simsoptpp/surfacerzfourier.cpp
index e0ca3ac41..98e158923 100644
--- a/src/simsoptpp/surfacerzfourier.cpp
+++ b/src/simsoptpp/surfacerzfourier.cpp
@@ -535,7 +535,7 @@ Array SurfaceRZFourier<Array>::dgamma_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];
@@ -558,7 +558,7 @@ Array SurfaceRZFourier<Array>::dgamma_by_dcoeff_vjp(Array& v) {
         for (int i = 0; i < num_dofs(); ++i) {
             resptr_private[i] = 0.;
         }
-#pragma omp for ordered
+#pragma omp for
         for (int k1 = 0; k1 < numquadpoints_phi; ++k1) {
             double phi  = 2*M_PI*quadpoints_phi[k1];
             double sinphi = sin(phi);
@@ -613,7 +613,7 @@ Array SurfaceRZFourier<Array>::dgamma_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];
@@ -752,7 +752,7 @@ Array SurfaceRZFourier<Array>::dgammadash1_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];
@@ -834,7 +834,7 @@ Array SurfaceRZFourier<Array>::dgammadash1_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];
@@ -1084,7 +1084,7 @@ Array SurfaceRZFourier<Array>::dgammadash2_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];
@@ -1166,7 +1166,7 @@ Array SurfaceRZFourier<Array>::dgammadash2_by_dcoeff_vjp(Array& v) {
 #pragma omp for ordered
     for( int t=0; t<omp_get_num_threads(); ++t )
             {
-                #pragma omp ordered critical
+                #pragma omp ordered
                 {
                     for(int i=0; i<num_dofs(); ++i) {
                         resptr[i] += resptr_private[i];