# 1 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
;hcchdigp3d
;avance-version (14/02/11)
;HCCH-TOCSY
;3D sequence with
;   inverse correlation using multiple inept transfer and
;      C-C DIPSI3 spinlock
;
;      F1(H,t1) -> F2(C,t2) -> F2(C') -> F1(H',t3)
;
;off resonance C=O pulse using shaped pulse
;phase sensitive (t1)
;phase sensitive (t2)
;spinlock during z-filter
;(use parameterset HCCHDIGP3D)
;
;(L.E. Kay, G.Y. Xu, A.U. Singer, D.R. Muhandiram & J. D. Forman-Kay
;   J. Magn. Reson. B 101, 333 - 337 (1993))
;
;$CLASS=HighRes
;$DIM=3D
;$TYPE=
;$SUBTYPE=
;$COMMENT=


prosol relations=<triple>
# 1 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Avance.incl" 1
;Avance3.incl
;   for AV III
;
;avance-version (13/06/17)
;
;$CLASS=HighRes Incl
;$COMMENT=


# 169 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Avance.incl"
;$Id: Avance3.incl,v 1.9.8.3 2013/08/30 09:44:49 ber Exp $
# 29 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d" 2

# 1 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Grad.incl" 1
;Grad2.incl  -  include file for Gradient Spectroscopy
;   for TCU3
;
;avance-version (07/01/17)
;
;$CLASS=HighRes Incl
;$COMMENT=

# 27 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Grad.incl"
define list<gradient> EA=<EA>


# 31 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Grad.incl"
;$Id: Grad2.incl,v 1.13.8.1 2012/01/31 17:56:17 ber Exp $
# 30 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d" 2

# 1 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl" 1
;Delay.incl  -  include file for commonly used delays
;
;version (13/08/07)
;
;$CLASS=HighRes Incl
;$COMMENT=

# 9 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl"
;general delays

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 DELTA9
define delay DELTA10
define delay DELTA11
define delay DELTA12
define delay DELTA13
define delay DELTA14
define delay DELTA15
define delay DELTA16

define delay TAU
define delay TAU1
define delay TAU2
define delay TAU3
define delay TAU4
define delay TAU5
define delay TAU6
define delay TAU7
define delay TAU8
define delay TAU9
# 40 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl"
define delay INCR1
define delay INCR2
define delay INCR3
define delay INCR4
define delay INCR5
define delay INCR6


;delays for centering pulses
# 50 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl"
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


;loop counters
# 69 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl"
define loopcounter COUNTER
define loopcounter SCALEF
define loopcounter FACTOR1
define loopcounter FACTOR2
define loopcounter FACTOR3
define loopcounter FACTOR4
define loopcounter FACTOR5
define loopcounter FACTOR6


# 80 "/opt/topspin3.2/exp/stan/nmr/lists/pp/Delay.incl"
;$Id: Delay.incl,v 1.13.8.2 2013/08/30 09:44:49 ber Exp $
# 31 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d" 2


# 34 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
"p2=p1*2"
"p4=p3*2"
"p22=p21*2"
"d11=30m"
"d12=20u"

"d4=1.6m"			;tau a
"d21=1.1m"			;tau c
"d23=475u"			;tau b

# 45 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
"p16=500u"
"p19=2m"
"p29=300u"
"p30=5m"
"p31=4.4m"

"d16=150u"

# 54 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
"d0=3u"
"d10=3u"

"in0=inf1/2"
"in10=inf2/2"

# 61 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
"DELTA1=d4-p16-d16-4u"
"DELTA2=d4-p16-d16-4u+d0*2+p4"
"DELTA3=d23-p29-d16"
"DELTA4=p22+p2+d10*2+4u"
"DELTA5=d21-p16-d16-4u"
"DELTA6=d4-p16-d16-p3*2-15u"


"d31=(p9*54.33*4)*l1"
# 72 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
"spoff5=bf2*(cnst21/1000000)-o2"


"acqt0=-p1*2/PI"
# 78 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
aqseq 312


# 1 "mc_line 81 file /opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d exp. def. part of mc cmd. before ze"
; dimension 3 aq-mode (F2) States-TPPI (F1) States-TPPI  F2->F1
define delay MCWRK
define delay MCREST
define delay d0orig
"d0orig=d0"
define delay d10orig
"d10orig=d10"
define loopcounter t2loop
"t2loop=0"
define loopcounter t1loop
"t1loop=0"
define loopcounter ph1loop
"ph1loop=0"
define loopcounter ph2loop
"ph2loop=0"
define loopcounter ST1CNT
"ST1CNT = td2 / 2"
define loopcounter ST2CNT
"ST2CNT = td1 / 2"
"MCREST = d11 - d11"
"MCWRK = 0.071429*d11"

    dccorr
# 81 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
1 d11 ze
# 1 "mc_line 81 file /opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d exp. def. of mc cmd. after ze"
      "d0=d0orig + t1loop * in0 "
      "phval3 = (ph1loop % 4) * 90"
      "phval3 = (phval3 + t1loop * -90)%360"
      MCWRK ip3 + phval3
      "d10=d10orig + t2loop * in10 "
      "phval4 = (ph2loop % 4) * 90"
      MCWRK ip4 + phval4
# 82 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
  d31 pl12:f2
# 1 "mc_line 83 file /opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d exp. start label for mc cmd."
2 MCWRK  * 3 do:f2
LBLF2, MCWRK  * 3
LBLST2, MCWRK  * 4
LBLF1, MCWRK  * 3
LBLST1, MCWRK 
  MCREST
# 84 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
3 d1
  50u setnmr3|0 setnmr0|34|32|33 ctrlgrad 0
  d12 pl1:f1
# 88 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
  (p1 ph3)
  4u
  p16:gp1
  d16
  DELTA1 pl2:f2
  d0
  (p4 ph1):f2
  d0
  (p2 ph1)
  4u
  p16:gp1
  d16
  DELTA2 pl3:f3
  (p1 ph2)

  p19:gp3
  d16

  (p3 ph4):f2
  d10
  (p22 ph1):f3
  4u
  p29:gp2
  d16
  DELTA3 pl0:f2
  (p14:sp5 ph1):f2
  4u
  p29:gp2
  d16
  DELTA3 pl2:f2
  p2 ph1
  d10
  (p4 ph1):f2
  DELTA4
  p29:gp2
  d16
  DELTA3
  (p14:sp5 ph1):f2
  4u
  p29:gp2
  d16
  DELTA3 pl2:f2
  (p3 ph2):f2
  4u
  d12 pl15:f2
						;begin DIPSI3
9 (p9*2.722 ph7):f2
  (p9*4.389 ph9):f2
  (p9*2.778 ph7):f2
  (p9*3.056 ph9):f2
  (p9*0.333 ph7):f2
  (p9*2.556 ph9):f2
  (p9*4.000 ph7):f2
  (p9*2.722 ph9):f2
  (p9*4.111 ph7):f2
  (p9*3.778 ph9):f2
  (p9*3.889 ph7):f2
  (p9*2.889 ph9):f2
  (p9*3.000 ph7):f2
  (p9*0.333 ph9):f2
  (p9*2.500 ph7):f2
  (p9*4.050 ph9):f2
  (p9*2.830 ph7):f2
  (p9*4.389 ph9):f2
  (p9*2.722 ph9):f2
  (p9*4.389 ph7):f2
  (p9*2.778 ph9):f2
  (p9*3.056 ph7):f2
  (p9*0.333 ph9):f2
  (p9*2.556 ph7):f2
  (p9*4.000 ph9):f2
  (p9*2.722 ph7):f2
  (p9*4.111 ph9):f2
  (p9*3.778 ph7):f2
  (p9*3.889 ph9):f2
  (p9*2.889 ph7):f2
  (p9*3.000 ph9):f2
  (p9*0.333 ph7):f2
  (p9*2.500 ph9):f2
  (p9*4.050 ph7):f2
  (p9*2.830 ph9):f2
  (p9*4.389 ph7):f2
  (p9*2.722 ph9):f2
  (p9*4.389 ph7):f2
  (p9*2.778 ph9):f2
  (p9*3.056 ph7):f2
  (p9*0.333 ph9):f2
  (p9*2.556 ph7):f2
  (p9*4.000 ph9):f2
  (p9*2.722 ph7):f2
  (p9*4.111 ph9):f2
  (p9*3.778 ph7):f2
  (p9*3.889 ph9):f2
  (p9*2.889 ph7):f2
  (p9*3.000 ph9):f2
  (p9*0.333 ph7):f2
  (p9*2.500 ph9):f2
  (p9*4.050 ph7):f2
  (p9*2.830 ph9):f2
  (p9*4.389 ph7):f2
  (p9*2.722 ph7):f2
  (p9*4.389 ph9):f2
  (p9*2.778 ph7):f2
  (p9*3.056 ph9):f2
  (p9*0.333 ph7):f2
  (p9*2.556 ph9):f2
  (p9*4.000 ph7):f2
  (p9*2.722 ph9):f2
  (p9*4.111 ph7):f2
  (p9*3.778 ph9):f2
  (p9*3.889 ph7):f2
  (p9*2.889 ph9):f2
  (p9*3.000 ph7):f2
  (p9*0.333 ph9):f2
  (p9*2.500 ph7):f2
  (p9*4.050 ph9):f2
  (p9*2.830 ph7):f2
  (p9*4.389 ph9):f2
  lo to 9 times l1
						;end DIPSI3
  d12 pl10:f1
  (p17 ph1)
  (p17*2 ph2)
  4u
  p30:gp4
  d16 pl1:f1
  (p1 ph1)
  4u
  p31:gp4
  d16 pl2:f2
# 219 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
  (p3 ph2):f2
  4u
  p16:gp1
  d16
  DELTA5
  (center (p2 ph1) (p4 ph1):f2 )
  4u
  p16:gp1
  d16
  DELTA5
  (p3 ph1):f2

  (p1 ph1)
  4u
  p16:gp1
  d16
  DELTA1
  (center (p2 ph1) (p4 ph1):f2 )
  4u
  p16:gp1
  d16
  DELTA6 
  4u setnmr3^0 setnmr0^34^32^33 ctrlgrad 7
  (p3 ph1 3u p3 ph5):f2
  4u pl12:f2 
  (p1 ph1)
  go=2 ph31 cpd2:f2
# 1 "mc_line 246 file /opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d exp. mc cmd. in line"
  MCWRK  do:f2 wr #0 if #0 zd 

  "ph1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval3 = (ph1loop % 4) * 90"
      "phval3 = (phval3 + t1loop * -90)%360"
      MCWRK ip3 + phval3
      "d10=d10orig + t2loop * in10 "
      "phval4 = (ph2loop % 4) * 90"
      MCWRK ip4 + phval4
  lo to LBLF2 times 2
  MCWRK 

  "t1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval3 = (ph1loop % 4) * 90"
      "phval3 = (phval3 + t1loop * -90)%360"
      MCWRK ip3 + phval3
      "d10=d10orig + t2loop * in10 "
      "phval4 = (ph2loop % 4) * 90"
      MCWRK ip4 + phval4
  lo to LBLST2 times ST2CNT
  MCWRK 
  "t1loop=0"
  "ph1loop=0"
  MCWRK 

  "ph2loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval3 = (ph1loop % 4) * 90"
      "phval3 = (phval3 + t1loop * -90)%360"
      MCWRK ip3 + phval3
      "d10=d10orig + t2loop * in10 "
      "phval4 = (ph2loop % 4) * 90"
      MCWRK ip4 + phval4
  lo to LBLF1 times 2
  MCWRK 

  "t2loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval3 = (ph1loop % 4) * 90"
      "phval3 = (phval3 + t1loop * -90)%360"
      MCWRK ip3 + phval3
      "d10=d10orig + t2loop * in10 "
      "phval4 = (ph2loop % 4) * 90"
      MCWRK ip4 + phval4
  lo to LBLST1 times ST1CNT
  MCWRK 
  "t2loop=0"
  "ph2loop=0"
  MCWRK
# 249 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
exit

# 252 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
ph1=0
ph2=1
ph3=0 2
ph4=0 0 2 2
ph5=0 0 0 0 2 2 2 2
ph7=1
ph9=3
ph31=0 2 2 0


;pl0 : 0W
;pl1 : f1 channel - power level for pulse (default)
;pl2 : f2 channel - power level for pulse (default)
;pl3 : f3 channel - power level for pulse (default)
;pl10: f1 channel - power level for TOCSY-spinlock (trim pulse)
;pl12: f2 channel - power level for CPD/BB decoupling
;pl15: f2 channel - power level for TOCSY-spinlock
;sp5: f2 channel - shaped pulse 180 degree   (C=O off resonance)
;p1 : f1 channel -  90 degree high power pulse
;p2 : f1 channel - 180 degree high power pulse
;p3 : f2 channel -  90 degree high power pulse
;p4 : f2 channel - 180 degree high power pulse
;p9 : f2 channel -  90 degree low power pulse
;p14: f2 channel - 180 degree shaped pulse
;p16: homospoil/gradient pulse                         [500 usec]
;p17: f1 channel - trim pulse                          [2.5 msec]
;p19: gradient pulse 2                                 [2 msec]
;p22: f3 channel - 180 degree high power pulse
;p29: gradient pulse 3                                 [300 usec]
;p30: gradient pulse 4                                 [5 msec]
;p31: gradient pulse 5                                 [4.4 msec]
;d0 : incremented delay (F1 in 3D)                     [3 usec]
;d1 : relaxation delay; 1-5 * T1
;d4 : 1/(4J(CH)) - tau a                               [1.6 msec]
;d10: incremented delay (F2 in 3D)                     [3 usec]
;d11: delay for disk I/O                               [30 msec]
;d12: delay for power switching                        [20 usec]
;d16: delay for homospoil/gradient recovery
;d21: 1/(6J'(CH)) - tau c                              [1.1 msec]
;d23: tau b                                            [475 usec]
;d31: length of DIPSI-3 cycle as executed = (p9*54.33*4)*l1
;cnst21: CO chemical shift (offset, in ppm)
;cnst23: Caliphatic chemical shift (offset, in ppm)
;o2p: Caliphatic chemical shift (cnst23)
;l1: loop for DIPSI cycle:
;       mixing time = ((p9*54.33*4) * l1)              [12 msec]
;inf1: 1/SW(Hali) = 2 * DW(Hali)
;inf2: 1/SW(C) = 2 * DW(C)
;in0: 1/(2 * SW(Hali)) = DW(Hali)
;nd0: 2
;in10: 1/(2 * SW(C)) = DW(C)
;nd10: 2
;ns: 8 * n
;ds: 32
;td1: number of experiments in F1
;td2: number of experiments in F2
;FnMODE: States-TPPI in F1 (not suitable for TPPI)
;FnMODE: States-TPPI (or TPPI) in F2
;cpd2: decoupling according to sequence defined by cpdprg2
;pcpd2: f2 channel - 90 degree pulse for decoupling sequence
# 314 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
;use gradient ratio:    gp 1 : gp 2 : gp 3 : gp 4
;                         16 :   16 :   30 :   60

;for z-only gradients:
;gpz1: 16%
;gpz2: 16%
;gpz3: 30%
;gpz4: 60%

;use gradient files:
;gpnam1: SMSQ10.50
;gpnam2: SMSQ10.50
;gpnam3: SMSQ10.100
;gpnam4: SMSQ10.100
# 331 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
;Processing

;SR(F1): 1/4 SWH(F1)

# 337 "/opt/topspin3.2/exp/stan/nmr/lists/pp/hcchdigp3d"
;$Id: hcchdigp3d,v 1.18.6.2 2014/02/11 11:32:18 ber Exp $
