# composite180.in (SIMPSON 3.1)

#composite pulse X(90°)Y(180°)X(90°), 
#Mz versus offset varying from 0 to np kHz in steps of 1 kHz

spinsys {
  channels 29Si
  nuclei   29Si
  shift 1 0 0 0 0 0 0
}

par {
  proton_frequency 400e6
  variable tsw     0.5
  sw               1.0e6/tsw
  np               71
  crystal_file     alpha0beta0
  gamma_angles     1
  start_operator   I1z
  detect_operator  2*I1z
  verbose          1101
  variable rf      50000
  variable p1      5
  variable p2      10
}

proc pulseq {} {
  global par

  for {set i 0} {$i < $par(np)} {incr i} {
    
    offset [expr $i*1000]

    pulse $par(p1) $par(rf) x
    pulse $par(p2) $par(rf) y
    pulse $par(p1) $par(rf) x

    acq

    reset
  }
}

proc main {} {
  global par

  fsave [fsimpson] $par(name).fid

  puts "Larmor frequency (Hz) of 29Si: "
  puts [resfreq 29Si $par(proton_frequency)]
}

#SIMP
#NP=71
#SW=2000000
#TYPE=FID
#DATA
#-1 0
#-0.999201142 0
#-0.996818244 0
#-0.992892115 0
#-0.98749005 0
#-0.980704759 0
#-0.972652891 0
#-0.963473165 0
#-0.953324156 0
#-0.942381735 0
#-0.930836233 0
#-0.918889346 0
#-0.906750829 0
#-0.894635045 0
#-0.882757383 0
#-0.87133064 0
#-0.860561377 0
#-0.85064633 0
#-0.841768924 0
#-0.834095927 0
#-0.827774307 0
#-0.822928342 0
#-0.819657007 0
#-0.8180317 0
#-0.818094331 0
#-0.819855805 0
#-0.823294932 0
#-0.828357769 0
#-0.834957433 0
#-0.84297436 0
#-0.852257047 0
#-0.862623239 0
#-0.873861583 0
#-0.885733699 0
#-0.897976674 0
#-0.910305935 0
#-0.922418464 0
#-0.933996333 0
#-0.944710501 0
#-0.954224826 0
#-0.962200259 0
#-0.968299147 0
#-0.972189597 0
#-0.973549855 0
#-0.972072634 0
#-0.967469333 0
#-0.9594741 0
#-0.947847683 0
#-0.932381017 0
#-0.912898499 0
#-0.88926091 0
#-0.861367939 0
#-0.829160281 0
#-0.792621271 0
#-0.751778042 0
#-0.706702172 0
#-0.657509822 0
#-0.604361359 0
#-0.547460449 0
#-0.487052651 0
#-0.423423507 0
#-0.356896151 0
#-0.287828476 0
#-0.216609867 0
#-0.143657563 0
#-0.0694126623 0
#0.00566415856 0
#0.0810971812 0
#0.156400376 0
#0.231082274 0
#0.304650925 0
#END