# ------------------  INPUTS TO MAIN PROGRAM  -------------------
max_step = 0

#if AMREX_SPACEDIM == 2

geometry.coord_sys   =  1
geometry.prob_lo     =  0.0   -1.6e9
geometry.prob_hi     =  1.6e9  1.6e9
amr.n_cell           =  32     64
castro.lo_bc         =  3      2
castro.hi_bc         =  2      2

#elif AMREX_SPACEDIM == 3

geometry.coord_sys   =  0
geometry.prob_lo     = -1.6e9 -1.6e9 -1.6e9
geometry.prob_hi     =  1.6e9  1.6e9  1.6e9
amr.n_cell           =  64     64     64
castro.lo_bc         =  2      2      2
castro.hi_bc         =  2      2      2

#endif

amr.max_level        = 1
amr.ref_ratio        = 2 2 2 2 2 2 2 2 2 2 2
amr.n_error_buf      = 2 2 2 2 2 
amr.blocking_factor  = 16
amr.max_grid_size    = 32

castro.max_tagging_radius = 0.6

# WHICH PHYSICS
castro.do_hydro = 1
castro.do_grav  = 1
castro.do_rotation = 1
castro.rotational_period = 100.0

castro.do_scf_initial_model = 1
castro.scf_maximum_density = 1.3519e6
castro.scf_equatorial_radius = 1.0e9
castro.scf_polar_radius = 1.0e9
castro.scf_relax_tol = 1.e-4

# GRAVITY
gravity.gravity_type = PoissonGrav # Full self-gravity with the Poisson equation
gravity.max_multipole_order = 6    # Multipole expansion includes terms up to r**(-max_multipole_order)

# DIAGNOSTICS & VERBOSITY
castro.sum_interval   = 1       # timesteps between computing integrals
amr.data_log          = grid_diag.out

castro.v = 1
amr.v = 1

# CHECKPOINT FILES
amr.checkpoint_files_output = 1
amr.check_file        = chk      # root name of checkpoint file
amr.check_int         = 1        # timesteps between checkpoints

# PLOTFILES
amr.plot_files_output = 1
amr.plot_file         = plt      # root name of plotfile
amr.plot_per          = 1        # timesteps between plotfiles
amr.derive_plot_vars  = ALL

# refinement

amr.refinement_indicators = denerr

amr.refine.denerr.max_level = 14
amr.refine.denerr.value_greater = 1.e6
amr.refine.denerr.field_name = density