!                              -*-f90-*-     (for emacs)
!
! Initialisation parameters for MRI turbulence
!
&init_pars
  cvsid='$Id: start.in,v 1.1 2016/02/22 05:22:33 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)
  unit_magnetic=1e0             ! (= 1 G)
  xyz0=-2.,-6.28,-3.
  xyz1= 2., 6.28,+3.
  lperi=T,T,F
  bcz = 's','s','a','a2','a2','s','s','a'
  ip=14, lwrite_ic=F, lwrite_aux=T
/
&eos_init_pars
  cs0=16., rho0=1e-6, mu=.6
/
&hydro_init_pars
  inituu='gaussian-noise', ampluu=1e-3
  Omega=12.
/
&density_init_pars
  initlnrho='stratification', strati_type='lnrho_ss'
/
&magnetic_init_pars
  initaa='Alfvenz-rot-shear','coswave-Ay-kx'
  amplaa=1e-6,2e-7
  amplaa=1e-3,2e-4
  kz_aa=6.,0.
  B_ext=0.,0.,200.
/
&grav_init_pars
  gravz_profile='linear',nu_epicycle=12.
/
&entropy_init_pars
/
&shear_init_pars
  qshear=1.5
/
&radiation_init_pars
  radx=1,rady=1,radz=1,rad2max=1
  bc_rad='p:p','p:p','0: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
/
