plotname = "cylinder2d_single_refinement_n20_re100_16s"

load 'common/velocity_fluid.gnuplot'

plot 'data/cylinder2d_single_refinement_n20_re100_16s_1.csv' using ($6*10):($7*10):(sqrt($2*$2+$3*$3)) palette pt 5 ps 1.3 notitle,\
     'data/cylinder2d_single_refinement_n20_re100_16s_0.csv' using ($6*10):($7*10):(sqrt($2*$2+$3*$3)) palette pt 5 ps 0.7 notitle

load 'common/velocity_wrapper.gnuplot'

set obj 1 rect from 1.0,0.55 to 9,3.55
set obj 1 fillstyle empty border -1 front lw 1

load 'common/include_fluid.gnuplot'