# onepowderMAS-delay.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
    delay 1.0
    acq -y
  }
}

proc main {} {
  global par

  fsave [fsimpson] $par(name).fid

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

# SIMP
# NP=21
# SW=1000000
# TYPE=FID
# DATA
# 0 0
# 0.182423558 -0.0356466173
# 0.0989424049 -0.0578138302
# -0.0905746502 -0.0303358849
# -0.0995905549 -0.00300984106
# 0.0637974991 -0.0142674177
# 0.129704616 -0.049368046
# 0.0234087454 -0.0532752705
# -0.0800400754 -0.029275804
# -0.0394933071 -0.0137163123
# 0.0702655428 -0.0250494197
# 0.0709753475 -0.0479183045
# -0.0126538606 -0.0414017975
# -0.0487570241 -0.0256550845
# 0.00684655207 -0.028157438
# 0.0564105692 -0.0392794355
# 0.0369406466 -0.0352583419
# -0.0167652653 -0.0294893819
# -0.018874834 -0.0204315996
# 0.0129901036 -0.0262808214
# 0.034967544 -0.0365599687
# END