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


# x = xh/zbr
# xw = xbrw/zbr

set xlabel 'Barrier Width / Barrier Height'
set ylabel 'Effective Barrier Porosity'

# increase effective porosity with barrier width
pb(xw,pbr) = pbr + (1 - exp(-0.5*xw))*0.3*(1-pbr)

set xrange [0:15]
set yrange [0:1.05]
set key bottom

plot pb(x,1.0) ls 1 with lines title 'optical porosity = 1.0' \
   , pb(x,0.8) ls 2 with lines title 'optical porosity = 0.8' \
   , pb(x,0.6) ls 3 with lines title 'optical porosity = 0.6' \
   , pb(x,0.4) ls 4 with lines title 'optical porosity = 0.4' \
   , pb(x,0.2) ls 5 with lines title 'optical porosity = 0.2' \
   , pb(x,0.0) ls 6 with lines title 'optical porosity = 0.0'

#pause -1