!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  abiabatically stratified polytropic atmosphere
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2012/02/09 16:59:13 pkapyla Exp $'
  nt=1000000, it1=50, isave=50, itorder=3, ialive=10, cdt=.6
  nt=1000000, it1=20, isave=50, itorder=3, ialive=10, cdt=.6
  dsnap=1000., dvid=100.
  bcz = 's','s','a','a2','a','s','s','a'  !(vertical field condition)
  bcz = 's','s','a','a2','a','a','a','s'  !(perfect conductor)
  lpencil_check_small=T
!  ip=13
/
&eos_run_pars
/
&hydro_run_pars
  lupw_uu=T
/
&density_run_pars
  diffrho_hyper3=5e-11, idiff='hyper3'
/
&grav_run_pars
/
&entropy_run_pars
  tau_cool_ss=.316
!  tau_cool_ss=1e-2
  iheatcond='chi-const'
  chi_t=3e-2
  chi=1e-9
  lupw_ss=T
/
&magnetic_run_pars
  eta=5e-3, B_ext=.0, .15, .0
  eta=3e-2, B_ext=.0, .1, .0
  luse_Bext_in_b2=T
/
&magn_mf_run_pars
  lmeanfield_jxb=T
!  meanfield_Beq_profile='exp(z/H)', meanfield_Beq_height=-2., meanfield_Beq=1.
  meanfield_Beq=1.
  meanfield_qs=0., meanfield_qp=40., meanfield_qe=0.
  meanfield_Bs=.1, meanfield_Bp=.05
  qp_model='rational'
/
&viscosity_run_pars
  nu=3e-2
/
