set term svg size 800,600 dynamic enhanced

set termoption dashed

set style line 1 lt 1 lc rgb "black" lw 1
set style line 2 lt 2 lc rgb "black" lw 1
set style line 3 lt 3 lc rgb "black" lw 1
set style line 4 lt 4 lc rgb "black" lw 1
set style line 5 lt 5 lc rgb "black" lw 1
set style line 6 lt 1 lc rgb "black" lw 2


ismin(x,y) = (x>y)?y:x
ismax(x,y) = (x>y)?x:y

set xlabel 'Ridge Height / Ridge Spacing'
set ylabel 'Z0 / Ridge Height'

# Hagen and Armbrust, 1992
#0.0399192,0.0143622
#0.0549222,0.0367874
#0.0699034,0.0415748
#0.10954,0.0665197
#0.109747,0.0745827
#0.147615,0.103559
#0.175998,0.110362
#0.204366,0.105323
#0.210072,0.0982677

wzorgh(hl) = 1/(-64.1+135.5*hl+(20.84/sqrt(hl)))

set xrange [0:0.22]
#set yrange [0:1]
set key bottom

plot wzorgh(x) ls 1 with lines title 'Predicted' \
   , '-' using 1:2 ls 3 title 'Measured' with points
0.0399192 0.0143622
0.0549222 0.0367874
0.0699034 0.0415748
0.10954 0.0665197
0.109747 0.0745827
0.147615 0.103559
0.175998 0.110362
0.204366 0.105323
0.210072 0.0982677
EOF
#pause -1