; hmqc_ipap_noesy.jfy
; constant time in both indirect dimensions
; for measuring coulpings between 2'-OH and C1'/C2'
; Jinfa Ying, 5/19/05

#include "bits.nt"

#define CARBON
#define PROTON

;p1 = 90 deg (10us) 1H pulse @pl1
;p3 = 90 deg 13C pulse @pl3
;p7 = 90 deg 15N pulse @pl7

"d2=3.0m"
"d3=d2*0.5-10u-p22-p3"
"d5=25m"

"d9=3u"
"d6=d5-p7*2-d9"

"d21=p3*2.22-p1*2.11"
"d22=p3*2"

"d11=50m"
"d12=1m"

"d28=3u"
"d29=d2-p3*4.44-d28-4u"
"d0=d5*0.5-p1*2.11-2u+d28+p3*2.22+2u+p1*0.637"   ;d0+d30+p1*4.22+4u=25m
"d30=d5*0.5-p1*2.11-2u-d28-p3*2.22-2u-p1*0.637"  ;d2-d28*2-p3*4.44-4u+d0-p1*0.637=d30+d2+p1*0.637

"d10=d5*0.5-p3"
"d20=d5*0.5-p3"

"d7=d8-p24-2.5m"   ; d8 = noesy mixing time

"in10=in20"
"in6=in20"      ; dd6 needed to move 15N dec pulse
"in9=in20"
"in30=in0"
"in28=in0*0.25"
"in29=in28"

1    ze
     1m RESET
     1m
2    d12
     d11
     2m
20   d12
3    d12
4    d12
5    d12
6    d12 do:C1
7    1m LOCK_ON
     d1 do:N do:H
     50u UNBLKGRAD
     1m LOCK_OFF
     20u pl1:H
     20u pl3:C1
     20u pl7:N
     (p3 ph0):C1  ;kill Boltzman
     50m

;**************** MQ preparation, 13C & 1H evolution
    (p1 ph10):H
     d28
    (d21 p1 ph0 2u p1*2.22 ph1 2u p1 ph0 d21):H (p3 ph0 2u p3*2.44 ph1 2u p3 ph0):C1
     d29
    (p3 ph11):C1

;???*****1u has to be added in the following line in either C1 or H pulse train
;???*****otherwise, error message "duration out of range" occurs on our drx800
    (1u d10 p3*2 ph5 d20 1u):C1 (d0 p1 ph0 3u p1*2.22 ph1 3u p1 ph0 d30):H (1u d6 p7*2 ph0 d9 1u):N

    (p3 ph0):C1
     10u
     p22:gp22    ; 1m @ 3~5%
     d3
if "l1==2" goto 32
     d22               ;(on i.e. p3*2)/(off i.e. d22) to get in-phase/anti-phase
     10u
     p22:gp23    ; 1m @ -3~-5%
     d3
    (p1 ph9):H
goto 33
32  (p3*2 ph0):C1
     10u
     p22:gp23    ; 1m @ -3~-5%
     d3
    (p1 ph8):H

33   d7                  ;d8 = mixing time
     p24:gp24
     2.5m pl1:H
    (p1 ph2:r):H         ; jump-return
     d4                  ; delay to optimize jump-return for 2'-OH
    (p0 ph0):H

;***************Excitation Sculpting
     10u
     p25:gp25
     100u pl2:H
    (p2 ph14:r):H
     3u
     3u pl1:H
    (p1*2 ph4):H
     3u
     3u pl2:H
    (p2 ph14:r):H
     10u
     p25:gp25
     100u BLKGRAD
     go=2 ph31
     1m do:C1
     1m LOCK_ON
     d11 wr #0 if #0 zd

;*******interleave loop
     d12 iu1
  lo to 20 times 2
     d12*0.5 ru1

#ifdef CARBON
     d12*0.5 ip11
  lo to 3 times 2
     d12*0.2 ip31*2
     d12*0.2 dd10
     d12*0.2 id20
     d12*0.2 dd6
     d12*0.2 id9
  lo to 4 times l3          ; 13C loop
     d12*0.2 rd10
     d12*0.2 rd20
     d12*0.2 rd6
     d12*0.2 rd9
#endif

#ifdef PROTON
     d12*0.2 ip10
  lo to 5 times 2
     d12*0.25 id28
     d12*0.25 dd29
     d12*0.125 dd0
     d12*0.125 id30
     d12*0.25 ip31*2
  lo to 6 times l4         ; 1H loop
#endif


1m RESET
1m LOCK_ON

exit


ph0=0
ph1=1
ph2=2 
ph4=1
ph5=0 1 2 3
ph8=1
ph9=2
ph10=0
ph11=0
ph14=3
ph31=0 2