!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Convection in vertically stratified atmosphere/solar convection zone
! Run parameters
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2013/11/12 09:11:13 brandenb Exp $'
  unit_length=1e8, unit_velocity=1e5, unit_density=1e-6, unit_magnetic=1e0
  nt=500000, it1=100, isave=1000, itorder=3, cdt=.7
  nt=10, it1=1, isave=1000, itorder=3, cdt=.7
  dsnap=0.1, dvid=0.005
  dsnap=.1, dvid=500.
  dsnap=5., dvid=10.
  iz=100
  bcz = 'a:s','a:s','a','StS','StS','fB:pot','fB:nil','s:nil','s'
  lpencil_check_small=F
  lwrite_aux=T
/
&eos_run_pars
  xHe=0.1
  yMetals=1e-4
/
&hydro_run_pars
  tdamp=1.0
  dampu=10.0
  ldamp_fade=T
/
&density_run_pars
  lupw_lnrho=T
/
&forcing_run_pars
  iforce='helical', force=1e-36, relhel=0.
  iforce_profile='surface_z', r_ff=0., width_ff=.1
/
&grav_run_pars
  gravz_profile='reduced_top'
  gravz=-275.
  reduced_top=0.5
  zgrav=0.1
/
&entropy_run_pars
  chi=5e-3
  lupw_lnTT=T
  TT_cor=4e3,tau_heat_cor=1e-3,zcor=0.1
/
&magnetic_run_pars
  iresistivity='eta-const','shock-perp'
  lweyl_gauge=T
  eta=5e-3, eta_shock=5.
  B_ext=0., 0., 10.
  va2max_jxb=2500., va2power_jxb=4
/
&radiation_run_pars
  bc_rad='p:p','p:p','S:c'
  radx=1, rady=1, radz=1, rad2max=1
  TT_top=4e3,tau_top=0.2
  lradflux=T
/
&viscosity_run_pars
  ivisc='nu-const','nu-shock'
  nu=5e-3, nu_shock=5.
/
&shock_run_pars
  lshock_first=F
/
