# Box of 1000 water molecules.
variable run  string tip3p_aa_02
echo            log
log             ${run}.log

# Neighbour list.
neighbor        2.0 bin
neigh_modify    check yes

# Force field definition.         
boundary p p p
units           real
atom_style      full
pair_style      lj/charmm/coul/long 8.0 10.0
bond_style      harmonic
angle_style     charmm
dihedral_style  charmm
kspace_style    pppm 0.0001

# Atom Definition
read_data "tip3p_box.data"

# Settings
mass 1  15.9994
mass 2   1.0079

pair_coeff  1 1 0.102 3.188
pair_coeff  2 2 0.0 0.0

bond_coeff  1  450.0 0.9572
angle_coeff 1  55.0 104.52 0.0 0.0

# Physical observables.
thermo          1
thermo_style custom step etotal evdwl ecoul elong ke pe temp
timestep        2.
velocity        all create 300 391624 dist gaussian
fix             1 all shake .0001 10 0 m 1.0079 a 1
fix             2 all nvt temp 300 300 200.
run 10 
run 0  # Same value
run 0  # Here it changes!
run 0 
run 0 
run 0