# Point dipoles in a 2d box

echo                screen
units               lj
neighbor            0.3 bin
atom_style          hybrid molecular dipole
dimension           2
boundary            p p p

read_data           data.test
#read_restart        equ.1

pair_style          soft 1.12246 
pair_coeff          * * 5.0


reset_timestep      0

neigh_modify        delay 0 exclude molecule all

velocity            all create 1.0 14001

#fix                  11 all rigid molecule langevin 1.2 1.2 0.01 11145514
fix                  11 all rigid/nve molecule langevin 1.2 1.2 0.01 11145514
fix                  12 all  momentum 1000 linear 1 1 1 angular
fix                  13 all enforce2d 

compute              1 all erotate/rigid 11
compute              2 all ke/rigid 11

dump                 1 all custom 1000 test.lammpstrj id type x y
timestep            0.001
thermo              100
thermo_style        custom step  f_11 c_1 c_2 
restart             1000 equ.1 equ.2 

run                 10000