# Settings

units lj
dimension 3
atom_style molecular
pair_style lj/cut 1.12246
bond_style fene
angle_style cosine
special_bonds lj 1 1 1

read_data test.dat

mass * 1.0

pair_coeff * * 1.0 1.0
bond_coeff 1 30.0 1.50 0.0 0.0
angle_coeff 1 0.0


# change neighbor list parameters to avoid dangerous builds
neighbor        0.9 bin
neigh_modify    delay 0 every 1 check yes 

timestep 0.01

velocity all create 1.0 1234

compute c1 all temp

compute c2 all gyration

compute c3 all com

thermo_style custom step c_c1 c_c2 c_c3[*]

fix f2 all lb/fluid 5 1 2.80 1.0 dx 1.0 noise 1.0 1234 &
 bodyforce 0.0 0.0 0.0 setGamma 20.0 D3Q19 trilinear

 fix f3 all lb/viscous

 fix f6 all nve

thermo 10

run 10

quit