!                           -*-f90-*-     (for Emacs)
!
!  Run parameters
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2013/03/24 05:05:15 wlyra Exp $',
  ip=14,
  nt=100, it1=10, isave=1000, itorder=3
  cdt=0.4,
  cdtv=0.4,
  cdtv3=0.1
  dtmin=1e-6,
  dsnap=6.283 ,dvid=1.0
!  tmax=5.
!
  wborder_int=0.05
  wborder_ext=0.075
!
  lpencil_check=T
  lpencil_check_small=T
  
  !tmax=5.0389660475700824
  !tmax=6.
!
  tborder=0.1
!  tmax=13.
!
!  lisotropic_advection=T
  max_walltime=93600.
!  tmax=629.
  tmax=1260.
  lwrite_dvar=T
  lread_aux=T
/
&eos_run_pars
/
&hydro_run_pars
  borderuu='initial-condition'
  lpressuregradient_gas=F
/
&density_run_pars
  borderlnrho='initial-condition'
  idiff='hyper3-cyl' ,'shock'
  diffrho_hyper3=5e-4 ,diffrho_shock=10.
  density_floor=1e-3
/
&grav_run_pars
/
&viscosity_run_pars
  ivisc='hyper3-cyl' ,'nu-shock','simplified'
  nu_hyper3=5e-4 ,nu_shock=10.,nu=1e-4
/
&special_run_pars
  ldust_pressureforce=T
  factor_photoelectric=1.0    ! photo heating
!  Sentropy=0.0                ! polytropic 
  factor_localiso=0.5         ! isothermal
/
&shock_run_pars
/
&particles_run_pars
  ldragforce_dust_par=T
  ldragforce_gas_par=T
  lmigration_redo=T
  lcheck_exact_frontier=F
  lparticlemesh_tsc=T
  ldraglaw_variable=T
  gravr_profile='newtonian'
  lcylindrical_gravity_par=T
  lcommunicate_rhop=T
/