# al-zfilter-amplitude-24p2p3-filtering-p3b.in

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

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

spinsys {
  channels 27Al
  nuclei   27Al
  quadrupole 1 2 3e6 0 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   24*I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     9300
  variable p1      4
  variable p2      4 
}

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

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

  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 27Al: "
  puts [resfreq 27Al $par(proton_frequency)]
}

# SIMP
# NP=17
# SW=4000000
# TYPE=FID
# DATA
# 0 0
# -3.19256558e-09 3.56061129e-12
# -9.12431989e-08 8.10299998e-11
# -5.86790254e-07 6.59804346e-11
# -1.98784469e-06 -2.86328025e-09
# -4.6598091e-06 -1.93929209e-08
# -8.5994695e-06 -7.42522734e-08
# -1.35059427e-05 -2.09736708e-07
# -1.91312035e-05 -4.83468957e-07
# -2.56372492e-05 -9.62280164e-07
# -3.36773433e-05 -1.71732414e-06
# -4.43261121e-05 -2.82092226e-06
# -5.92455393e-05 -4.34302255e-06
# -8.1077853e-05 -6.35047152e-06
# -0.000113643046 -8.91362065e-06
# -0.000161651159 -1.21188049e-05
# -0.000230031391 -1.6079536e-05
# END