Use nucleus displaced centre in vertex move

Samo Penic

2016-05-31 5c3143e7b7d0979d6eed07b6dfbe9f214159bdbb

 src/main.c

@@ -42,10 +42,24 @@

    if(command_line_args.dump_from_vtk[0]!=0){
      ts_fprintf(stdout,"************************************************\n");
      ts_fprintf(stdout,"***** Dumping vesicle from VTK points list *****\n");
      ts_fprintf(stdout,"**** Restoring vesicle from VTK points list ****\n");
      ts_fprintf(stdout,"************************************************\n\n");
      vesicle = parseDump(command_line_args.dump_from_vtk);
      tape = vesicle->tape;
      int arguments_no;
      FILE *fd=fopen(".status","r");
      if(fd!=NULL){
         arguments_no=fscanf(fd,"%u", &start_iteration);
         if(arguments_no==0){
            ts_fprintf(stdout,"No information of start iteration in .status file");
            }
         fclose(fd);
         start_iteration++;
         printf("nucleus coords: %.17e %.17e %.17e\n",vesicle->nucleus_center[0], vesicle->nucleus_center[1], vesicle->nucleus_center[2]);
      }
      else
         ts_fprintf(stdout,"No .status file. The iteration count will start from 0");
/* Here you should read new tape file, reassign some values in vertex from the tape and assign read tape to vesicle->tape */
//        tape=parsetape(command_line_args.tape_fullfilename);
  //      vesicle=vtk2vesicle(command_line_args.dump_from_vtk,tape);
    }