# na-two-pulse-24-phase-cycling-p1.in

# approach: 24-phase cycling the two pulses and the receiver

# spin-3/2 central-transition echo amplitude optimization
# versus the first-pulse duration
# in two-pulse phase-modulated -3Q-MAS sequence,
# the second-pulse duration p2 = 5 micro seconds

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

par {
  proton_frequency 400e6
  spin_rate        10000
  variable tsw     0.5
  sw               1.0e6/tsw
  np               41
  crystal_file     rep100
  gamma_angles     10
  start_operator   I1z
  detect_operator  I1c
  verbose          1101
  variable rf      100000
  variable p2      5
}

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

  acq $par(ph31)
  for {set i 1} {$i < $par(np)} {incr i} {
    pulse $par(tsw) $par(rf) $par(ph1)
    store 1
    pulse $par(p2)  $par(rf) $par(ph2)
    acq                      $par(ph31)
    reset
    prop 1    
  }
}

proc main {} {
  global par
  
  foreach p {{   0    0  270}\
             {  60    0   90}\
             { 120    0  270}\
             { 180    0   90}\
             { 240    0  270}\
             { 300    0   90}\
             {   0   90   90}\
             {  60   90  270}\
             { 120   90   90}\
             { 180   90  270}\
             { 240   90   90}\
             { 300   90  270}\
             {   0  180  270}\
             {  60  180   90}\
             { 120  180  270}\
             { 180  180   90}\
             { 240  180  270}\
             { 300  180   90}\
             {   0  270   90}\
             {  60  270  270}\
             { 120  270   90}\
             { 180  270  270}\
             { 240  270   90}\
             { 300  270  270}} {
    set par(ph1)  [lindex $p 0]
    set par(ph2)  [lindex $p 1]
    set par(ph31) [lindex $p 2]

    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 23Na: "
  puts [resfreq 23Na $par(proton_frequency)]
}

# SIMP
# NP=41
# SW=2000000
# TYPE=FID
# DATA
# 0 0
# 0.00125679417 -2.77555756e-17
# 0.0339282577 -2.35922393e-16
# 0.195744737 -1.27675648e-15
# 0.566434533 -4.4408921e-16
# 1.09245834 -7.77156117e-16
# 1.63751862 -1.11022302e-16
# 2.10867748 -1.66533454e-16
# 2.46906347 -8.32667268e-17
# 2.71331348 3.33066907e-16
# 2.86908614 -2.49800181e-16
# 2.95538752 -2.77555756e-17
# 2.97951815 -1.11022302e-16
# 2.96661858 -1.94289029e-16
# 2.90320646 -7.21644966e-16
# 2.73287153 -2.77555756e-16
# 2.45520283 -4.99600361e-16
# 2.12660292 -8.32667268e-16
# 1.78857522 -6.66133815e-16
# 1.46644566 2.01227923e-16
# 1.18593134 -2.15105711e-16
# 0.968441965 -2.08166817e-16
# 0.845548698 -2.77555756e-16
# 0.859731422 -6.9388939e-16
# 1.03120914 -4.56015434e-16
# 1.33527906 -3.88578059e-16
# 1.69308901 -3.46944695e-16
# 1.99106038 -1.11022302e-16
# 2.14858566 -1.66533454e-16
# 2.14846382 -8.32667268e-17
# 2.00810777 -2.53269627e-16
# 1.76322478 -1.80411242e-16
# 1.45576869 -4.16333634e-17
# 1.13837145 -4.99600361e-16
# 0.884901173 -4.71844785e-16
# 0.72046018 -7.21644966e-16
# 0.574821363 -9.71445147e-16
# 0.377560716 -7.91033905e-16
# 0.140554076 -2.68882139e-16
# -0.0765976066 -2.91433544e-16
# -0.222412628 -2.08166817e-16
# END