Added a funciton that calculates everything needed for spherical harmonics

Samo Penic

2018-07-02 f977e8fe32b6a1a2ae85662ab35fd7c9f647cc3e

 src/io.c

@@ -980,8 +980,10 @@

   if(vesicle->tape->stretchswitch==1){
      fprintf(fh,"<DataArray type=\"Float64\" Name=\"stretching_energy\" format=\"ascii\">");
      for(i=0;i<blist->n+monono*polyno+filno*(fonono-1)
;i++){
      for(i=0;i<blist->n;i++){
         fprintf(fh, "0.0 ");
      }
      for(i=0;i<monono*polyno+filno*(fonono-1);i++){
         fprintf(fh,"0.0 ");
      }
      for(i=0;i<vesicle->tlist->n;i++){