!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Initialisation parameters for interstellar
!
&init_pars
  cvsid='$Id: start.in,v 1.4 2020/09/04 08:21:02 brandenb Exp $'
  lperi=T,T,T
  xyz_units='pi','',''
  xyz0=-0.5,-0.5,-0.5
  xyz1=+1.5,+1.5,+1.5
  lwrite_aux=T,lread_aux=F
/
&eos_init_pars
  cs0=1.,gamma=1. 
/
&hydro_init_pars
/
&density_init_pars
  initlnrho='blob', ampllnrho=10., radius_lnrho=.5
/
&particles_init_pars
  initxxp='gaussian-x+random', gaussian_x_frac=.3
  initvvp='constant'
  vpx0=0.
  vpy0=0.
  vpz0=0.
  bcpx='p'
  bcpy='p'
  bcpz='p'
  xp0 =.35 !(/sqrt(2), R^2 = 2*sigma^2, sigma=R/sqrt(2)
  !plot_io,var.x,var.np_ap1 & oplot,col=122,var.x,275.+850.*exp(-.5*var.x^2/.35^2)
  interp_pol_rho='cic'
  interp_pol_TT='cic'
  interp_pol_uu='cic'
  lparticlemesh_tsc=T
  lnp_ap_as_aux=T
/
&particles_radius_init_pars
  initap='constant'
  !ap0=5e-6, 1e-5, 2e-5, 5e-5, 1e-4, 2e-4, 5e-4, 1e-3
  !ap0=1e-4, 2e-4, 5e-4, 1e-3, 2e-3, 5e-3, 1e-2
  ap0=1e-4
  lsweepup_par=F
  rhopmat=1e3
/
