diff --git a/lattice/example/soleil/Example_input.prm b/lattice/example/soleil/Example_input.prm
new file mode 100644
index 0000000000000000000000000000000000000000..c03d692bf0dbaf5b831dbde9458905b36e87be51
--- /dev/null
+++ b/lattice/example/soleil/Example_input.prm
@@ -0,0 +1,205 @@
+#*******************************************************************
+# 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 /home/zhang/codes/TracyIII/lattice/   
+
+# 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 
+
+# Deactivate quadrupole Fringefield
+  QuadFringeOffFlag
+
+# 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
+
+#======================================
+#specific for SOLEIL lattice
+#======================================
+
+# read the multipole field error, and add 
+# the corresponding components of the element
+  multipole_file multipole_error_fullquad.dat
+# files for looking for the multipole field errors of 
+# correctors and skew quadrupoles, which are integrated in the sextupoles
+  fic_hcorr corh.txt
+  fic_vcorr corv.txt
+  fic_skew  corqt.txt
+
+  ReadMultipoleFlag  
+
+
+# file of source of coupling; only works for SOLEIL lattice
+  virtualskewquad_file   virtual_skew_quad_currents.txt 
+
+#### Virtual Skew quadrupole FACTORY
+# read the virtual setting of skew quadrupoles from a file
+  ReadVirtualSkewquadFlag 
+
+# add random rotation error to the two half quadrupoles
+# to introduce coupling into the lattice
+  ErrorCoupling2Flag   0  0.0007
+
+# fit tune for the lattice with each quadrupole is cut into two pieces
+# name, flag, tuned quadrople,tuned quadrople,targetnudx    targetnudz
+  FitTune4Flag     qp7a qp7b  qp9a qp9b  18.202 10.317
+
+#=================================================
+## general command (NOT machine dependent)
+#=================================================  
+
+
+
+
+# Read the multipole field errors, and replacement the corresponding components
+# of the element.
+ReadfefileFlag   dip.fe
+
+# read the misalignment error from a file 
+  ReadaefileFlag    soleil1.ae
+
+#####################
+#  orbit correction
+  CODCorrectFlag
+
+#parameters used to do orbit correction
+# file to specify which correctors are used for H orbit correction  
+  hcorr_file  hcorr_56nom.state
+# file to specify which correctors are used for V orbit correction  
+  vcorr_file  vcorr_56nom.state
+
+#number of iterations for orbit correction, should NOT small than 1
+n_orbit       5     
+#number of singular values in H-plane, must be not larger than the number of correctors used for orbit correction
+nwh           60    
+#number of singular values in V-plane, must be not larger  than the number of correctors used for orbit correction
+nwv           60    
+
+
+
+
+
+
+
+#  print the twiss paramters to a file
+  PrintTwissFlag     twisstest.out
+
+# print close orbit to a file
+  PrintCODFlag testcod.out
+
+#print  coordinates tracking around COD at each element to a file 
+ PrintTrackFlag  tracktest.out  0.001 0.000  0.0  0.0  0.0  0.0   50
+
+
+# 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
+
+# tune shift with amplitude 
+  AmplitudeTuneShiftFlag nudxtest.out nudztest.out   50 30 516 0.035 0.02 0.0
+
+# tune shift with energy
+  EnergyTuneShiftFlag nudptest.out  31 1026 0.06 
+
+# frequency map analysis for on momentum particle
+ FmapFlag  fmaptest.out 31 21 516 0.025 0.005 0.0 true
+
+# frequency map analysis for off momentum particle
+ FmapdpFlag  fmapdptest.out 101 121 1026 25e-3 0.06 0.3e-3 true
+
+
+# add random rotation error to the full quadrupole
+# to introduce coupling into the lattice
+  ErrorCouplingFlag   0  0.0007
+
+# calcualte the coupling factor in the lattice 
+  CouplingFlag  
+
+# calculate momentum acceptance 
+  MomentumAccFlag momacctest.out  6D 1 189 0.01 0.05 100 -0.01 -0.05 100 1026  0.0002
+
+# calculate phase space
+  PhaseSpaceFlag phasetest.out 6D 1e-6 0.0 1e-6 0.0 0.012 0.0 1000 false
+
+
+#===============================================
+### OTHER (for debug ......)
+#===============================================
+# calculate induced amplitude
+  InducedAmplitudeFlag  
+
+# specific for test
+  CodeComparaisonFlag  
+
+# calculate Eta function
+  EtaFlag  
+
+
+#================================================
+#        TO BE UPDATED
+#===============================================
+
+### Touschek lifetime
+# calculate Touschek lifetime
+#  TouschekFlag 
+# Intra Beam Scattering
+#  IBSFlag  
+# Track momentum acceptance and then get Touschek lifetime
+#  TousTrackFlag 
+
+# print girder
+#PrintGirderFlag testgirder.out
+
+
+
+