# na-shifted-echo-antiecho-phase-filtering-p3.in

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

# spin-3/2 echo amplitude optimization
# versus the third-pulse duration
# in three-pulse shifted-echo-antiecho 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   96*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 {1}

  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.153390697 6.07406384e-17
# -0.597967623 1.60551805e-16
# -1.27569429 2.37145228e-16
# -2.09327349 4.12873671e-16
# -2.94341567 5.70292215e-16
# -3.73053353 6.78092392e-16
# -4.39273843 7.74382581e-16
# -4.91200417 8.34053485e-16
# -5.30850506 8.45163523e-16
# -5.6218286 8.20621556e-16
# -5.88758683 7.84709724e-16
# -6.11973632 6.91779484e-16
# -6.30568774 5.98592275e-16
# -6.4148143 5.1899559e-16
# -6.41470629 5.18491169e-16
# -6.28662083 6.33325876e-16
# END