Re: [ESPResSo-users] particle does not exist

[Xikai Jiang]

Thank you, Stefan. I check the content of confout, it shows correct particle information as

{particles {id pos v}

 { 0 1 2 3 0 0 0 }

 ...

}

Just find that it's a simple mistake I made that I'm accessing a particle that was already

deleted.. So it's important to detect which particles are deleted.

When I retrieve the deleted particle position, Espresso will return a value of "na", to

find which atom is deleted, I find that one way is to compare it with a number such as:

if { [part $id print pos] >=100 } {

   puts "it's deleted particle"

}

> Subject: Re: [ESPResSo-users] particle does not exist
> From: address@hidden
> Date: Tue, 17 Jun 2014 19:00:49 +0200
> CC: address@hidden
> To: address@hidden
>
> Dear Xikai,
> please check the content of confout for particle 4633. There must be a section like
> Are you sure that there is something like { particle 4633 pos 1 2 3 } in there?
> Just after reading the configuration you should then be able to e.g. print its position,
> puts [ part 4633 print pos ]
> and also to delete it. It is very likely that this is a simple error in your particle indexing or something similar.
> Cheers and good luck
> Stefan
>
> On Jun 17, 2014, at 6:54 PM, Xikai Jiang <address@hidden> wrote:
>
> > Dear all:
> >
> > I've got a question about deleting particles from a simulation continued
> > from previous one.
> >
> > In the older simulation, I wrote blockfiles to store particle information such as pid,
> > positions, velocities.
> >
> > When the old simulation finished, I continued to run it by reading saved blockfiles by
> > the following command:
> >
> > set last [open "confout" "r"]
> > while { [ blockfile $last read auto] != "eof" } { }
> > close $last
> >
> > The simulation can continue running without problem, however, when I attempt to
> > delete a particle, it said "particle 4633 does not exist and cannot be removed".
> > I'm sure that particle 4633 exists in the saved configuration file.
> >
> > Does anyone encounter this problem before?
> >
> > Thanks for your time!
> >
> > XIkai
>