!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
!  run parameters for interstellar
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2008/07/25 12:21:56 mkorpi Exp $',
  nt=80000, it1=10, isave=50, itorder=3, ialive=1
  dspec=2e-2, ou_spec=T, ro_spec=T, ss_spec=T !, fft_switch='fftpack'
  dsnap=1.0, dvid=1e-2, dtmin=1e-19
  ip=20
  bcx = 'p','p','p', 'p',    'p'
  bcy = 'p','p','p', 'p',    'p'
  bcz = 's','s','ouf:a','f:a2','a2:cT'
  lwrite_aux=T
/
&eos_run_pars
/
&hydro_run_pars
/
&density_run_pars
  lupw_lnrho=T
  !cs2bot=100.0    ! =0.01
  !cs2top=100000.0  !cs2top=10000.0 ! cs2top is the right paramerter 
  cs2top=10000.
/
&forcing_run_pars
  !iforce='gaussianpot', force=50.
/
&entropy_run_pars
  lupw_ss=T
  !iheatcond='K-const',hcond0=1e-8
/
&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=0.05,nu_shock=1.0
/
&shock_run_pars
  lshock_first=T
  lshock_max5=F
  lforce_periodic_shockviscosity=T
/