ASSIGN, PRINT="DCO4.print.txt" !, echo="error.out"
  OPTION, -ECHO, DOUBLE, VERIFY
  SETPLOT, XSIZE=25.4, YSIZE=20.32
  SETPLOT, LWIDTH=5, LSCALE=1.5, SSCALE=1.5, RSCALE=1.5  
                   
! ==============================================================================
! Load XSIF file
! ------------------------------------------------------------------------------

  CALL, FILENAME="kvals.DCO4.72.xsif"
!  CALL, FILENAME="kvals.DCO4.90.xsif"
!  CALL, FILENAME="kvals.DCO4.100.xsif"

  CALL, FILENAME="DCO4.xsif"

! ==============================================================================
! SUBROUTINEs ---- 
! ------------------------------------------------------------------------------

comment

  USE, ARC
  CELL
  
    VARY, K1.QARCF, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QARCD, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, MUX=ARCMUX
    CONSTRAINT, #E, MUY=ARCMUY

    SIMPLEX, CALLS=5000
    MIGRAD,  CALLS=5000

  ENDMATCH

USE, ARC
SAVEBETA, LABEL=ENDARC, PLACE=#E
TWISS,SAVE

ARC22 :  LINE=(-ARC2)

USE, ARC2
SAVEBETA, LABEL=ENDARC2, PLACE=#E
TWISS,SAVE

USE, ARC22
SAVEBETA, LABEL=ENDARC22, PLACE=#E
TWISS,SAVE


! ------------------------------------------------------------------------------
!         MUX & MUY 
! ------------------------------------------------------------------------------

  USE, PHT
  SAVEBETA, LABEL=PHTW, PLACE=#E
  TWISS, SAVE

!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
!        VMIN=0,-1, VMAX=50,1                               &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Phase Trombone"
!
! ------------------------------------------------------------------------------
! ------------------------------------------------------------------------------

  IDSPPHT :  LINE=(-DSP)
  
  DSPPHT :  LINE=(DSP, PHT)


  USE, IDSPPHT
  MATCH, BETX=PHTW[BETX], ALFX=-PHTW[ALFX], &
         BETY=PHTW[BETY], ALFY=-PHTW[ALFY], DX=0.0, DPX=0.0
  
    VARY, K1.QDSP4, STEP=1E-5, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QDSP5, STEP=1E-5, LOWER=-0.9, UPPER=+1.2
    VARY, K1.QDSP6, STEP=1E-5, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QDSP7, STEP=1E-5, LOWER=-0.7, UPPER=+0.7

    VARY, K1.QDSP0, STEP=1E-5, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QDSP1, STEP=1E-5, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QDSP2, STEP=1E-5, LOWER=-0.9, UPPER=+0.7
    VARY, K1.QDSP3, STEP=1E-5, LOWER=-0.7, UPPER=+0.9

    WEIGHT, DX=50, DPX=700

    CONSTRAINT, #E, BETX=ENDARC2[BETX], ALFX=-ENDARC2[ALFX]
    CONSTRAINT, #E, BETY=ENDARC2[BETY], ALFY=-ENDARC2[ALFY]
    CONSTRAINT, #E, DX=ENDARC2[DX],     DPX=-ENDARC2[DPX]

    CONSTRAINT, QDSP2[1], BETY<40, BETX<50
    CONSTRAINT, QDSP3[1], BETY<40, BETX<50

    CONSTRAINT, QDSP4[1], BETY<40, BETX<50
    CONSTRAINT, QDSP5[1], BETY<40, BETX<50
    CONSTRAINT, QDSP6[1], BETY<40, BETX<50

    SIMPLEX, CALLS=100000
    MIGRAD,  CALLS=100000

  ENDMATCH

! USE, DSPPHT
! TWISS, BETA0=ENDARC2, SAVE
!
! PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
!        VMIN=0,-1, VMAX=100,1                              &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="STRSEC SECTION"
!
! ------------------------------------------------------------------------------
! ------------------------------------------------------------------------------

  USE, PSD
  MATCH, LINE= WGC
  
    VARY, K1.QPSD0, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD1, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD2, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD3, STEP=1E-5, LOWER=-0.9, UPPER=+0.9

    VARY, K1.QPSD4, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD5, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD6, STEP=1E-5, LOWER=-0.9, UPPER=+0.9
    VARY, K1.QPSD7, STEP=1E-5, LOWER=-0.9, UPPER=+0.9

    WEIGHT, DX=50, DPX=500

    CONSTRAINT, QPSD4[1],BETY<50,BETX<50
    CONSTRAINT, QPSD5[1],BETY<50,BETX<50
    CONSTRAINT, QPSD6[1],BETY<50,BETX<50
    CONSTRAINT, #E, LINE=ARC2

    CONSTRAINT, QPSD1[1],BETY<50,BETX<50
!    CONSTRAINT, QPSD2[1],BETY<50,BETX<50
    CONSTRAINT, QPSD3[1],BETY<50,BETX<50

    SIMPLEX, CALLS=40000
    MIGRAD,  CALLS=40000

  ENDMATCH

!  USE, PSD
!  TWISS,LINE=WGC,SAVE
!
!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
!        VMIN=0,-1, VMAX=100,1                              &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="STRSEC SECTION"
!
!-------------------------------------------------------------------------------
! ------------------------------------------------------------------------------

  USE, CCH
  CELL
  
    VARY, K1.QCCH1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QCCH2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, BETX=40
    CONSTRAINT, #E, BETY=16

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH

! ------------------------------------------------------------------------------
!
!  USE, WIG
!  CELL
!  
!    VARY, K1.QWGC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
!    VARY, K1.QWGC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
!
!    SIMPLEX, CALLS=1000
!    MIGRAD,  CALLS=1000
!
!  ENDMATCH
!
! ------------------------------------------------------------------------------
! ------------------------------------------------------------------------------

  USE, RFC
  CELL
  
    VARY, K1.QRFC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QRFC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, MUX=0.2
    CONSTRAINT, #E, MUY=0.2

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH

! ------------------------------------------------------------------------------


  USE, IJ0
  CELL
  
    VARY, K1.QINJ1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QINJ2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

!    CONSTRAINT, DINJ3[1], BETX=71.0  !70.0
    CONSTRAINT, #E, MUX=0.25, MUY=0.4093382

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH

! ------------------------------------------------------------------------------


  USE, IFD
  CELL
  
    VARY, K1.QIFD1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QIFD2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, MUX=0.13849+0.000427
    CONSTRAINT, #E, MUY=0.21393-0.000106

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH

! ------------------------------------------------------------------------------
! ------------------------------------------------------------------------------

  USE, MFC
  MATCH, LINE=IFD
  
    VARY, K1.QMFC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMFC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMFC3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMFC4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, LINE=CCH

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH


  USE, MFC
  TWISS, LINE=IFD, SAVE

!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
!        VMIN=0,-1, VMAX=50,1                               &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching MFC"
!
! ------------------------------------------------------------------------------


  USE, MPI
  MATCH, LINE=PHT !CCH
  
    VARY, K1.QMPI1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7  !
    VARY, K1.QMPI2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMPI3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7  !
    VARY, K1.QMPI4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7


    CONSTRAINT, #E, LINE=IJ0

    SIMPLEX, CALLS=10000
    MIGRAD,  CALLS=10000

  ENDMATCH

  USE, MPI
  TWISS, LINE=PHT, SAVE

!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
!        VMIN=0, VMAX=100                                    &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching Chicane to Injection"
!
! ------------------------------------------------------------------------------


  USE, MIF
  MATCH, LINE=IJ0
  
    VARY, K1.QMIF1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMIF2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMIF3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMIF4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, LINE=IFD

    SIMPLEX, CALLS=10000
    MIGRAD,  CALLS=10000

  ENDMATCH

  USE, MIF
  TWISS, LINE=IJ0, SAVE

!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
!        VMIN=0, VMAX=100                                    &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching Chicane to Injection"
!
! ------------------------------------------------------------------------------


  USE, MRC
  MATCH, LINE=RFC
  
    VARY, K1.QMRC1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMRC2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMRC3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMRC4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, LINE=CCH

    SIMPLEX, CALLS=1000
    MIGRAD,  CALLS=1000

  ENDMATCH

! ------------------------------------------------------------------------------


  USE, MPR
  MATCH, LINE=PHT
  
    VARY, K1.QMPR1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMPR2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMPR3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMPR4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
 !   VARY, DMPR1[L], STEP=1E-4
 !   VARY, DMPR2[L], STEP=1E-4

    CONSTRAINT, #E, LINE=RFC

    SIMPLEX, CALLS=50000
    MIGRAD,  CALLS=10000

  ENDMATCH

!------------------------------------------------------------------------------


  USE, MDW
  MATCH, LINE=CCH
  
    VARY, K1.QMDW1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMDW2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMDW3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QMDW4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, LINE=WGC

    SIMPLEX, CALLS=20000
    MIGRAD,  CALLS=20000

  ENDMATCH

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!           fine matching of the straights to even # of Pi
!-------------------------------------------------------------------------------

 USE, STRSECR
 SAVEBETA, LABEL=ENDSTRSEC, PLACE=#E
 TWISS,SAVE


 MUXPHT:=PHTW[MUX]-(ENDSTRSEC[MUX]-13.0-ARCMUX)/6.0
 MUYPHT:=PHTW[MUY]-(ENDSTRSEC[MUY]-12.0-ARCMUY)/6.0


  USE, PHT
  CELL
  
    VARY, K1.QPHT1, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QPHT2, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QPHT3, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QPHT4, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QPHT5, STEP=1E-6, LOWER=-0.7, UPPER=+0.7
    VARY, K1.QPHT6, STEP=1E-6, LOWER=-0.7, UPPER=+0.7

    CONSTRAINT, #E, BETX=PHTW[BETX], ALFX=PHTW[ALFX], MUX=MUXPHT
    CONSTRAINT, #E, BETY=PHTW[BETY], ALFY=PHTW[ALFY], MUY=MUYPHT

    CONSTRAINT, #S, BETX=PHTW[BETX], ALFX=PHTW[ALFX]
    CONSTRAINT, #S, BETY=PHTW[BETY], ALFY=PHTW[ALFY]

 !   WEIGHT, MUX=400

 !   CONSTRAINT, QPHT1[1],BETX<50
 !   CONSTRAINT, QPHT2[1],BETX<50
 !   CONSTRAINT, QPHT3[1],BETX<50
 !   CONSTRAINT, QPHT4[1],BETX<50
 !   CONSTRAINT, QPHT5[1],BETX<50
 !   CONSTRAINT, QPHT6[1],BETX<50

    SIMPLEX, CALLS=50000
    MIGRAD,  CALLS=20000

    ENDMATCH

!-------------------------------------------------------------------------------
! Instrumentation and correction components
!-------------------------------------------------------------------------------
   LDBPM := 0.25 !0.25  !Drift from quadrupole to beam position monitor
   LDCOR := 0.1  !0.10  !Drift from orbit corrector to quadrupole
   LCOR  := 0.15 !0.15  !Orbit corrector length
!-------------------------------------------------------------------------------

  USE, RING
  HARMON

    HVARY, K2.SARC1, STEP=1E-6, LOWER=   0.0, UPPER=+100.0
    HVARY, K2.SARC2, STEP=1E-6, LOWER=-100.0, UPPER=   0.0

    HTUNE

  ENDHARM

endcomment

! ==============================================================================
! COMMANDs
! ------------------------------------------------------------------------------

 RF_ON:=0
 SXT_ON := 1

! ------------------------------------------------------------------------------
! Produce some plots
! ------------------------------------------------------------------------------

  USE, ARC
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Arc Cell"


  USE, DSP
  TWISS, LINE=ARC2, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Dispersion Suppressor"
 
       
  USE, PHT
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Phase Trombone"
        
  
  USE, MPR
  TWISS, LINE=PHT, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching PHT to RF"


  USE, RFC
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="RF Cavities"
        

  USE, MRC
  TWISS, LINE=RFC, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching RF to Chicane"
        
  USE, CCH
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Circumference Chicane"

        
  USE, MDW
  TWISS, LINE=CCH, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,             &
        VMIN=0, VMAX=50                                     &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching Chicane to Wiggler"
        

  USE, WGC
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0, VMAX=50                                    &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Wiggler"

        
  USE, PSD
  TWISS, LINE=WGC, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=50,1                                    &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Reverse Dispersion Suppressor"


  USE, STRSECR
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="STRSECR with RF and WIG"


  USE, HRINGR
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Half Ring (STRSECR+ARCs)"


!*************************  INJ  *******************************

  USE, MPI
  TWISS, LINE=PHT, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0, VMAX=100                                    &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching PHT to Injection"


  USE, INJ
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0, VMAX=100                                   &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Injection"

  USE, EXT
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0, VMAX=100                                   &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Extraction"


  USE, MIF
  TWISS, LINE=IJ0, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0, VMAX=100                                    &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching Injection to IFD"

      
  USE, IFD
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="IFD"


  USE, MFC
  TWISS, LINE=IFD, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,            &
        VMIN=0,-1, VMAX=50,1                               &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Matching IFD to Chicane"

               
  USE, PSD
  TWISS, LINE=WGC, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=50,1                                    &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Reverse Dispersion Suppressor"


  USE, STRSECI
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="STRSECI"


  USE, HRINGI
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="Half Ring (STRSECI+ARC)"


!************************** RING *********************************

  USE, RING
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="DR Ring"


  USE, RNGI
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        VMIN=0,-1, VMAX=100,1                              &
        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="DR starts from IP"

!************************** SUPP ***********************************
!
!  USE, A0LD
!  TWISS, LINE=ARC, SAVE
!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
!        VMIN=0,-1, VMAX=60,1                               &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="ARC+last a.cell+d.supp to PHT"
!
!  USE, PRA0
!  TWISS, LINE=WGC, SAVE
!  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY,VAXIS2=DX,  &
!        VMIN=0,-1, VMAX=60,1                               &
!        STYLE=1, COLOUR=100, SPLINE=.T., TITLE="d.supp+first a.cell+ARC"
!

!-----------------------------------------------------------------------
!
! ------------------------------------------------------------------------------
! Output the lattice functions
! ------------------------------------------------------------------------------

  USE, RNGI
  TWISS, SAVE
  SURVEY, TAPE="DCO4.survey.txt"

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E

  OPTICS, COLUMNS= NAME, S, L,               &
                   BETX, ALFX, MUX, DX, DPX, &
                   BETY, ALFY, MUY,          &
                   K0L, E1, E2, K1L, K2L,    &
                   FILENAME= "DCO4.lfunctions.txt"

  OPTICS, COLUMNS= NAME, KEYWORD, S, L,                    &
                   K0L, E1, E2, K1L, KS1L, K2L, K3L, TILT, &
                   FILENAME= "DCO4.lattice.txt"


! ------------------------------------------------------------------------------
! Calculate the beam parameters
! ------------------------------------------------------------------------------

  ASSIGN, PRINT="DCO4.BMPM.txt"

  RF_ON  := 1
  SXT_ON := 1
  USE, RNGI
  TWISS, SAVE
  BEAM, PARTICLE=POSITRON, ENERGY=E0, RADIATE
  BMPM


  stop