;3D J-modulated [13C,1H] ct HSQC
;
;Reference:
;M. Ottiger, F. Delaglio, J.L. Marquardt, N. Tjandra and A. Bax
;J. Magn. Reson. 1998, 134 (1998)
;
;NOTE need ds to be (2n+1)*ns !!!!!!!!
;16 plains, modulated by JCH; needs REAL FT in F2
;constant time [13C,1H] in t1
;constant time tuned to 1/Jcc , modified for hsec refocusing pulse

#include "Gradnt_new.incl"
;#define PULSECHECK		

;p1   high power 1H
;p18  low power 1H @ pl2
;p3   high power 13C
;p4   low power 13C null at CO
;p6   hyperbolic secant pulse
;p16  180CO sinc1.0 
;p31  low power dec on 13C
;p30  low power dec on 15N

#ifdef PULSECHECK
"p13=p3*2"			; 180 degree, zero signal
#else
"p13=p3"			; 90 degree
#endif

;defining all of the gradients to be used
#define GRAD20(ampl) p20:ngrad:c34
#define GRAD21(ampl) p21:ngrad:c34
#define GRAD22(ampl) p22:ngrad:c34
#define GRAD23(ampl) p23:ngrad:c34
#define GRAD24(ampl) p24:ngrad:c34

;Gradient Pulses
"p20=3.5m"
"p21=3.978m"
"p22=0.7m"
"p23=0.3m"
"p24=4.0m"

;Set increments as follows:
;in22=in19=in23=in16=912.5u

"d2=1.3m"
"d11=50u"
"d12=10u"
"d15=p16-p1*2"
"d19=5u"
"d20=10u"
"d21=14.005m"  
"d16=d21+p3*1.26-in22*8-p1*2-d20"	; adapt number of increments
"d22=in22*8-p21-p16"			; if changing in22,19,23,16
"d23=d21-d20-p16-d19-p21+p3*1.26"
"d24=d2-p22"
"d26=p3-p1"
"d25=0.5*d26"
"d27=d2-p23"



     1m RESET
1    ze 
     1m pl2:H
     1m ru5
     (p18 ph0):H
     1m ru1
2    d11 do:C1 do:N
3    d12*3
4    d12*3
5    10u
     10u pl1:H
     10u pl3:C1
     10u LOCK_ON
     d1 pl30:N
     10u LOCK_OFF
     (p3 ph5):C1
     GRAD24(cnst20)
     2m
     (p1 ph0):H
     d2
     (d26 p1*2 ph0):H (p3*2 ph1):C1
     d2  
     (p1 ph11):H
     2u
     GRAD20(cnst20)		;need huge gradient
     1m
   if "l1 == 1" goto 10
     (p3 ph2):C1
     goto 12
10     (p3 ph12):C1     
;start CT
12   5u
     5u pl16:C1
	5u
     (p6:sp1 ph0):C1           ;Hsec compensation refoc. pulse
     10u
     5u pl17:C1
     if "l5 == 9" goto 20
     d20 cpd2:N
     (d19 p1*2 ph0 d15):H (p16:sp2 ph0 d19):C1
     5u
     GRAD21(cnst21)
     d23 pl16:C1
     goto 30
20   d20 cpd2:N
     (p16:sp2 ph0):C1
     5u
     GRAD21(cnst21)
     d22
     (p1*2 ph8):H
     d16 pl16:C1
30   5u
     (p6:sp1 ph13):C1
     8u
     d21 pl3:C1
     2u do:N
;end CT 
     (d25 p1 ph0):H (p13 ph3):C1
     2u
     GRAD22(cnst20)
     d24 
     (d26 p1*2 ph4):H (p3*2 ph5):C1
     2u 
     GRAD23(cnst21)
     d27 pl31:C1
     go=2 ph31 cpd4:C1
     5u iu1
     d11 do:C1 wr #0 if #0 zd
     lo to 3 times 2 
	d12 ru1 
	d12 id20
	d12 dd21
     lo to 4 times l3  
	10u rd20
	10u rd21
	if "l5 == 9" goto 88
	10u iu5
	10u id19
	10u dd23
     goto 99
88	10u
	10u id22
	10u dd16
99      1m
	lo to 5 times 16	; if changing in22,19,23,16: adjust d22,16!!

    d12 LOCK_ON
    d12 do:N
    d12 do:C1

    exit 


ph0=0 
ph1=1
ph11=1 3
ph2=0 0 1 1 2 2 3 3
ph12=0 0 3 3 2 2 1 1
ph3=0
ph4=0  
ph5=0 
ph8=0 2
ph13=(360) 1
ph31=0 2 1 3 2 0 3 1