diff --git a/CMakeLists.txt b/CMakeLists.txt
index 82025a83eed4fc2a43b9de3a1b918eea9902a461..33a4ba7ac41d7d55a53a805c777a15d777a22c33 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -373,6 +373,7 @@ file(COPY test/decoupled/2p2c/test_dec2p2c.input DESTINATION parameters)
 file(COPY test/decoupled/2p2c/test_multiphysics2p2c-reference.vtu DESTINATION references)
 file(COPY test/common/generalproblem/generallens_box-reference.vtu DESTINATION references)
 file(COPY test/common/generalproblem/generallens_decoupled-reference.vtu DESTINATION references)
+file(COPY tutorial/tutorial_coupled.input DESTINATION parameters)
 if(SUPERLU_FOUND)
   file(COPY test/freeflow/stokes/stokes-reference.vtu DESTINATION references)
   file(COPY test/freeflow/stokes/test_stokes.input DESTINATION parameters)
@@ -413,7 +414,7 @@ add_test(test_impes_adaptive     bin/runTest.sh references/output2padaptive-refe
 add_test(test_dec2p2c   bin/runTest.sh references/test_dec2p2c-reference.vtu   test_dec2p2c-00021.vtu     test/decoupled/2p2c/test_dec2p2c --parameter-file=parameters/test_dec2p2c.input)
 add_test(test_multiphysics2p2c   bin/runTest.sh references/test_multiphysics2p2c-reference.vtu   test_multiphysics2p2c-00021.vtu     test/decoupled/2p2c/test_multiphysics2p2c)
 
-add_test(tutorial_coupled   tutorial/tutorial_coupled   1 1)
+add_test(tutorial_coupled   tutorial/tutorial_coupled --parameter-file=parameters/tutorial_coupled.input)
 add_test(tutorial_decoupled tutorial/tutorial_decoupled 1)
 
 if(SUPERLU_FOUND)
diff --git a/test/boxmodels/1p2c/outflow-reference.vtu b/test/boxmodels/1p2c/outflow-reference.vtu
index 686831dde646503c5cfcda0b71d1f91c91743a65..4298be9b9a2cb2a297149a8e38e8b4ac632e5ca0 100644
--- a/test/boxmodels/1p2c/outflow-reference.vtu
+++ b/test/boxmodels/1p2c/outflow-reference.vtu
@@ -28,44 +28,44 @@
           0.00997054 0.00997056 0.00997057
         </DataArray>
         <DataArray type="Float32" Name="Vy" NumberOfComponents="1" format="ascii">
-          2.9018e-18 -7.25448e-19 -7.2545e-19 -1.08817e-18 6.52903e-18 1.08817e-18 2.90178e-18 3.26451e-18 -5.07812e-18 2.53906e-18 0 3.62722e-19
-          -7.25443e-18 1.45089e-18 -3.62721e-18 2.53905e-18 -7.97985e-18 -1.08816e-18 1.45088e-18 4.35264e-18 7.25438e-19 7.25438e-19 -1.45087e-18 -1.45087e-18
-          0 -3.62718e-19 7.25435e-19 -2.1763e-18 3.62717e-19 1.81358e-19 0 -1.63222e-18 7.25432e-19 -9.06789e-19 3.62715e-19 -1.2695e-18
-          -1.81357e-18 -3.26443e-18 0 0 0 0 -4.3527e-18 -1.4509e-18 -4.35268e-18 3.62723e-18 1.01562e-17 7.25444e-19
-          1.01562e-17 8.7053e-18 5.80352e-18 7.25439e-18 7.25438e-19 -1.45087e-18 -7.25436e-19 -5.07804e-18 0 -3.26445e-18 -2.53901e-18 -2.90172e-18
-          -4.71529e-18 0 0
+          -2.9018e-18 -1.4509e-18 -3.62725e-18 -1.4509e-18 -5.07813e-18 -3.62724e-18 -5.80357e-18 -1.45089e-18 -7.97989e-18 -3.62723e-18 7.25444e-19 3.62722e-19
+          5.0781e-18 -7.25443e-19 0 0 -7.2544e-19 -3.6272e-19 -1.45088e-18 2.90176e-18 -2.17631e-18 -7.25438e-19 -1.45087e-18 -1.45087e-18
+          -7.25436e-19 -7.25436e-19 -5.07804e-18 -2.1763e-18 0 1.81358e-19 -3.26445e-18 -7.25433e-19 -2.17629e-18 1.45086e-18 3.62715e-19 3.62715e-19
+          -1.81357e-18 -1.81357e-18 0 0 0 0 -4.3527e-18 -1.4509e-18 -2.17634e-18 2.90178e-18 7.25445e-19 0
+          -6.52898e-18 0 0 7.25439e-18 7.25438e-19 -1.45087e-18 -7.25436e-19 7.25435e-19 3.62717e-19 1.81358e-18 5.07802e-18 3.62715e-19
+          -1.81357e-18 0 0
         </DataArray>
         <DataArray type="Float32" Name="x_H2O" NumberOfComponents="1" format="ascii">
-          0.783938 0.783939 0.783938 0.783939 0.78394 0.78394 0.783946 0.783946 0.783961 0.783961 0.783994 0.783994
-          0.784056 0.784056 0.784162 0.784162 0.78433 0.78433 0.784579 0.784579 0.78493 0.78493 0.785405 0.785405
-          0.786025 0.786025 0.786808 0.786808 0.787773 0.787773 0.788935 0.788935 0.790306 0.790306 0.791894 0.791894
-          0.793707 0.793707 0.795747 0.795747 0.796822 0.796822 0.783938 0.783939 0.78394 0.783946 0.783961 0.783994
-          0.784056 0.784162 0.78433 0.784579 0.78493 0.785405 0.786025 0.786808 0.787773 0.788935 0.790306 0.791894
-          0.793707 0.795747 0.796822
+          0.783938 0.783938 0.783938 0.783938 0.78394 0.78394 0.783945 0.783945 0.783958 0.783958 0.783988 0.783988
+          0.784043 0.784043 0.784139 0.784139 0.784292 0.784292 0.784519 0.784519 0.784841 0.784841 0.78528 0.78528
+          0.785854 0.785854 0.786584 0.786584 0.787488 0.787488 0.788581 0.788581 0.789877 0.789877 0.791386 0.791386
+          0.793114 0.793114 0.795067 0.795067 0.796098 0.796098 0.783938 0.783938 0.78394 0.783945 0.783958 0.783988
+          0.784043 0.784139 0.784292 0.784519 0.784841 0.78528 0.785854 0.786584 0.787488 0.788581 0.789877 0.791386
+          0.793114 0.795067 0.796098
         </DataArray>
         <DataArray type="Float32" Name="x_N2" NumberOfComponents="1" format="ascii">
-          0.216062 0.216061 0.216062 0.216061 0.21606 0.21606 0.216054 0.216054 0.216039 0.216039 0.216006 0.216006
-          0.215944 0.215944 0.215838 0.215838 0.21567 0.21567 0.215421 0.215421 0.21507 0.21507 0.214595 0.214595
-          0.213975 0.213975 0.213192 0.213192 0.212227 0.212227 0.211065 0.211065 0.209694 0.209694 0.208106 0.208106
-          0.206293 0.206293 0.204253 0.204253 0.203178 0.203178 0.216062 0.216061 0.21606 0.216054 0.216039 0.216006
-          0.215944 0.215838 0.21567 0.215421 0.21507 0.214595 0.213975 0.213192 0.212227 0.211065 0.209694 0.208106
-          0.206293 0.204253 0.203178
+          0.783938 0.783938 0.783938 0.783938 0.78394 0.78394 0.783945 0.783945 0.783958 0.783958 0.783988 0.783988
+          0.784043 0.784043 0.784139 0.784139 0.784292 0.784292 0.784519 0.784519 0.784841 0.784841 0.78528 0.78528
+          0.785854 0.785854 0.786584 0.786584 0.787488 0.787488 0.788581 0.788581 0.789877 0.789877 0.791386 0.791386
+          0.793114 0.793114 0.795067 0.795067 0.796098 0.796098 0.783938 0.783938 0.78394 0.783945 0.783958 0.783988
+          0.784043 0.784139 0.784292 0.784519 0.784841 0.78528 0.785854 0.786584 0.787488 0.788581 0.789877 0.791386
+          0.793114 0.795067 0.796098
         </DataArray>
         <DataArray type="Float32" Name="X_H2O" NumberOfComponents="1" format="ascii">
-          0.7 0.7 0.7 0.7 0.700002 0.700002 0.700009 0.700009 0.700028 0.700028 0.700069 0.700069
-          0.700146 0.700146 0.700278 0.700278 0.700486 0.700486 0.700795 0.700795 0.701231 0.701231 0.70182 0.70182
-          0.702589 0.702589 0.703563 0.703563 0.704763 0.704763 0.70621 0.70621 0.707919 0.707919 0.709903 0.709903
-          0.712171 0.712171 0.714727 0.714727 0.716076 0.716076 0.7 0.7 0.700002 0.700009 0.700028 0.700069
-          0.700146 0.700278 0.700486 0.700795 0.701231 0.70182 0.702589 0.703563 0.704763 0.70621 0.707919 0.709903
-          0.712171 0.714727 0.716076
+          0.7 0.7 0.7 0.7 0.700002 0.700002 0.700008 0.700008 0.700025 0.700025 0.700061 0.700061
+          0.70013 0.70013 0.700249 0.700249 0.700438 0.700438 0.70072 0.70072 0.70112 0.70112 0.701664 0.701664
+          0.702377 0.702377 0.703285 0.703285 0.704409 0.704409 0.70577 0.70577 0.707385 0.707385 0.709268 0.709268
+          0.711429 0.711429 0.713874 0.713874 0.715167 0.715167 0.7 0.7 0.700002 0.700008 0.700025 0.700061
+          0.70013 0.700249 0.700438 0.70072 0.70112 0.701664 0.702377 0.703285 0.704409 0.70577 0.707385 0.709268
+          0.711429 0.713874 0.715167
         </DataArray>
         <DataArray type="Float32" Name="X_N2" NumberOfComponents="1" format="ascii">
-          0.3 0.3 0.3 0.3 0.299998 0.299998 0.299991 0.299991 0.299972 0.299972 0.299931 0.299931
-          0.299854 0.299854 0.299722 0.299722 0.299514 0.299514 0.299205 0.299205 0.298769 0.298769 0.29818 0.29818
-          0.297411 0.297411 0.296437 0.296437 0.295237 0.295237 0.29379 0.29379 0.292081 0.292081 0.290097 0.290097
-          0.287829 0.287829 0.285273 0.285273 0.283924 0.283924 0.3 0.3 0.299998 0.299991 0.299972 0.299931
-          0.299854 0.299722 0.299514 0.299205 0.298769 0.29818 0.297411 0.296437 0.295237 0.29379 0.292081 0.290097
-          0.287829 0.285273 0.283924
+          0.3 0.3 0.3 0.3 0.299998 0.299998 0.299992 0.299992 0.299975 0.299975 0.299939 0.299939
+          0.29987 0.29987 0.299751 0.299751 0.299562 0.299562 0.29928 0.29928 0.29888 0.29888 0.298336 0.298336
+          0.297623 0.297623 0.296715 0.296715 0.295591 0.295591 0.29423 0.29423 0.292615 0.292615 0.290732 0.290732
+          0.288571 0.288571 0.286126 0.286126 0.284833 0.284833 0.3 0.3 0.299998 0.299992 0.299975 0.299939
+          0.29987 0.299751 0.299562 0.29928 0.29888 0.298336 0.297623 0.296715 0.295591 0.29423 0.292615 0.290732
+          0.288571 0.286126 0.284833
         </DataArray>
         <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
           998.237 998.234 998.237 998.234 998.232 998.232 998.23 998.23 998.227 998.227 998.225 998.225