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'
set ylabel 'Aerodynamic Roughness'

# calc. for ridge aerodynamic roughness
hl(szrgh,sxprg) = ismin(0.20,szrgh/sxprg)
wzorg(szrgh,sxprg) = (szrgh>5)?szrgh * 1/(-64.1+135.5*hl(szrgh,sxprg)+(20.84/sqrt(hl(szrgh,sxprg)))):0

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

plot wzorg(x,400) ls 1 with lines title '400mm spacing' \
   , wzorg(x,1000) ls 2 with lines title '1000mm spacing' \
   , wzorg(x,2000) ls 3 with lines title '2000mm spacing' \
   , wzorg(x,4000) ls 4 with lines title '4000mm spacing' \
   , wzorg(x,8000) ls 5 with lines title '8000mm spacing'

#pause -1