;real.sel2hmbc_gs.kaf>
; sel2hmbc_gs.kaf   selective hmbc with gradients  
; this pulse sequence is used when the selective 180 is  
; greater than approx 1.5 ms
; set o1 on proton peak of interest 
; measure distance in hz to homonuclear coupling partner
; selective 180= (1/(distance in hz/1.73205))/2
; calibrate power of HL1 to give calculated 180
; Bax, Farley, and Walker  JMR  Series A 119, 134-138 (1996)
; for more information about this pulse sequence, see:
; sel1hmbc_gs.kaf
;BMRB Pulse Sequence Accession Number: 7

p2=p1*2
d0=3u
d11=1m
d13=4u
d16=125u
d27=p16/l21-2u
d2=1s/(cnst2*2)   ; one bond J filter cnst2=145 for 1 bond Jch

1 ze
  1m dl1
2 d1 setf2|3^9
  100u hl1
3 p1 ph1
  d2 ; one bond J filter
  p3:d ph5
  d6 ; evolution of long range couplings
  p3:d ph3
  d0
  50u setf2^3|9 
4 2u:ngrad
  d27
  lo to 4 times l21
  d16:ngrad setf2|3
  d5
  p3:d ph2
  d13
  p2 ph0
  50u setf2^3
5 2u:ngrad
  d27
  lo to 5 times l21
  d16:ngrad setf2|3
  d0
  go=2 ph31 
  d1 do wr #0 if #0 id0 zd
  d11 setf2^9|3
  lo to 3 times td1
  d11 setf2^9|3
exit 
  
ph0=0
ph1=0 
ph2=0 0 2 2
ph3=0 2
ph5=0 0 0 0 2 2 2 2
ph31=0 2 2 0

; Remember: NDO=1 !!!!  (second DO is outside of evolution period)
; p1 90 1H pulse
; p3 90 13C pulse
; cnst2 should=145 hz
; l[21]= 32 or 64 (value from gradpulse program)
; p16= 1ms
; hl1=1dB
; DL[1]=16db for a 13uS 90 pulse (Hard X channel)
;d6 : delay for evolution of long range couplings 60-80 ms
;d16: delay for homospoil/gradient recovery 1ms
;cnst2: 1 bond filter should be approximately 145
;d27: delay for shaped gradient
;L21: loop for shaped gradients
;in0: 1/(2 * SW(X)) = DW(X)
;NS: 2 * n
;DS: 16
;td1: number of experiments
;MC2: QF 
;use gradient ratio 	8:10
;use sine for wave form