;pdsdexch.mh 
;2D 13C exchange spectroscopy
;Szeverenyi NM, Sullivan MJ, Maciel GE.
;J. Magn. Reson., volume 47, 462, 1982

;p3 = 1H 90 degree pulse
;p15 = CP contact pulse
;pl2 = 1H CP power 
;pl12 = 1H dec. power during acquisition
;d0 = t1 evolution time
;d2 = exchange mixing

;pl1: 13C CP power
;pl2: 1H excitation power
;pl3: 1H CP power
;pl11: 13C hard pulse power
;pl12: 1H dec. power
;p1: 13C 90
;p2: 13C 180
;p3: 1H 90
;p15: 1H/13C contact time
;p8: TPPM dec. pulse, <= 180 deg

;d3: Hahn echo delay
;d1: recycle delay
;d2: spin diffusion mixing 

  "l1=aq/(p8*2)"
  "p2=2*p1"
  "d3=1s/l31-p2/2"
  "d13=d3-10u"
define loopcounter td1half
  "td1half=td1/2"

1 ze
2 d1  
  10u pl1:f1
  10u pl2:f2
  2u:f1 ph2 		
  2u:f2 ph1 
 
3 p3:f2 ph1			;1H excitation pulse
  0.5u pl3:f2			;1H power for CP                   	
  (p15 ph2):f1 (p15 ph3):f2		;spin lock
  1u cw:f2 pl22:f2 pl11:f1

  d0				;evolution, does not have to be rotor synch.

4 0.5u:f1 ph4			;hard 90 pulse
  p1:f1 ph4
  1u do:f2

  d2				;spin diffusion

  1u cw:f2 pl12:f2
6 0.5u:f1 ph6			;hard 90 flip back pulse
  p1:f1 ph6

  d3:f1 ph9 
  p2:f1 ph9
  d13

  1u:f1 ph0			;dead time, reset RF phase to detect
  1u adc ph31			;start ADC with ph31 signal routing
7 p8:f2 ph13			;p8 < 180 degree
  p8:f2 ph14			
lo to 7 times l1
  aq do:f2			;data sampling
  rcyc=2 			;do NS scans

;2D				;States' detection
  50m wr #0 if #0 ip4 zd	;write file to disk, increment, zero data
  lo to 2 times 2		;cosine/sine detection
  10m ip4 id0 	 		;increments loop by a defined number
  10m ip4
  lo to 1 times td1half

exit

ph0= 0
ph1= +x +x +x +x -x -x -x -x	;1H excitation
ph3= +y  			;1H spin lock
ph2= +x +x +x +x +x +x +x +x 	;13C spin lock
     +y +y +y +y +y +y +y +y
     -x -x -x -x -x -x -x -x 
     -y -y -y -y -y -y -y -y	
ph8= +x +x +x +x -x -x -x -x 
     +y +y +y +y -y -y -y -y 
ph4= +y +y +y +y +y +y +y +y	;storage, cosine detection 
     -x -x -x -x -x -x -x -x 
     -y -y -y -y -y -y -y -y	
     +x +x +x +x +x +x +x +x
ph6= +x +y -x -y +x +y -x -y	 ;read out, full CYLCOPS
     +y -x -y +x +y -x -y +x 
     -x -y +x +y -x -y +x +y
     -y +x +y -x -y +x +y -x 
ph9= +y -x -y +x +y -x -y +x	;Hahn echo phase
     -x -y +x +y -x -y +x +y

ph13= (360) 14
ph14= (360) 0
ph31= 1 2 3 0 3 0 1 2		;SQ detection
      2 3 0 1 0 1 2 3
      3 0 1 2 1 2 3 0
      0 1 2 3 2 3 0 1