!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  sample file for 1-D sound wave
!
&run_pars
  cvsid="$Id: run.in,v 1.2 2026/08/24 18:30:42 brandenb Exp $"
  nt=99000000, it1=2000, isave=2000, itorder=3
  dvid=-1e3, dsnap=-1e5, dtmin=-1e-12, cdtv=0.10
  dvid=1e8, dsnap=1e8, dtmin=-1e-12, cdtv=0.10
  dspec=.05, lspec_start=F, ou_spec=T, ab_spec=T, ele_spec=T, sp_spec=T, mu_spec=T
  Lor_spec=T, OmU_spec=T, ro_spec=T, gph_spec=T, aBE_spec=T, ABE2_spec=T, uBE_spec=T
  lwrite_aux=T, lpencil_check=F
  lpencil_check_small=F
  ip=16
  lit1_logspacing=T, tmax_logspacing=1e8
/
&eos_run_pars
/
&hydro_run_pars
  lupw_uu=T
/
&density_run_pars
  lupw_lnrho=T
  lrelativistic_eos=T
  idiff='kap-tdep', diffrho=2e4, Sc=1., tdep_kap_type='Sc'
/
&viscosity_run_pars
  ivisc='nu-tdep', nu=2e4, PrM=10., tdep_nu_type='PrM_sigEm'
/
&magnetic_run_pars
  ldensity_add_je_heating=F
  llorentzforce=F
/
&power_spectrum_run_pars
  lcorrect_integer_kcalc=T
  trigger_spec='lnascale'
/
&backreact_infl_run_pars
  lbackreact_infl=T, ldt_backreact_infl=T, Ndiv=0, cdt_phi=.1, lold_ldt_phi=F
  echarge_type='erun'
  lrho_chi_inhom=T
  Gamma_phi0=1e-9, lrho_rad_apply2=F
  lheating_keep_on=T, lwstate_crit=T
  solve_phi_criterion='rhophi', a4rhophim_crit=1e-30
  lold_lrho_chi_dtconstraint=F
  lrho_chi=T
  lswitch_toMHD_when_heating=F, lswitch_toMHD_when_nophi=F
  lswitch_toMHD_at_lna=T, lna_switch_toMHD=4.17
  lsmooth_Gamma_phi=T, dlnascale_reheating=.5, lg_Gamma_phi_fraction_firststep=-8.
  linclude_rho_EB_in_wstate=T
/
&disp_current_run_pars
  luse_scale_factor_in_sigma=T
  lcollinear_EB_aver=T, sigB_prefactor=1., sigE_prefactor=1.
  lallow_bprime_zero=T
  lphi_hom=F, lphi_linear_regime=F
  ldensity_add_je_heating=T, lohmic_heating_ee=T, lohmic_heating_justee=T
  lapply_Gamma_corr=F
  alpf=90.
  llorentzforce_ee=T
  cdt_sigE=1.
  lresistive_gauge_ee=T
  llate_reset_el_pencil=T
  linclude_dphiB_in_MHD=T, linclude_gphixE_in_MHD=T
  lcorrect_sign_adphiB_term=T, lignore_adphiB_term_in_MHD_current=F
/