!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for forced helical MHD turbulence
!
&run_pars
  cvsid='$Id: run.in,v 1.2 2019/08/28 02:51:10 brandenb Exp $',
  nt=40000, it1=1, isave=1000, itorder=3
  nt=40, it1=1, isave=1000, itorder=3
  nt=400, it1=1, isave=1000, itorder=3
  dsnap=10., dvid=.02, ialive=10
  dspec=.1, ou_spec=T, ab_spec=F, GWs_spec=T, GWh_spec=T
  lpoint=3
  dt=.02
  dt=.01
  lpencil_check=F
  lpencil_check_small=F
  lkinflow_as_aux=T
/
&hydro_run_pars
  !kinematic_flow='delta_correlated'
  !kinematic_flow='potential', kx_uukin=1., ky_uukin=1., kz_uukin=1.
  !kinematic_flow='Beltrami-x', kx_uukin=1., relhel_uukin=1., omega_kinflow=0.
  !kinematic_flow='poshel-roberts', relhel_uukin=0.
  !kinematic_flow='Galloway-Proctor-nohel', omega_kinflow=1., eps_kinflow=1.
  !kinematic_flow='potential', kx_uukin=1., ky_uukin=1., kz_uukin=1., cx_uukin=1., cy_uukin=1., cz_uukin=1., lrandom_ampl=T, ampl_kinflow=10.
  kinematic_flow='potential_random', kx_uukin=1., ky_uukin=1., kz_uukin=1.
  lrandom_ampl=T, dtforce=.1, lwrite_random_ampl=T, ampl_kinflow=3.
  lrandom_location=T, lwrite_random_location=T
  lrandom_wavenumber=T, lwrite_random_wavenumber=T
  !kinematic_flow='Beltrami-x', kx_uukin=2., relhel_uukin=0., omega_kinflow=0.
  !ampl_kinflow=1.4142136  !(so that .5*<u^2>=1
/
&special_run_pars
  !cstress_prefactor='1'
  cstress_prefactor='16pi'
  !ctrace_factor='1/3'
  !nscale_factor=.0
  ldebug_print=F, lkinGW=T
/
&power_spectrum_run_pars
/
