To generate the black squares that appear in Figure 4, comment out line 247 in flow_shear.f90,

246            if (shift_state(iky) > shift_times(iky)) then
247               g0k(iky, shift_start) = 0.0
248            end if

recompile, and re-run the input files given in the global/ directory.
