!                              -*-f90-*-     (for emacs)
!
! Initialisation parameters for MRI turbulence
!
&init_pars
  cvsid='$Id: start.in,v 1.2 2019/01/02 07:50:07 brandenb Exp $',
  unit_length=1e8               ! (= 1 Mm)
  unit_velocity=1e5             ! (= 1 km/s)
  unit_density=1e0              ! (= 1 g/cm^3)
  unit_temperature=1e0          ! (= 1 K)
  xyz0=-2.0,-2.0,-0.0
  xyz1= 2.0, 2.0,+1.0
  lperi=T,T,F
  bcz = 's','s','a','s:a2','s:a2','s'
  ip=14, lwrite_ic=F, lwrite_aux=T
/
&eos_init_pars
  cs0=40., rho0=1.276e-7, mu=.6
  cs0=30., rho0=1.520e-7, mu=.6
  cs0=16., rho0=2.75e-7, mu=.6  !(gives Sigma=7.1e-7)
  cs0=10., rho0=5.00e-7, mu=.6  !(gives Sigma=7.1e-7)
  cs0=5., rho0=8.78e-7, mu=.6  !(gives Sigma=7.1e-7)
/
&hydro_init_pars
  inituu='gaussian-noise', ampluu=0e-3
  Omega=12.
/
&density_init_pars
  !initlnrho='isothermal'
  initlnrho='stratification', strati_type='lnrho_ss'
/
&grav_init_pars
  gravz_profile='linear',nu_epicycle=12.
/
&entropy_init_pars
/
&shear_init_pars
  qshear=1.5
/
&radiation_init_pars
  radx=0,rady=0,radz=1,rad2max=1
  bc_rad='p:p','p:p','S:0', lradflux=T
  expo_rho_opa=1., expo_temp_opa=4.  !(Hminus)
  expo2_rho_opa=1., expo2_temp_opa=-3.5  !(Kramers)
  opacity_type='kappa_double_power_law', ref_rho_opa=1e-5, ref_temp_opa=2e4
  kappa_cst=2e12
/