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 pt 9

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

set xlabel 'Effective Biomass Drag Coefficient'
set ylabel 'U* / U*v'

wus_ratio(brcd)=(brcd>0.0001)?ismin(1,0.86*exp(-brcd/0.0298)+0.25*exp(-brcd/0.356)):1

#set log x
#set log y
set xrange [0.00001:0.25]
set yrange [0:]
set key top right

# using decimal in numbers to ensure real division, not integer divistion
plot wus_ratio(x) ls 1 with lines title 'predicted' \
   , '-' using 1:2 ls 6 title 'Lyles and Allison (1976)' with points \
   , '-' using 1:2 ls 4 title 'van de Ven et al. (1989)' with points
0.0131341 0.870564
0.0172101 0.741127
0.0276268 0.605428
0.0280797 0.528184
0.0321558 0.674322
0.0326087 0.59499
0.0430254 0.359081
0.0475543 0.450939
0.0484601 0.302714
0.048913 0.321503
0.0647645 0.354906
0.0683877 0.250522
0.0692935 0.26096
0.0738225 0.23382
0.0819746 0.246347
0.0973732 0.300626
0.113225 0.202505
0.116395 0.175365
0.118659 0.160752
0.162138 0.146138
0.243659 0.102296
0.0176311 0.753506
0.028481 0.507069
0.0488246 0.202097
0.0646474 0.354366
0.0687161 0.18523
0.162297 0.196975
0.243671 0.13784
0.00769857 0.970772
0.00770328 0.960334
0.0685585 0.348643
0.0686529 0.139875
0.0740353 0.2401
EOF
0.00271625 0.991644
0.00868279 0.793678
0.0145683 0.774879
0.0146531 0.587378
0.0445639 0.424629
0.0472538 0.474606
0.0594449 0.507838
0.0626123 0.501562
0.078519 0.316012
EOF

#pause -1