diff --git a/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/Allclean b/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/Allclean
index 8fe461efbe..29e49eb825 100755
--- a/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/Allclean
+++ b/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/Allclean
@@ -5,6 +5,6 @@ cd "${0%/*}" || exit                                # Run from this directory
 
 cleanCase0
 
-rm -f validation/*.eps
+rm -f validation/*.png
 
 #------------------------------------------------------------------------------
diff --git a/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/validation/plot b/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/validation/plot
index d4cc841205..2cd0026da9 100755
--- a/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/validation/plot
+++ b/tutorials/heatTransfer/buoyantSimpleFoam/buoyantCavity/validation/plot
@@ -31,65 +31,46 @@ cd "${0%/*}" || exit                                # Run from this directory
 
 #------------------------------------------------------------------------------
 
-plot_T()
+plot_x_vs_T_U()
 {
     index="$1"
     sampleFile="$2"
-    benchmarkFile="$3"
+    benchmarkDir="$3"
 
-    image="OF_vs_EXPT_T$index.png"
-
-    gnuplot<<PLT_T
+    gnuplot<<PLT_T_U
     set terminal pngcairo font "helvetica,20" size 1000, 1000
     set grid
     set key left top
     set xrange [0:0.08]
+    set xlabel "Channel width, x [m]"
+
+    # Benchmark - Experimental
+        benchmarkT="$benchmarkDir/mt_z0_${index}0_lo.dat"
+        benchmarkU="$benchmarkDir/mv_z0_${index}0_lo.dat"
+
+    # OpenFOAM
+        samples="$sampleFile"
+
+
+    # plot temperature profiles
     set yrange [285:310]
-    set xlabel "Channel width, x / [m]"
-    set ylabel "Temperature / [K]"
-    set output "$image"
-
-    # Benchmark - Experimental
-        benchmark="$benchmarkFile"
-
-    # OpenFOAM
-        samples="$sampleFile"
+    set ylabel "Temperature [K]"
+    set output "OF_vs_EXPT_T$index.png"
 
     plot \
-        benchmark u (\$1/1000):(\$2+273.15) t "Expt 0.$index" w p lt 1 pt 6, \
-        samples t "OpenFOAM 0.$index" w l lt -1
-PLT_T
-}
+        benchmarkT u (\$1/1000):(\$2+273.15) t "Expt 0.$index" w p lt 1 pt 6, \
+        samples u 1:2 t "OpenFOAM 0.$index" w l lt -1
 
 
-plot_U()
-{
-    index=$1
-    sampleFile="$2"
-    benchmarkFile="$3"
-
-    image="OF_vs_EXPT_U$index.png"
-
-    gnuplot<<PLT_U
-    set terminal pngcairo font "helvetica,20" size 1000, 1000
-    set grid
-    set key left top
-    set xrange [0:0.08]
+    # plot velocity profiles
     set yrange [-0.2:0.2]
-    set xlabel "Channel width, x / [m]"
-    set ylabel "Vertical velocity component, Uy / [m/s]"
-    set output "$image"
-
-    # Benchmark - Experimental
-        benchmark="$benchmarkFile"
-
-    # OpenFOAM
-        samples="$sampleFile"
-
+    set ylabel "Vertical velocity component, Uy [m/s]"
+    set output "OF_vs_EXPT_U$index.png"
     plot \
-        benchmark u (\$1/1000):(\$2) t "Expt 0.$index" w p lt 1 pt 6, \
-        samples u 1:3 title "OpenFOAM 0.$index" w l lt -1
-PLT_U
+        benchmarkU u (\$1/1000):(\$2) t "Expt 0.$index" w p lt 1 pt 6, \
+        samples u 1:4 title "OpenFOAM 0.$index" w l lt -1
+
+PLT_T_U
 }
 
 
@@ -121,34 +102,16 @@ OFDATAROOT="$SETSDIR/$LATESTTIME"
 
 EXPTDATAROOT=./exptData
 
-# generate temperature profiles
-TSets="1 3 4 5 6 7 9"
-for i in $TSets
+# generate temperature and velocity profiles
+set="1 3 4 5 6 7 9"
+for i in $set
 do
 
-    echo "    processing temperature profile at y/yMax of 0.$i"
+    echo "    processing temperature and velocity profiles at y/yMax of 0.$i"
 
-    OF="$OFDATAROOT/y0.${i}_T.xy"
+    OF="$OFDATAROOT/y0.${i}_T_U.xy"
 
-    EXPT="$EXPTDATAROOT/mt_z0_${i}0_lo.dat"
-
-    plot_T "$i" "$OF" "$EXPT"
-
-done
-
-
-# generate velocity profiles
-USets="1 3 4 5 6 7 9"
-for i in $USets
-do
-
-    echo "    processing velocity profile at y/yMax of 0.$i"
-
-    OF="$OFDATAROOT/y0.${i}_U.xy"
-
-    EXPT="$EXPTDATAROOT/mv_z0_${i}0_lo.dat"
-
-    plot_U "$i" "$OF" "$EXPT"
+    plot_x_vs_T_U "$i" "$OF" "$EXPTDATAROOT"
 
 done
 
diff --git a/tutorials/incompressible/boundaryFoam/steadyBoundaryLayer/plot b/tutorials/incompressible/boundaryFoam/steadyBoundaryLayer/plot
index 5421eb9765..ffe9e2b65b 100755
--- a/tutorials/incompressible/boundaryFoam/steadyBoundaryLayer/plot
+++ b/tutorials/incompressible/boundaryFoam/steadyBoundaryLayer/plot
@@ -28,9 +28,8 @@ plot_yPlus_vs_uPlus() {
     set key spacing 0.75
     set xlabel "y^+"
     set ylabel "u^+"
-    set offset .05, .05
     set parametric
-    set trange [0:1.5]
+    set trange [0:30]
     set logscale x
     set format x "10^{%T}"
     set output "$image"
diff --git a/tutorials/incompressible/simpleFoam/T3A/Allclean b/tutorials/incompressible/simpleFoam/T3A/Allclean
index 8fe461efbe..29e49eb825 100755
--- a/tutorials/incompressible/simpleFoam/T3A/Allclean
+++ b/tutorials/incompressible/simpleFoam/T3A/Allclean
@@ -5,6 +5,6 @@ cd "${0%/*}" || exit                                # Run from this directory
 
 cleanCase0
 
-rm -f validation/*.eps
+rm -f validation/*.png
 
 #------------------------------------------------------------------------------
diff --git a/tutorials/incompressible/simpleFoam/T3A/validation/plot b/tutorials/incompressible/simpleFoam/T3A/validation/plot
index 45d6298697..b936020214 100755
--- a/tutorials/incompressible/simpleFoam/T3A/validation/plot
+++ b/tutorials/incompressible/simpleFoam/T3A/validation/plot
@@ -2,6 +2,48 @@
 cd "${0%/*}" || exit                                # Run from this directory
 #------------------------------------------------------------------------------
 
+plot()
+{
+    gnuplot<<PLT
+    set terminal pngcairo font "helvetica,20" size 1000, 1000
+    set encoding utf8
+    set termoption dash
+    set style increment user
+    set style line 1 lt 1 linecolor rgb "blue"  linewidth 1.5
+    set style line 11 lt 2 linecolor rgb "black" linewidth 1.5
+
+    # time = system("foamListTimes -case .. -latestTime")
+    time=$(foamListTimes -case .. -latestTime)
+
+    # benchmark
+    benchmark="exptData/T3A.dat"
+
+    # plot turbulent intensity profile
+    samples="../postProcessing/kGraph/".time."/line_k.xy"
+    set xlabel "x"
+    set ylabel "u'"
+    set output "x_vs_u.png"
+    set title "T3A - Flat plate - Turbulent intensity"
+    plot [:1.5][:0.05] \
+        samples u (\$1-0.04):(1./5.4*sqrt(2./3.*\$2)) t "kOmegaSSTLM" w l ls 1, \
+        benchmark u (\$1/1000):(\$3/100) t "Experiment" w p ls 11
+
+
+    # plot skin-friction coefficient profile
+    samples="../postProcessing/wallShearStressGraph/".time."/line_wallShearStress.xy"
+    set xlabel "Re_x"
+    set ylabel "c_f"
+    set output "Rex_vs_cf.png"
+    set title "T3A - Flat plate - C_f"
+    plot [:6e+5][0:0.01] \
+        samples u ((\$1-0.04)*5.4/1.5e-05):(-\$2/0.5/5.4**2) t "kOmegaSSTLM" w l, \
+        benchmark u (\$1/1000*5.4/1.51e-05):2 t "Experiment" w p ls 11
+PLT
+}
+
+
+#------------------------------------------------------------------------------
+
 # Require gnuplot
 command -v gnuplot >/dev/null || {
     echo "FOAM FATAL ERROR: gnuplot not found - skipping graph creation" 1>&2
@@ -11,32 +53,7 @@ command -v gnuplot >/dev/null || {
 
 #------------------------------------------------------------------------------
 
-gnuplot<<GNUPLOT
-    set term post enhanced color solid linewidth 2.0 20
-    set out "graphs.eps"
-    set encoding utf8
-    set termoption dash
-    set style increment user
-    set style line 1 lt 1 linecolor rgb "blue"  linewidth 1.5
-    set style line 11 lt 2 linecolor rgb "black" linewidth 1.5
+plot
 
-    time = system("foamListTimes -case .. -latestTime")
-
-    set xlabel "x"
-    set ylabel "u'"
-    set title "T3A - Flat Plate - turbulent intensity"
-    plot [:1.5][:0.05] \
-        "../postProcessing/kGraph/".time."/line_k.xy" \
-        u (\$1-0.04):(1./5.4*sqrt(2./3.*\$2))title "kOmegaSSTLM" w l ls 1, \
-        "exptData/T3A.dat" u (\$1/1000):(\$3/100) title "Exp T3A" w p ls 11
-
-    set xlabel "Re_x"
-    set ylabel "c_f"
-    set title "T3A - Flat Plate - C_f"
-    plot [:6e+5][0:0.01] \
-        "../postProcessing/wallShearStressGraph/".time."/line_wallShearStress.xy" \
-        u ((\$1-0.04)*5.4/1.5e-05):(-\$2/0.5/5.4**2) title "kOmegaSSTLM" w l, \
-        "exptData/T3A.dat" u (\$1/1000*5.4/1.51e-05):2 title "Exp" w p ls 11
-GNUPLOT
 
 #------------------------------------------------------------------------------
diff --git a/tutorials/incompressible/simpleFoam/backwardFacingStep2D/plot b/tutorials/incompressible/simpleFoam/backwardFacingStep2D/plot
index 62159831df..06f6c420db 100755
--- a/tutorials/incompressible/simpleFoam/backwardFacingStep2D/plot
+++ b/tutorials/incompressible/simpleFoam/backwardFacingStep2D/plot
@@ -14,16 +14,12 @@ plot_U() {
     gnuplot<<EOF1
     set terminal pngcairo font "helvetica,20" size 1000, 1000
     set grid
-    set key top left
+    set key top right
     set xrange [-0.4:1.2]
     set yrange [0:3]
     set xlabel "U/U_0"
     set ylabel "y/h"
     set output "$image"
-    set lmargin 7.5
-    set rmargin 1.5
-    set tmargin 0.1
-    set bmargin 3.2
     set format x "%.1f"
     set format y "%.1f"
 
@@ -31,12 +27,13 @@ plot_U() {
         samples="$sampleDir"
         Uref="$Uref"
         href="$href"
+        fields="p_U_turbulenceProperties:devReff.xy"
 
     plot \
-        "$sampleDir/x_by_h_01_U.xy" u (\$2/Uref):(\$1/href) w l lw 2 lc rgb "red" t "x/h = 1", \
-        "$sampleDir/x_by_h_04_U.xy" u (\$2/Uref):(\$1/href) w l lw 2 lc rgb "green" t "x/h = 4", \
-        "$sampleDir/x_by_h_06_U.xy" u (\$2/Uref):(\$1/href) w l lw 2 lc rgb "blue" t "x/h = 6", \
-        "$sampleDir/x_by_h_10_U.xy" u (\$2/Uref):(\$1/href) w l lw 2 lc rgb "black" t "x/h = 10"
+        samples."/x_by_h_01_".fields u (\$3/Uref):(\$1/href) w l lw 2 lc rgb "red" t "x/h = 1", \
+        samples."/x_by_h_04_".fields u (\$3/Uref):(\$1/href) w l lw 2 lc rgb "green" t "x/h = 4", \
+        samples."/x_by_h_06_".fields u (\$3/Uref):(\$1/href) w l lw 2 lc rgb "blue" t "x/h = 6", \
+        samples."/x_by_h_10_".fields u (\$3/Uref):(\$1/href) w l lw 2 lc rgb "black" t "x/h = 10"
 EOF1
 }
 
@@ -66,10 +63,6 @@ plot_tau() {
     set xlabel "x/h"
     set ylabel "C_f"
     set output "$image"
-    set lmargin 10
-    set rmargin 1.5
-    set tmargin 0.1
-    set bmargin 3.2
 
     # OpenFOAM
         Uref="$Uref"
@@ -107,7 +100,7 @@ sampleDir=postProcessing/sample/"$timeDir"
 
 #------------------------------------------------------------------------------
 
-Uref=$(awk '{print $2}' $sampleDir/Uref_U.xy)
+Uref=$(awk '{print $2}' $sampleDir/Uref_p_U_turbulenceProperties:devReff.xy)
 
 echo ""
 echo "# Plots the U profiles"
diff --git a/tutorials/incompressible/simpleFoam/rotatingCylinders/Allclean b/tutorials/incompressible/simpleFoam/rotatingCylinders/Allclean
new file mode 100755
index 0000000000..ad9bdeaa8d
--- /dev/null
+++ b/tutorials/incompressible/simpleFoam/rotatingCylinders/Allclean
@@ -0,0 +1,10 @@
+#!/bin/sh
+cd "${0%/*}" || exit                                # Run from this directory
+. ${WM_PROJECT_DIR:?}/bin/tools/CleanFunctions      # Tutorial clean functions
+#------------------------------------------------------------------------------
+
+cleanCase
+
+rm -f *.png
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/incompressible/simpleFoam/rotatingCylinders/plot b/tutorials/incompressible/simpleFoam/rotatingCylinders/plot
index a5dd1278c7..2fb9f14dc9 100755
--- a/tutorials/incompressible/simpleFoam/rotatingCylinders/plot
+++ b/tutorials/incompressible/simpleFoam/rotatingCylinders/plot
@@ -26,17 +26,19 @@ gnuplot<<EOF
     C = 0.5*A*A*R1*R1 + 2.*A*B*log(R1) - 0.5*B*B/R1/R1
     analytical_p(r) = 0.5*A*A*r*r + 2.*A*B*log(r) - 0.5*B*B/r/r - C
 
+    samples="postProcessing/sample1/1000/centreLine_p_U_U:Transformed.xy"
+
     set ylabel "Pressure, p"
     set output "p.png"
     set samples 20
     set key bottom right
-    plot "postProcessing/sample1/1000/centreLine_p.xy" w lines lw 2 lc "black" t "OpenFOAM", \
+    plot samples u 1:2 w lines lw 2 lc "black" t "OpenFOAM", \
         analytical_p(x) w p ps 2 pt 6 lw 2 lc rgb "red" t "analytical"
 
     set ylabel "U_{/Symbol q}, rad/s"
     set output "Utheta.png"
     set samples 20
     set key top right
-    plot "postProcessing/sample1/1000/centreLine_U_U:Transformed.xy" u 1:6 w lines lw 2 lc "black" t "OpenFOAM", \
+    plot samples u 1:7 w lines lw 2 lc "black" t "OpenFOAM", \
         analytical_utheta(x) w p ps 2 pt 6 lw 2 lc rgb "red" t "analytical"
 EOF
diff --git a/tutorials/incompressible/simpleFoam/turbulentFlatPlate/Allrun b/tutorials/incompressible/simpleFoam/turbulentFlatPlate/Allrun
index 2467e39b35..6da543b847 100755
--- a/tutorials/incompressible/simpleFoam/turbulentFlatPlate/Allrun
+++ b/tutorials/incompressible/simpleFoam/turbulentFlatPlate/Allrun
@@ -1,4 +1,4 @@
-#!/bin/sh
+#!/bin/bash
 cd "${0%/*}" || exit                                # Run from this directory
 . ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions        # Tutorial run functions
 #------------------------------------------------------------------------------
diff --git a/tutorials/verificationAndValidation/atmosphericModels/atmDownstreamDevelopment/plot b/tutorials/verificationAndValidation/atmosphericModels/atmDownstreamDevelopment/plot
index b27704976e..cf95a9505f 100755
--- a/tutorials/verificationAndValidation/atmosphericModels/atmDownstreamDevelopment/plot
+++ b/tutorials/verificationAndValidation/atmosphericModels/atmDownstreamDevelopment/plot
@@ -375,7 +375,7 @@ do
 
     endTime=$( \
         foamDictionary results/$setup/system/controlDict \
-            -disableFunctionEntries -entry endTime -value \
+            -entry endTime -value \
         )