diff --git a/tracy/tools/demo/demo1/Input_demo1.prm b/tracy/tools/demo/demo1/Input_demo1.prm
new file mode 100644
index 0000000000000000000000000000000000000000..d3d16dba3aa8505c09a3b7cc1853c1adc7d1a8a5
--- /dev/null
+++ b/tracy/tools/demo/demo1/Input_demo1.prm
@@ -0,0 +1,73 @@
+#*******************************************************************
+# 15 March 2010
+# File for checking TRacy is running well
+#
+# input script to set the file and bool flag
+#
+# comments:
+# 1) The rules to define commands are described in the tracy mannual.
+#
+# 2) Test verison for the development of Tracy 3, always the nearest
+# input file verison.
+#
+# 12/2011 Jianfeng Zhang
+#
+#******************************************************************
+
+#================================================
+# mandatory parts
+#================================================
+
+# relative directory of the files in the script
+ in_dir ../
+
+# lattice file must be wihout .lat extension
+ lat_file soleil_120cor
+
+# H-corrector name in lattice file
+ h_corr HCM
+# V-corrector name in lattice file
+ v_corr VCM
+# girder start symbol
+ gs GS
+# girder end symbol
+ ge GE
+# BPM name in lattice file
+ bpm_name BPM
+# skew quadrupole
+ qt QT
+
+#======================================
+# GLOBAL PARAMETER PART
+#======================================
+
+# Activate quadrupole Fringefield
+ QuadFringeOnFlag
+
+# VACCUUM CHAMBER PART
+# read vacuum chamber from chamber_file
+ ReadChamberFlag example_Chamber.dat
+
+# set RF voltage, if false, then use the definition value
+# of RF cavity in the lattice
+ RFvoltageFlag 3000000
+
+#=================================================
+## general command (NOT machine dependent)
+#=================================================
+
+# print the twiss paramters to a file
+ PrintTwissFlag twisstest.out
+
+# tune obtained by tracking
+ TuneTracFlag
+
+# chromaticity obtained by tracking
+ ChromTracFlag
+
+# fit tune for the lattice with full quadrupoles
+ FitTuneFlag q7 q9 18.202 10.317
+
+# fit chromaticity
+ FitChromFlag sx9 sx10 2.0 2.6
+
diff --git a/tracy/tools/demo/demo2/Input_demo2.prm b/tracy/tools/demo/demo2/Input_demo2.prm
new file mode 100644
index 0000000000000000000000000000000000000000..d7a8aa02c741d3938cb7bf805504ec424ec12f73
--- /dev/null
+++ b/tracy/tools/demo/demo2/Input_demo2.prm
@@ -0,0 +1,75 @@
+#*******************************************************************
+# 15 March 2010
+# File for checking TRacy is running well
+#
+# input script to set the file and bool flag
+#
+# comments:
+# 1) The rules to define commands are described in the tracy mannual.
+#
+# 2) Test verison for the development of Tracy 3, always the nearest
+# input file verison.
+#
+# 12/2011 Jianfeng Zhang
+#
+#******************************************************************
+
+#================================================
+# mandatory parts
+#================================================
+
+# relative directory of the files in the script
+ in_dir ../
+
+# lattice file must be wihout .lat extension
+ lat_file soleil_120cor
+
+# H-corrector name in lattice file
+ h_corr HCM
+# V-corrector name in lattice file
+ v_corr VCM
+# girder start symbol
+ gs GS
+# girder end symbol
+ ge GE
+# BPM name in lattice file
+ bpm_name BPM
+# skew quadrupole
+ qt QT
+
+#======================================
+# GLOBAL PARAMETER PART
+#======================================
+
+# Activate quadrupole Fringefield
+ QuadFringeOnFlag
+
+# VACCUUM CHAMBER PART
+# read vacuum chamber from chamber_file
+ ReadChamberFlag example_Chamber.dat
+
+# set RF voltage, if false, then use the definition value
+# of RF cavity in the lattice
+ RFvoltageFlag 3000000
+
+#=================================================
+## general command (NOT machine dependent)
+#=================================================
+
+# print the twiss paramters to a file
+ PrintTwissFlag twisstest.out
+
+# tune obtained by tracking
+ TuneTracFlag
+
+# chromaticity obtained by tracking
+ ChromTracFlag
+
+# fit tune for the lattice with full quadrupoles
+ FitTuneFlag q7 q9 18.202 10.317
+
+# fit chromaticity
+ FitChromFlag sx9 sx10 2.0 2.6
+
+# frequency map analysis for on momentum particle
+ FmapFlag fmaptest.out 3 2 516 0.025 0.005 0.0 true
diff --git a/tracy/tools/demo/example_Chamber.dat b/tracy/tools/demo/example_Chamber.dat
new file mode 100644
index 0000000000000000000000000000000000000000..1d5463343ffe61fa618d2b5bdb034901061fc793
--- /dev/null
+++ b/tracy/tools/demo/example_Chamber.dat
@@ -0,0 +1,26 @@
+#**************************************************************************
+# Script to set the vacuum chamber
+#
+# Comments:
+# 1) Line start with "#" is comment
+# 2) The format of the vacuum chamber definition is
+# start element, end element, minimum x, maximum x,
+# minimum y, maximum y.
+# 3) The first line is to define the global
+# vacuum chamber limit around the ring, and
+# the key words should be "Start","All".
+# 4) To set the vacuum chamber, it is needed to add two markers
+# in the lattice, such as MK1 and MK2, MK1 is before the first element,
+# and MK2 is after the end element for the vacuum chamber region.
+# 5) the units are meters for minimum x, maximum x, minimum y, maximum y
+#****************************************************************************
+#
+# MK1 MK2 dxmin dxmax dymin dymax (Apertures in meters)
+Start All -35e-3 35e-3 -12.5e-3 12.5e-3
+#sdm1 esdm -21e-3 21e-3 -5e-3 5e-3
+debut ehu600 -35e-3 35e-3 -7e-3 7e-3
+ssep esep -20e-3 35e-3 -7e-3 7e-3
+ssdm esdm -21e-3 21e-3 -5e-3 5e-3
+ssdac esdac -35e-3 25e-3 -2.5e-3 2.5e-3
+
+
diff --git a/tracy/tools/demo/soleil_120cor.lat b/tracy/tools/demo/soleil_120cor.lat
new file mode 100644
index 0000000000000000000000000000000000000000..c77787675251d00117261ad39d3f57913b26186e
--- /dev/null
+++ b/tracy/tools/demo/soleil_120cor.lat
@@ -0,0 +1,320 @@
+{*******************************}
+{* *}
+{* SOLEIL Storage Ring *}
+{* The IDEAL LATTICE updated *}
+{* *}
+{* *}
+{* 06/06/09 Laurent Nadolski *}
+{* lattice for orbit correction*}
+{* 56 girder *}
+{* Do not split element *}
+{* no quad FF in Tracy 3.2 *}
+{* *}
+{*******************************}
+
+define lattice;
+
+intmeth= 4;
+driftkick = 2;
+
+{***** System parameters *****}
+Energy= 2.739; { GeV }
+dP = 1.0d-10;
+CODeps= 1.0d-15;
+
+{***** drift space *****}
+
+ SD1 : drift, L= 6.000000-0.0201;
+ SD1a : drift, L= 0.900000;
+ SD1b : drift, L= 4.300000;
+ SD1c : drift, L= 3.700000;
+ SD1d : drift, L= 0.800000-0.0200-0.291000;
+ SD1e : drift, L= 5.658900;
+ SEP : drift, L= 0.200000;
+ SD2 : drift, L= 0.390000-0.08-0.0201;
+ SD3 : drift, L= 0.200000-0.08-0.0181;
+ SD4 : drift, L= 0.390000-0.08-0.0181;
+ SD41 : drift, L= 0.252100-0.08;
+ SD42 : drift, L= 0.205000-0.08;
+ SD43 : drift, L= 0.205100-0.08;
+ SD5 : drift, L= 0.200000-0.08-0.0201;
+ SD6 : drift, L= 0.790000-0.08;
+ SD7 : drift, L= 0.440000-0.0201;
+ SD8 : drift, L= 0.200000-0.08-0.0201;
+ SD9 : drift, L= 0.476340-0.08-0.0201;
+ SD9a : drift, L= 0.204200;
+ SD91 : drift, L= 0.251240;
+ SD92 : drift, L= 0.204300;
+ SD93 : drift, L= 0.251300;
+ SD10 : drift, L= 0.470000-0.08-0.0201;
+ SD10a: drift, L= 0.172300;
+ SD12 : drift, L= 0.470000-0.08-0.0201;
+ SD13 : drift, L= 3.482550-0.0201-0.320900;
+ SD13a: drift, L= 3.141452;
+ SD13i: drift, L= 0.200000;
+ SD13j: drift, L= 0.150000;
+ SD13k: drift, L= 0.732550-0.0202-0.3209;
+ SD14 : drift, L= 0.590000-0.0181-0.0201;
+ SD14a: drift, L= 0.385000;
+ SDAC : drift, L= 1.904680;
+ SDAC1: drift, L= 1.48428;
+ SDAC3: drift, L= 1.904680-0.08;
+ SDAC5: drift, L= 1.004680-0.160500;
+ SDAC4: drift, L= 1.004680-0.3404-0.08;
+ SDAC2: drift, L= 0.200000;
+ SDACm: drift, L= 1.904680-1.80/2;
+ SDB1i : drift, L= 0.291000;
+ SDB1 : drift, L= 0.320900;
+ SDB2 : drift, L= 0.166800;
+ SDB3 : drift, L= 0.252000-0.08;
+ SDB4 : drift, L= 0.277600-0.08;
+ SDB5 : drift, L= 0.205000-0.08;
+ SDB6 : drift, L= 0.119800;
+ SDB7 : drift, L= 0.166900;
+ SDB8 : drift, L= 0.252000-0.08;
+ SDB9 : drift, L= 0.119800;
+ SDB10: drift, L= 0.166900;
+ SDB11: drift, L= 0.251900-0.08;
+ SDB12: drift, L= 0.204900-0.08;
+ SDB13: drift, L= 0.119800;
+ SDB14: drift, L= 0.166800;
+ SDB15: drift, L= 0.252000-0.08;
+ SDB17: drift, L= 0.205000-0.08;
+ SDB18: drift, L= 0.119900;
+ SDC12: drift, L= 0.320900;
+ SDC13: drift, L= 0.320350;
+ SDC44: drift, L= 0.420400-0.08;
+ SDC45: drift, L= 0.421000-0.08;
+ SDC2 : drift, L= 0.079000-0.08;
+ SD89 : drift, L= 0.340400;
+ SD88 : drift, L= 0.341000;
+ SD87 : drift, L= 0.3209;
+ SD86 : drift, L= 0.3204;
+ SD85 : drift, L= 0.320350;
+ SDC30: drift, L= 0.320900;
+ SD141: drift, L= 0.431900;
+
+{** Quadrupole **}
+Nq=8/2; {Number of slices}
+LQC=0.3602;
+LQL=0.4962;
+
+Q1 : quadrupole, L=LQC,K= -1.073038, method=intmeth,N=Nq;
+Q2 : quadrupole, L=LQL,K= 1.603455, method=intmeth,N=Nq;
+Q3 : quadrupole, L=LQC,K= -0.649160, method=intmeth,N=Nq;
+Q4 : quadrupole, L=LQC,K= -1.259864, method=intmeth,N=Nq;
+Q5 : quadrupole, L=LQC,K= 1.696440, method=intmeth,N=Nq;
+Q6 : quadrupole, L=LQC,K= -1.007805, method=intmeth,N=Nq;
+Q7 : quadrupole, L=LQL,K= 2.084913, method=intmeth,N=Nq;
+Q8 : quadrupole, L=LQC,K= -1.551946, method=intmeth,N=Nq;
+Q9 : quadrupole, L=LQC,K= -1.510915, method=intmeth,N=Nq;
+Q10: quadrupole, L=LQC,K= 1.756880, method=intmeth,N=Nq;
+
+{** bending **}
+beta_gap=37e-3;
+tracy_gap=beta_gap*2*0.724;
+
+BEND0 : bending, L= 1.05243, T= 11.25, T1=5.5906, T2=5.67658, K=0.00204, N=4,
+ method=intmeth,gap=tracy_gap;
+
+{** Sextupole **}
+NqSx=1; {Number of slices}
+Lsext = 0.16;
+coef=1.0/Lsext;
+method4sextu = 4;
+sextfringe = 0;
+
+SX1 : sextupole, L=Lsext, K = 1.719190*coef, method=method4sextu, N = NqSx;
+SX2 : sextupole, L=Lsext, K = -4.104560*coef, method=method4sextu, N = NqSx;
+SX3 : sextupole, L=Lsext, K = -2.169380*coef, method=method4sextu, N = NqSx;
+SX4 : sextupole, L=Lsext, K = 3.604650*coef, method=method4sextu, N = NqSx;
+SX5 : sextupole, L=Lsext, K = -3.698210*coef, method=method4sextu, N = NqSx;
+SX6 : sextupole, L=Lsext, K = 3.246670*coef, method=method4sextu, N = NqSx;
+SX7 : sextupole, L=Lsext, K = -5.003520*coef, method=method4sextu, N = NqSx;
+SX8 : sextupole, L=Lsext, K = 4.193720*coef, method=method4sextu, N = NqSx;
+SX9 : sextupole, L=Lsext, K = -3.152295*coef, method=method4sextu, N = NqSx;
+SX10: sextupole, L=Lsext, K = 1.936432*coef, method=method4sextu, N = NqSx;
+
+{** Cavities **}
+CAV: Cavity, harnum=416, Frequency=352.202d6, Voltage=2.8d6;
+
+{** Horizontal correctors **}
+HCM : corrector, horizontal, method=intmeth;
+
+{* Dipole for modeling dipole errors in bending magnets*}
+DIP : corrector, horizontal, method=intmeth;
+
+BEND : DIP, BEND0;
+
+{** Vertical correctors **}
+VCM : corrector, vertical, method=intmeth;
+
+
+
+{** Skew quads **}
+Nsq = 4;
+
+{**QT : quadrupole, tilt=45.0, K= 0.0, method=intmeth,N=Nsq; **}
+QT : multipole, N = 1, method = intmeth;
+
+{** BPM **}
+BPM: Beam Position Monitor ;
+
+{** girder start (gs) and end (ge) *}
+
+gs : multipole, N = 1, method = intmeth; {girder start }
+ge : multipole, N = 1, method = intmeth; {girder end}
+
+{** Marker **}
+ SECT1:marker;
+ FIN : marker;
+ DEBUT : marker;
+ ssep : marker;
+ esep : marker;
+ eesep : marker;
+ sSDM : marker;
+ eSDM : marker;
+ eeSDM : marker;
+ sSDAC : marker;
+ eSDAC : marker;
+ eHU600: marker;
+
+{**Superperiods**}
+SINJ: SD1a,ssep,3*SEP,esep,SD1c,eHU600,SD1d;
+
+SD13L: SD13k,sSDM, 13*SD13i,2*SD13j,13*SD13i,
+ eSDM,eeSDM,SD13k;
+SDACL: SDAC5,sSDAC, 9*SDAC2, eSDAC,SDAC5;
+SDACL1: SDAC4,sSDAC, 9*SDAC2, eSDAC,SDAC4;
+
+
+SUP1:SINJ, BPM, SDB1i, gs, Q1, SD2, SX1, QT, HCM, VCM, SD3, Q2, SDB2,
+ BPM, SD14a, Q3, SD5, SX2, HCM, VCM, ge, SD6, BEND,
+ SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB3, BPM, SD9a, Q5, SD12,
+ SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD91, BPM, SDB5, SX3, HCM, VCM, SD8,
+ Q4, ge, SD7, BEND, SD7, gs, Q6, SD5, SX5, HCM, VCM, SD41, BPM, SDB6, Q7,
+ SD3, SX6, HCM, VCM, SD2, Q8, ge, SDC12, BPM, SD13L, BPM, SDC13, gs,
+ Q8, SD2 ,SX8, QT, HCM, VCM, SD3, Q7, SDB7, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB8, BPM, SD9a, Q10, SD8, SX10,
+ QT, HCM, VCM, ge, SDC44, BPM, SDACL1, BPM, SDC45, gs,
+ SX10, HCM, VCM, SD8, Q10, SD91, BPM, SD42, SX9, HCM, VCM, SD8, Q9, ge, SD7,
+ BEND, SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB9, Q7, SD3, SX8,
+ QT, HCM, VCM, SD2, Q8, ge, SD87, BPM, 2*SD13a, BPM, SD86, gs,
+ Q8, SD2, SX8, HCM, VCM, SD3, Q7, SDB10, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB11, BPM, SD92, Q10,
+ SD8, SX10, HCM, VCM, ge, SD89, BPM, 2*SDAC1, BPM, SD88, gs,
+ SX10, QT, HCM, VCM, SD8, Q10, SD93, BPM, SDB12, SX9, HCM, VCM, SD8,
+ Q9, ge, SD7, BEND,
+ SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB13, Q7, SD3,
+ SX8, HCM, VCM, SD2, Q8, ge, SD87, BPM,
+ 2*SD13a, BPM, SD85, gs,
+ Q8, SD2, SX6, HCM, VCM, SD3, Q7, SDB14, BPM,
+ SD43, SX5, HCM, VCM, SD5, Q6, ge, SD7,
+ BEND, SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB15, BPM,
+ SD9a, Q5, SD10, SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD93, BPM, SDB17,
+ SX3, HCM, VCM, SD8, Q4, ge, SD7,
+ BEND, SD6, gs, SX2, HCM, VCM, SD5, Q3, SD141, BPM, SDB18, Q2,
+ SD3, SX1, QT, HCM, VCM, SD2, Q1, ge, SDC12, BPM, SD1e;
+
+SUP2: SD1e, BPM, SDB1, gs,
+ Q1, SD2, SX1, QT, HCM, VCM, SD3, Q2, SDB7,
+ BPM, SD14a, Q3, SD5, SX2, HCM, VCM, ge, SD6, BEND,
+ SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB3, BPM, SD9a, Q5, SD12,
+ SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD91, BPM, SDB5, SX3, HCM, VCM, SD8,
+ Q4, ge, SD7, BEND, SD7, gs, Q6, SD5, SX5, HCM, VCM, SD41, BPM, SDB6, Q7,
+ SD3, SX6, HCM, VCM, SD2, Q8, ge, SDC12, BPM, SD13L, BPM, SD87, gs,
+ Q8, SD2 ,SX8, QT, HCM, VCM, SD3, Q7, SDB7, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB8, BPM, SD9a, Q10, SD8, SX10,
+ QT, HCM, VCM, ge, SDC44, BPM, SDACL1, BPM, SDC45, gs,
+ SX10, HCM, VCM, SD8, Q10, SD91, BPM, SD42, SX9, HCM, VCM, SD8, Q9, ge, SD7,
+ BEND, SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB9, Q7, SD3, SX8,
+ QT, HCM, VCM, SD2, Q8, ge, SD87, BPM, 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX8, HCM, VCM, SD3, Q7, SDB10, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB11, BPM, SD92, Q10,
+ SD8, SX10, HCM, VCM, ge, SD89, BPM, 2*SDAC1, BPM, SD88, gs,
+ SX10, QT, HCM, VCM, SD8, Q10, SD93, BPM, SDB12, SX9, HCM, VCM, SD8,
+ Q9, ge, SD7, BEND,
+ SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB13, Q7, SD3,
+ SX8, HCM, VCM, SD2, Q8, ge, SD87, BPM,
+ 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX6, HCM, VCM, SD3, Q7, SDB14, BPM,
+ SD43, SX5, HCM, VCM, SD5, Q6, ge, SD7,
+ BEND, SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB15, BPM,
+ SD9a, Q5, SD10, SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD93, BPM, SDB17,
+ SX3, HCM, VCM, SD8, Q4, ge, SD7,
+ BEND, SD6, gs, SX2, HCM, VCM, SD5, Q3, SD141, BPM, SDB18, Q2,
+ SD3, SX1, QT, HCM, VCM, SD2, Q1, ge, SDC12, BPM, SD1e;
+
+SUP3: SD1e, BPM, SDB1, gs,
+ Q1, SD2, SX1, QT, HCM, VCM, SD3, Q2, SDB2,
+ BPM, SD14a, Q3, SD5, SX2, HCM, VCM, ge, SD6, BEND,
+ SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB3, BPM, SD9a, Q5, SD12,
+ SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD91, BPM, SDB5, SX3, HCM, VCM, SD8,
+ Q4, ge, SD7, BEND, SD7, gs, Q6, SD5, SX5, HCM, VCM, SD41, BPM, SDB6, Q7,
+ SD3, SX6, HCM, VCM, SD2, Q8, ge, SDC12, BPM, SD13L, BPM, SD87, gs,
+ Q8, SD2 ,SX8, QT, HCM, VCM, SD3, Q7, SDB7, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB8, BPM, SD9a, Q10, SD8, SX10,
+ QT, HCM, VCM, ge, SDC44, BPM, SDACL1, BPM, SDC45, gs,
+ SX10, HCM, VCM, SD8, Q10, SD91, BPM, SD42, SX9, HCM, VCM, SD8, Q9, ge, SD7,
+ BEND, SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB9, Q7, SD3, SX8,
+ QT, HCM, VCM, SD2, Q8, ge, SD87, BPM, 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX8, HCM, VCM, SD3, Q7, SDB10, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB11, BPM, SD92, Q10,
+ SD8, SX10, HCM, VCM, ge, SD89, BPM, 2*SDAC1, BPM, SD88, gs,
+ SX10, QT, HCM, VCM, SD8, Q10, SD93, BPM, SDB12, SX9, HCM, VCM, SD8,
+ Q9, ge, SD7, BEND,
+ SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB13, Q7, SD3,
+ SX8, HCM, VCM, SD2, Q8, ge, SD87, BPM,
+ 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX6, HCM, VCM, SD3, Q7, SDB14, BPM,
+ SD43, SX5, HCM, VCM, SD5, Q6, ge, SD7,
+ BEND, SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB15, BPM,
+ SD9a, Q5, SD10, SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD93, BPM, SDB17,
+ SX3, HCM, VCM, SD8, Q4, ge, SD7,
+ BEND, SD6, gs, SX2, HCM, VCM, SD5, Q3, SD141, BPM, SDB18, Q2,
+ SD3, SX1, QT, HCM, VCM, SD2, Q1, ge, SDC12, BPM, SD1e;
+
+
+SUP4: SD1e, BPM, SDB1, gs,
+Q1, SD2, SX1, QT, HCM, VCM, SD3, Q2, SDB2,
+ BPM, SD14a, Q3, SD5, SX2, HCM, VCM, ge, SD6, BEND,
+ SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB3, BPM, SD9a, Q5, SD12,
+ SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD91, BPM, SDB5, SX3, HCM, VCM, SD8,
+ Q4, ge, SD7, BEND, SD7, gs, Q6, SD5, SX5, HCM, VCM, SD41, BPM, SDB6, Q7,
+ SD3, SX6, HCM, VCM, SD2, Q8, ge, SDC12, BPM, SD13L, BPM, SD87, gs,
+ Q8, SD2 ,SX8, QT, HCM, VCM, SD3, Q7, SDB7, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB8, BPM, SD9a, Q10, SD8, SX10,
+ QT, HCM, VCM, ge, SDC44, BPM, SDACL1, BPM, SDC45, gs,
+ SX10, HCM, VCM, SD8, Q10, SD91, BPM, SD42, SX9, HCM, VCM, SD8, Q9, ge, SD7,
+ BEND, SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB9, Q7, SD3, SX8,
+ QT, HCM, VCM, SD2, Q8, ge, SD87, BPM, 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX8, HCM, VCM, SD3, Q7, SDB10, BPM, SD42, SX7, HCM, VCM,
+ SD5, Q6, ge, SD7, BEND,
+ SD7, gs, Q9, SD8, SX9, HCM, VCM, SDB11, BPM, SD92, Q10,
+ SD8, SX10, HCM, VCM, ge, SD89, BPM, 2*SDAC1, BPM, SD88, gs,
+ SX10, QT, HCM, VCM, SD8, Q10, SD93, BPM, SDB12, SX9, HCM, VCM, SD8,
+ Q9, ge, SD7, BEND,
+ SD7, gs, Q6, SD5, SX7, HCM, VCM, SD41, BPM, SDB13, Q7, SD3,
+ SX8, HCM, VCM, SD2, Q8, ge, SD87, BPM,
+ 2*SD13a, BPM, SD87, gs,
+ Q8, SD2, SX6, HCM, VCM, SD3, Q7, SDB14, BPM,
+ SD43, SX5, HCM, VCM, SD5, Q6, ge, SD7,
+ BEND, SD7, gs, Q4, SD8, SX3, HCM, VCM, SDB15, BPM,
+ SD9a, Q5, SD10, SX4, QT, HCM, VCM, SDB4, BPM, SD10a, Q5, SD93, BPM, SDB17,
+ SX3, HCM, VCM, SD8, Q4, ge, SD7,
+ BEND, SD6, gs, SX2, HCM, VCM, SD5, Q3, SD141, BPM, SDB18, Q2,
+ SD3, SX1, QT, HCM, VCM, SD2, Q1, ge, SDC12, BPM, SD1e;
+
+RING: DEBUT,SUP1,SUP2,SUP3,SUP4,CAV,FIN;
+
+{***** Define CELL structure ****}
+ CELL: RING, symmetry = 1;
+
+end;