;lgcp2d.mh
;Measure X-H dipolar coupling as the indirection dimension of a 2D
;magic-angle spin lock on 1H for CP without spin diffusion
;wr generally higher than 10 kHz
;Hong M, Yao XL, Jakes K, Huster D.
;J. Phys. Chem. B, volume 106, 7355, 2002

;p3: 1H excitation pulse
;p21: 1H 55 deg. pulse
;p15: CP contact time
;p8: 1H TPPM dec. pulse
;pl2: 1H excitation power
;pl3: 1H LG-CP power
;pl1: 13C CP power, w1H,eff (+/-) wr
;FQ2: LG-CP
;FQ3: on resonance
;d1: recycle delay
;l31: spinning speed in Hz


  "p21=p3*353/900"
  "l1=aq/(p8*2)"
  "d3=1s/l31"

1 ze

2 d1 
  10u pl1:f1
  10u pl2:f2		
  2u:f2 ph1
  2u:f1 ph2

3 p3:f2 ph1			;90 deg pulse
  p21:f2 ph3			;35 deg pulse
  0.5u fq2:f2 			;put 1H off resonance
  0.5u pl3:f2
  (p15 ph2):f1 (p15 ph10):f2	;CP, t1	

  1u fq3:f2		;put 1H back on res for decoupling
  1u cw:f2 pl12:f2 	;pl11:f1

  1u:f1 ph0		;dead time delay, reset RF phase to detect
  2u adc ph31		;start ADC with ph31 signal routing
6 p8:f2 ph13		;TPPM decoupling
  p8:f2 ph14
lo to 6 times l1
  aq do:f2
  d3
  rcyc=2
;  wr #0

;2D				;States' detection
  200m wr #0 if #0 zd		;write file to disk, increment, zero data
  10m ipu15
  lo to 1 times td1

exit

ph0= +x
ph1= -y	+y 		;1H excitation
ph3= +y			;35 deg pulse
ph10= +x 		;1H LG-CP
ph2= +x +x +y +y -x -x -y -y 		;X CP spin lock
ph4= -x -x -y -y +x +x +y +y
ph13=(360) 0
ph14=(360) 10
ph31= 2 0 3 1 0 2 1 3