# onepowderMAS-deadtime.in 
# Spin-3/2 central-line intensity calculation
# for a powder rotating at the magic angle,
# submitted to the first- and the second-order
# quadrupole interactions.

spinsys {
  channels 23Na
  nuclei   23Na
  quadrupole 1 2 8e6 1 0 0 0
}

par {
  proton_frequency 400e6
  spin_rate        15000
  variable tsw     1
  sw               1.0e6/tsw
  np               21
  crystal_file     rep100
  gamma_angles     3
  start_operator   0.4*I1z
  verbose          1101
  variable rf      100000
  variable rf2     1
}

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

  matrix set detect elements {{2 3}}

  acq

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

proc main {} {
  global par

  fsave [fsimpson] $par(name).fid

  puts "Larmor frequency (Hz) of 23Na: "
  puts [resfreq 23Na $par(proton_frequency)]
}