#include "bits.nt"
;2D CT-TROSY-JNC'
;James Chou and Ad Bax   July 12, 2000 
;Interleaved 2D CT-TROSY with and without 99.3 ms of J(NC')+D(NC') dephasing
;for quantitative measurement of J(NC')+D(NC').

;--RF Pulses---------------
;p1 = 90 deg (10us) 1H pulse @pl1
;p2 = 90 deg @pl2 1H; should be ~1ms
;p6 = sp0 90 sinc1.0 on C' (~100 us at 800 MHz) 
;p7 = 90 deg (50us) 15N pulse @pl7
;p8 = sp3 sinc1.0 on CA (~100 us at 800 MHz), centerlobe sinx/x 
;p16= sp2 180 hsec.3 on C' (~1 ms at 800 MHz), hyperbolic secant 

;--Gradients---------------
;p21=1m, gp1=5%x       sine.50     ;sine bell shaped gradient
;p21=1m, gp3=30%z      sine.50

;--Options-----------------
#define INTERLEAVE

;--Evolutions--------------
define loopcounter NLOOP
"NLOOP=184"                        ;max = 184 for 500u increment in 15N dim
;15N evolution    
;in20=2*in21=2*in22=2*in23=2*in24
;for the current experiment,
;in20 = 250u
;in21 = in22 = 125u
;in23 = in24 = 125u
;Total = 500u

;--Delays------------------
"d3=2.5m-p21"
"d4=2.6m-p21-p2"
"d11=50m"
"d16=p16"         
"d20=5u"          
"d21=24m"
"d22=24m"
"d23=24m" 
"d24=24m-p2-16u"
"d26=p7-p1"
;
1       ze
        1m RESET
2       d11 do:N
3       3m
7       5u do:C1
        d1 LOCK_ON
        1m LOCK_OFF
        10u pl1:H
        10u pl7:N
        10u pl3:C1
;------- start 90-degree on hn -----------------------
        (p1 ph0):H                     ;the first 1H 90
	3u
	p21:gp1                        ;1m @ 5%x
        d3
        (d26 p1*2 ph0):H (p7*2 ph0):N
	p21:gp1                        ;1m @ 5%x
	d3
        (p1 ph1):H             
	10u pl2:H
        (p2 ph10):H                    ;selective water 90 flip back ~1m, -x
        5u  
;------- INEPT to 15N -----------------------
        if "l1==1" goto 78             ;select the appropriate TROSY phase 
	(p7 ph17):N                    ;ph17 = 1 2 3 0  --|
        goto 79                        ;                  |
78      (p7 ph7):N                     ;ph7 = 1 0 3 2   --|
;------- Start CT 15N evolution--------------
79      d20                            ;increment d20 = 5u
	if "l2==2" goto 20             ;interleaving expt with 0 and 99.3 msec NC' J-evolution
	(p16:sp2 ph0):C1                
        d21                            ;decrement d21 = 24m
        (p8:sp3 ph0):C1                
        d22                            ;increment d22 = 24m
        d16                             
        5u
        goto 21
20      d16                            
        d21
        (p8:sp3 ph0):C1               
        d22                            ;increment d22 = 24m
        5u
        (p16:sp2 ph0):C1
21	(p7*2 ph0):N
        5u
        d16
        d23                            ;decrement d23 = 24m
        (p8:sp3 ph0):C1
        5u
        d16
        d24                            ;decrement d24 = 24m-p2-16u
        5u pl2:H
        (p2 ph0):H                     ;selective water 90, ~1m, x
        6u
        5u pl1:H
;------- End CT 15N evolution--------------
	if "l1==1" goto 88             ;select the appropriate TROSY phase
	(p1 ph12):H                    ;ph12 = 1  --|
	goto 89                        ;            |
88	(p1 ph2):H                     ;ph2 = 3   --|
89	5u
	p21:gp1                        ;1m @ 5%x
	d3
        (d26 p1*2 ph0):H (p7*2 ph0):N
	5u
	p21:gp1                        ;1m @ 5%x
	d3
        if "l1==1" goto 98             ;select appropriate TROSY phase cycle
        (p1 ph0):H (p7 ph14):N         ;DOUBLE 90  ph14 = 1  --|
	goto 99                        ;                       |
98      (p1 ph0):H (p7 ph4):N          ;DOUBLE 90  ph4 = 3   --|
99      5u
	p21:gp3                        ;1m @ 30%z
	d4 pl2:H
	(p2 ph15):H
	3u
	5u pl1:H
        (d26 p1*2 ph0):H (p7*2 ph0):N
	8u pl2:H
	(p2 ph15):H
	5u
	p21:gp3                        ;1m @ 30%z
	d4 pl30:C1
	(p7 ph0):N
        go=2 ph31 cpd1:C1              ;CPD for C' during acquisition
        10u do:C1
        1m LOCK_ON
        d11 wr #0 if #0 zd
;
#ifdef INTERLEAVE
        10u iu2                        ;increment loopcounter l2
        lo to 3 times 2         
        10u ru2                        ;reset loopcounter l2
#endif
;
0.1m iu1                               ;increment loopcounter l1
lo to 3 times 2
10u id20
10u dd21
10u id22
10u dd23
10u dd24
10u ru2
10u ru1                                ;reset loopcounter l1
lo to 3 times NLOOP 
10u rd20
10u rd21
10u rd22
10u rd23
10u rd24
1m do:C1
1m do:N
1m RESET
1m
exit    
        
ph0=0
ph1=1
ph2=3
ph12=1
ph4=3
ph14=1
ph15=(360)178
ph10=2
ph7=1 0 3 2
ph17=1 2 3 0  
ph31=1 2 3 0             ;receiver phase