;hchetcorfil2.mh
;13C-1H HETCOR with C-H dipolar filter
;to remove directly bonded 13C-1H cross peaks.
;Yao XL, Schmidt-Rohr K, Hong M. J. Magn. Reson., volume 149, 139, 2001

;p2: 13C 180 pulse
;p3: 1H excitation
;p6: 1H 180 pulse
;p8: 1H TPPM pulse
;p9: one tr for CH filtering
;p10: t360*960, filter 
;p11: first t1 time=0.5us
;p12: t360*960, t1
;p13: t1 dwell time for 1H CS evolution
;p15: CP contact
;p21: 55 deg pulse
;p22: 35 deg pulse
;d1: recycle delay
;d10: filter time

;pl1: 13C spin lock
;pl2: 1H hard pulses
;sp1: FSLG during filter
;sp2: FSLG during t1
;sp3: ramp CP
;pl3: 1H spin lock
;pl11: 13C hard pulse
;pl12: TPPM decoupling during acq

define loopcounter td1half
  "td1half=(td1/2)"		;for STATES-TPPI
  "p22=(p3*353)/900"
  "p21=(p3*547)/900"		;magic-angle pulse
  "l1=aq/(p8*2)"
  "d11=1s/l31-p2-d10"

  100m
  "p10=p9;"       
  "p12=p11;"			;initial t1 value, should be 0.5u

1 ze                             ;zero data

12 d1 do:f2 			;recycle time, decoupler off
  10u pl11:f1 			;set power level of f1
  p1:f1 ph20			;dstroy 13C magnetization
  5m
  p1:f1 ph21
  5m
  p1:f1 ph20			;dstroy 13C magnetization
  5m
  p1:f1 ph21
  5m

  10u reset:f2 			;reset F2 DDS
  1u fq2:f2 			;offset during evolution
  10u pl2:f2 			;set power level of f2
  10u pl11:f1 			;set power level of f1

14 p3:f2 ph1			;excitation pulse
  p21:f2 ph3       		;55 deg pulse

16 (p10:sp1 ph4):f2  (d10 p2 ph5 d11):f1 	;tr
  (1u ph6):f2 (1u pl2):f2
  (p6 ph6):f2 		;(p2 ph9):f1  		;180 pulse
17 (p10:sp1 ph4):f2 (d10 p2 ph15 d11):f1 	;tr

  (p12:sp2 ph4):f2				;t1 
 
  0.5u pl2:f2
  (p22 ph7):f2 			;35 deg pulse. to -z
  (p3 ph8):f2			;90 deg pulse, to +/-x
  (p22 ph10):f2			;35 deg pulse, +y

  1u fq3:f2 			;LG-CP
  1u pl3:f2 pl1:f1
18 (p15 ph2):f1 (p15 ph11):f2	;CP

  1u fq2:f2
  1u cw:f2 pl12:f2
  1u:f1 ph0			;dead time delay, reset RF phase to detect
  2u adc ph31			;start ADC with ph31 signal routing
19 p8:f2 ph13			;TPPM decoupling
  p8:f2 ph14
lo to 19 times l1
  aq do:f2			
  rcyc=12 			;do NS scans

;2D cosine detection
22  100m wr #0 if #0 zd		;write  and increment file
 "p12=p12+p13;"			;calculate next t1 period
  1m
lo to 12 times td1		;number of rows

exit                       

ph0= 0
;ph1= 0			;1D sine slice, or 2D States-TPPI
ph1= 1 			;cosine det. for TPPI version
ph3= 3			;1H 55 deg. pulse
ph4= 0 			;1H FSLG decoupling phase
ph5= 0 2		;13C 180 
ph15= 2 0
ph9= 0
ph6= 1			;1H 180 
ph7= 3 			;1H 35 deg. pulse
ph8= 1 3		;1H 90 deg. pulse
ph10 = 1 
ph11= 0			;1H CP
ph2= 0 0 2 2 1 1 3 3	;13C CP
ph13=(360) 0
ph14=(360) 10
ph20=0			;destroy 13C magnetization
ph21=1			;destroy 13C magnetization
ph31=0 2 2 0 1 3 3 1	;detection