!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for forced helical MHD turbulence
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2019/02/17 07:19:16 brandenb Exp $',
  lpoint=14, mpoint=4, npoint=4, lpoint2=4, mpoint2=4, npoint2=4
  nt=20000000, it1=50, isave=1000, itorder=3
  nt=10, it1=1, isave=1000, itorder=3
  nt=1, it1=1, isave=1000, itorder=1 !(test)
  nt=5000, it1=10, isave=1000, itorder=1
  nt=30, it1=1, isave=1000, itorder=3
  nt=30, it1=1, isave=1000, itorder=3
  nt=18, it1=1, isave=1000, itorder=3
  nt=170, it1=5, isave=1000, itorder=3
  nt=17, it1=1, isave=1000, itorder=3
  nt=4, it1=1, isave=1000, itorder=3
  nt=40, it1=1, isave=1000, itorder=3
  nt=4, it1=1, isave=1000, itorder=3
  nt=400, it1=1, isave=1000, itorder=3
  nt=1, it1=1, isave=1000, itorder=3
  dsnap=10., dvid=6., ialive=10
  dsnap=10., dvid=.02, ialive=10
  dspec=.01, ou_spec=T, ab_spec=T, GWs_spec=T, GWh_spec=T
  !lpencil_check=F
  !lpencil_check_small=F
/
&eos_run_pars
/
&hydro_run_pars
/
&density_run_pars
/
&special_run_pars
  nscale_factor_conformal=0.
  cstress_prefactor='16pi'
  cstress_prefactor='6'
  !nscale_factor=.0
  !ldebug_print=F, lkinGW=T
  !diffhh=0e-12
/
&forcing_run_pars
  !iforce='gaussianpot'
  !width_ff=1.
  !width_ff=.2
  !force=.05
  !dtforce=200.
/
&magnetic_run_pars
  eta=1e-1
/
&viscosity_run_pars
  nu=1e-2
/
&power_spectrum_run_pars
  !lhalf_factor_in_GW=T
/
