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

# approach: 12-phase cycling the first pulse and the receiver

# spin-3/2 central-transition echo amplitude optimization
# versus the second-pulse duration
# in two-pulse phase-modulated -3Q-MAS sequence,
# the first-pulse duration p1 = 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 p1      5
}

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

  acq $par(ph31)
  pulse $par(p1) $par(rf) $par(ph1)

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

 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
# 2.09922822 -1.66533454e-16
# 6.34304518 1.38777878e-17
# 9.00910329 2.22044605e-15
# 8.66050705 2.77555756e-15
# 6.79038071 3.66373598e-15
# 5.45470498 3.55271368e-15
# 5.19311644 2.47024623e-15
# 5.00222595 1.47104551e-15
# 4.01503049 1.49880108e-15
# 2.86908614 4.99600361e-16
# 2.84886138 5.55111512e-16
# 4.09893333 2.22044605e-16
# 5.56975162 9.99200722e-16
# 6.32589372 1.33226763e-15
# 6.18966891 2.10942375e-15
# 5.44609664 1.55431223e-15
# 4.43025526 1.94289029e-15
# 3.43797585 1.44328993e-15
# 2.90616202 1.44328993e-15
# 3.35147662 6.66133815e-16
# 4.8430532 8.04911693e-16
# 6.63926553 1.47104551e-15
# 7.63088381 1.56819002e-15
# 7.34253752 1.34614542e-15
# 6.37997143 1.99840144e-15
# 5.62709452 1.88737914e-15
# 5.26425859 1.99840144e-15
# 4.95258841 1.11022302e-15
# 4.60580358 1.63757896e-15
# 4.57054021 1.31838984e-15
# 5.20369812 1.63410951e-15
# 6.34406795 1.27675648e-15
# 7.29706122 1.52655666e-15
# 7.57284428 1.44328993e-15
# 7.30935259 1.22124533e-15
# 6.79580111 9.15933995e-16
# 6.23946044 1.71390679e-15
# 6.00356474 1.08246745e-15
# 6.32093294 1.66533454e-15
# 6.98926573 1.80411242e-15
# END