# 3d flow around a sphere

seed	    	    12345
dimension   	    3

global              gridcut 0.1

boundary	    o r r

create_box  	    -2 2 -2 2 -2 2

create_grid         20 20 20

balance_grid        rcb cell

global		    nrho 1.0 fnum 0.0005

species		    ../../data/air.species N O
mixture		    air N O vstream 100.0 0 0 

read_surf           1 sdata.sphere
surf_collide	    1 diffuse 300.0 0.0
surf_modify         collide 1 1

collide		    vss air ../../data/air.vss

fix		    in inflow air xlo 

timestep 	    0.0001

compute             2 surf all n press ke
fix                 save ave/surf 1 50 50 c_2 ave running
region              slab block INF INF INF INF -0.1 0.1
dump                2 image 50 image.*.ppm type type pdiam 0.03 &
		    view 70 120 size 512 512 axes yes 0.9 0.02 &
                    gridz -0.8 proc gline yes 0.005 &
                    surf f_save[2] 0.0
dump_modify	    2 pad 4 region slab
dump_modify         2 cmap surf min max cf 0.0 2 min orange max green

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