# particles emitted from a 2nd circle used as boundary

seed	    	    12345
dimension   	    2
global              gridcut 0.0 comm/sort yes

boundary	    o r p

create_box  	    0 10 0 10 -0.5 0.5
create_grid 	    10 10 1 
balance_grid        rcb cell

global		    nrho 1.0 fnum 0.001

species		    air.species N O
mixture		    air N O vstream 10.0 0 0 temp 10.0

# 12 would be fine, 12.1 is simply to make arc visible in images

read_surf           data.circle
read_surf           data.circle group boundary invert origin 5 5 0 &
                    atrans 12.1 5 0 scale 4 4 1 clip
surf_collide	    1 diffuse 300.0 0.0
surf_modify         all collide 1

collide             vss air air.vss

region              slit block INF INF 4.5 5.5 INF INF

#fix		    in emit/surf air all normal yes perspecies no region slit
#fix		    in emit/surf air all perspecies yes #region slit
#fix		    in emit/face air xlo region slit

fix		    in emit/surf air boundary perspecies yes

timestep 	    0.0001

#dump                2 image all 10 image.*.ppm type type pdiam 0.1 &
#                    surf proc 0.01 size 512 512 axes yes 0.9 0.02
#dump_modify	    2 pad 4

fix                 1 balance 10 1.0 rcb part

stats		    100
stats_style	    step cpu np nattempt ncoll nscoll nscheck
run 		    300