# composite123y.in (SIMPSON 3.1)

#composite pulse X(90°)-X(180°)X(270°), positive gyromagnetic ratio nucleus
#My 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*I1y
  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
#6.47106067e-016 0
#0.00097013321 0
#0.00386198213 0
#0.00862010888 0
#0.0151528389 0
#0.0233333698 0
#0.0330013731 0
#0.0439651323 0
#0.0560042636 0
#0.0688730651 0
#0.0823045314 0
#0.0960150554 0
#0.109709819 0
#0.123088848 0
#0.135853672 0
#0.147714506 0
#0.158397807 0
#0.167654064 0
#0.175265602 0
#0.181054179 0
#0.184888115 0
#0.18668868 0
#0.186435466 0
#0.184170456 0
#0.180000534 0
#0.174098196 0
#0.166700275 0
#0.158104527 0
#0.148663995 0
#0.138779152 0
#0.12888788 0
#0.119453434 0
#0.110950629 0
#0.103850566 0
#0.0986042864 0
#0.095625817 0
#0.0952751264 0
#0.0978415546 0
#0.103528304 0
#0.11243859 0
#0.124564027 0
#0.139775807 0
#0.157819155 0
#0.17831148 0
#0.200744533 0
#0.224490793 0
#0.248814139 0
#0.272884771 0
#0.295798159 0
#0.316597687 0
#0.334300495 0
#0.34792591 0
#0.356525706 0
#0.359215375 0
#0.355205464 0
#0.343832028 0
#0.324585173 0
#0.29713473 0
#0.261352078 0
#0.21732726 0
#0.165380608 0
#0.106068226 0
#0.0401808702 0
#-0.0312640983 0
#-0.107037771 0
#-0.185721663 0
#-0.265736259 0
#-0.345375612 0
#-0.422847226 0
#-0.496316308 0
#-0.563953263 0
#END