;2D IPAP-H(N)CO (1H-coupled in F1(C'))
;
;Reference:
;Ottiger, M., Delaglio, F. & Bax, A. 
;J. Magn. Reson. 131, 373-378 (1998)
;
;needed files:
;Gradnt_new.incl        standard Bruker include file

#include "Gradnt_new.incl"

#define HYPERBOLIC	; hyperbolic sec. CA-decoupling 

;Pulse definitions:
;p1	proton 90 at pl1
;p2	proton 90 for WATERGATE at pl2 (~42 dB)
;p3	carbon 90 at pl3
;p4	carbon 90 at pl4 (~23us @ 600MHz, see paper)
;p14	sel. carbon 180 hseq.3  at sp2
;p5	sel. carbon  90 sinc1.0 at sp1
;	pl5 = 120 dB to preset shape on C2 
;p7	nitrogen 90 at pl7
;p31	nitrogen 90 at pl31 for decoupling (cpd2)

"p6=p5*2"


;Gradient pulses (may have to be set manually):
"p20=5.0m" 	; gp0 = -z-x gradient 50%
"p21=1.1m"	; gp1 = +xyz gradient 50%
"p22=0.5m"	; gp2 = -y   gradient 50%
"p23=0.6m"	; gp3 = +x   gradient 50%
"p24=0.7m"	; gp4 = -x   gradient 50%
"p25=0.4m"	; gp5 = +z   gradient 50%


;Delay calculations
#ifdef HYPERBOLIC
   "d0=in0-p4*0.6366-p14/2+10u"	; C' evolution
#else
   "d0=in0/2-p4*0.6366-p7"	; C' evolution
#endif

"d2=2.3m"		; 1/4JNH * 0.85
"d3=13.5m-5.3m"		; 1/4JNC' * 0.8 
"d4=d2-p25-p2-12u"
"d11=50u"
"d12=10u"
"d25=p1*4.4 + 4u"
"d26=p7-p1"
"d27=p5-p7"
"d28=p7-p4"
"d29=p14/2-p7"



	10m RESET
1	d11 ze
        1m ru2
2	d11
        d12
3	20u
4	d12*3
5	10u               
	5u do:N
	d1 LOCK_ON
	10u pl1:H
	10u pl3:C1
	10u pl5:C2
	10u pl7:N
	1m LOCK_OFF
	(p7 ph0):N (p3 ph0):C1
	5u
	p20:gp0	; GRAD20(cnst20)
	5m
;------------- start INEPT to 15N -----------------------
	(p1 ph0):H
	d2
	(d26 p1*2 ph0):H (p7*2 ph0):N
	d2
	(p1 ph2):H
        5u
	p21:gp1	;GRAD21(cnst20)
	1m pl5:C2
;--------------------- INEPT to C' ----------------------
	(p7 ph0):N
	p26:gp6
        1.65m
	if "l2==2" goto 20
	d25                 ;p1*2 comp delay
	goto 21
20	(p1 ph0 2u p1*2.4 ph1 2u p1 ph0):H
21	2.65m
	d3 
	(p6:sp1 ph0):C2 (d27 p7*2 ph9):N
	d3 
	d25
	p26:gp6
	4.3m
	if "l2==2" goto 30
        (p7 ph30):N 
	goto 31
30	(p7 ph8):N 
31      (p1*2 ph0):H
	p22:gp2	
	0.5m pl4:C2
;--------------- start evolution on C' ------------------
	(p4 ph10):C2 
#ifdef HYPERBOLIC
	d0 pl5:C1
        (d29 p7*2 ph11):N (p14:sp2 ph11):C1
#else
	d0 pl4:C1
        (p7*2 ph11):N (d28 p4*2 ph11):C1
#endif
	d0 ; pl5:C2
	(p4 ph0):C2 
	5u
	p23:gp3	; GRAD27(cnst20)
	0.5m
;----------------- INEPT back to 15N --------------------
	(p7 ph1):N
	5.3m
	d3 pl5:C2
	(p6:sp1 ph0):C2 (d27 p7*2 ph0):N
	d3 
	5.3m 
	(p7 ph0):N
	5u
	p24:gp4	; GRAD27(cnst20)
	0.5m pl1:H
;------------- INEPT back to H, WATERGATE ---------------
	(p1 ph0):H
	5u
	p25:gp5	; GRAD22(cnst21)
	50u
	d4 pl2:H
	(p2 ph14):H
	2u
	5u pl1:H
	(p1*2 ph15):H
	2u
	5u pl2:H
	(p2 ph14):H  (p7*2 ph0):N
	5u
	p25:gp5	; GRAD22(cnst21)
	d4 pl31:N
	go=2 ph31 cpd2:N
	d11 do:N wr #0 if #0 zd
	d12 iu2
     lo to 3 times 2
	10u ip10
	10u ru2
     lo to 4 times 2
	d12 id0
	d12 ip31
	d12 ip31
     lo to 5 times l3
	d12 do:H
	d12 do:N
	d12 do:C1
	d12 do:C2
	1m LOCK_ON
	100u RESET
	exit	
	
	
	
ph0=0
ph1=1
ph2=1 3
ph8=0 0 0 0 2 2 2 2 
ph30=1 1 1 1 3 3 3 3 
ph10=0 0 2 2
ph11=0
ph14=2
ph15=0
ph21=0
ph9= 0 0 0 0 0 0 0 0 ;1 1 1 1 1 1 1 1 
;     2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3
ph31=0 2 2 0 2 0 0 2 ; 2 0 0 2 0 2 2 0