!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Convection in a spherical wedge
!
&run_pars
  cvsid='$Id: run.in,v 1.1 2014/03/28 18:27:28 mkorpi Exp $',
  nt=5000000, it1=200, isave=1000, itorder=3, tmax=50000
  dsnap=200., dvid=1000., ialive=0,
  cdtv=0.30, d2davg=10., cdt=0.9
  slice_position='s', xtop_slice=0.98, ytop_slice=1.57
!  lread_oldsnap_nomag=T
/
&eos_run_pars
!  sigmaSBt=4.8070364E-03 ! drho=11.3
!  sigmaSBt=0.01 ! drho=13.9, Hp/dr=13
!  sigmaSBt=0.015 ! drho=16.1, Hp/dr=11.3
  sigmaSBt=0.020 ! drho=16.1, Hp/dr=11.3
/
&hydro_run_pars
  lupw_uu=T
  Omega=0.69
  lremove_mean_angmom=T
/
&density_run_pars
  lupw_lnrho=T
  cs2top=0.162371
/
&grav_run_pars
  gravx_profile='kepler',gravx=3.,
/
&entropy_run_pars
  Fbot=3.0701169E-05
  iheatcond='K-profile',hcond0=1e-8 ! not used but needed to get 
				    ! into calc_heatcond
  lupw_ss=T
  lhcond_global=T
!  chi_t=2e-5
!  chit_prof1=0., chit_prof2=12.5
!  chi_t=5e-5
!  chit_prof1=0., chit_prof2=5
!  chi_t=1e-4
!  chit_prof1=0., chit_prof2=2.5
  chi_t=2e-4
  chit_prof1=0., chit_prof2=1.35
!  xbot=0.75, xtop=0.98, widthss=0.01
  xbot=0.75, xtop=0.95, widthss=0.01
/
&viscosity_run_pars
  ivisc='nu-const'
  nu=5e-5
/
!&magnetic_run_pars
!  lupw_aa=T
!  eta=5e-5
!  lrun_initaa=T
!/