# composite123z.in (SIMPSON 3.1)

#composite pulse X(90°)-X(180°)X(270°), 
#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
  variable p3      15
}

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) -x
    pulse $par(p3) $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.999202405 0
#-0.996838339 0
#-0.992992912 0
#-0.987804515 0
#-0.981459773 0
#-0.974186731 0
#-0.966246527 0
#-0.957923868 0
#-0.949516636 0
#-0.941325026 0
#-0.933640597 0
#-0.926735662 0
#-0.920853377 0
#-0.916198929 0
#-0.912932129 0
#-0.911161688 0
#-0.910941412 0
#-0.912268431 0
#-0.915083563 0
#-0.919273783 0
#-0.924676725 0
#-0.931087048 0
#-0.938264427 0
#-0.945942889 0
#-0.953841123 0
#-0.961673384 0
#-0.969160578 0
#-0.976041079 0
#-0.982080874 0
#-0.987082609 0
#-0.990893175 0
#-0.993409494 0
#-0.994582265 0
#-0.994417445 0
#-0.99297538 0
#-0.990367545 0
#-0.986750956 0
#-0.982320403 0
#-0.977298747 0
#-0.971925578 0
#-0.96644463 0
#-0.961090382 0
#-0.956074334 0
#-0.951571465 0
#-0.947707386 0
#-0.944546718 0
#-0.942083161 0
#-0.940231716 0
#-0.938823435 0
#-0.937603017 0
#-0.936229475 0
#-0.934279996 0
#-0.931257036 0
#-0.926598552 0
#-0.919691195 0
#-0.909886168 0
#-0.896517357 0
#-0.878921268 0
#-0.856458196 0
#-0.828534046 0
#-0.794622126 0
#-0.754284253 0
#-0.707190496 0
#-0.653136891 0
#-0.592060527 0
#-0.524051436 0
#-0.44936082 0
#-0.368405234 0
#-0.281766447 0
#-0.190186839 0
#END