ZGOUBI USERS' GUIDE - HEP

4.5 Output Procedures 123¨ ¢ < ¤qq¥wFAISCEAU, FAISCNL, FAISTORE: Print/Store particle coordinatesFAISCEAU can be introduced anywhere in a structure. It produces a print of initial and actual coordinates of the particlesat the location where it stands, together with their tagging indices and letters, following the same format as for FAISCNL(except for SORT(I) which is not printed) .FAISCNL has a similar effect, except that the information is stored in a dedicated file FNAME (standard name is FNAME= ‘zgoubi.fai’ for post-processing with zpop). This file may further on be read by means of OBJET, option KOBJ = 3, orused for other purposes such as graphics (see Part D of the Guide). The data written to that file are formatted and orderedaccording to the FORTRAN sequence belowOPEN (UNIT = NL, FILE = FNAME, STATUS = ‘NEW’)DO 1 I=1, IMAXWRITE(NL,110) LET(I),IEX(I),(FO(J,I),J=1,7),(F(J,I),J=1,7),KineticE,> I(I),IREP(I),SORT(I),Mass,Charge,G-Factor,com-Life-time,unused,RET(I), DPR(I),> BORO,IPASS,KLEY,LBL1,LBL2,NOEL110 FORMAT(1X,A1,I2,1P,7E16.8,2 /,3E24.16,3 /,4E24.16,E16.8,4 /,2I6,8E16.8,5 /,E16.8, I6, A8, 2A10, I5)The meaning of main data is the following (see the keyword OBJET)¢ : one-character string, for tagging particle number ¤¤, ¤ , ¤1|¢¤ : flag, particle number, index! w G §: coordinates O , , , , ¥ and path length at the origin of the structure< w G "$¤: idem, at the current position·=! 1 "$¤: path length at which the particle has possibly been stopped¤¤ , O ¥1|¢(see CHAMBR or COLLIMA)DPR = momentum dispersion (MeV/c) (see CAVITE)IPASS: turn number (see REBELOTE)etc. :¤ : synchrotron phase space coordinates; RET =phase (radian),1FAISTORE has an effect similar to FAISCNL, with two more features. On the first data line, FNAME may be followedby a series of up to ¨¦´¡W¢¨ 10 ’s proper to the elements of the data file at the exit of which the print should occur; if there isno label, the print occurs by default at the location of FAISTORE; if there are labels the print occurs right downstream ofall optical elements wearing those labels (and no longer at the FAISTORE location). The next data line gives a parameter. For instance the data listqFAISTOREzgoubi.fai HPCKUP VPCKUP12¤¦¥ : printing will occur every ¤¥ other pass, if using REBELOTE with NPASS Ž ¤¦¥will result in output prints into zgoubi.fai, every 12 other pass, each time elements of the zgoubi.dat data list labeled eitherHPCKUP or VPCKUP are encountered.NoteBinary storage can be obtained from FAISCNL and FAISTORE. This for the sake of compactness and access speed, forinstance in case voluminous amounts of data would have to be manipulated.This is achieved by giving the storage file a name of the form b FNAME or B FNAME (e.g., ‘b zgoubi.fai’). TheFORTRAN WRITE list is the same as in the FORMATTED case above.This is compatible with the READ statements in zpop that will recognize binary storage from that very radical’b ’ or ’B ’.