I'm using fix adapt to fill a volume to the correct packing fraction after an insert/pack operation:
compute rad all property/atom radius
variable diam atom c_rad*2
variable size atom ramp(1,4)*v_diam
fix grow all adapt 1 atom diameter v_size scale no
However when I run it, the particles have quadrupled in size on tilmestep 40, not timestep 1000. I think the problem is that it's evaluating the radii at each time step and plugging it back into the 'size' variable, but I'm not sure how to stop it from doing that. Does anyone have any ideas how to do this? I'm pretty sure I had it working a couple years ago, but can't seem to figure out what I did.
Thanks,
Thomas