!  -*-f90-*-  (for Emacs)    vim:set filetype=fortran:  (for vim)
!
! Spreading layer on the surface of NS
! Initialisation parameters
!
&init_pars
  cvsid='$Id: start.in,v 1.1 2013/10/17 21:31:08 brandenb Exp $',
  ip=14,                        ! debugging parameter
  unit_system='cgs'
  unit_length=1e3,  unit_velocity=1.e9
  unit_density=2e-1 !hot_star 1e-1
  ! unit_temperature=30.e6
  pretend_lnTT=F
  xyz0  = 0., 0., 1000.,        ! first corner of box
  Lxyz  =  12.,   2,  15.,   ! box size  y-> 10/32=0.3125
  lperi =  F ,   T ,   F
  bcx = 'stp','stp','stp', 'stp','stp'
  bcz = 'stp','stp','stp', 'stp','stp'
! bot
    fbcx2(1)=0.,   fbcx2(2)=0., fbcx2(3)=0.,
    fbcz1(1)=0.,
    fbcz1_1(2)=0., 
    fbcz1_1(3)=0.,
    fbcz1_1(4)=100.,
 
! top
    fbcz2(1)=0.
  ! fbcz2_1(2)=0.,    
   fbcz2_1(3)=0., 
   fbcz2_1(4)=20.,
 /
&eos_init_pars
  cs0=1.
  mu=.6
/
&hydro_init_pars
  ! inituu='step_xz', uu_left=0., uy_left=0., uy_right=0.
/
&density_init_pars
 /
&entropy_init_pars
   lviscosity_heat=T
/
&radiation_init_pars
  bc_rad='p:p','p:p','S+F:S-F'
  radx=1, rady=0, radz=1, rad2max=1
  opacity_type='kappa_es'
  Frad_boundary_ref=0.
!  lnatalia=T
/
&viscosity_init_pars
  ivisc='nu-const'
  nu=1e-2  !1e-3
 /
&neutron_star_init_pars
 
  uu_init=0. !everywhere

  laccelerat_zone=T
  ldecelerat_zone=T
  lsurface_zone=F
 
  leffective_gravity=T

  rho_star=100., rho_disk=20.

  lmass_source_NS=F
 
  L_disk_point=46

  lnstar_entropy=T
  lnstar_T_const=F
  lnstar_1D=F 
 
  R_star=1000.
  M_star=1.9e5
  T_star=0.03  !0.01
  T_disk=0.03  !0.01
!  accretion_flux=1e-3

  l1D_cooling=F
  l1D_heating=F
 

  ltop_velocity_kep=T
  lextrapolate_bot_density=F
  lgrav_x_mdf=F

  nu_for_1D=1e-4
  star_rot=1.
  hot_star=T
  T_disk_const=T
  !  rotat=1
/