!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for kinematic dynamo run with Glen Roberts flow
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2013/05/27 05:52:58 brandenb Exp $',
  nt=10000000, it1=500, isave=20000, itorder=3
  dsnap=.0540, dvid=.05, dtmin=1e-19
  lpoint=26, mpoint=8, npoint=1, iproc_point=0
  bcx='s:fix','a:spd','a:spd'
  slice_position='w'
  lkinflow_as_aux=T
/
&hydro_run_pars
  kinematic_flow='spoke-like-NSSL', uphi_rbot=.7, uphi_rtop=.97
  uphi_step_width=.03, ampl_kinflow=5e4
/
&magnetic_run_pars
  iresistivity='eta-xdep', xdep_profile='step', eta=1e-2
  eta_width=-.01, eta_jump=100., eta_x0=1.
/
&magn_mf_run_pars
  x_surface2=.7, x_surface=1.
  meanfield_kf=1., meanfield_kf_profile='sqrt_surface_x1'
  meanfield_kf_width2=.02, meanfield_kf_width=.001
  meanfield_etat=1.
  meanfield_Beq_profile='sqrt(kf_x)'
  alpha_effect=0e-4, alpha_quenching=.01
  alpha_profile='surface_x2*cosy'
  alpha_profile='sqrt_surface_x2*cosy'
  alpha_width2=.03, alpha_width=.001
  meanfield_etat_profile='sqrt_surface_x1', meanfield_etat_width=.05
  lalpha_Omega_approx=T
/
&magn_mf_demfdt_run_pars
  tau1_emf=1., eta_emf_over_etat=1.
/
&special_run_pars
  kf_alpm=1., alpmdiff=1.
/
