# al-split-t1-echo-96-phase-cycling-p3.in

# oQ -> 3Q -> 1Q -> -1Q
# approach: 96-phase cycling the first pulse, the third pulse and the receiver

# spin-5/2 echo amplitude optimization
# versus the third-pulse duration
# in three-pulse split-t1 phase-modulated 3Q-MAS sequence,
# the first-pulse duration p1 = 9 micro seconds
# the second-pulse duration p2 = 1 micro seconds

spinsys {
  channels 27Al
  nuclei   27Al
  quadrupole 1 1 3e6 1 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   I1z
  detect_operator  I1c
  verbose          1101
  variable rf      90000
  variable rf3     90000
  variable p1      9 
  variable p2      1
  variable N1      12
  variable N3      8
  variable deltap1 -3
  variable deltap3 2
}

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

  acq $par(ph31)
  pulse  $par(p1) $par(rf) $par(ph1)
  pulse  $par(p2) $par(rf) x

  for {set i 1} {$i < $par(np)} {incr i} {
    pulse $par(tsw) $par(rf3) $par(ph3) 
    acq       [expr $par(ph31) - 90]

  }
}

proc main {} {
  global par

  for {set j1 0} {$j1 < $par(N1)} {incr j1} {
    set par(ph1)   [expr $j1*360./$par(N1)]
    for {set j3 0} {$j3 < $par(N3)} {incr j3} {
      set par(ph3)   [expr $j3*360./$par(N3)]
      set par(ph31)  [expr $par(deltap1)*$par(ph1) + $par(deltap3)*$par(ph3) ]

      set g [fsimpson]
      if [info exists f] {
        fadd $f $g
        funload $g
      } else {
        set f $g
      }
    }
  }

  fsave $f $par(name).fid
  funload $f
  puts "Larmor frequency (Hz) of 27Al: "
  puts [resfreq 27Al $par(proton_frequency)]
}

# SIMP
# NP=17
# SW=4000000
# TYPE=FID
# DATA
# 0 0
# 0.154920405 -8.32667268e-16
# 0.743556466 -1.33226763e-15
# 1.28550259 -2.99760217e-15
# 0.664035368 -3.60822483e-15
# -2.13624432 -1.22124533e-15
# -7.20995796 -2.9559688e-15
# -13.3320502 1.42941214e-15
# -18.3668067 3.10862447e-15
# -20.3272547 4.77395901e-15
# -18.4306837 3.27515792e-15
# -13.4899886 5.88418203e-15
# -7.40362317 2.10942375e-15
# -2.11851943 2.44249065e-15
# 1.26525627 4.4408921e-15
# 2.76741745 -3.33066907e-16
# 3.13594975 6.9388939e-16
# 3.17180728 -3.74700271e-16
# 3.22446638 -3.05311332e-16
# 3.07790859 -2.60902411e-15
# 2.18781702 -5.21804822e-15
# 0.0522594512 -5.77315973e-15
# -3.48536442 -4.6629367e-15
# -8.07003457 -5.8009153e-15
# -12.8393641 -2.08166817e-15
# -16.6359978 9.43689571e-16
# -18.4074521 -1.50866732e-15
# -17.6495193 6.59194921e-16
# -14.6597964 1.79370407e-15
# -10.4325231 -7.21644966e-16
# -6.19924696 3.88578059e-16
# -2.85693548 -5.27355937e-16
# -0.658782996 7.21644966e-16
# 0.650006115 2.22044605e-15
# END