# al-split-t1-antiecho-phase-filtering-p1.in

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

# spin-5/2 antiecho amplitude optimization
# versus the first-pulse duration
# in three-pulse split-t1 phase-modulated -3Q-MAS sequence,
# the second-pulse duration p2 = 4 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.5
  sw               1.0e6/tsw
  np               34
  crystal_file     rep100
  gamma_angles     10
  proton_frequency 800e6
  start_operator   96*I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     90000
  variable p2      4 
  variable p3      2 
}

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

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

  acq -y        
  for {set i 1} {$i < $par(np)} {incr i} {
    pulse $par(tsw) $par(rf) x
    store 2
    filter 1
    pulse $par(p2)  $par(rf) x
    filter 2
    pulse $par(p3) $par(rf3) x
    acq -y
    reset
    prop 2
  }
}

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
# 0.00639701201 -2.38465192e-16
# 0.162297846 -2.53589315e-16
# 0.628573069 -1.71903999e-16
# 1.05159054 -2.78452812e-16
# 0.963883756 3.94701703e-16
# 0.499743815 4.58085208e-16
# 0.185551782 -4.27888495e-16
# 0.209956954 -2.37444557e-16
# 0.619058221 2.65902291e-16
# 1.60568476 -1.60515581e-16
# 3.0948265 -5.13556488e-18
# 4.72771476 -3.40554638e-16
# 6.32090024 -9.05198513e-16
# 7.86276471 -6.8583263e-16
# 9.23715297 -8.10545441e-16
# 10.1778291 -1.49777898e-15
# 10.4966992 -1.7626464e-15
# 10.115117 -4.0983924e-16
# 9.06593061 5.85251577e-16
# 7.54371316 -3.80803365e-16
# 5.7139584 -7.05682301e-16
# 3.81403449 -8.5184156e-16
# 2.16729874 -1.03484308e-15
# 0.998499488 1.25326672e-16
# 0.32900288 4.42553825e-16
# -0.0584787644 -5.928312e-16
# -0.298674881 -1.3873031e-15
# -0.386008848 -3.25999762e-16
# -0.277161016 5.59699904e-16
# 0.0352134653 -5.73617955e-16
# 0.569081452 -1.18446125e-15
# 1.42393358 -1.21442071e-15
# 2.4867949 -6.91560549e-16
# END