# 1 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
# 1 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy" 1
;Delay.incl  -  include file for commonly used delays
;
;version 00/02/07
;
;$CLASS=HighRes Incl
;$COMMENT=


# 15 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
define pulse PG1
define pulse PG2
define pulse PG3
define pulse PG4
define pulse PG5
define pulse PG6
define pulse PG7
define pulse PG8
define pulse PG9
define pulse PG10
define pulse PG11
define pulse PG12
define pulse PG13
define pulse PG14
define pulse PG15
define pulse PG16
define pulse PG17
define pulse PG18
define pulse PG20
define pulse PG21
define pulse PG22
define pulse PG23
define pulse PG24
define pulse PG25
define pulse PG26
define pulse PG27
define pulse PG28
define pulse PG29
define pulse PG30
define pulse PG31


;general delays
# 49 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
define delay DELTA
define delay DELTA1
define delay DELTA2
define delay DELTA3
define delay DELTA4
define delay DELTA5
define delay DELTA6
define delay DELTA7
define delay DELTA8

define delay TAU
define delay TAU1
define delay TAU2
define delay TAU3
define delay TAU4
define delay TAU5

# 67 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;delays for centering pulses

define delay CEN_HN1
define delay CEN_HN2
define delay CEN_HN3
define delay CEN_HC1
define delay CEN_HC2
define delay CEN_HC3
define delay CEN_HC4
define delay CEN_HP1
define delay CEN_HP2
define delay CEN_CN1
define delay CEN_CN2
define delay CEN_CN3
define delay CEN_CN4
define delay CEN_CP1
define delay CEN_CP2

# 86 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;loop counters

define loopcounter COUNTER
define loopcounter SCALEF
define loopcounter FACTOR1
define loopcounter FACTOR2
define loopcounter FACTOR3

# 96 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;$Id: Delay.incl,v 1.13 2009/07/02 16:40:43 ber Exp $
;Avance3.incl
;   for AV III
;
;avance-version (07/12/14)
;
;$CLASS=HighRes Incl
;$COMMENT=


# 261 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;$Id: Avance3.incl,v 1.9 2009/07/02 16:40:43 ber Exp $
;Grad2.incl  -  include file for Gradient Spectroscopy
;   for TCU3
;
;avance-version (07/01/17)
;
;$CLASS=HighRes Incl
;$COMMENT=


# 288 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
define list<gradient> EA=<EA>


# 292 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;$Id: Grad2.incl,v 1.13 2009/07/02 16:40:43 ber Exp $
;De2.incl
;   for AVi
;
;avance-version (08/02/15)
;
;$CLASS=HighRes Incl
;$COMMENT=


define delay rde1
# 305 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
"rde1=de-de1"


# 316 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;for multi receive


# 325 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;$Id: De2.incl,v 1.7 2009/07/02 16:40:43 ber Exp $
;Daz.incl
;
;Dante-Z
;version  00/02/07
;
;$CLASS=HighRes Incl
;$COMMENT=


define delay DZ_TAU
define delay DZ_TAU2
define delay DZ_TAU3
define delay DZ_DELT			; 1/Bandwidth (cnst8)
# 340 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
define pulse DZ_THET
define pulse DZ_THE2
define pulse DZ_THE3
define pulse DZ_180
define pulse DZ_180a
define pulse DZ_180b
define pulse DZ_180c

define loopcounter DZ_LC1
define loopcounter DZ_LC2
define loopcounter DZ_LC3
define loopcounter DZ_LC4
define loopcounter DZ_LC5

define pulse P_TOTAL
# 357 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
"DZ_THET=10u"
"DZ_THE2=10u*0.5"
"DZ_THE3=10u*0.25"

"DZ_DELT=1s/cnst8"

"DZ_LC1=DZ_DELT/(dw)"
"DZ_LC2=(DZ_DELT/(dw))/3"
"DZ_LC3=(DZ_DELT/(dw))*2/3"
"DZ_LC4=(DZ_DELT/(dw))*2/9"
"DZ_LC5=(DZ_DELT/(dw))*6/9"
# 369 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
"DZ_TAU=dw*0.5"
"DZ_TAU2=dw*0.5*(DZ_LC1/(DZ_LC2*2+DZ_LC3))"
"DZ_TAU3=dw*0.5*(DZ_LC1/(DZ_LC4*4+DZ_LC5))"

"DZ_180=DZ_THET*DZ_LC1*2"
"DZ_180a=DZ_THET*DZ_LC1*2"
"DZ_180b=DZ_THET*DZ_LC2*2+DZ_THET*DZ_LC3*2" 	
				;=DZ_THE2*DZ_LC2*4+DZ_THET*DZ_LC3*2
"DZ_180c=DZ_THET*DZ_LC4+DZ_THET*DZ_LC4*2+DZ_THET*DZ_LC5*2"
				;=DZ_THE3*DZ_LC4*4+DZ_THE2*DZ_LC4*4+DZ_THET*DZ_LC5*2

# 390 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;#define DAZ363(phase1,phase2)                  \n;        1011 (DZ_TAU2) (DZ_THE2 phase1):f1     \n;             (DZ_TAU2) (DZ_THE2 phase2):f1     \n;             lo to 1011 times DZ_LC2           \n\ 


# 394 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;        1012 (DZ_TAU2) (DZ_THET phase1):f1     \n;             (DZ_TAU2) (DZ_THET phase2):f1     \n;             lo to 1012 times DZ_LC3           \n\ 


;        1013 (DZ_TAU2) (DZ_THE2 phase1):f1     \n;             (DZ_TAU2) (DZ_THE2 phase2):f1     \n;             lo to 1013 times DZ_LC2            
# 402 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;#define DAZ22622(phase1,phase2)                \n;        1021 (DZ_TAU3) (DZ_THE3 phase1):f1     \n;             (DZ_TAU3) (DZ_THE3 phase2):f1     \n;             lo to 1021 times DZ_LC4           \n\ 


# 406 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;        1022 (DZ_TAU3) (DZ_THE2 phase1):f1     \n;             (DZ_TAU3) (DZ_THE2 phase2):f1     \n;             lo to 1022 times DZ_LC4           \n\ 


;        1023 (DZ_TAU3) (DZ_THET phase1):f1     \n;             (DZ_TAU3) (DZ_THET phase2):f1     \n;             lo to 1023 times DZ_LC5           \n\ 
# 412 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;        1024 (DZ_TAU3) (DZ_THE2 phase1):f1     \n;             (DZ_TAU3) (DZ_THE2 phase2):f1     \n;             lo to 1024 times DZ_LC4           \n\ 


;        1025 (DZ_TAU3) (DZ_THE3 phase1):f1     \n;             (DZ_TAU3) (DZ_THE3 phase2):f1     \n;             lo to 1025 times DZ_LC4            
# 418 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;


# 422 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/bits.jfy"
;$Id: Daz.incl,v 1.13 2009/07/02 16:40:43 ber Exp $
# 1 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy" 2


;written by abx 7/28/98
;taken from trosy-ge.dw by Jordan 10/11/04
;HN(CA)CB trosy for full transfer to Cb
# 7 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
;p1 = 90 deg (10us) 1H pulse @pl1
;p7 = 90 deg (50us) 15N pulse @pl7
;p4 = selective 180 deg (23.7*2us) 13CA 
;     pulse @pl4 (f4) and @pl67 (f5)

;Gradient Pulses
;p20    2m 
;p21    300u
;
;gp0    (50,50,0) sine.50
;gp1    (50,0,0)  sine.50
;gp2    (50,0,0)  sine.50
;gp3    (0,50,0)  sine.50

# 22 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
"in0=inf3*0.5"
"in10=inf2*0.5"
"in20=inf1*0.5"


"d25=p22"
"d26=p22+100u"
# 30 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
;"d10=in10*0.5-p1*1.27"
;"d20=in20*0.5-p1*1.27"

;"d10=in10*0.5"
"d10=3u"
"d16=d10*2+6u+p7*4.4"

;"d20=3u"
;"d15=d20*2+18u+p7*4.4"
"d20=in20*0.5-10u"
# 41 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
"d0=105u+p7*0.635-p4*2-2u"


# 45 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
"d11=50m"
"d21=2.1m-p21"
"d22=2.5m-p21-p2"
"d23=2.1m-p2-14u"

;"d9=p7*2.2-d10*0.5-p1"
;"d19=p7*2.2-d20*0.5-p1"

;"d7=d8-p7-p27-1m-d9"
;"d17=d18-p7-p27-1m-d19-d9"
"d7=d8-p7-p27-1m"
"d17=d18-p7-p27-26m"
# 1 "mc_line 59 file /opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy exp. def. part of mc cmd. before ze"
; dimension 4 aq-mode (F3) Echo-Antiecho (F2) States-TPPI (F1) States-TPPI  F2->F1
define delay MCWRK
define delay MCREST
define delay d20orig
"d20orig=d20"
define delay d10orig
"d10orig=d10"
define delay d0orig
"d0orig=d0"
define list<loopcounter> SPARSELIST = <$NUSLIST>
define loopcounter t3loop
"t3loop=0"
define loopcounter t2loop
"t2loop=0"
define loopcounter t1loop
"t1loop=0"
define loopcounter ph1loop
"ph1loop=0"
define loopcounter ph2loop
"ph2loop=0"
define loopcounter ph3loop
"ph3loop=0"
define loopcounter STNCNT
"STNCNT=SPARSELIST.len/3"
"MCWRK = 0.016949*10m"
"MCREST = 10m - 10m"

    dccorr
# 59 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
# 59 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
1       ze
# 1 "mc_line 59 file /opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy exp. def. of mc cmd. after ze"
  MCWRK
  "t3loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph3loop = 0"
  MCWRK
  "t2loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph2loop = 2 * (t2loop % 2)"
  MCWRK
  "t1loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph1loop = 2 * (t1loop % 2)"
      "d20=d20orig + t1loop * in20 "
      "phval5 = (ph1loop % 4) * -90"
      MCWRK ip5 + phval5
      "d10=d10orig + t2loop * in10 "
      "phval6 = (ph2loop % 4) * 90"
      MCWRK ip6 + phval6
      "d0=d0orig + t3loop * in0 "
      "phval2 = (ph3loop % 2) * 180"
      MCWRK ip2 + phval2
      "phval12 = (ph3loop % 2) * 180"
      MCWRK ip12 + phval12
      "phval4 = (ph3loop % 2) * 180"
      MCWRK ip4 + phval4
      "phval7 = (t3loop * 180)%360"
      MCWRK ip7 + phval7
      "phval31 = (t3loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph3loop * 1"
      MCWRK setgrad EA
# 60 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
        1m
# 1 "mc_line 61 file /opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy exp. start label for mc cmd."
2 MCWRK  * 11 do:f2  do:f3 
LBLF3, MCWRK  * 13
LBLF2, MCWRK  * 13
LBLF1, MCWRK  * 21
LBLSTN, MCWRK setnmr3^0
  MCREST
# 62 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
        d11
        0.1m
3       3m
4       3m
5       3m
6       1m do:f3
7       5u do:f2
        10u pl4:f2
        d1 setnmr3^0 setnmr0^34^32^33 ctrlgrad 7
        1m setnmr3|0 setnmr0|34|32|33 ctrlgrad 0
        10u pl0:f1 
        10u pl7:f3
        10u fq=cnst8(bf ppm):f1
;------- start 90-degree on hn -----------------------
;       (center (d19 p1 ph0 d20 p1 ph5:r d19):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)
;       (center (p15:sp5 ph0 d20 p15:sp6 ph5:r):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)
;       (center (p15:sp5 ph0 5u d20 gron15 5u groff 5u d20 gron15*-1 5u groff p15:sp6 ph5:r):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)
       (center (p15:sp5 ph0 5u d20 5u 5u d20 5u p15:sp6 ph5:r):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)

;        10u pl0:f1
;       (p15:sp5 ph0):f1
;        3u
;        d20 gron15
;        3u groff 
;       (p7 ph0 p7*2.4 ph1 p7 ph0):f3
;        3u pl1:f1
;        d20 gron15*-1
;        3u groff
;       (p1 ph0 p1*2.3 ph1 p1 ph0):f1
;        d15
;       (p1 ph0 p1*2.3 ph1 p1 ph0):f1
;        3u
;        3u pl0:f1
;       (p15:sp6 ph5):f1
        20u
        12.5m ;gron17
        40u ;groff
        12.5m ;gron17*-1
        40u ;groff
        d17
       (p7 ph0):f3
        100u 
        p27:gp8
        800u
;       (center (d9 p1 ph0 d10 p1 ph6:r d9):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)
;       (center (p15:sp5 ph6:r d10 d10 p15:sp6 ph0):f1 (p7 ph0 p7*2.4 ph1 p7 ph0):f3)
       (p15:sp5 ph6:r):f1
        d10
       (p7 ph0 p7*2.4 ph1 p7 ph0):f3
        d10 pl1:f1
       (p1 ph0 p1*2.3 ph1 p1 ph0):f1
        d16
       (p1 ph0 p1*2.3 ph1 p1 ph0):f1
        3u
        3u pl0:f1
       (p15:sp6 ph0):f1

;--------NOE mixing starts--------------
        d7
       (p7 ph0):f3
        100u fq=0:f1
        p27:gp7  ;1.7m @44%
        900u pl1:f1
;------- start 90-degree on hn -----------------------
       (p1 ph0):f1
        3u
        p21:gp1
        d21
        300u
       (center (p1*2 ph0):f1 (p7*2 ph0):f3)
        d21
        p21:gp1
        303u
       (p1 ph1):f1                ;INEPT to 15N
        3u
        3u pl12:f1
       (p2 ph11:r):f1
;goto 999
        6u
        4u pl1:f1
        p20:gp0
        200u pl6:f2
        10u pl4:f2
;*****************************************
       (p7 ph7):f3
        10u
        p23:gp5*EA
        200u
       (p7*2 ph7):f3 
        p23:gp6*EA
        d0 
       (p4*2 ph0):f2
        4u
       (p4*2:sp8 ph0):f2
        d0
       (p1 ph12):f1
        3u
        3u pl2:f1
       (p2:sp0 ph2:r):f1
;goto 999
        6u gron2
        d23 pl1:f1
        2u groff
       (center (p1*2 ph0):f1 (p7*2 ph0):f3)
        3u gron2
        d23 pl2:f1
        3u groff
       (p2:sp1 ph0:r):f1
        6u
        2u pl1:f1
       (p1 ph0):f1 (p7 ph4):f3    ;DOUBLE 90
        5u
;goto 999
        de
        p21:gp3           
        d22 pl2:f1
       (center (d26 p2 ph10:r 3u 2u pl1 p1*2 ph0 3u 2u pl2 p2 ph10:r):f1 (p7*2 ph0):f3)
;goto 999
        5u pl1:f1
        p21:gp3
        d22
       (p7 ph0):f3
        p22:gp4 
999     95u pl30:f2
        5u setnmr0^34^32^33 ctrlgrad 7
       go=2 ph31 cpd2:f2
# 1 "mc_line 188 file /opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy exp. mc cmd. in line"
  MCWRK  do:f2 do:f3 setnmr3^0 wr #0 if #0 zd 

  MCWRK
  "ph3loop+=1"
      "d20=d20orig + t1loop * in20 "
      "phval5 = (ph1loop % 4) * -90"
      MCWRK ip5 + phval5
      "d10=d10orig + t2loop * in10 "
      "phval6 = (ph2loop % 4) * 90"
      MCWRK ip6 + phval6
      "d0=d0orig + t3loop * in0 "
      "phval2 = (ph3loop % 2) * 180"
      MCWRK ip2 + phval2
      "phval12 = (ph3loop % 2) * 180"
      MCWRK ip12 + phval12
      "phval4 = (ph3loop % 2) * 180"
      MCWRK ip4 + phval4
      "phval7 = (t3loop * 180)%360"
      MCWRK ip7 + phval7
      "phval31 = (t3loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph3loop * 1"
      MCWRK setgrad EA
  lo to LBLF3 times 2
  MCWRK 
  MCWRK
  "ph3loop-=2"
  MCWRK 

  MCWRK
  "ph2loop+=1"
      "d20=d20orig + t1loop * in20 "
      "phval5 = (ph1loop % 4) * -90"
      MCWRK ip5 + phval5
      "d10=d10orig + t2loop * in10 "
      "phval6 = (ph2loop % 4) * 90"
      MCWRK ip6 + phval6
      "d0=d0orig + t3loop * in0 "
      "phval2 = (ph3loop % 2) * 180"
      MCWRK ip2 + phval2
      "phval12 = (ph3loop % 2) * 180"
      MCWRK ip12 + phval12
      "phval4 = (ph3loop % 2) * 180"
      MCWRK ip4 + phval4
      "phval7 = (t3loop * 180)%360"
      MCWRK ip7 + phval7
      "phval31 = (t3loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph3loop * 1"
      MCWRK setgrad EA
  lo to LBLF2 times 2
  MCWRK 
  MCWRK
  "ph2loop-=2"
  MCWRK 

  MCWRK
  "ph1loop+=1"
      "d20=d20orig + t1loop * in20 "
      "phval5 = (ph1loop % 4) * -90"
      MCWRK ip5 + phval5
      "d10=d10orig + t2loop * in10 "
      "phval6 = (ph2loop % 4) * 90"
      MCWRK ip6 + phval6
      "d0=d0orig + t3loop * in0 "
      "phval2 = (ph3loop % 2) * 180"
      MCWRK ip2 + phval2
      "phval12 = (ph3loop % 2) * 180"
      MCWRK ip12 + phval12
      "phval4 = (ph3loop % 2) * 180"
      MCWRK ip4 + phval4
      "phval7 = (t3loop * 180)%360"
      MCWRK ip7 + phval7
      "phval31 = (t3loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph3loop * 1"
      MCWRK setgrad EA
  lo to LBLF1 times 2
  MCWRK 
  MCWRK
  "ph1loop-=2"
  MCWRK 

  MCWRK
  "t3loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph3loop = 0"


  MCWRK
  "t2loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph2loop = 2 * (t2loop % 2)"


  MCWRK
  "t1loop = SPARSELIST"
  MCWRK  SPARSELIST.inc
  MCWRK
  "ph1loop = 2 * (t1loop % 2)"
      "d20=d20orig + t1loop * in20 "
      "phval5 = (ph1loop % 4) * -90"
      MCWRK ip5 + phval5
      "d10=d10orig + t2loop * in10 "
      "phval6 = (ph2loop % 4) * 90"
      MCWRK ip6 + phval6
      "d0=d0orig + t3loop * in0 "
      "phval2 = (ph3loop % 2) * 180"
      MCWRK ip2 + phval2
      "phval12 = (ph3loop % 2) * 180"
      MCWRK ip12 + phval12
      "phval4 = (ph3loop % 2) * 180"
      MCWRK ip4 + phval4
      "phval7 = (t3loop * 180)%360"
      MCWRK ip7 + phval7
      "phval31 = (t3loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph3loop * 1"
      MCWRK setgrad EA
  lo to LBLSTN times STNCNT
  MCWRK
# 192 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
# 193 "/opt/topspin3.1/exp/stan/nmr/lists/pp/user/doubleNOESY_trosy.jfy"
1m do:f2
1m do:f3
exit    
        
ph0=0
ph1=1          ;original setting1     ;check right phase for Boltzmann !!!!!
ph2=3          ;original setting
ph4=1
ph5=0 0 2 2
ph6=0
ph7=1 3
ph10=2
ph11=0 
ph12=1         ;original seting
ph31=1 3 3 1


