# na-zfilter-phase-filtering-p3.in

# approach: filtering the desired coherence (-3Q) after the first pulse
#           and (0Q) after the second pulse

# spin-3/2 echo amplitude optimization
# versus the third-pulse duration
# in three-pulse z-filtered phase-modulated -3Q-MAS sequence,
# the first-pulse duration p1 = 4 micro seconds
# the second-pulse duration p2 = 4 micro seconds

spinsys {
  channels 23Na
  nuclei   23Na
  quadrupole 1 1 1e6 1 0 0 0
}

par {
  spin_rate        5000
  variable tsw     0.25
  sw               1.0e6/tsw
  np               17
  crystal_file     rep100
  gamma_angles     10
  proton_frequency 800e6
  start_operator   48*I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     93000
  variable p1      4 
  variable p2      4 
}

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

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

  acq -y        
  pulse  $par(p1) $par(rf) x
  filter 1
  pulse  $par(p2) $par(rf) x
  filter 2

  for {set i 1} {$i < $par(np)} {incr i} {
    pulse $par(tsw) $par(rf3) x 
    acq -y
  }
}

proc main {} {
  global par

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

# SIMP
# NP=17
# SW=4000000
# TYPE=FID
# DATA
# 0 0
# 0.127795698 8.03159512e-17
# 0.153849443 1.5656893e-16
# 0.0071825619 2.20119894e-16
# -0.326891787 3.16838941e-16
# -0.802336705 4.00302421e-16
# -1.33025188 5.05346851e-16
# -1.81051246 5.48786499e-16
# -2.16410041 5.76803086e-16
# -2.35439796 5.36857831e-16
# -2.39152474 5.51956568e-16
# -2.3208468 5.36745579e-16
# -2.20201985 5.62221981e-16
# -2.08685205 5.84366829e-16
# -2.00313094 6.33485335e-16
# -1.94858498 6.42783431e-16
# -1.89563253 6.29921118e-16
# END