!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for interstellar
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2008/08/21 11:39:41 mkorpi Exp $',
  nt=100000, it1=10, isave=50, itorder=3, ialive=1
  dspec=2e-2, ou_spec=T, ro_spec=T, ss_spec=T !, fft_switch='fftpack'
  dsnap=0.07, dvid=1e-2, dtmin=1e-19
  ip=20
  !bcx = 'p','p','p', 'p',    'p'
  !bcy = 'p','p','p', 'p',    'p'
  !bcz = 's','ouf:s','ouf:a','a2','a2:cT'
  lwrite_aux=T
/
&eos_run_pars
/
&hydro_run_pars
/
&density_run_pars
  lupw_lnrho=T
  !cs2bot=100.0    ! =0.01
  cs2top=500000.0  !cs2top=10000.0 ! cs2top is the right paramerter 
/
&forcing_run_pars
  !iforce='gaussianpot', force=50.
/
&entropy_run_pars
  lupw_ss=T
  !iheatcond='K-const',hcond0=1e-8
/
&magnetic_run_pars
  eta=5e-1
/
&interstellar_run_pars
  lSNI=F, lSNII=F, t_next_SNI=100.0
  !heating_select='cst'
  !heating_rate=0.015
  !cooling_select='SSr'
  cooling_select='off'
  heating_select='off'
/
&viscosity_run_pars
  ivisc='nu-const','nu-shock'
  nu=5e-1, nu_shock=1.0
/
&shock_run_pars
  lshock_first=T
  lshock_max5=F
  lforce_periodic_shockviscosity=T
/

