来自密俊霞的邮件

Dear LAMMPS developers:
thermo 100
variable pg3z5 equal {pg3z} variable a loop 100000 label loop variable pg3z2 equal 0.1*sin(0.000628*(v_a+1)) run 1 variable pg3z4 equal {pg3z5}+{pg3z2}-{pg3z}
displace_atoms g3 move 0 0 ${pg3z4} units box
run 0
next a
jump in.Layer1c loop

I’m trying to output thermo info every 100 steps,but it thermo info every 1 step.What should I do?
Best wishes
Sincerely yours
Mi junxia

Dear LAMMPS developers:
thermo 100
variable pg3z5 equal \{pg3z\} variable a loop 100000 label loop variable pg3z2 equal 0\.1\*sin\(0\.000628\*\(v\_a\+1\)\) run 1 variable pg3z4 equal {pg3z5}+\{pg3z2\}\-{pg3z}
displace_atoms g3 move 0 0 ${pg3z4} units box
run 0
next a
jump in.Layer1c loop

I'm trying to output thermo info every 100 steps,but it thermo info every 1
step.What should I do?

don't use "run 1" and use fix move instead of displace_atoms.

axel.

Dear LAMMPS developers:
I try to add a vibration in Z direction.

region rt block INF INF INF INF 3.4 INF units box
group gt region rt
velocity gt create 300 4928439 rot yes dist gaussian units box
fix 3 gt langevin 300 300 0.1 699483 zero yes
timestep 0.001
thermo 100
run 200
variable v0 equal 1.0-v_vgtx
velocity gt set {v0} 0.0 0.0 sum yes units box variable a loop 1000 label loop variable pgtz2 equal 0.1*sin(0.000628*(v_a+1)) run 1 variable pgtz4 equal {pgtz5}+{pgtz2}-{pgtz}
displace_atoms gt move 0 0 ${pgtz4} units box
run 0
next a
jump in.Layer1c loop

It works well but output thermo info every 1 step,I just want it to output thermo info every 100 steps.

I try to use fix move command instead of displace_atoms

region rt block INF INF INF INF 3.4 INF units box
group gt region rt
velocity gt create 300 4928439 rot yes dist gaussian units box
fix 3 gt langevin 300 300 0.1 699483 zero yes
timestep 0.001
thermo 100
run 200
variable v0 equal 1.0-v_vgtx
velocity gt set ${v0} 0.0 0.0 sum yes units box
fix vib gt move wiggle NULL NULL 0.1 10 units box
run 1000

after 100 steps,it warn:One or more atoms are time integrated more than once
after 300 steps,it report error:
ERROR: Lost atoms: original 13776 current 13752

Do you have any idea how to correct this behaviour?Thank you very much.

Best wishes
Sincerely yours
Mi junxia