;ha_hbhn_caconhtrwg.fl

;$OWNER=doetsch
#include <Avance.incl>
#include <Grad.incl>

"p2=p1*2"
"p22=p21*2"

"p11=3.873*1000000/((cnst22-cnst21)*bf2*4)"
"p12=1.732*1000000/((cnst22-cnst21)*bf2*2)"

"in0=inf1"
"d0=in0/2-p1*1.273"

;d9:DELTA
"d4=1.75m-p17"
"d28=d9/2-p17"
"d29=d9/2-1.75m-p12" ; = d28-d4-p12 = d9/2-p17-1.75m+p17-p12
"d24=d9/2-p17-p11-3u-p14" ; = d28-p11-3u-p14
;d2:(zeta-DELTA)/2
"d3=d9/2-p17-d2-p15" ; = d28-d2-p15
"d26=d2*2-p11-p14+p15-p1-d7-4.6m" ; = d24+d2+3u-d3-p1-d7-4.6m
;d7:adjust such that d7+p14+4.6m approximates d26+3.5m, i.e. d7=d26-1.1m-p14

"d14=4.6m-p17"
"d18=d17+3u-d14"

"d25=3u"
"d15=d30-p14-d25"
"d16=d30-p18-500u-p1+p21*1.273"
"d6=600u+d5*2-p7*2+p22-d13*2-d12*2-p2" ;;=600u+d5-(p7-d5-p22)-d13*2-d12*2-p2-p7
"d19=p7*2-d5*2-p22+d13*2+d12*2+p2-p21-p19-8u" ;;=600u-d6-p21-p19-8u

"in10=inf2/4"
"in30=(d30-4u)/(td2/2-1)"
"if ( in30 > in10) { in25=0; } else { in25=in10*(1-in30/in10)*2; }"
"if ( in30 > in10) { in15=in10; } else { in15=in30; }"
"if ( in30 > in10) { in16=in10; } else { in16=in10*(2-in30/in10); }"
"if ( in30 > in10) { in30=in10; }"

"spoff1=bf2*((cnst22)/1000000)-o2"
"spoff2=bf2*((cnst22)/1000000)-o2"
"spoff3=bf2*((cnst21)/1000000)-o2"
"spoff4=bf2*((cnst21)/1000000)-o2"
"spoff5=0"
"spoff6=0"

aqseq 321

1 d11 ze
  d11 pl0:f2 pl3:f3
2 d11 
3 d1 pl1:f1
  50u UNBLKGRAD
  (p11:sp4 ph1):f2
  p17:gp1
  1m

  p1 ph10
  p17:gp2
  d28
  p2 ph10
  d29
  (p12:sp3 ph1):f2
  p17:gp2
  d4
  (center (p1 ph11 d0 p1 ph12):f1 (p4:sp6 ph1):f2 (p22 ph1):f3 )

  (p11:sp4 ph3):f2
  3u
  (p14:sp1 ph1):f2
  p17:gp3
  d24
  p2 ph1
  d2
  (p15:sp5 ph4):f2
  p17:gp3
  d3
  p1 ph2
  d26
  p2 ph1
  d7
  (p14:sp1 ph1):f2
  4.6m
  (p11:sp4 ph1):f2
  p17:gp4
  200u

  (p13:sp2 ph5):f2
  p17:gp5
  d14
  (p12:sp3 ph1):f2
  d18
  (center (p14:sp1 ph1):f2 (p22 ph1):f3 )
  p17:gp5
  d17
  (p12:sp3 ph1):f2
  3u
  (p13:sp2 ph1):f2
  p17:gp6
  200u

  (p21 ph6):f3
  d30 
  (p12:sp3 ph1):f2
  d30
  (p22 ph1):f3
  d25
  (p14:sp1 ph1):f2
  d15
  (p12:sp3 ph1):f2
  d16
  p18:gp7*EA
  500u

  p1 ph8
  d13
  d12 pl0:f1
  (p7:sp7 ph17:r):f1
  p20:gp8
  d8 pl1:f1
  (center (p2 ph1):f1 (p14:sp1 ph1):f2 (p22 ph1):f3 )
  p20:gp8
  d8 pl0:f1
  (p7:sp7 ph7:r):f1
  d13
  d12 pl1:f1 
  (p1 ph1):f1 (p21 ph9):f3
  p16:gp9
  600u pl0:f1
  (p7:sp7 ph7:r):f1 (d5 p22 ph1):f3
  d13
  d12 pl1:f1
  p2 ph1
  d13
  d12 pl0:f1
  (p7:sp7 ph7:r):f1
  p16:gp9
  d6 pl1:f1
  (p21 ph18):f3
  p19:gp10
  d19 
  4u BLKGRAD
  go=2 ph31
  d11 mc #0 to 2
	F1PH(rd30 & rd25 & rd15 & rd16 & ip10 & ip11, id0)
	F2EA(igrad EA & ip8*2 & ip18*2 & ip17*2, dd30 & id25 & id15 & id16 & ip6*2 & ip31*2)		
exit

ph1 = 0  
ph2 = 1
ph3 = 0 2
ph4 = 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1
      2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3
ph5 = 0 0 0 0 2 2 2 2
ph6 = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
      2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
ph7 = 2
ph8 = 1
ph9 = 1
ph10= 0 0 2 2
ph11= 0 0 0 0 2 2 2 2
ph12= 0 0 0 0 0 0 0 0 2 2 2 2 2 2 2 2
ph17= 1
ph18= 0

ph31= 0 2 2 0 2 0 0 2 2 0 0 2 0 2 2 0
      2 0 0 2 0 2 2 0 0 2 2 0 2 0 0 2

;cnst21: Ca chemical shift (offset, in ppm)
;cnst22: CO chemical shift (offset, in ppm)
;o2p: Cb chemical shift

;p16: 800u
;p17: 500u
;p18: 1.4m
;p19: 175u
;p20: 300u

;for z-only gradients:
;gpz1: -9%
;gpz2: 10%
;gpz3: -12%
;gpz4: -7%
;gpz5: 13%
;gpz6: 14%
;gpz7: 78.9%
;gpz8: 8%
;gpz9: 25%
;gpz10: 64%