# script for plotting raw data experiment I 
# the whole plots will be in a unique file :psuml1.ps
# 
# then each plot is saved individually in encapsulated Postscript format  
#    names: ml1DDDDm.eps   for mean velocity
#           ml1DDDDf.eps   for fluctuating part
set term postscript
set output "../../ps/psuml1.ps"
set data style linespoint
set yrange [20:45]
set xrange [-60:60]
set title "PLANE MIXING LAYER: Mean Velocity Profile"
set key 50,25
set ylabel "<U> [m/s]"
set xlabel "y [mm]"
set grid
plot "../../data/0030.ml1" using 2:3 title "X=30 mm"
plot "../../data/0040.ml1" using 2:3 title "X=40 mm"
plot "../../data/0050.ml1" using 2:3 title "X=50 mm"
plot "../../data/0070.ml1" using 2:3 title "X=70 mm"
plot "../../data/0100.ml1" using 2:3 title "X=100 mm"
plot "../../data/0200.ml1" using 2:3 title "X=200 mm"
plot "../../data/0250.ml1" using 2:3 title "X=250 mm"
plot "../../data/0300.ml1" using 2:3 title "X=300 mm"
plot "../../data/0350.ml1" using 2:3 title "X=350 mm"
plot "../../data/0400.ml1" using 2:3 title "X=400 mm"
plot "../../data/0450.ml1" using 2:3 title "X=450 mm"
plot "../../data/0500.ml1" using 2:3 title "X=500 mm"
plot "../../data/0550.ml1" using 2:3 title "X=550 mm"
plot "../../data/0600.ml1" using 2:3 title "X=600 mm"
plot "../../data/0650.ml1" using 2:3 title "X=650 mm"
plot "../../data/0725.ml1" using 2:3 title "X=725 mm"
plot "../../data/0750.ml1" using 2:3 title "X=750 mm"
plot "../../data/0800.ml1" using 2:3 title "X=800 mm"
plot "../../data/0850.ml1" using 2:3 title "X=850 mm"
plot "../../data/0900.ml1" using 2:3 title "X=900 mm"
plot "../../data/0950.ml1" using 2:3 title "X=950 mm"
plot "../../data/1000.ml1" using 2:3 title "X=1000 mm"
plot "../../data/1050.ml1" using 2:3 title "X=1050 mm"

set title "PLANE MIXING LAYER: <u'u'> Velocity Profile"
set key
set ylabel "<u'u'> [(m/s)**2]"
set xlabel "y [mm]"
set yrange [-1:12]
plot "../../data/0030.ml1" using 2:4 title "X=30 mm"
plot "../../data/0040.ml1" using 2:4 title "X=40 mm"
plot "../../data/0050.ml1" using 2:4 title "X=50 mm"
plot "../../data/0070.ml1" using 2:4 title "X=70 mm"
plot "../../data/0100.ml1" using 2:4 title "X=100 mm"
plot "../../data/0200.ml1" using 2:4 title "X=200 mm"
plot "../../data/0250.ml1" using 2:4 title "X=250 mm"
plot "../../data/0300.ml1" using 2:4 title "X=300 mm"
plot "../../data/0350.ml1" using 2:4 title "X=350 mm"
plot "../../data/0400.ml1" using 2:4 title "X=400 mm"
plot "../../data/0450.ml1" using 2:4 title "X=450 mm"
plot "../../data/0500.ml1" using 2:4 title "X=500 mm"
plot "../../data/0550.ml1" using 2:4 title "X=550 mm"
plot "../../data/0600.ml1" using 2:4 title "X=600 mm"
plot "../../data/0650.ml1" using 2:4 title "X=650 mm"
plot "../../data/0725.ml1" using 2:4 title "X=725 mm"
plot "../../data/0750.ml1" using 2:4 title "X=750 mm"
plot "../../data/0800.ml1" using 2:4 title "X=800 mm"
plot "../../data/0850.ml1" using 2:4 title "X=850 mm"
plot "../../data/0900.ml1" using 2:4 title "X=900 mm"
plot "../../data/0950.ml1" using 2:4 title "X=950 mm"
plot "../../data/1000.ml1" using 2:4 title "X=1000 mm"
plot "../../data/1050.ml1" using 2:4 title "X=1050 mm"


set term postscript eps
set data style linespoint
set yrange [20:45]
set xrange [-60:60]
set title "PLANE MIXING LAYER: Mean Velocity Profile"
set key 50,25
set ylabel "<U> [m/s]"
set xlabel "y [mm]"
set grid
set output "../../ps/ml10030m.eps"
plot "../../data/0030.ml1" using 2:3 title "X=30 mm"

set output "../../ps/ml10040m.eps"
plot "../../data/0040.ml1" using 2:3 title "X=40 mm"

set output "../../ps/ml10050m.eps"
plot "../../data/0050.ml1" using 2:3 title "X=50 mm"

set output "../../ps/ml10070m.eps"
plot "../../data/0070.ml1" using 2:3 title "X=70 mm"

set output "../../ps/ml10100m.eps"
plot "../../data/0100.ml1" using 2:3 title "X=100 mm"

set output "../../ps/ml10200m.eps"
plot "../../data/0200.ml1" using 2:3 title "X=200 mm"

set output "../../ps/ml10250m.eps"
plot "../../data/0250.ml1" using 2:3 title "X=250 mm"

set output "../../ps/ml10300m.eps"
plot "../../data/0300.ml1" using 2:3 title "X=300 mm"

set output "../../ps/ml10350m.eps"
plot "../../data/0350.ml1" using 2:3 title "X=350 mm"

set output "../../ps/ml10400m.eps"
plot "../../data/0400.ml1" using 2:3 title "X=400 mm"

set output "../../ps/ml10450m.eps"
plot "../../data/0450.ml1" using 2:3 title "X=450 mm"

set output "../../ps/ml10500m.eps"
plot "../../data/0500.ml1" using 2:3 title "X=500 mm"

set output "../../ps/ml10550m.eps"
plot "../../data/0550.ml1" using 2:3 title "X=550 mm"

set output "../../ps/ml10600m.eps"
plot "../../data/0600.ml1" using 2:3 title "X=600 mm"

set output "../../ps/ml10650m.eps"
plot "../../data/0650.ml1" using 2:3 title "X=650 mm"

set output "../../ps/ml10725m.eps"
plot "../../data/0725.ml1" using 2:3 title "X=725 mm"

set output "../../ps/ml10750m.eps"
plot "../../data/0750.ml1" using 2:3 title "X=750 mm"

set output "../../ps/ml10800m.eps"
plot "../../data/0800.ml1" using 2:3 title "X=800 mm"

set output "../../ps/ml10850m.eps"
plot "../../data/0850.ml1" using 2:3 title "X=850 mm"

set output "../../ps/ml10900m.eps"
plot "../../data/0900.ml1" using 2:3 title "X=900 mm"

set output "../../ps/ml10950m.eps"
plot "../../data/0950.ml1" using 2:3 title "X=950 mm"

set output "../../ps/ml11000m.eps"
plot "../../data/1000.ml1" using 2:3 title "X=1000 mm"


set output "../../ps/ml11050m.eps"
plot "../../data/1050.ml1" using 2:3 title "X=1050 mm"


set title "PLANE MIXING LAYER: <u'u'> Velocity Profile"
set key
set ylabel "<u'u'> [(m/s)**2]"
set xlabel "y [mm]"
set yrange [-1:12]


set output "../../ps/ml10030f.eps"
plot "../../data/0030.ml1" using 2:4 title "X=30 mm"

set output "../../ps/ml10040f.eps"
plot "../../data/0040.ml1" using 2:4 title "X=40 mm"

set output "../../ps/ml10050f.eps"
plot "../../data/0050.ml1" using 2:4 title "X=50 mm"

set output "../../ps/ml10070f.eps"
plot "../../data/0070.ml1" using 2:4 title "X=70 mm"

set output "../../ps/ml10100f.eps"
plot "../../data/0100.ml1" using 2:4 title "X=100 mm"

set output "../../ps/ml10200f.eps"
plot "../../data/0200.ml1" using 2:4 title "X=200 mm"

set output "../../ps/ml10250f.eps"
plot "../../data/0250.ml1" using 2:4 title "X=250 mm"

set output "../../ps/ml10300f.eps"
plot "../../data/0300.ml1" using 2:4 title "X=300 mm"

set output "../../ps/ml10350f.eps"
plot "../../data/0350.ml1" using 2:4 title "X=350 mm"

set output "../../ps/ml10400f.eps"
plot "../../data/0400.ml1" using 2:4 title "X=400 mm"

set output "../../ps/ml10450f.eps"
plot "../../data/0450.ml1" using 2:4 title "X=450 mm"

set output "../../ps/ml10500f.eps"
plot "../../data/0500.ml1" using 2:4 title "X=500 mm"

set output "../../ps/ml10550f.eps"
plot "../../data/0550.ml1" using 2:4 title "X=550 mm"

set output "../../ps/ml10600f.eps"
plot "../../data/0600.ml1" using 2:4 title "X=600 mm"

set output "../../ps/ml10650f.eps"
plot "../../data/0650.ml1" using 2:4 title "X=650 mm"

set output "../../ps/ml10725f.eps"
plot "../../data/0725.ml1" using 2:4 title "X=725 mm"

set output "../../ps/ml10750f.eps"
plot "../../data/0750.ml1" using 2:4 title "X=750 mm"

set output "../../ps/ml10800f.eps"
plot "../../data/0800.ml1" using 2:4 title "X=800 mm"

set output "../../ps/ml10850f.eps"
plot "../../data/0850.ml1" using 2:4 title "X=850 mm"

set output "../../ps/ml10900f.eps"
plot "../../data/0900.ml1" using 2:4 title "X=900 mm"

set output "../../ps/ml10950f.eps"
plot "../../data/0950.ml1" using 2:4 title "X=950 mm"

set output "../../ps/ml11000f.eps"
plot "../../data/1000.ml1" using 2:4 title "X=1000 mm"


set output "../../ps/ml11050f.eps"
plot "../../data/1050.ml1" using 2:4 title "X=1050 mm"