#include <bits.gc>
;#include <bits_750.gc>

; NOTE !!!
; uncomment use bits_750.gc above only for the 750!!! 

; NOTE !!!  set ds to be (2n+1)*ns !!! (i.e. odd multiple of ns) 
;DS = (2n+1)*NS (cannot be 0!)

;#define PULSECHECK		

; used to meas. ch RDCs
;3D version, modulated in 2d dimension by JCH; needs REAL FT in F2
;bodenhausen-ruben  constant time in t1
;constant time tuned to 1/Jcc , modified for hsec refocusing pulse (AB 11/97)
;uses extra 80 W amplifier for dlo:
;GV & ABX 14-1-91 JMR, in press
;  bodenrubenct.gv

;l3=128, l5=1
;p1 = high power 1H
;p18 = low power 1H @pl2 (turned off, pl2=120dB)

;o3p=56ppm for CAHA, 18 ppm for Me 
;p3 = high power 13C @pl3
;p13 = p3*2 only for pulsecheck 
;p6 spnam1 hsec.3 @ sp1 (~-2dB; p180=500us@600 Mhz), pl16=120dB 
;p16 = 180CO spnam2=fasinc.256 @sp2, spoffs2=121ppm, pl17=120dB

;pcpd2 low power dec on 15N (p90=230u) @pl30 cpdprg2 waltz16
;pcpd3 low power dec on 13C (p90=90u) garp @pl31 

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

; NOTE gradfile brct3djz.r, cnst20=80, cnst21=-80

;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

"d26=p3-p1"
"d25=0.5*d26"
"d2=1.3m"
"d11=50m"
"d12=1m"
"d15=p16-p1*2"
"d19=5u"
"d22=in22*10-p21-5u-p16"
"d20=10u"
"d21=14.005m"  

;Gradient Pulses
"p20=3.5m"
"p21=3.975m" ; optimize length for max signal
"p22=0.7m"
"p23=0.3m"
"p24=4.0m"

"d23=d21-d20-p16-d19-p21+p3*1.26"

"d16= d21 + p3*1.26 - in22*10 -p1*2 +5u"
"d24=d2-p22"
"d27=d2-p23"

;in19=in23=in16=in22=730u
;ysw=1/(in20+in21)=1/220u
;zsw=sw(J)=1/(in22+in16)=1/1460 (intensity is J-modulated in 20 real planes)

     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 pl1:H
     10u pl3:C1
     10u LOCK_ON
     d1 pl30:N
     10u LOCK_OFF
     (p3 ph5):C1
     5u
     GRAD24(cnst20)
     2m
     (p1 ph0):H
     d2
     (d26 p1*2 ph0):H (p3*2 ph1):C1
     d2  
     (p1 ph11):H
     5u
     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 == 11" 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 ph0):H
     d16 pl16:C1
30   5u
     (p6:sp1 ph13):C1
     8u
     d21 pl3:C1
     5u do:N
;end CT 
     (d25 p1 ph0):H (p13 ph3):C1
;     (d25 p1 ph0):H (p3 ph3):C1
     5u
     GRAD22(cnst20)
     d24 
     (d26 p1*2 ph4):H (p3*2 ph5):C1
     5u 
     GRAD23(cnst21)
     d27 pl31:C1
     go=2 ph31 cpd3: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 == 11" goto 88
	10u iu5
	10u id19
	10u dd23
     goto 99
88	10u id22
	10u dd16
99      1m
	lo to 5 times 10 ; change back to 10 !!! 

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 
ph6=0  ; adjust for power level dl1 180 pulse
;ph13=(360) 1
ph13=(360) 0
ph31=0 2 1 3 2 0 3 1