+ if(!strncmp(word[1],"ins_atoms",9)) {
+ iap.ins_steps=atoi(word[2]);
+ iap.ins_atoms=atoi(word[3]);
+ iap.element=atoi(word[4]);
+ iap.element=atoi(word[4]);
+ iap.brand=atoi(word[5]);
+ for(i=0;i<strlen(word[6]);i++) {
+ switch(word[6][i]) {
+ case 'b':
+ iap.attr|=ATOM_ATTR_VB;
+ break;
+ case 'h':
+ iap.attr|=ATOM_ATTR_HB;
+ break;
+ case 'v':
+ iap.attr|=ATOM_ATTR_VA;
+ break;
+ case 'f':
+ iap.attr|=ATOM_ATTR_FP;
+ break;
+ case '1':
+ iap.attr|=ATOM_ATTR_1BP;
+ break;
+ case '2':
+ iap.attr|=ATOM_ATTR_2BP;
+ break;
+ case '3':
+ iap.attr|=ATOM_ATTR_3BP;
+ break;
+ default:
+ break;
+ }
+ }
+ // only rand mode by now
+ if(word[8][0]=='t')
+ iap.type=INS_TOTAL;
+ else {
+ iap.type=INS_REGION;
+ iap.x0=atof(word[8]);
+ iap.y0=atof(word[9]);
+ iap.z0=atof(word[10]);
+ iap.x1=atof(word[11]);
+ iap.y1=atof(word[12]);
+ iap.z1=atof(word[13]);
+ }
+ add_stage(mdrun,STAGE_INSERT_ATOMS,&iap);
+ }
+ else if(!strncmp(word[1],"continue",8)) {
+ cp.runs=atoi(word[2]);
+ add_stage(mdrun,STAGE_CONTINUE,&cp);
+ }
+ else if(!strncmp(word[1],"anneal",6)) {
+ ap.count=0;
+ ap.runs=atoi(word[2]);
+ ap.dt=atof(word[3]);
+ add_stage(mdrun,STAGE_ANNEAL,&ap);
+ }
+ else {
+ printf("%s unknown stage type: %s\n",
+ ME,word[1]);
+ return -1;
+ }