!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Convection in vertically stratified atmosphere/solar convection zone
! Run parameters
!
&run_pars
  cvsid='$Id$'
  nt=41, it1=10, isave=200, itorder=3
  dsnap=20., dvid=0.2
  bcx = 'p','p','p', 'p',    'p'
  bcy = 'p','p','p', 'p',    'p'
  bcz = 's','s','a','a2','c1:cT'
  cdtv=0.30, cdt=0.4
  lpencil_check=T
/
&eos_run_pars
/
&hydro_run_pars
  Omega=0.1 ! angular velocity
  theta=0.  ! direction of angular velocity vector (0 : vertical)
/
&density_run_pars
  cdiffrho=0.
  !lupw_lnrho=T
  ! upwinding cannot be used togehter with non-equidistant
  ! grid because it involves higher derivatives
/
&grav_run_pars
  gravz_profile='const',gravz=-1.
/
&entropy_run_pars
  Luminosity=0., wheat=0.1, cool=15., wcool=0.2,
  iheatcond='K-profile', hcond0=8.e-3
  !lupw_ss=T ! use upwinding for entropy to avoid wiggles
/
&viscosity_run_pars
  ivisc='nu-const'
  nu=4e-3
/
! Nomenclature for boundary conditions:
!   Variables are (currently): ux,uy,uz,lam,ent,ax,ay,az
!   p      : periodic
!   s/a/a2 : symmetric/antisymmetric/antisymmetric w.r.t. boundary value
!   s:a2   : bc `s' at lower (left) and `a2' at upper (right) boundary

! tinit    -- initial time for smoothly switching on heating, etc