;nchetcoreven.mh
;2D 15N-13C chemical shift correlation
;Hong M, Griffin RG. J. Am. Chem. Soc., 120, 7113, 1998

;p3: proton excitation
;p15: contact time
;p2: 13C 180 pulse
;p1: 13C 90 pulse
;p12: 15N 180 pulse
;p11: 15N 90 pulse
;p8: TPPM 180 pulse length
;p18: 1H TPPM 180 pulse during REDOR
;d0: 15N evolution time
;pl1: 13C CP power
;pl2: 1H excitation power
;pl12: dec. during acq 
;pl22: dec. during REDOR
;pl3: 1H CP power
;pl11: 13C hard power for REDOR
;pl21: 13C hard power for 90 pulses
;pl4: 15N hard power
;l0: REDOR time, l0*tr = exc
;l1: TPPM cycles
;l31: spinning speed

define loopcounter td1half
   "td1half=td1/2"
   "d3=2s/l31-p1*0.5-p11*1.5-1.5u"	;z-filter
   "d21=(p11-p1)/2"
   "d23=0.25s/l31-(p1/2)"
   "d24=0.25s/l31-(p11/2)-0.5u"		;before 15N pulse
   "d25=0.25s/l31-1u"			;REDOR
   "d26=0.25s/l31-p12/2"		;REDOR
   "d27=0.25s/l31-p2/2"			;REDOR
   "d2=0.5s/l31"			;make up one rotor period
   "l1=aq/(p8*2)"			;TPPM


1 ze
2 d1 		
  10u pl1:f1
  10u pl2:f2			;set decoupler power,1H
  10u pl4:f3			;3rd nucleus power,often 15N
  2u:f2 ph1
  2u:f1 ph2

  1m rpp5			;reset phase lists for redor

  p3:f2 ph1			;90 degree proton pulse
  0.5u pl3:f2
  (p15 ph2):f1 (p15 ph3):f2	
  1u cw:f2 pl22:f2

  d25 pl11:f1 				;tr/4			
3 d26:f3 ph5
  p12:f3 ph5^       		;15N 180 deg. pulse
  d26
lo to 3 times l0		;an even number

  d27:f1 ph6			;tr/4
  d2				;tr/2
  p2:f1 ph6			;13C 180 deg. pulse
  d2
  d27

4 d26:f3 ph5
  p12:f3 ph5^      		;15N 180 deg. pulse
  d26
lo to 4 times l0
  0.5u pl21:f1
  (d24 ph9):f1 (d24 ph10):f3	;phase switch

  (d21 p1 d21 ph9):f1 (p11 ph10):f3	;simultaneous 90 pulses
  
  d0				;t1

  p11:f3 ph4			;15N 90
  d3 do:f2			;z-filter
  1u cw:f2 pl22:f2
  p1:f1 ph11

  d23 pl11:f1 				;tr/4			
6 d26:f3 ph5
  p12:f3 ph5^       		;15N 180 deg. pulse
  d26
lo to 6 times l0		;an even number

  d27:f1 ph8			;tr/4
  d2
  p2:f1 ph8			;13C refocus
  d2
  d27

7 d26:f3 ph5
  p12:f3 ph5^      		;15N 180 deg. pulse
  d26
lo to 7 times l0
  d25

  1u pl12:f2
  1u:f1 ph0		;dead time delay, reset RF phase to detect
  1u adc ph31		;start ADC with ph31 signal routing
8 p8:f2 ph13		;1H TPPM dec. 
  p8:f2 ph14
lo to 8 times l1
  aq do:f2		;allow time for data sampling
  rcyc=2 		;do NS scans;
;  wr #0

;2D				;States' detection
  100m wr #0 if #0 ip4 zd
  lo to 2 times 2
  10m ip4 id0
  10m ip4
  lo to 1 times td1half

exit

ph0= 0
ph1= +x -x			;1H excitation
ph2= +x +x +y +y -x -x -y -y  	;13C CP
ph3= +y				;1H CP

ph5= +x +y +x +y +y +x +y +x 	;15N REDOR pulses  
     -x -y -x -y -y -x -y -x	
ph6= +x +x +y +y -x -x -y -y	;13C refocus 	
ph8= +y -y -x +x -y +y +x -x	;13C refocus,

ph9= +x +x +y +y -x -x -y -y	;13C read in 
ph10= +x +x +y +y -x -x -y -y 	;15N read-in, 13C read out
ph11= -x -x -y -y +x +x +y +y 	;13C flip back
ph4= -x -x -y -y +x +x +y +y	;15N read-out
     +x +x +y +y -x -x -y -y

ph13= (360) 14
ph14= (360) 0
ph31= 0 2 1 3 2 0 3 1  		;does not accept #4!!
      2 0 3 1 0 2 1 3