# Program parameters
sim_type = MD
backend = CPU
use_edge = true
seed = 1
cells_auto_optimisation = false
max_density_multiplier = 10
# verlet_skin = 0.5
verlet_skin = 0.1
interaction_type = DNA2_nomesh
dh_half_charged_ends = 1
list_type = no

# Simulation parameters
steps = 20000
newtonian_steps = 103
diff_coeff = 2.5
thermostat = john
dt = 0.003
ensemble = nvt
# T = 20C
T = 296.15K
salt_concentration = 0.5

use_average_seq = true

# Input/output
topology = sys.top
conf_file = init.conf
lastconf_file = last_conf.dat
refresh_vel = 1
log_file = sim.log
restart_step_counter = 0
energy_file = energy.dat
print_conf_interval = 500
print_conf_ppc = 51
print_energy_every = 500
trajectory_file = output.dat
time_scale = linear


# oxDNA observables
data_output_1 = {
  name = split_energy.dat
  print_every = 500
  col_1 = {
          type = step
          units = MD
  }
  col_2 = {
          type = potential_energy
          # type = pair_energy
          split = true
  }
}
data_output_2 = {
  name = potential.dat
  print_every = 500
  col_1 = {
          type = step
          units = MD
  }
  col_2 = {
          type = potential_energy
          split = false
  }
}
data_output_3 = {
  name = pair.dat
  print_every = 500
  col_1 = {
          type = step
          units = MD
  }
  col_2 = {
          type = pair_energy
          split = false
  }
}