!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Convection in vertically stratified atmosphere/solar convection zone
! Initialisation parameters
!
&init_pars
  cvsid='$Id: start.in,v 1.1 2014/08/03 03:27:38 brandenb Exp $'
  unit_length=1e8, unit_velocity=1e5, unit_density=1e0, unit_magnetic=1e0
  xyz0=-2.,-2., 0.0
  xyz1= 2., 2., 8.0
  lperi= T, T, F
  lwrite_aux=T
  lwrite_ic=F
/
&eos_init_pars
  lss_as_aux=T, lpp_as_aux=T, lcp_as_aux=T, lgamma_as_aux=T, lcs_as_aux=T
  lconst_yH=T, yH_const=1.
  xHe=0.1
/
&hydro_init_pars
  !inituu='gaussian-noise', ampluu=0.1
/
&density_init_pars
  !initlnrho='stratification', strati_type='lnrho_lnTT'
  initlnrho='exp_zbot', rho_left=1e-3, Hrho=5.
  lreduced_sound_speed=T, reduce_cs2=.02
/
&grav_init_pars
  gravz_profile='const', gravz=-274.
/
&entropy_init_pars
  initlnTT='const_TT', TT_const=5e4
/
!&magnetic_init_pars
!/
&radiation_init_pars
  radx=0,rady=0,radz=1,rad2max=1     ! use only vertical rays
  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=1e3, kappa2_cst=1e3
/
