# 1 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"
;c13noe3d.fa
;3D (1H,13C,1H) NOESY-HSQC-SE with optional 15N decoupling
;Bruker Avance/Xwin-nmr version
;Written up by F. Abildgaard, NMRFAM (abild@nmrfam.wisc.edu)
;
;  $Id: pulseprogram,v 1.1 2021/02/11 13:54:21 bmrbsvc Exp $
;
;  Disclaimer: This pulse program is provided "as is" for your
;  information. Support for the use of this pulse program is only
;  provided to users of the National Magnetic Resonance Facility
;  at Madison (NMRFAM). Users of this pulse program employ it at
;  their own risk. Neither NMRFAM nor University of Wisconsin-Madison
;  are liable for any physical or other damage incurred during the
;  use of this pulse program.
;
;f1: 1H, f2: 15N, f3: 13C (channel assignments may be changed below)
;o1p: 4.7ppm
;o2p: 118ppm
;o3p: 43ppm
;
;p1 90 H at pl1
;p3 90 C at pl3
;d8 mixing time
;pl0 120dB
;
;C13 GARP decoupling (cpdprg3), pcpd 90 dgr pulse at pl13
; or WURST decoupling (cpdprg3), pcpd 90 dgr pulse at sp2 (define WURST below)
; For side-band suppression, define SBSUPR below (8 step phase cycle)
;N15 Waltz-16 (cpdprg2), pcpd 90 dgr pulse at pl12
;C evolution:
;  in0, SW(C)=1/(2*in0), typically 24ppm
;  l4 complex points
;  Define C_90180 below for DW/2 initial delay (90,-180 phase)
;  Process as Echo/Anti-Echo
;H1 evolution:
;  in10, SW(H)=1/(2*in10), typically 8ppm
;  set cnst10 to 0 (preferably) or 1 to make d10 the smallest possible
;  positive delay. cnst10=0 gives (90,-180) phase distortion in F1.
;  cnst10=1 gives (270,-540) phase distortion in F1 (use LP to correct).
;  Processing: extend the FID backwards cnst10 points 
;  and apply (90,-180) phase correction
;  In addition a -45dgr phase correction is required (ph0).
;  l6 complex points
;  Process as States (although data is acquired as States-TPPI)
;
;ns=4*n (or 8*n if SBSUPR defined), ds=4, 8, ...
;
;Recommendations for gradients:
;gpz1:  12%
;gpz2:   7%
;gpz3:  24%
;gpz4:  29%
;gpx5:  54% adjust for magic angle
;gpz5:  30%
;gpz6:  24%
;gpz7:  15%
;gpx8:  54% adjust for magic angle
;gpz8:  30%
;gpnam1: sine.50
;gpnam2: sine.100
;gpnam3: sine.32
;gpnam4: sine.100
;gpnam5: sine.100
;gpnam6: sine.32
;gpnam7: sine.32
;gpnam8: sine.20
;
;#define ONE_D		; uncomment for 1D experiment





;#define SBSUPR
;#define WURST
;#define MESSERLE
;#define OPTIM_P19	; uncomment if you want to optimize p19 (GRAD8)

;			; the correct expt time (works with XWIN-NMR 2.x)
;
;Define channel assignments:



;
;You shouldn't have to worry about anything beyond this point :-)
;
;sanity checks
;
# 93 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

# 96 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

;
define delay TAUA3
define delay TAUB6
define delay TAUA7
define delay TAUC
define delay DELTA
define delay CEN_HC1
define delay CEN_HC2
define pulse H1_90
define pulse H1_180
define pulse C13_90
define pulse C13_180
define pulse C13_59
define pulse C13_298
define pulse C13_240
define pulse GRAD1
define pulse GRAD2
define pulse GRAD3
define pulse GRAD4
define pulse GRAD5
define pulse GRAD6
define pulse GRAD7
define pulse GRAD8
;
"d11=100m"                      ;disk i/o
"d12=10u"                       ;power switching etc.
"d13=5u"                        ;a short delay
"d14=60u"			;ip,id etc.
"d16=300u"                      ;gradient recovery
"H1_90=p1"
"H1_180=H1_90*2"
"C13_90=p3"
"C13_180=C13_90*2"
"C13_59=C13_90*59.4/90.0"
"C13_240=C13_90*240.0/90.0"
"C13_298=C13_90*298.0/90.0"

"GRAD1=500u"
"GRAD2=1.0m"
"GRAD3=250u"
"GRAD4=1.0m"
"GRAD5=720u"
"GRAD6=250u"
"GRAD7=300u"

  "p19=181u"

"GRAD8=p19"


  "d0=d13+in0/2"
# 150 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

;"d2=1.866m"		;1/4J
;"d3=0.933m"		;1/8J
"d2=1.6m"		;1/4J
"d3=0.8m"		;1/8J
"d6=H1_90"
"d7=C13_90"
"TAUA3=d2-GRAD3-d13"
"TAUB6=d3-GRAD6-d13"
"TAUA7=d2-GRAD7-d13"
"TAUC=GRAD5+d16+d13-d13*2-H1_180"
# 163 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  "DELTA=GRAD8+d16+d13*2+d12"

"d18=d8-C13_90-d13-GRAD2-d16"

  "d10=((cnst10*2+1)*in10-(H1_90*1.273+C13_90*4+d13*2))/2"



"CEN_HC1=(C13_90*2+C13_180+d13*2-H1_180)/2"
"CEN_HC2=(C13_59*2+C13_298+d13*2-H1_180)/2"



# 182 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"



  "d31=2*(TAUA3+GRAD6+TAUB6+GRAD7+TAUA7)+GRAD1+GRAD2+GRAD3+GRAD4+TAUC+GRAD5+DELTA+GRAD8"




# 1 "/u/exp/stan/nmr/lists/pp//Avance.incl"
;Avance.incl
;
;version  99/08/02





























# 190 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

# 1 "/u/exp/stan/nmr/lists/pp//Grad.incl"
;Grad.incl  -  include file for Gradient Spectroscopy
;avance-version
;version 99/02/19
















define list<gradient> EA=<EA>
# 191 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

1 ze
2 d13
  d14 do:f3
  d11 setnmr2^0
  3m
  d14
3 d14
  d14
  d14
  d14
  d14
4 d14
  d14
  d14
5 d14
  d14
  d14
6 d13


# 1 "/u/exp/stan/nmr/lists/pp//faexptcorr.incl"
; This is here to trick "expt" to calculate the experiment time
; correctly.
  if "1 == 1" goto nod31
  d31
nod31, d13
# 212 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  d1 pl1:f1 pl12:f2 pl3:f3
  d13 setnmr2|0
  d13 setnmr0|34
  (C13_90 ph0):f3
  d13
  GRAD1:gp1		; 500u, 8G/cm, sine.50
  d16

  d13 cpd2:f2


  (H1_90 ph5):f1
  d10
  (C13_90 ph0 d13 C13_180 ph1 d13 C13_90 ph0):f3
  d10
# 231 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  (H1_90 ph0):f1

  d13 do:f2

  d18
  (C13_90 ph0):f3
  d13
  GRAD2:gp2		;1m, 5G/cm, sine.100
  d16
  (H1_90 ph0):f1
  d13
  GRAD3:gp3		;250u, 16G/cm, sine.32
  TAUA3
  (CEN_HC1 H1_180 ph0):f1 (C13_90 ph0 d13 C13_180 ph1 d13 C13_90 ph0):f3
  d13
  GRAD3:gp3		;250u, 16G/cm, sine.32
  TAUA3
# 252 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  (H1_90 ph1):f1
  d13
  GRAD4:gp4		;1m, 20G/cm, sine.100
;#ifdef 	; uncomment for N15 decoupling
;  d16 cpd2:f2
;#else
  d16
;#endif
  (C13_90 ph7):f3
  d0 
  (H1_180 ph0):f1
  d0 
;#ifdef 	; uncomment for N15 decoupling
;  TAUC do:f2
;#else
  TAUC
;#endif
  (C13_180 ph0):f3
  d13
  GRAD5:gp5*EA*-1	;720u, -/+40G/cm, sine.100
  d16
  (C13_90 ph8):f3
  d6
  (H1_90 ph0):f1
  d13
  GRAD6:gp6		;250u, 16G/cm, sine.32
  TAUB6
  (C13_180 ph0):f3
  (H1_180 ph0):f1 
  d13
  GRAD6:gp6		;250u, 16G/cm, sine.32
  TAUB6
  d7
  (C13_90 ph1):f3
  (H1_90 ph1):f1 
  d13
  GRAD7:gp7		;300u, 10G/cm, sine.32
  TAUA7                                                       
  (CEN_HC1 H1_180 ph0):f1 (C13_90 ph0 d13 C13_180 ph1 d13 C13_90 ph0):f3
  d13
  GRAD7:gp7		;300u, 10G/cm, sine.32
  TAUA7 
  (H1_90 ph0):f1
  DELTA
  (H1_180 ph0):f1
  d13
  GRAD8:gp8		;181u, 40G/cm, sine.20
  d13
  d16 setnmr0^34

# 305 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

# 308 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  d12 pl13:f3


  go=2 ph31 cpds3:f3

# 316 "/u/exp/stan/nmr/lists/pp/c13noe3d.gc"

  d11 do:f3 wr #0 if #0 zd

  d13 setnmr2^0

  3m igrad EA
  d14 ip8
  d14 ip8
  lo to 3 times 2
  d14 id0
  d14 ip7
  d14 ip7
  d14 ip31
  d14 ip31
  lo to 4 times l4
  d14 rd0


  d14 ip5
  d14 ip5
  lo to 5 times 2
  d14 id10
  d14 ip31
  d14 ip31
  lo to 6 times l6


exit
  
ph0=0
ph1=1
ph2=2
ph3=3
ph4=0 2
ph5=(8) 1 5
ph7=0 0 2 2
ph8=0
ph13=0 0 0 0 2 2 2 2
ph31=0 2 2 0