!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for forced helical MHD turbulence
!
&run_pars
  cvsid='$Id: run.in,v 1.3 2017/06/20 17:06:58 brandenb Exp $',
  nt=1900000, it1=100, isave=5000, itorder=3
  dsnap=50., dvid=10., ialive=100, d2davg=100.
  tmax=14800
  slice_position='s'
  xtop_slice=1.
  ytop_slice=1.30
/
&eos_run_pars
/
&hydro_run_pars
  lupw_uu=T 
/
&density_run_pars
lupw_rho=T
!  diffrho_hyper3=6e-11, idiff='hyper3'
  ldensity_profile_masscons=T
/
&grav_run_pars
  gravx_profile='kepler',gravx=3.,
/
&forcing_run_pars
  iforce='helical', force=0.2, relhel=1.
  iforce_profile='surface_x_stepy'
  !iforce_profile='surface_stepx_cosy*siny**n_hel_sin_pow', n_hel_sin_pow=6.
  !iforce_profile='surface_helx_cosy*siny**n',r_ff=.8,width_ff=0.01
  width_ff=.02
  lscale_kvector_fac=T
  scale_kvectorx=20.9440
  scale_kvectory=20.9440
  scale_kvectorz=20.9440
! iforce_profile='surface_x' 
  r_ff=1.,lhelical_test=F
/
&magnetic_run_pars
  !eta=5e-5 !(good for 128 points in radius)
  eta=2e-4
/
&viscosity_run_pars
  !nu=5e-5
  nu=2e-4
/
