# al-zfilter-amplitude-12p2p3-filtering-p1b.in

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

# spin-5/2 echo and antiecho amplitude optimization
# versus the first-pulse duration
# in three-pulse z-filtered amplitude-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 2 3e6 0 0 0 0
}

par {
  spin_rate        5000
  variable tsw     0.25
  sw               1.0e6/tsw
  np               34
  crystal_file     rep100
  gamma_angles     10
  proton_frequency 800e6
  start_operator   12*I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     9300
  variable p2      4 
  variable p3      2 
}

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

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

  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.254006376 -0.000213894106
# 0.522012845 -0.000713438891
# 0.759390044 -0.0015038028
# 0.93578504 -0.00252388104
# 1.04488836 -0.00367498893
# 1.09926023 -0.00484861799
# 1.11642039 -0.00595674449
# 1.10676482 -0.0069382218
# 1.07018609 -0.00772715186
# 1.0008938 -0.00822791281
# 0.895256615 -0.00833943224
# 0.757941001 -0.00800996702
# 0.603604373 -0.00727089403
# 0.453456301 -0.00622351469
# 0.327748853 -0.00498998996
# 0.238054371 -0.00365614188
# 0.184689924 -0.00224254003
# 0.160121801 -0.000730395545
# 0.154416862 0.000874988161
# 0.159051485 0.00250085035
# 0.168084453 0.00402064624
# 0.178423662 0.00529471852
# 0.1897555 0.00619675016
# 0.202010199 0.00663577825
# 0.211018986 0.00658971585
# 0.206866219 0.00613271101
# 0.177806907 0.00543042968
# 0.11751184 0.00469083414
# 0.0300481058 0.00408937939
# -0.0706740363 0.00370880405
# -0.165750939 0.00351778522
# -0.237739561 0.00339609282
# -0.276042132 0.00319745858
# END