Input file for FRIB (LS1+FS1)
Closed this issue · 1 comments
I'm trying to debug main.cpp with FRIB lattice file.
The input lattice file already exist at FLAME/python/uscsi/test/latticeout_IMP_withPV_consolidate.lat .
And executable compiled file 'main' also requires detailed data from TLM_Lite source codes.
So I have executed "main" with 2 arguments.
$ ./main /${path}/latticeout_IMP_withPV_consolidate.lat /${path}/tlm_lite
This returns
Parsing succeeds
Main exception: NCharge
Accordingly, I add
NCharge = [10111.0, 10531.0];
in latticeout_IMP_withPV_consolidate.lat .
And then, I was able to run the "main" file.
But, the calculated result of the beam parameters (eg. BeamRMS.out etc.) is differ greatly from TLM Java code (tlm_lite/example/runtlm).
Do you have any idea? or Is there any wrong command in above?
(Now I'm checking physics variables in the cpp source code also.)
Add information about the result.
After I add
NCharge = [10111.0, 10531.0];
,
$ ./main /${path}/latticeout_IMP_withPV_consolidate.lat /${path}/tlm_lite
returns below
Parsing succeeds
Ion charge states:
0.13866 0.14286
Ion charge amount:
10111.0 10531.0
Barycenter:
-----------
omission
-----------
Beam envelope:
-----------
omission
-----------
InitLong:
IonZ = 0.13866 IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432
InitLong: 0.01832 sec
InitLattice:
IonZ = 0.13866 IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432
IonZ = 0.14286 IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432
*** GetTransitFac: beta out of range 0.0327619
and
tail -2 BeamRMS.out
9.855276000000000e+00 7.023517629255475e+08 2.697435992480285e+05 7.325989069535193e+08 2.812918864926295e+05 3.013622197197078e-01 8.764218641711942e-04
9.855276000000000e+00
The main code stopped near s~10.0 [m] and the x RMS reached to 7.0e8 [mm].