# spinlockz.in (SIMPSON 3.1)

#spinlock pulse X(90°)Y(90°)-detection, for positive gyromagnetic ratio nuclei, 
#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      5
}

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

    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
#2.77555756e-016 0
#0.039988575 0
#0.0799086882 0
#0.119692187 0
#0.15927139 0
#0.198579231 0
#0.237549404 0
#0.276116505 0
#0.314216173 0
#0.351785226 0
#0.388761797 0
#0.425085465 0
#0.460697381 0
#0.495540396 0
#0.52955918 0
#0.562700335 0
#0.594912509 0
#0.6261465 0
#0.656355358 0
#0.68549448 0
#0.713521701 0
#0.740397373 0
#0.766084449 0
#0.790548551 0
#0.813758033 0
#0.835684044 0
#0.856300575 0
#0.875584507 0
#0.893515649 0
#0.910076764 0
#0.925253599 0
#0.939034897 0
#0.95141241 0
#0.962380899 0
#0.971938129 0
#0.980084858 0
#0.98682482 0
#0.992164698 0
#0.996114087 0
#0.998685463 0
#0.999894129 0
#0.99975817 0
#0.998298386 0
#0.995538233 0
#0.991503752 0
#0.986223488 0
#0.979728412 0
#0.972051832 0
#0.963229299 0
#0.953298512 0
#0.942299216 0
#0.930273093 0
#0.917263655 0
#0.90331613 0
#0.888477341 0
#0.872795591 0
#0.856320538 0
#0.83910307 0
#0.821195174 0
#0.802649815 0
#0.783520798 0
#0.76386264 0
#0.743730438 0
#0.723179734 0
#0.702266386 0
#0.681046433 0
#0.659575965 0
#0.637910991 0
#0.61610731 0
#0.594220387 0
#0.57230522 0
#END