dimension 2
boundary p p p
units lj

# Define box and place one particle
region total block -10 10 -10 10 -10 10
create_box 1 total
create_atoms 1 single 0 0 0
mass 1 1.0

# Set up Langevin dynamics in 2D
fix thermo all langevin ${T} ${T} 1 9876
fix keepinplace all enforce2d
fix timeintegrator all nve

# Run
timestep 0.01
run 1000 
#dump hdump all h5md 100 trajectory.h5 position image velocity
run 20000