;h2bc3d.ab
;avance-version (2006-10-18)
;
;H2BC3D:
;Heteronuclear two-bond correlation
;Double constant time experiment, using a single constant time period.
;phase sensitive using Echo/Antiecho gradient selection in F1,
;phase sensitive using States-TPPI in F2,
;uses a third order low-pass J-filter for supression of 1J(XH)
;3D version
;can also be used to setup projection reconstruction experiments
;
;Andrew J. Benie & Ole W. Sørensson
;
;Pulse Program Library
;Carlsberg Research Center
;
; DISCLAIMER
; This pulse sequence is tested and working on our spectrometers, but
; adjustments in pulse program syntax and/or parameter sets may be necessary on
; other hardware configurations or systems. Please pay attention to the
; requirements of your system, when implementing pulse sequences from this
; library. Carlsberg Laboratory or the Instrument Center does not make any
; warranties with respect to the software. In no event will Carlsberg Laboratory
; or the Instrument Center be liable for any loss, damages of any nature arising
; or relating to the use or performance of the software.
;
;Copyright (c), 2006, Carlsberg Research Center
;All rights reserved

prosol relations=<triple>

#include <Avance.incl>
#include <Grad.incl>
#include <Delay.incl>


define list<gradient> EA1 = {1.0000 -1.0000}
define delay LP1
define delay LP2
define delay LP3
define delay LP4
define delay TAUC1
define delay TAUC2
define delay PI


"p2  = p1*2"
"p4  = p3*2"

"d0  = 3u"
"d10 = 1u"
"d11 = 30m"
"d12 = 50u"
"d31 = d1 - d0"

"DELTA = p16 + d16"
"DELTA1= 2*d0 + p2"
"TAU   = (1s/(cnst6 + cnst8))"
"TAUC1  = TAU - d12"
"TAUC2  = TAU - d12*2"
"TAU1  = (1s/((cnst6 + 0.07*(cnst8 - cnst6))*2))"
"TAU3  = (1s/((cnst8 - 0.07*(cnst8 - cnst6))*2))"
"LP1   = TAU3 - DELTA"
"LP2   = TAU - DELTA - p3 - 0.5*larger(p2,p3)"
"LP3   = TAU1 - DELTA - 0.5*larger(p2,p3) - de"
"LP4   = TAU + TAU3 - TAU1 - de"
"in0=inf1/2"
"in10=inf2/2"

"d27   = cnst15/2 - (TAU + p3 + p2/2 + d0)"
"d29   = cnst15/2 - (TAU + 2*DELTA + p4 + p3 + p2 + p2/2 + 3*d0 + 2*d12)"
"d30 = d29 + d10 - d0"
"d28 = d27 - d10 - d0"

aqseq 321

1 ze
2 d11 do:f2
  3m
3 6m
4 6m
5 1u
  "d31 = d1 - d0*2"
  d1 BLKGRAD
  d12 pl12:f2
  d12 cpd2:f2
  (p1 ph12):f1
  d28 
  TAUC1 do:f2 
  d12 pl2:f2
  (p3 ph2):f2
  d0 UNBLKGRAD
  (p2 ph10):f1
  d0
  p16:gp1*EA1
  d16
  (p4 ph3):f2
  p16:gp1*EA1*-1
  d16
  DELTA1
  (p3 ph1):f2
  TAUC2 pl12:f2
  d30 cpd2:f2
  d12 do:f2
  d12 pl2:f2
  (p1 ph11):f1 (p3 ph4):f2
  p16:gp1*0.07545
  d16
  LP1
  (p3 ph5):f2
  p16:gp1*0.17605
  d16
  LP2
  (center (p2 ph10):f1 (p3 ph5):f2)
  p16:gp1*0.7545
  d16 BLKGRAD
  LP3 pl14:f2 
  LP4 cpd2:f2
  go=2 ph31
  d11 do:f2 mc #0 to 2
  F1EA(igrad EA1 & rd10, id0)
  F2PH(ip12, id10)
exit
	
; --- H phases ------------------------
ph10 = 0
ph11 = 1
ph12 = 0
; --- X phases ------------------------
ph1  = 0 2 2 0
ph2  = 0 0 2 2 2 2 0 0
ph3  = 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3
ph4  = {0}*16 {2}*16
ph5  = 0
ph6  = 1
; --- Receiver phases -----------------
ph31 = 0 2


;pl1     :f1 power level for pulse (default)
;pl2     :f2 power level for pulse (default)
;pl12    :f2 power level for cpd/BB decoupling
;p1      :f1 90 degree high power pulse
;p2      :f1 180 degree high power pulse
;p3      :f2 90 degree high power pulse
;p4      :f2 180 degree high power pulse
;p16     :gradient pulse
;
;d0      :increment delay (3usec) F1
;d1      :relaxation delay: 1-5*T1
;d10	 :incremented delay (3D)                  
;d11     :delay for disk I/O			[30 msec]
;d12     :delay for power switching		[20 usec]
;d16     :delay for homospoil/gradient recovery	[200 msec]
;d27     :First decoupling period
;d31     :Actual initial decoupling period (calculated)
;d28     :Second decoupling period
;d24     :Actual second decoupling period (calculated)
;d30     :Adjusted d1 to include heating compensation
;
;TAU     :1/(J(min) + J(max))
;TAU1    :1/2(J(min)+k(J(max)-J(min))
;TAU3    :1/2(J(max)-k(J(max)-J(min))
;DELTA   :Gradient duration, p16 + d16 + 50u
;DELTA1  :Correction for initial t1 period
;
;cnst6   :1J(min), ~125 Hz
;cnst8   :1J(max), ~165 Hz
;cnst15  :T, ~0.018 sec
;
;
;nd0     :2
;nd10	 :2
;td1: number of experiments in F1
;td2: number of experiments in F2
;FnMODE: Echo-AnitEcho in F1
;FnMODE: States-TPPI (or TPPI) in F2
;EA1     :Echo-Antiecho gradient program, [+1 -1]
;
;pcpd2   :Pulse length for cpd on f2		[65 usec]
;cpdprg2 :cpd-program				[garp]
;
;gpnam1  :Gradient pulse 1 file name		[SINE.100]
;gpz1    :Gradient strength 1			[80%]
;
;ns      :2, 4, 8, 16, 32, ....
;ds      :128

;Processing:
;
; F1     : SSB = 2, WDW = SINE, PH_mod = no
; F2     : SSB = 2, WDW = SINE, PH_mod = no
; F3     : SSB = 4, WDW = SINE, PH_mod = pk