Changeset 11 in TRACY3 for trunk/tracy/tools

Timestamp:

Oct 21, 2013, 10:40:43 AM (11 years ago)

Author:

zhangj

Message:

 

File:

• trunk/tracy/tools/soltracy.cc (modified) (7 diffs)

trunk/tracy/tools/soltracy.cc

@@ -41 +41 @@
   globval.H_exact = false;
 
+  
+  //output files
+  char fic_twiss[S_SIZE + 4] = "";   //twiss file
+  char fic_summary[S_SIZE + 4]="";   //summary file
+  
+  
   /* parameters to read the user input script .prm */
   long i=0L; //initialize the for loop to read command string
@@ -99 +105 @@
   else
     globval.ge = ElemIndex(ge_name);
-  
+ 
+ 
 //  globval.g = ElemIndex("g");  /* get family index of  girder*/
-   
-  /* print the summary of the element in lattice */  
+ strcpy(fic_twiss,lat_file);
+  strcpy(fic_summary,lat_file);
+  strcat(fic_twiss,".twi"); 
+  strcat(fic_summary,".sum");
+  
+    //generate the twiss and summary files with proper appendix  
+  /* print the summary of the element in lattice to the screen and an external file*/  
   printglob();   
+  printglob2file(fic_summary);
+ //print the twiss file
+   /* print out lattice functions, with all the optical information for the lattice with design values */
+  printlatt(fic_twiss);  
+  
   /************************************************************************
     print files, very important file for debug
@@ -246 +263 @@
    }
    
+   else if(strcmp(CommandStr,"PrintTrackElemFlag") == 0) {
+      cout << "\n";
+      cout << "print the tracked coordinates to file: "
+           << UserCommandFlag[i]._PrintTrackElem_track_file << "\n";
+        
+      PrintTrackElem(UserCommandFlag[i]._PrintTrackElem_track_file,
+                 UserCommandFlag[i]._PrintTrackElem_x, UserCommandFlag[i]._PrintTrackElem_px, 
+                 UserCommandFlag[i]._PrintTrackElem_y, UserCommandFlag[i]._PrintTrackElem_py, 
+                 UserCommandFlag[i]._PrintTrackElem_delta,UserCommandFlag[i]._PrintTrackElem_ctau,
+                 UserCommandFlag[i]._PrintTrackElem_nelem1,UserCommandFlag[i]._PrintTrackElem_nelem2);  
+   }
     //print the girder
   // else if(strcmp(CommandStr,"PrintGirderFlag") == 0) {
@@ -462 +490 @@
  //        UserCommandFlag[i]._FmapFlag_delta, 
  //        UserCommandFlag[i]._FmapFlag_diffusion,
+ //        UserCommandFlag[i]._FmapFlag_printloss,
  //        numprocs,myid);
 
@@ -514 +543 @@
            UserCommandFlag[i]._FmapFlag_ymax, 
            UserCommandFlag[i]._FmapFlag_delta, 
-           UserCommandFlag[i]._FmapFlag_diffusion);
+           UserCommandFlag[i]._FmapFlag_diffusion,
+           UserCommandFlag[i]._FmapFlag_printloss);
       //#endif
   }
@@ -542 +572 @@
  //          UserCommandFlag[i]._FmapdpFlag_z,
  //          UserCommandFlag[i]._FmapdpFlag_diffusion,
+ //            UserCommandFlag[i]._FmapdpFlag_printloss,
  //          numprocs,myid);
 
@@ -594 +625 @@
               UserCommandFlag[i]._FmapdpFlag_emax, 
               UserCommandFlag[i]._FmapdpFlag_z,
-              UserCommandFlag[i]._FmapdpFlag_diffusion);
+              UserCommandFlag[i]._FmapdpFlag_diffusion,
+               UserCommandFlag[i]._FmapdpFlag_printloss);
        //  #endif
   }