;IPAP-[15N,1H]-HSQC (1H-coupled in F1(N))
;
;Reference:
;Ottiger, M., Delaglio, F. & Bax, A. 
;J. Magn. Reson. 131, 373-378 (1998)
;
;needed files:
;Gradnt_new.incl        standard Bruker include file
;FQ2LIST                frequency list (CA - CO)

#include "Gradnt_new.incl"

;#define CARBON_LABEL	; uncomment if also 13C labelled

;Pulse definitions:
;p1	proton 90 at pl1
;p2	proton 90 for WATERGATE at pl2 (~42 dB)
;p5	carbon 180 at pl5 (~46us @ 600MHz, see paper)
;p7	nitrogen 90 at pl7
;p31	nitrogen 90 at pl31 for decoupling (cpd2)

;Gradient pulses (may have to be set manually):
"p20=2.0m"	
"p11=1.0m"
"p12=0.4m"

#ifdef CARBON_LABEL
"d0=in0*0.5-4u-p5*2.0-p7*0.63"	;nitrogen incremental delay
#else
"d0=in0*0.5 - p7*0.63 -p1"
#endif

"d3=2.3m"		;1/4JNH  * 0.85
"d4=d3-p12-2u"	
"d5=d3-p12-p2-14u"
"d11=50u"		;write delay
"d12=10u"		;increment delay
"d13=25m"
"d22=p2"
"d23=p3"
"d24=p17*1.26 -p1"
"d25=p1*2"
"d26=p7-p1"




1	ze
	10u ru2
2	d11 do:N
	d12
3	d12*3
4	d12*3
5	10u do:C1
	10u pl5:C1 
	10u pl7:N
#ifdef CARBON_LABEL
	10u fq2:C1	;jump to 56ppm
#endif
	d1
	1m LOCK_OFF
	10u pl1:H 
	10u pl5:C1
	(p7 ph0):N         ;eliminate Boltzmann
	p20:gp3
	3m
	(p1 ph0)
	2u
	p12:gp5
	d4
	(p7*2 ph6):N (d26 p1*2 ph0)
	2u
	p12:gp5
	d4
;------------- start INEPT to 15N -----------------------
	(p1 ph6)
	2u
	4u pl2:H
	p20:gp0
	1m pl1:H
	(p17 ph3):N
	if "l2==1" goto 20
	1m  ;p23:gp3
	1.7m
	(p7*2 ph5):N (d26 p1*2 ph0):H
	1m  ;p23:gp3
	1.7m
	d24
	(p1 ph4):H
20	d0 pl5:C1
#ifdef CARBON_LABEL
	2u
	(p5*2 ph10):C1	;(p1*2 ph0):H
	2u
	4u fq2:C1 	;jump to 177ppm
	(p5*2 ph10):C1
#else
 	d25  ;	(p1*2 ph0)   ;no 180 for F1-coupled spectrum
#endif
	d0
	if "l2==2"goto 40
	(p7 ph7):N
	goto 41
40	(p7 ph17):N
41	2u
	p11:gp1
	2m pl1:H
;------------- INEPT back to H, WATERGATE ---------------
	(p1 ph0)
	5u
	p12:gp0
	d5 pl2:H
	2u
	(p2 ph14)
	2u
	5u pl1:H
	(p1*2 ph15)
	2u
	5u pl2:H
	(p7*2 ph10):N (p2 ph14)
	2u
	p12:gp0
	d5 ph0
        5u pl31:N
	go=2 ph31 cpd2:N
	d11 wr #0 if #0 zd
	1m LOCK_ON
	10u do:N
	d12 iu2
     lo to 3 times 2
	d12 ru2
	d12 dp7
	d12 dp17
     lo to 4 times 2
	d12 id0
	d12 ip31
	d12 ip31
	lo to 5 times l3
	d12 do:C1
	d12 do:N
	1m LOCK_ON
	100u RESET
	exit

ph0=0
ph1=0
ph3=0 0 2 2
ph4=0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 2
ph5=0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3
ph6=1 3
ph7=0
ph17=1 1 1 1 3 3 3 3
ph10=0
ph14=2
ph15=0
ph31=0 2 2 0