;NOESY-P selective CT-NOESY interleaved for measuring H-P couplings
;added dephasing during shaped pulse
;Justin Wu and Ad bax, JBNMR, 19, 367-370 (2001)

#include "Gradnt_new.incl"
; SFO2 = 242.940 for nucleic acids
; This experiment works correctly

;"d20=35.1m"
"d0 = 5u"
"d3=p3"
"d4=p6*0.5-p3"
"d18 = d20-d0" 
"d19=d18-p3*2"

1  ze
2  50m
3  10u ;fq1:f1
   10u pl18:f1
   10u cw:f1
   d1 pl3:f2  
   10u do:f1
;   10u fq1:f1
   10u pl1:f1
   p1 ph1
   d20
   5u pl18:f1
   if "l1==1" goto 20
   (p6:sp0 ph2):H (d4 p3*2 ph0):f2
   10u
   goto 21
20 (p6:sp0 ph2):H 
   10u
21 10u
   d18
   3u pl1:f1
   (p1 ph11 3u p1*2 ph12 3u p1 ph11):f1 (p3*2 ph0):f2
   if "l1==1" goto 10
   3u
   (p6:sp0 ph15):H (d4 p3*2 ph0):f2
   6u
   goto 11
10 3u
   (p6:sp0 ph15):H
   6u
11 d0 pl1:f1
   p1 ph0
   d5   ;NOE mixing
   p1 ph3
   go=2 ph31
   50m wr #0 if #0 zd
   10u iu1
   lo to 3 times 2
   10u ru1
   10u ip1
   lo to 3 times 2
   10u id0
   10u dd20
   10u ip31
   10u ip31
   lo to 3 times l3
exit

ph0=0
ph1=0 2 
ph2=0 0 1 1 2 2 3 3
ph3=0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3
ph11=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
     1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
     2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
     3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
ph12=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
     2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
     3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
     0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
ph15=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
     1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
     2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
     3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
ph31=0 2 2 0 1 3 3 1 2 0 0 2 3 1 1 3