# 1 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
# 1 "<built-in>"
# 1 "<command-line>"
# 1 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
;hsqcnoef3gpsi
;avance-version (12/01/11)
;2D H-1/X correlation via double inept transfer
;   using sensitivity improvement
;for measuring H1-N15 NOEs
;phase sensitive using Echo/Antiecho-TPPI gradient selection
;with decoupling during acquisition
;using f3 - channel
;recording NONOE and NOE interleaved
;(use parameterset HSQCNOEF3GPSI)
;
;$CLASS=HighRes
;$DIM=2D
;$TYPE=
;$SUBTYPE=
;$COMMENT=


prosol relations=<triple>
# 1 "/opt/topspin3.6.4/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.6.4/exp/stan/nmr/lists/pp/Avance.incl"
;$Id: Avance3.incl,v 1.10.2.2 2013/08/30 09:43:33 ber Exp $
# 23 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi" 2

# 1 "/opt/topspin3.6.4/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.6.4/exp/stan/nmr/lists/pp/Grad.incl"
define list<gradient> EA=<EA>


# 31 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/Grad.incl"
;$Id: Grad2.incl,v 1.14 2012/01/31 17:49:21 ber Exp $
# 24 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi" 2

# 1 "/opt/topspin3.6.4/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.6.4/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.6.4/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.6.4/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.6.4/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.6.4/exp/stan/nmr/lists/pp/Delay.incl"
;$Id: Delay.incl,v 1.14.2.1 2013/08/30 09:43:33 ber Exp $
# 25 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi" 2


# 28 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
"p2=p1*2"
"p22=p21*2"
"d11=30m"
"d12=20u"
"d24=1s/(cnst4*cnst11)"
"d25=1s/(cnst4*cnst12)"
"d26=1s/(cnst4*4)"

"p0=p1*4/3"

# 39 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
"d0=3u"

"in0=inf1/2"

		
"DELTA2=p16+d16+8u"
"DELTA3=d25-p16-d16-4u"
# 50 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
"DELTA1=d25-p16-d16-p2-d0*2"


# 54 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
"l0=1"
"l4=d1/(p0+5m)"


"spoff1=0"
"spoff13=bf2*((cnst21+cnst22)/2000000)-o2"
# 61 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
		
# 1 "mc_line 62 file /opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi exp. def. part of mc cmd. before ze"
; dimension 2D; AQ_mode  (F1) Echo-Antiecho
define delay MCWRK
define delay MCREST
define delay d0orig
"d0orig=d0"
define loopcounter t1loop
"t1loop=0"
define loopcounter ph1loop
"ph1loop=0"
define loopcounter ST1CNT
"ST1CNT = trunc(td1 / ( 2 * 2 ))"
"MCREST = d11 - d11"
"MCWRK = 0.062500*d11"

    dccorr
# 62 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
1 ze
# 1 "mc_line 62 file /opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi exp. def. of mc cmd. after ze"
LBLAV, MCWRK
      "d0=d0orig + t1loop * in0 "
      "phval5 = (ph1loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval6 = (t1loop * 180)%360"
      MCWRK ip6 + phval6
      "phval8 = (t1loop * 180)%360"
      MCWRK ip8 + phval8
      "phval31 = (t1loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph1loop * 1"
      MCWRK setgrad EA
# 63 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
  d11 pl16:f3
# 1 "mc_line 64 file /opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi exp. start label for mc cmd."
2 MCWRK  do:f3
LBLF1I1, MCWRK  * 7
LBLF1, MCWRK  * 7
LBLST1, MCWRK 
  MCREST
# 65 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
3 d12

  if "l0 %2 == 1"
     {
     d1
     50u setnmr3|0 setnmr0|34|32|33 ctrlgrad 0
     4u pl0:f1
     (p11:sp1 ph1:r):f1
     4u
     4u pl1:f1
     }
  else
     {
4    (p0 ph1)
     5m
     lo to 4 times l4
     4u
     50u setnmr3|0 setnmr0|34|32|33 ctrlgrad 0
     }

  (p1 ph1)
  p16:gp1
  d16 pl3:f3
# 89 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
  (p21 ph8):f3
  4u
  p16:gp2*-1*EA
  d16
  DELTA3
  (p22 ph6):f3
  d0 


# 100 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
  (p2 ph9)


  d0
  p16:gp2*EA
  d16
  DELTA1
# 108 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
  (center (p1 ph1) (p21 ph4):f3 )
  d24
  (center (p2 ph1) (p22 ph1):f3 )
  d24
  (center (p1 ph2) (p21 ph5):f3 )
  d26
  (center (p2 ph1) (p22 ph1):f3 )
  d26
  (p1 ph1)
  DELTA2
  (p2 ph1)
  4u
  p16:gp3
  d16 pl16:f3
  4u  setnmr3^0 setnmr0^34^32^33 ctrlgrad 7
  go=2 ph31 cpd3:f3
# 1 "mc_line 124 file /opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi exp. mc cmd. in line"
  MCWRK  do:f3 wr #0 if #0 zd iu0
  lo to LBLF1I1 times 2
  MCWRK 

  "ph1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval5 = (ph1loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval6 = (t1loop * 180)%360"
      MCWRK ip6 + phval6
      "phval8 = (t1loop * 180)%360"
      MCWRK ip8 + phval8
      "phval31 = (t1loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph1loop * 1"
      MCWRK setgrad EA
  lo to LBLF1 times 2
  MCWRK 

  "t1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval5 = (ph1loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval6 = (t1loop * 180)%360"
      MCWRK ip6 + phval6
      "phval8 = (t1loop * 180)%360"
      MCWRK ip8 + phval8
      "phval31 = (t1loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph1loop * 1"
      MCWRK setgrad EA
  lo to LBLST1 times ST1CNT
  MCWRK 
  "t1loop=0"
  "ph1loop=0"
  MCWRK rf #0
  lo to LBLAV times tdav
# 127 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
exit
   

ph1=0 
ph2=1
ph4=0 0 2 2
ph5=3 3 1 1
ph6=0
ph7=3
ph8=1 3 
ph9=0 0 2 2
ph31=0 2 2 0 
  

;pl0 : 0W
;pl1 : f1 channel - power level for pulse (default)
;pl3 : f3 channel - power level for pulse (default)
;pl16: f3 channel - power level for CPD/BB decoupling
;sp1 : f1 channel - shaped pulse  90 degree  (H2O on resonance)
;sp13: f2 channel - shaped pulse 180 degree  (Ca and C=O, adiabatic)
;p0 : f1 channel - 120 degree high power pulse
;p1 : f1 channel -  90 degree high power pulse
;p2 : f1 channel - 180 degree high power pulse
;p8 : f2 channel - 180 degree shaped pulse for inversion (adiabatic)
;p11: f1 channel -  90 degree shaped pulse           [1 msec]
;p16: homospoil/gradient pulse
;p21: f3 channel -  90 degree high power pulse
;p22: f3 channel - 180 degree high power pulse
;d0 : incremented delay (2D)                         [3 usec]
;d1 : relaxation delay; 1-5 * T1
;d11: delay for disk I/O                             [30 msec]
;d12: delay for power switching                      [20 usec]
;d16: delay for homospoil/gradient recovery
;d24: 1/(4J)YH for YH
;     1/(8J)YH for all multiplicities
;d25: 1/(4J)YH for YH
;     1/(8J)YH for all multiplicities
;d26: 1/(4J(YH))
;cnst4: = J(YH)
;cnst11: for multiplicity selection = 4 for NH, 8 for all multiplicities
;cnst12: for multiplicity selection = 4 for NH, 8 for all multiplicities
;cnst21: CO chemical shift (offset, in ppm)
;cnst22: Calpha chemical shift (offset, in ppm)
;l0: flag to switch between NOE and NONOE
;l4: loop for saturation sequence: = d1 / (p0 + 5m)
;inf1: 1/SW(X) = 2 * DW(X)
;in0: 1/(2 * SW(X)) = DW(X)
;nd0: 2
;ns: 1 * n
;ds: >= 16
;td1: total number of experiments
;        = number of experiments for each 2D * 2
;FnMODE: echo-antiecho
;cpd3: decoupling according to sequence defined by cpdprg3
;pcpd3: f3 channel - 90 degree pulse for decoupling sequence
# 184 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
;use gradient ratio:	gp 1 : gp 2 : gp 3
;			  50 :   80 : 16.2      for N-15

;for z-only gradients:
;gpz1: 50%
;gpz2: 80%
;gpz3: 16.2%

;use gradient files:
;gpnam1: SMSQ10.100
;gpnam2: SMSQ10.100
;gpnam3: SMSQ10.100
# 198 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
                                          ;preprocessor-flags-start
;LABEL_CN: for C-13 and N-15 labeled samples start experiment with 
;             option -DLABEL_CN (eda: ZGOPTNS)
                                          ;preprocessor-flags-end


;use AU-program split [2] to separate NOE and NONOE data into
;   different datasets
# 209 "/opt/topspin3.6.4/exp/stan/nmr/lists/pp/hsqcnoef3gpsi"
;$Id: hsqcnoef3gpsi,v 1.11 2012/01/31 17:49:26 ber Exp $
