+ fp.lx=atoi(word[2]);
+ fp.ly=atoi(word[3]);
+ fp.lz=atoi(word[4]);
+ fp.lc=atof(word[5]);
+ mdrun->lc=fp.lc;
+ if(!strncmp(word[1],"lc",2)) {
+ if(wcnt==8) {
+ fp.fill_element=atoi(word[6]);
+ fp.fill_brand=atoi(word[7]);
+ }
+ else {
+ fp.fill_element=mdrun->element1;
+ fp.fill_brand=0;
+ }
+ }
+ else {
+ switch(word[6][0]) {
+ case 'i':
+ if(word[6][1]=='r')
+ fp.p_type=PART_INSIDE_R;
+ else
+ fp.p_type=PART_INSIDE_D;
+ break;
+ case 'o':
+ if(word[6][1]=='r')
+ fp.p_type=PART_OUTSIDE_R;
+ else
+ fp.p_type=PART_OUTSIDE_D;
+ break;
+ default:
+ break;
+ }
+ }
+ if((fp.p_type==PART_INSIDE_R)||
+ (fp.p_type==PART_OUTSIDE_R)) {
+ fp.p_vals.r=atof(word[7]);
+ fp.p_vals.p.x=atof(word[8]);
+ fp.p_vals.p.y=atof(word[9]);
+ fp.p_vals.p.z=atof(word[10]);
+ }
+ if((fp.p_type==PART_INSIDE_D)||
+ (fp.p_type==PART_OUTSIDE_D)) {
+ fp.p_vals.p.x=atof(word[7]);
+ fp.p_vals.p.y=atof(word[8]);
+ fp.p_vals.p.z=atof(word[9]);
+ fp.p_vals.d.x=atof(word[10]);
+ fp.p_vals.d.y=atof(word[11]);
+ fp.p_vals.d.z=atof(word[12]);
+ }
+ fp.lattice=mdrun->lattice;
+ add_stage(mdrun,STAGE_FILL,&fp);
+ }
+ else if(!strncmp(word[0],"thermal_init",12)) {
+ add_stage(mdrun,STAGE_THERMAL_INIT,NULL);