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

#Raupach 1992 data
#brcd	wz0/h
#0.0110812818	0.0085952794 # 6mm high cylinders
#0.022403713	0.0223604222 # 6mm high cylinders
#0.0452949725	0.0335698963 # 6mm high cylinders
#0.0898028252	0.0541452047 # 6mm high cylinders
#0.1815611814	0.0624956803 # 6mm high cylinders
#0.2341137923	0.1442679963
#0.4641592396	0.0832568547
#0.8025250549	0.1007998765
#1.6872366479	0.0503987773
#2.8053238244	0.0503987773
#3.0934706506	0.0297885936
#4.313402642	0.0610196904
#4.7564397512	0.0305091424

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

set xlabel 'Effective Biomass Drag Coefficient'
set ylabel 'Aerodynamic Roughness (mm) / Standing Biomass Height'

wzzovdht(brcd)=1./(17.27-(1.254*log(brcd)/brcd)-(3.714/brcd))

set log x
set log y
set xrange [0.01:10]
set yrange [0.001:1]
set key box top left opaque

# using decimal in numbers to ensure real division, not integer divistion
plot wzzovdht(x) ls 1 with lines title 'Predicted' \
   , '-' using 1:($2) ls 5 title 'Raupach (1992) 6mm height' with points \
   , '-' using 1:($2) ls 3 title 'Raupach (1992) 200mm height' with points
0.0110812818	0.0085952794 # 6mm high cylinders
0.022403713	0.0223604222 # 6mm high cylinders
0.0452949725	0.0335698963 # 6mm high cylinders
0.0898028252	0.0541452047 # 6mm high cylinders
0.1815611814	0.0624956803 # 6mm high cylinders
EOF
 0.2341137923	0.1442679963
 0.4641592396	0.0832568547
 0.8025250549	0.1007998765
 1.6872366479	0.0503987773
 2.8053238244	0.0503987773
 3.0934706506	0.0297885936
 4.313402642	0.0610196904
 4.7564397512	0.0305091424
EOF

#pause -1