Dear LAMMPS authors, I write to report a bug in the file velocity.cpp line 451 on the current distribution (29Apr10):
double v_lo,v_hi;
if (v_dim == 0) {
v_lo = xscale*atof(arg[1]);
v_hi = xscale*atof(arg[2]);
} else if (v_dim == 1) {
v_lo = yscale*atof(arg[1]);
v_hi = yscale*atof(arg[2]);
} else if (v_dim == 0) {
v_lo = zscale*atof(arg[1]);
v_hi = zscale*atof(arg[2]);
}
It should be:
double v_lo,v_hi;
if (v_dim == 0) {
v_lo = xscale*atof(arg[1]);
v_hi = xscale*atof(arg[2]);
} else if (v_dim == 1) {
v_lo = yscale*atof(arg[1]);
v_hi = yscale*atof(arg[2]);
} else if (v_dim == 2) {
v_lo = zscale*atof(arg[1]);
v_hi = zscale*atof(arg[2]);
}
It was preventing from use ramped velocity distribution in the z direction.
Best regards,
Fabio Pavia
Brown University
Engineering Division, Mechanics of Solids