Recently, i simulate a case needing pull some atoms back to the domain. I tried boundary command of "f","s","m","p" , but they are not right. Because the domain lower boundary is fixed and upper boundary is open. Some atoms escape from upper bounday, and needed to pull back after each step. It seems the "shrink wrap" is not right. Thanks!