# 1 "/nmrdata/nmrlib/pp/B_HNCO"
;IBS_Best_HNCO
;BH-HNCO
;using broadband 15N pulses
;ZGOPTION -DINEPT
;BB 10/01/14
;
;$CLASS=BB-Assign
;$DIM=3D
;$TYPE=BEST-HSQC
;$SUBTYPE=
;$COMMENT=


prosol relations=<IBS>
# 1 "/opt/topspin3.5pl7/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.5pl7/exp/stan/nmr/lists/pp/Avance.incl"
;$Id: Avance3.incl,v 1.10.2.2 2013/08/30 09:43:33 ber Exp $
# 17 "/nmrdata/nmrlib/pp/B_HNCO" 2

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


# 31 "/opt/topspin3.5pl7/exp/stan/nmr/lists/pp/Grad.incl"
;$Id: Grad2.incl,v 1.14 2012/01/31 17:49:21 ber Exp $
# 18 "/nmrdata/nmrlib/pp/B_HNCO" 2

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


# 25 "/nmrdata/nmrlib/pp/B_HNCO"
"p41=7.2/(cnst2*bf1/1000000)" 
"spw25=plw1*(pow((p1*1.01/p41)/0.125,2))" 
"spoff25=bf1*(cnst1/1000000)-o1"  
;"spoal25=0.5"

"p42=4.875/(cnst2*bf1/1000000)" 
"spw26=plw1*(pow((p1*1.97/p42)/0.0798,2))"   
"spoff26=bf1*(cnst1/1000000)-o1" 
;"spoal26=0.5"

"p43=4.6/(cnst2*bf1/1000000)" 
"spw28=plw1*(pow((p1*1.04/p43)/0.06103,2))"   
"spoff28=bf1*(cnst1/1000000)-o1" 
# 39 "/nmrdata/nmrlib/pp/B_HNCO"
"spw29=plw1*(pow((p1*1.04/p43)/0.06103,2))"   

;"spoal28=0"
"spoff29=bf1*(cnst1/1000000)-o1" 
;"spoal29=1.0"

"p44 =p1*8.0"    
"spoff30=0.0" 
"spw30=plw1"  
# 53 "/nmrdata/nmrlib/pp/B_HNCO"
"p13=4.6/(60.0*bf2/1000000)" 
"spw2=plw2*(pow((p3*1.04/p13)/0.06103,2))"   
"spw8=plw2*(pow((p3*1.04/p13)/0.06103,2))"   
;"spoal2=0"
;"spoal8=1"

"p14=4.875/(70.0*bf2/1000000)" 
"spw3=plw2*(pow((p3*1.97/p14)/0.0798,2))"   
"spw5=plw2*(pow((p3*1.97/p14)/0.0798,2))"   
;"spw7=plw2*(pow((p3*1.97/p14)/0.0798,2))"   

# 65 "/nmrdata/nmrlib/pp/B_HNCO"
"cnst21 = o2/bf2"    
"cnst22 = 39.0"     

"spoff2=0"
"spoff3=0"
"spoff5=bf2*((cnst22-cnst21)/1000000)" 
"spoff7=bf2*((cnst21-cnst22)/1000000)"   
"spoff8=0"
"spoff9=0"

# 78 "/nmrdata/nmrlib/pp/B_HNCO"
"p50 =500u"    
"spoff50=0.0" 
"spw50=plw3*(pow((p21*8/p50),2))"   

"p51=4.875/(40*bf3/1000000)" 
"spw51=plw3*(pow((p21*1.97/p51)/0.0798,2))"   
"spoff51=0.0"                
;"spoal51=0.5"

# 90 "/nmrdata/nmrlib/pp/B_HNCO"
"d11=30m"

"d25=2.5m"
"d26=2.7m"

# 98 "/nmrdata/nmrlib/pp/B_HNCO"
"d27=0"


# 106 "/nmrdata/nmrlib/pp/B_HNCO"
"d23=15m-p43"

"DELTA4=d23-d27-p16-d16-p14*2.0-p44+p21*4/PI"
"DELTA9=0"

# 119 "/nmrdata/nmrlib/pp/B_HNCO"
"DELTA1=d25-p41*0.5-p42*0.5+p50*0.5"
"DELTA2=d26-p16-d16-p42*0.5"
"DELTA3=d26-p17-d16-p43*0.5-p42*0.5"
"if (p50 > p14) DELTA5=p50-p14"
"if (p50 <= p14) DELTA5=0"
"DELTA6=d23-p51*0.5"
"DELTA7=d23-p51*0.5-d27-p44"
"DELTA8=p16+d16+de+8u"


# 133 "/nmrdata/nmrlib/pp/B_HNCO"
"d0=3u"
"in0=inf1/2"


# 140 "/nmrdata/nmrlib/pp/B_HNCO"
"d10=3u"
"d30=d23+3u"

"d29=p43"   

"in10=inf2/2"
# 147 "/nmrdata/nmrlib/pp/B_HNCO"
"FACTOR2=d30*10000000*2/td2"
"in30=FACTOR2/10000000"

"if ( in30 > in10 ) { in29 = 0; } else { in29=in10-in30; }"
"if ( in30 > in10 ) { in30 = in10; }"

# 155 "/nmrdata/nmrlib/pp/B_HNCO"
"acqt0=0"
baseopt_echo

aqseq 321

# 1 "mc_line 160 file /nmrdata/nmrlib/pp/B_HNCO exp. def. part of mc cmd. before ze"
; dimension 3D; AQ_mode  (F2) Echo-Antiecho (F1) States-TPPI
define delay MCWRK
define delay MCREST
define delay d0orig
"d0orig=d0"
define delay d10orig
"d10orig=d10"
define delay d29orig
"d29orig=d29"
define delay d30orig
"d30orig=d30"
define loopcounter t2loop
"t2loop=0"
define loopcounter t1loop
"t1loop=0"
define loopcounter ph1loop
"ph1loop=0"
define loopcounter ph2loop
"ph2loop=0"
define loopcounter ST1CNT
"ST1CNT = td1 / 2"
define loopcounter ST2CNT
"ST2CNT = td2 / 2"
"MCREST = d11 - d11"
"MCWRK = 0.033333*d11"

    dccorr
# 160 "/nmrdata/nmrlib/pp/B_HNCO"
1 d11 ze
# 1 "mc_line 160 file /nmrdata/nmrlib/pp/B_HNCO exp. def. of mc cmd. after ze"
LBLAV, MCWRK
      "d0=d0orig + t1loop * in0 "
      "phval14 = (ph1loop % 4) * 90"
      MCWRK ip14 + phval14
      "d10=d10orig + t2loop * in10 "
      "d29=d29orig + t2loop * in29 "
      "d30=d30orig - t2loop * in30 "
      "phval5 = (ph2loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval10 = (t2loop * 180)%360"
      MCWRK ip10 + phval10
      "phval31 = (t2loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph2loop * 1"
      MCWRK setgrad EA
# 161 "/nmrdata/nmrlib/pp/B_HNCO"
# 1 "mc_line 161 file /nmrdata/nmrlib/pp/B_HNCO exp. start label for mc cmd."
2 MCWRK  * 7 do:f3
LBLF2, MCWRK  * 7
LBLST2, MCWRK  * 8
LBLF1, MCWRK  * 7
LBLST1, MCWRK 
  MCREST
# 162 "/nmrdata/nmrlib/pp/B_HNCO"
3 5u
  d1 pl1:f1 pl2:f2 pl3:f3
  (p21 ph1):f3
  p16:gp7
  d16
  (p21 ph3):f3
  p16:gp7
  d16
  
  50u setnmr3|0 setnmr0|34|32|33 ctrlgrad 0
# 175 "/nmrdata/nmrlib/pp/B_HNCO"
  (p41:sp25 ph1)    
 
  DELTA1
  (center (p42:sp26 ph2) (p50:sp50 ph1):f3 )

  DELTA1
  (p41:sp25 ph2):f1   

  p16:gp3
  d16 pl3:f3
# 188 "/nmrdata/nmrlib/pp/B_HNCO"
 30u fq=cnst21(bf ppm):f2  

  (p21 ph1):f3
  DELTA6
  (center (p14:sp3 ph1):f2 (p51:sp51 ph1):f3 )  
  DELTA7 pl3:f3
  (p44:sp30 ph1)
  d27

# 199 "/nmrdata/nmrlib/pp/B_HNCO"
  (p21 ph2):f3


  (p44:sp30 ph1)
  p16:gp4
  d16
# 208 "/nmrdata/nmrlib/pp/B_HNCO"
  (p13:sp2 ph14):f2   
  d0
  (center (p14:sp5 ph1):f2 (p50:sp50 ph1):f3 )   
  d0
  (p14:sp3 ph1):f2
  4u
  (p14:sp5 ph1):f2   
  DELTA5
  10u
  (p13:sp8 ph1):f2   
  
  4u
  (p44:sp30 ph1)
  p17:gp3
  d16 pl3:f3


# 226 "/nmrdata/nmrlib/pp/B_HNCO"
  (p21 ph10):f3
  (p50:sp50 ph1):f3
  d10
  (p14:sp5 ph1):f2   
  DELTA9   
  (p44:sp30 ph1)
  DELTA4 


# 236 "/nmrdata/nmrlib/pp/B_HNCO"
  p16:gp1*EA

  d16
  (p14:sp3 ph1):f2   
  d29              
  (p50:sp50 ph1):f3
  d30  pl3:f3      
 (p43:sp28 ph1)   

# 249 "/nmrdata/nmrlib/pp/B_HNCO"
 (p21 ph4):f3
  p16:gp5
  d16
  DELTA2

  (center (p42:sp26 ph1) (p51:sp51 ph1):f3 )
  DELTA2
  p16:gp5
  d16 pl3:f3
  (p21 ph5):f3
  (p43:sp29 ph2)    

  p17:gp6
  d16
  DELTA3
  (center (p42:sp26 ph2) (p51:sp51 ph2):f3 )
  DELTA3
  p17:gp6
  d16
  (p43:sp28 ph1)   
# 270 "/nmrdata/nmrlib/pp/B_HNCO"
  DELTA8
  (p42:sp26 ph1)  
  p16:gp2
  d16  pl16:f3
  4u setnmr3^0 setnmr0^34^32^33 ctrlgrad 7


# 279 "/nmrdata/nmrlib/pp/B_HNCO"
  go=2 ph31 cpd3:f3
# 1 "mc_line 280 file /nmrdata/nmrlib/pp/B_HNCO exp. mc cmd. in line"
  MCWRK  do:f3 wr #0 if #0 zd 

  "ph2loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval14 = (ph1loop % 4) * 90"
      MCWRK ip14 + phval14
      "d10=d10orig + t2loop * in10 "
      "d29=d29orig + t2loop * in29 "
      "d30=d30orig - t2loop * in30 "
      "phval5 = (ph2loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval10 = (t2loop * 180)%360"
      MCWRK ip10 + phval10
      "phval31 = (t2loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph2loop * 1"
      MCWRK setgrad EA
  lo to LBLF2 times 2
  MCWRK 

  "t2loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval14 = (ph1loop % 4) * 90"
      MCWRK ip14 + phval14
      "d10=d10orig + t2loop * in10 "
      "d29=d29orig + t2loop * in29 "
      "d30=d30orig - t2loop * in30 "
      "phval5 = (ph2loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval10 = (t2loop * 180)%360"
      MCWRK ip10 + phval10
      "phval31 = (t2loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph2loop * 1"
      MCWRK setgrad EA
  lo to LBLST2 times ST2CNT
  MCWRK 
  "t2loop=0"
  "ph2loop=0"
  MCWRK 

  "ph1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval14 = (ph1loop % 4) * 90"
      MCWRK ip14 + phval14
      "d10=d10orig + t2loop * in10 "
      "d29=d29orig + t2loop * in29 "
      "d30=d30orig - t2loop * in30 "
      "phval5 = (ph2loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval10 = (t2loop * 180)%360"
      MCWRK ip10 + phval10
      "phval31 = (t2loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph2loop * 1"
      MCWRK setgrad EA
  lo to LBLF1 times 2
  MCWRK 

  "t1loop+=1"
      "d0=d0orig + t1loop * in0 "
      "phval14 = (ph1loop % 4) * 90"
      MCWRK ip14 + phval14
      "d10=d10orig + t2loop * in10 "
      "d29=d29orig + t2loop * in29 "
      "d30=d30orig - t2loop * in30 "
      "phval5 = (ph2loop % 2) * 180"
      MCWRK ip5 + phval5
      "phval10 = (t2loop * 180)%360"
      MCWRK ip10 + phval10
      "phval31 = (t2loop * 180)%360"
      MCWRK ip31 + phval31
      MCWRK
      "phval0 = ph2loop * 1"
      MCWRK setgrad EA
  lo to LBLST1 times ST1CNT
  MCWRK 
  "t1loop=0"
  "ph1loop=0"
  MCWRK rf #0
  lo to LBLAV times tdav
# 284 "/nmrdata/nmrlib/pp/B_HNCO"

exit

# 288 "/nmrdata/nmrlib/pp/B_HNCO"
ph1=0
ph2=3 
ph3=2
ph4=0
ph5=1
ph10=0 0 2 2
ph14=0 2
ph18=0 0 0 0 2 2 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)
;sp26: Reburp.1000
;sp25: Pc9_4_90.1000
;sp28: Eburp2.1000
;sp29: Eburp2tr.1000
;sp30 BIP
;p16: homospoil/gradient pulse                         [1 msec]
;p21: f3 channel -  90 degree high power pulse
;p22: f3 channel - 180 degree high power pulse
;p8 : f2 channel - 180 degree shaped pulse for inversion (adiabatic)
;p17: gradient pulse 3                                 [300 usec]
;p41: PC9
;p42: REBURP
;p43: EBURP
;d0 : incremented delay (F1)                           [3 usec]
;d1 : relaxation delay; 1-5 * T1
;d11: delay for disk I/O                               [30 msec]
;d16: delay for homospoil/gradient recovery
;d25: 1/(4J(NH)                                     
;d26: 1/(4J(NH)     
;d27: 1/(4J(NH)                            
;cnst1: H(N) excitation frequency (in ppm)
;cnst2: H(N) excitation band width (in ppm)
;cnst22: Calpha chemical shift offset (in ppm)
;cnst26: Call chemical shift (offset, in ppm)          [101 ppm]
;inf1: 1/SW(N) = 2 * DW(N)
;in0: 1/(2 * SW(N)) = DW(N)
;nd0: 2
;ns: 2 * n
;ds: >= 16
;td1: number of experiments in F1
# 334 "/nmrdata/nmrlib/pp/B_HNCO"
;for z-only gradients:

;use gradient files:   
;gpnam1: SMSQ10.32
;gpnam2: SMSQ10.32
;gpnam3: SMSQ10.32
;gpnam4: SMSQ10.32
;gpnam5: SMSQ10.32
;gpnam6: SMSQ10.32
;gpnam7: SMSQ10.32

# 347 "/nmrdata/nmrlib/pp/B_HNCO"
;Processing

;PHC0(F1): 45.0

# 353 "/nmrdata/nmrlib/pp/B_HNCO"
;$Id: b_trosyetf3gpsi,v 1.2.2.1.4.1 2012/01/31 17:56:18 ber Exp $
