processors 1 1 1 

units lj
atom_style bond

special_bonds lj/coul 0 1 1

read_data polymer_10c10m.dat

neighbor 0.4 bin
neigh_modify every 1 delay 1
comm_modify vel yes

bond_style fene
bond_coeff * 30.0 1.5 1.0 1.0

dump 	1 all custom 60000 *_Monomer_Coordinates id mol type x y z ix iy iz
compute b1 all property/local btype batom1 batom2
compute b2 all bond/local dist engpot force
dump    2 all local 60000 *_Bond_Analysis index c_b1[*] c_b2[*]

timestep 0.01
thermo 1000
thermo_style custom step temp pe etotal epair ebond density press pxx pyy pzz lx ly lz

pair_style dpd 1.0 1.0 122347
pair_coeff * * 25 4.5 1.0

velocity all create 1.0 17786140

fix 1 all nve
run 50000

pair_coeff * * 50.0 4.5 1.0
velocity all create 1.0 15086120
run 50
pair_coeff * * 100.0 4.5 1.0
velocity all create 1.0 15786120 
run 50
pair_coeff * * 150.0 4.5 1.0
velocity all create 1.0 15486120
run 50
pair_coeff * * 200.0 4.5 1.0
velocity all create 1.0 17986120
run 100
pair_coeff * * 250.0 4.5 1.0
velocity all create 1.0 15006120
run 100
pair_coeff * * 500.0 4.5 1.0
velocity all create 1.0 15087720
run 100
pair_coeff * * 1000.0 4.5 1.0
velocity all create 1.0 15086189
run 100

unfix 1
fix 2 all npt temp 1.0 1.0 2.0 iso 0.0 0.0 10.0

pair_style lj/cut 2.5
pair_coeff * * 1.0 1.0 2.5
run 6000000
write_data polymer_10x10.data
write_restart my_restart_lj