# al-split-t1-5Q-echo-phase-filtering-p2.in

# oQ -> -5Q -> 1Q -> -1Q
# approach: filtering the desired coherences (-5Q) after the first pulse
#           and (1Q) after the second pulse

# spin-5/2 echo amplitude optimization
# versus the second-pulse duration
# in three-pulse split-t1 phase-modulated 5Q-MAS sequence,
# the first-pulse duration p1 = 13.5 micro seconds
# the third-pulse duration p3 = 2 micro seconds

spinsys {
  channels 27Al
  nuclei   27Al
  quadrupole 1 1 3e6 1 0 0 0
}

par {
  spin_rate        5000
  variable tsw     0.25
  sw               1.0e6/tsw
  np               34
  crystal_file     rep10
  gamma_angles     10
  proton_frequency 800e6
  start_operator   160*I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     90000
  variable p1      13.5
  variable p3      2 
}

proc pulseq {} {
  global par
  maxdt $par(tsw)

  matrix set 1 coherence {-5}
  matrix set 2 coherence {1}

  pulse  $par(p1) $par(rf) x
  store 10

  reset
  acq -y        
  for {set i 1} {$i < $par(np)} {incr i} {
    reset
    prop 10
    filter 1
    pulse [expr ($i)*$par(tsw)] $par(rf) x 
    filter 2
    pulse $par(p3) $par(rf3) x
    acq -y
  }
}

proc main {} {
  global par

  fsave [fsimpson] $par(name).fid
  puts "Larmor frequency (Hz) of 27Al: "
  puts [resfreq 27Al $par(proton_frequency)]
}

# SIMP
# NP=17
# SW=4000000
# TYPE=FID
# DATA
# 0 0
# -2.11643964e-05 -2.3653399e-17
# -0.00134190369 -4.44440898e-18
# -0.0147437802 8.34201259e-18
# -0.0761366602 1.8713687e-17
# -0.252367237 2.900836e-17
# -0.619598191 1.63246694e-16
# -1.21789435 5.43288268e-16
# -2.00391062 6.82970855e-16
# -2.83040209 7.41896528e-16
# -3.4710145 1.25154936e-15
# -3.68601794 2.634691e-15
# -3.30807658 3.89621099e-15
# -2.31261356 2.08875499e-15
# -0.837142421 6.66779084e-16
# 0.856798331 -3.01114459e-16
# 2.46087516 -1.72080117e-15
# 3.70963781 -2.87469452e-15
# 4.44744628 -2.37153096e-15
# 4.64679627 -2.68735019e-15
# 4.38073857 -1.66111939e-15
# 3.77146244 -2.49429488e-15
# 2.94316173 -2.84117707e-15
# 2.00328576 -4.03923985e-15
# 1.05322133 -2.94571833e-15
# 0.202721973 -1.26969849e-15
# -0.434751777 3.56079608e-16
# -0.768014641 1.01819067e-15
# -0.76656563 2.94173444e-15
# -0.485770159 4.52331792e-15
# -0.0478113348 3.10984424e-15
# 0.409914876 2.94846561e-15
# 0.792598036 1.5991844e-15
# 1.04751369 -3.22535178e-16
# END