# Runtime options file for Phantom, written 16/02/2024 14:28:53.7
# Options not present assume their default values
# This file is updated automatically after a full dump

# job name
             logfile =   ext01.log    ! file to which output is directed
            dumpfile =   ext_00000    ! dump file to start from

# options controlling run time and input/output
                tmax =   1.168E+05    ! end time
               dtmax =  557.709602    ! time between dumps
                nmax =          -1    ! maximum number of timesteps (0=just get derivs and stop)
                nout =          -1    ! write dumpfile every n dtmax (-ve=ignore)
           nmaxdumps =          -1    ! stop after n full dumps (-ve=ignore)
            twallmax =      000:00    ! maximum wall time (hhh:mm, 000:00=ignore)
           dtwallmax =      024:00    ! maximum wall time between dumps (hhh:mm, 000:00=ignore)
           nfulldump =          10    ! full dump every n dumps
            iverbose =           0    ! verboseness of log (-1=quiet 0=default 1=allsteps 2=debug 5=max)

# options controlling accuracy
              C_cour =       0.300    ! Courant number
             C_force =       0.250    ! dt_force number
                tolv =   1.000E-02    ! tolerance on v iterations in timestepping
               hfact =       1.000    ! h in units of particle spacing [h = hfact(m/rho)^(1/3)]
                tolh =   1.000E-04    ! tolerance on h-rho iterations

# options controlling hydrodynamics, shock capturing
               alpha =  0.266669082  ! shock viscosity parameter
              alphau =       1.000    ! shock conductivity parameter
                beta =       2.000    ! beta viscosity

# options controlling damping
               idamp =           0    ! artificial damping of velocities (0=off, 1=constant, 2=star, 3=disc)

# options controlling equation of state
                ieos =           6    ! eqn of state (1=isoth;2=adiab;3=locally iso;8=barotropic)
                  mu =       2.381    ! mean molecular weight

# options controlling cooling
              C_cool =       0.050    ! factor controlling cooling timestep
            icooling =           0    ! cooling function (0=off, 1=library (step), 2=library (force),3=Gammie, 4=ISM, 5,6=KI02, 7=powerlaw)
          use_mcfost =           F    ! use the mcfost library
    use_mcfost_stars =           F    ! Fix the stellar parameters to mcfost values or update using sink mass
mcfost_computes_Lacc =           F    ! Should mcfost compute the accretion luminosity
     mcfost_uses_PdV =           T    ! Should mcfost use the PdV work and shock heating?
    mcfost_keep_part =       0.999    ! Fraction of particles to keep for MCFOST
                 ISM =           0    ! ISM heating : 0 -> no ISM radiation field, 1 -> ProDiMo, 2 -> Bate & Keto
    mcfost_dust_subl =           F    ! Should mcfost do dust sublimation (experimental!)

# options controlling sink particles
     h_soft_sinksink =       0.000    ! softening length between sink particles
               f_acc =       0.800    ! particles < f_acc*h_acc accreted without checks
      r_merge_uncond =       0.000    ! sinks will unconditionally merge within this separation
        r_merge_cond =       0.000    ! sinks will merge if bound within this radius

# options relating to external forces
      iexternalforce =           0    ! 1=star,2=coro,3=bina,4=prdr,5=toru,6=toys,7=exte,8=spir,9=Lens,10=dens,11=Eins,

# options controlling physical viscosity
           irealvisc =           0    ! physical viscosity type (0=none,1=const,2=Shakura/Sunyaev)
          shearparam =       0.100    ! magnitude of shear viscosity (irealvisc=1) or alpha_SS (irealvisc=2)
            bulkvisc =       0.000    ! magnitude of bulk viscosity

# options controlling dust
               idrag =           1    ! gas/dust drag (0=off,1=Epstein/Stokes,2=const K,3=const ts)
      ilimitdustflux =           T    ! limit the dust flux using Ballabio et al. (2018)
   icut_backreaction =           0    ! cut the drag on the gas phase (0=no, 1=yes)

# options for injecting/removing particles
               rkill =      -1.000    ! deactivate particles outside this radius (<0 is off)

# gravitational waves
                  gw =           F    ! calculate gravitational wave strain