!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!

&init_pars
  ip = 14
  xyz0 =-0.5, 0.0, -0.05
  xyz1 = 0.5, 2.0, 0.05
  fbcy1 = 0.0,1.0,0.,1.0
  lperi=T,F,T
  bcy = 'e2','e2','e2','e2'
/
&eos_init_pars
  cs0=10
  !cp=1.0
  gamma=1.0
/
&hydro_init_pars
/
&density_init_pars
  ldensity_nolog=T
  initlnrho='const_lnrho'
  lnrho_const=0.0
/
&solid_cells_init_pars
  ncylinders=1
  cylinder_radius=50.0e-3
  cylinder_xpos=0.0
  cylinder_ypos=1.0
  cylinder_zpos=0.0
  skin_depth_solid=4e-3
  initsolid_cells='cylinderstream_y'
  init_uu=1.
  grid_func_ogrid = 'sinh', , ! Grid stretching in radial direction
  xyz_star_ogrid=50.0e-3      ! Smallest cells at solid surface
  coeff_grid_o=18, ,          ! Streching paramater
  interpolation_method=1      ! Linear interoplation between grids
  lcheck_interpolation=.false.
  SBP=.true.                  ! Use summation-by-parts operators
  lfilter_solution=.true.     ! Use Padé filtering
  lfilter_rhoonly=.false.     ! Switch to true to only filter rho-component
  particle_interpolate=3      ! Quadratic interpolation of velocities to particle positions
  lparticle_uradonly=.false.  ! Quadratic interpolation for ur and uth
  lspecial_rad_int=.true.     ! Special particle interpolation near solid surface
/
&NSCBC_init_pars
  nscbc_bc = '','part_ref_inlet:part_ref_outlet',''
  nscbc_sigma_in=4.0
  nscbc_sigma_out=0.28
  p_infty=100
/
&particles_init_pars
  initxxp='random-box'
  initvvp='constant'
  vpx0=0.
  vpy0=1.0
  vpz0=0.
  bcpx='p'
  bcpy='rmv'
  bcpz='p'
  xp0=-50.0e-3
  yp0=0.01
  zp0=0
  Lx0=100.0e-3
  Ly0=0.05
  Lz0=0
  mean_free_path_gas=67e-9
  interp_pol_rho='cic'
  interp_pol_TT='cic'
  interp_pol_uu='cic'
  linsert_particles_continuously=T
  lnostore_uu=.false.
/
&particles_radius_init_pars
  initap='specify'    ! Specified distribution of particles
  radii_distribution=0.3,0.225,0.175,0.15,0.1,0.05 
  ap0 = 3.6742346142E-04, ! St = 0.6
        4.2426406871E-04, ! St = 0.8
        4.7434164903E-04, ! St = 1
        8.2158383626E-04, ! St = 3
        1.0606601718E-03, ! St = 5
        1.5000000000E-03, ! St = 10
  ! Note that more than 4 ap0 requires input in cparam.local
  ! See line about ndustrad contrubtion in cparam.local
  lsweepup_par=F
  rhopmat=1e3           ! Particle density / fluid density
/