!! TPS lattice 'd24p79h2' is used for CSR simulation
!! Revolution period= 1.72919627e-6 sec.
!! Nominal total beam current of 400 mA and 600 bunches are assumed
!! Single bunch current(ave.)= 2/3 mA for nominal operation
!! 2/3 mA/bunch== 1.1528e-9 C
!! PJ. Chou, April 3, 2008

tube     : maxamp, x_max=0.034, y_max=0.015, elliptical=1, exponent=2
DDD      : drift, l = 0.125000
DL1      : drift, l = 6.000000
dl1a     : drift, l=2.9
dl1b     : drift, l=2.9
DL2      : drift, l = 0.730000
DL2A     : drift, l = 0.355000
DL2B     : drift, l = 0.125000
DL3      : drift, l = 0.700000
DL4      : drift, l = 0.750000
DL4A     : drift, l = 0.125000
DL4B     : drift, l = 0.375000
DC5      : drift, l = 0.550000
DC6      : drift, l = 0.950000
DC6A     : drift, l = 0.275000
DC6B     : drift, l = 0.425000
DC7      : drift, l = 0.255000
DS1      : drift, l = 3.500000
DS2      : drift, l = 0.680000
DS2A     : drift, l = 0.305000
DS2B     : drift, l = 0.125000
DS3      : drift, l = 0.800000
DS3A     : drift, l = 0.425000
DS3B     : drift, l = 0.125000
DS4      : drift, l = 0.320000
      
QS5      : kquad, l = 0.300000, k1 =   1.601058, n_kicks=20
QS4      : kquad, l = 0.300000, k1 =  -0.975030, n_kicks=20
QS3      : kquad, l = 0.300000, k1 =  -1.353903, n_kicks=20
QS2      : kquad, l = 0.600000, k1 =   1.545859, n_kicks=20
QS1      : kquad, l = 0.300000, k1 =  -1.696006, n_kicks=20
QL3      : kquad, l = 0.300000, k1 =  -1.237476, n_kicks=20
QL2      : kquad, l = 0.600000, k1 =   1.297832, n_kicks=20
QL1      : kquad, l = 0.300000, k1 =  -1.078493, n_kicks=20

!!MB: Use edge_order=1 to ensure symplecticity
BEND: csbend, l=1.1, angle=0.1308996939, k1=0.0, e1=0.06544984695, e2=&
0.06544984695, integration_order=4, n_kicks=20, edge_order=1
bend8: csbend, l=1.1, angle=0.1308996939, k1=0.0, e1=0.06544984695, e2=&
0.06544984695, integration_order=4, n_kicks=20, edge_order=1

!Sextupole
!!MB: Usually don't need many kicks for sextupoles
S1       : ksext,l=0.125, k2 = 12.351504, n_kicks=2  
S2       : ksext,l=0.125, k2 = -28.418568, n_kicks=2 
S3       : ksext,l=0.125, k2 = 38.783568, n_kicks=2  
S4       : ksext,l=0.125, k2 = -42.464976, n_kicks=2 
S5       : ksext,l=0.125, k2 = 32.451576, n_kicks=2  
S6       : ksext,l=0.125, k2 = -40.431864, n_kicks=2 
SD       : ksext,l=0.125, k2 = -28.2858, n_kicks=2   
SF       : ksext,l=0.125, k2 = 38.693352, n_kicks=2  

BPM      : MONI, L=0.0
KHV      : KICKER, HKICK=0.000, VKICK = 0.000
!!MB: Use more exact frequency, set phase to 180 deg since we are not
!!MB: simulating radiation
rf       : rfca,l=0.2,volt=3.5e6,phase=180,freq="864 1.72919627e-6  /"

LLONG    :LINE=(DL1,BPM, QL1, DL2A, S1,KHV,S1, DL2B, QL2, DL3,&
               BPM,QL3, DL4A, S2,KHV,S2, DL4B)
LLONGrf  :LINE=(DL1A,rf,DL1B,BPM, QL1, DL2A, S1,KHV,S1, DL2B, QL2, DL3,&
               BPM,QL3, DL4A, S2,KHV,S2, DL4B)
RLONG    :LINE=(DS1,BPM, QS1, DS2A, S3,KHV,S3, DS2B, QS2, DS3A,&
               BPM,S4,KHV,S4, DS3B, QS3,DS4)
HLACH    :LINE=(DC5, QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7, &
               SF,KHV,SF,BPM)
SHORT    :LINE=(DS1,BPM, QS1, DS2A, S5,KHV,S5, DS2B, QS2, DS3A,&
               BPM,S6,KHV,S6, DS3B, QS3, DS4)
HSACH    :LINE=(DC5, QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7, &
               SF,KHV,SF,BPM)
HSACH1   :LINE=(DC5,QS4,BPM, DC6A, SD,KHV,SD, DC6B, QS5, DC7)

SCELL    :LINE=(SHORT, BEND, HSACH, -HSACH1, BEND, -SHORT)
LCELL    :LINE=(LLONG, BEND, HLACH, -HSACH1, BEND, -RLONG)
LCELLrf  :LINE=(LLONGrf, BEND, HLACH, -HSACH1, BEND, -RLONG)
ILCELL   :LINE=(RLONG, BEND, HSACH1,-HLACH, BEND8, -LLONG)

SUPERA   :LINE=(LCELL, SCELL, SCELL, ILCELL)
SUPERArf :LINE=(LCELLrf, SCELL, SCELL, ILCELL)                           
RING     :LINE=(tube, 5*SUPERA, superarf)  
return