!  -*-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.2 2013/10/29 23:23:33 brandenb Exp $'
  nt=500000, it1=100, 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:s','fB:s','s:a','s','s'
  !bcz = 'a:s','a:s','a','StS','StS','fB:c1','fB:nil','s:nil','s','s'
  !bcz = 'a:s','a:s','a','StS','StS','fB:pot','fB:nil','s:nil','s','s'
  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
/
&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,0.0,1.4104739
  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
  !ldivu_perp=T
/
