# twoPulseEcho_P2.in
# spin-5/2 central-line intensity optimization
# with P2 in two-pulse MAS sequence,
# echo transfer pathway of I = 5/2

spinsys {
  channels 27Al
  nuclei   27Al
  quadrupole 1 2 5e6 -1 0 0 0
}

par {
  spin_rate        15000
  variable tsw     0.25
  sw               1.0e6/tsw
  np               17
  crystal_file     rep100
  gamma_angles     10
  start_operator   I1z
  verbose          1101
  variable rf      100000
  proton_frequency 800e6
}

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

  matrix set 1 elements {{2 5}}
  matrix set detect elements {{3 4}}

  acq
  pulse 4 $par(rf) x
  filter 1

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

proc main {} {
  global par

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