!--------1---------2---------3---------4---------5---------6---------7--------8
  TITLE "ILC Damping Ring (OCS6) 20 June 2006"
!--------1---------2---------3---------4---------5---------6---------7--------8

  ASSIGN, PRINT="OCS6.print.txt"
  ASSIGN, ECHO="OCS6.echo.txt"

  SETPLOT, XSIZE=25.4, YSIZE=20.32
  SETPLOT, LWIDTH=5, LSCALE=1.5,SSCALE=1.5,RSCALE=1.5 

!--------1---------2---------3---------4---------5---------6---------7--------8

  OPTION, -ECHO, DOUBLE, VERIFY

  CALL, FILENAME="OCS6.xsif"

  OPTION, ECHO

!--------1---------2---------3---------4---------5---------6---------7--------8

  BEAM, PARTICLE=ELECTRON, ENERGY=E0

  RF_ON := 0

  USE, RING

  SELECT, OPTICS, CLEAR
  SELECT, OPTICS, #S/#E
  OPTICS, COLUMNS= NAME, KEYWORD, S, L, K0L, E1, E2, &
                   K1L, KS1L, K2L, K3L, TILT, TYPE, &
          FILENAME= "OCS6.lattice.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= "OCS6.lfunctions.txt"

!--------1---------2---------3---------4---------5---------6---------7--------8

  RF_ON  := 0
  SXT_ON := 1

  USE, ACELL
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS6.ACELL.ps"

  USE, CWIG
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS6.CWIG.ps"

  USE, SECTORW
  PRINT, #S/#E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS6.SECTORW.ps"

  USE, RING
  PRINT, #E
  TWISS, SAVE
  PLOT, TABLE=TWISS, HAXIS=S, VAXIS1=BETX,BETY, VAXIS2=DX, &
        STYLE=1, COLOUR=100, SPLINE=.T., FILE="OCS6.RING.ps"

!--------1---------2---------3---------4---------5---------6---------7--------8

  TWISS, TUNES, DELTAP=-0.01:0.01:0.001, SAVE
  PLOT, HAXIS=DELTAP, VAXIS1=QX, VAXIS2=QY, &
                      SPLINE, COLOUR=100, STYLE=1

!--------1---------2---------3---------4---------5---------6---------7--------8

  ASSIGN, PRINT="OCS6.bmpm.txt"

  USE, RING

  RF_ON  := 0
  SXT_ON := 0
  PRINT, #E
  TWISS, TUNES, SAVE

  RF_ON  := 1
  SXT_ON := 1
  PRINT, #E
  TWISS, TUNES, SAVE

! Note: BMPM does not work when skew quadrupoles are present in lattice

  BMPM

  BEAM, PARTICLE=electron, ENERGY=E0, NPART=2.0E10, RADIATE
  EMIT

  STOP