# one180.in (SIMPSON 3.1)

#one pulse 180°, 
#Mz versus offset varying from 0 to np kHz in steps of 1 kHz

spinsys {
  channels 29Si
  nuclei   29Si

  shift 1 0 0 0 0 0 0
}

par {
  proton_frequency 400e6
  variable tsw     0.5
  sw               1.0e6/tsw
  np               71
  crystal_file     alpha0beta0
  gamma_angles     1
  start_operator   I1z
  detect_operator  2*I1z
  verbose          1101
  variable rf      50000
  variable p1      10
}

proc pulseq {} {
  global par

  for {set i 0} {$i < $par(np)} {incr i} {
    
    offset [expr $i*1000]

    pulse $par(p1) $par(rf) x

    acq

    reset
  }
}

proc main {} {
  global par

  fsave [fsimpson] $par(name).fid

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


#SIMP
#NP=71
#SW=2000000
#TYPE=FID
#DATA
#-1 0
#-0.999200123 0
#-0.996801961 0
#-0.992809924 0
#-0.987231348 0
#-0.980076481 0
#-0.971358454 0
#-0.961093259 0
#-0.949299702 0
#-0.935999367 0
#-0.921216559 0
#-0.904978254 0
#-0.88731403 0
#-0.868256003 0
#-0.847838748 0
#-0.826099223 0
#-0.803076681 0
#-0.778812578 0
#-0.753350481 0
#-0.726735964 0
#-0.699016504 0
#-0.670241373 0
#-0.640461523 0
#-0.609729473 0
#-0.578099183 0
#-0.545625939 0
#-0.512366221 0
#-0.47837758 0
#-0.443718502 0
#-0.408448286 0
#-0.372626901 0
#-0.336314859 0
#-0.299573079 0
#-0.262462752 0
#-0.225045204 0
#-0.187381765 0
#-0.149533636 0
#-0.111561752 0
#-0.0735266546 0
#-0.0354883617 0
#0.00249376064 0
#0.0403611227 0
#0.078056034 0
#0.115521822 0
#0.15270295 0
#0.189545127 0
#0.225995421 0
#0.26200236 0
#0.297516033 0
#0.33248819 0
#0.366872329 0
#0.400623784 0
#0.433699808 0
#0.466059645 0
#0.497664607 0
#0.528478131 0
#0.558465846 0
#0.587595619 0
#0.61583761 0
#0.643164307 0
#0.669550564 0
#0.69497363 0
#0.719413172 0
#0.742851293 0
#0.765272538 0
#0.786663905 0
#0.807014838 0
#0.826317221 0
#0.844565365 0
#0.861755988 0
#0.877888188 0
#END