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

# approach: 12-phase cycling the first pulse 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   2*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) x
    acq                      $par(ph31)
    reset
    prop 1
    
  }
}

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

    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.22044605e-16
# 0.0339282577 8.32667268e-17
# 0.195744737 7.77156117e-16
# 0.566434533 -4.4408921e-16
# 1.09245834 -2.22044605e-16
# 1.63751862 3.33066907e-16
# 2.10867748 0
# 2.46906347 1.38777878e-16
# 2.71331348 4.4408921e-16
# 2.86908614 -1.11022302e-16
# 2.95538752 5.55111512e-17
# 2.97951815 1.94289029e-16
# 2.96661858 -5.55111512e-17
# 2.90320646 -4.99600361e-16
# 2.73287153 2.22044605e-16
# 2.45520283 3.33066907e-16
# 2.12660292 -9.43689571e-16
# 1.78857522 -1.22124533e-15
# 1.46644566 -6.38378239e-16
# 1.18593134 -5.41233725e-16
# 0.968441965 -8.32667268e-16
# 0.845548698 -9.43689571e-16
# 0.859731422 -1.36002321e-15
# 1.03120914 -6.83264209e-16
# 1.33527906 1.38777878e-16
# 1.69308901 -2.77555756e-17
# 1.99106038 6.66133815e-16
# 2.14858566 5.55111512e-17
# 2.14846382 -8.32667268e-17
# 2.00810777 -3.12250226e-16
# 1.76322478 -7.63278329e-16
# 1.45576869 -7.77156117e-16
# 1.13837145 -8.32667268e-16
# 0.884901173 -8.8817842e-16
# 0.72046018 -6.66133815e-16
# 0.574821363 -1.11022302e-15
# 0.377560716 -8.8817842e-16
# 0.140554076 9.36750677e-17
# -0.0765976066 -3.88578059e-16
# -0.222412628 -4.16333634e-16
# END