-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathcest_zz_H.cw
200 lines (164 loc) · 3.63 KB
/
cest_zz_H.cw
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
;1H CEST with saturation during zz period
;based on Sekhar et al. PNAS 113, E2794-801 (2016)
;
;Chris Waudby ([email protected])
prosol relations=<triple>
#include <Avance.incl>
#include <Grad.incl>
#include <Delay.incl>
;list of CEST saturation frequencies
;first line in file should be zero, indicating the reference plane
;give values in Hz relative to sfo1
define list<frequency> H1sat = <$FQ1LIST>
"p2=p1*2"
"p22=p21*2"
"d11=30m"
"d12=20u"
"d13=4u"
"d26=1s/(cnst4*4)"
"in0=inf2"
"DELTA=d26-larger(p21,p2)-p16-d16-4u"
"DELTA1=d26-p21-p16-d16-8u-p11"
"DELTA2=d26-p21-p16-d16-12u-p11"
"acqt0=0"
; initial delay for half dwell phase correction
"d0=in0/2-1.27324*p21-2*p2-8u"
"l2=1"
aqseq 312
1 ze
d11 pl16:f3
2 d1 do:f3
d12 pl1:f1 pl3:f3
50u UNBLKGRAD ; crush equilibrium Nz (p16, gp4)
(p21 ph10):f3
p16:gp1
d16*2
(p10:sp10 ph13:r):f1 ; H2O flip-down (sp10)
4u pl1:f1
(p1 ph10):f1 ; start INEPT
4u
p16:gp2
d16
DELTA
(center (p1 ph10 p2 ph11 p1 ph10):f1 (p22 ph10):f3 )
DELTA
4u
p16:gp2
d16
(p1 ph11):f1
; zz purge
4u
p16:gp3
d16
; 15N evolution
(p21 ph1):f3
2u
d0*0.5 gron4
2u groff
(p1 ph10 p2 ph11 p1 ph10):f1
2u
d0*0.5 gron4*-1
2u groff
(p21 ph2):f3
;goto 999 ;optimise sp10
; zz purge
4u
p16:gp5
d16
; CEST period
4u pl9:f1
4u H1sat:f1
if "l2==1" goto 77
4u LOCKH_ON
d18 cw:f1 ph10
4u do:f1
4u LOCKH_OFF
77 4u pl1:f1
4u fq=0:f1
; zz purge
4u
p16:gp6
d16
4u pl0:f1
; retro-INEPT
(p11:sp11 ph12):f1 ; H2O flip-down (sp11)
4u pl1:f1
(p1 ph10)
;goto 999 ; optimise sp11
4u
p16:gp7
d16
DELTA1
(p11:sp12 ph12):f1
4u pl1:f1
(center (p2 ph10):f1 (p22 ph10):f3)
4u pl0:f1
(p11:sp12 ph12):f1
DELTA2
p16:gp7
d16
999 4u BLKGRAD
4u pl16:f3
; acquisition
go=2 ph31 cpds3:f3
d1 do:f3 mc #0 to 2
F1QF(H1sat.inc & iu2)
F2PH(ru2 & ip2, id0)
exit
ph1= 0 2
ph2= 0
ph10= 0
ph11= 1
ph12= 2
ph13= 3
ph31= 0 2
;pl0 : 120 dB
;pl1 : f1 channel - power level for pulse (default)
;pl3 : f3 channel - power level for pulse (default)
;pl9 : f1 channel - power level for CEST saturation
;pl16: f3 channel - power level for CPD/BB decoupling
;sp10: 90 deg Squa1.1000 water flip-down pulse
;sp11: 90 deg Squa1.1000 water flip-down pulse
;sp12: 90 deg Squa1.1000 water flip-down pulse
;p1 : f1 channel - 90 degree high power pulse
;p2 : f1 channel - 180 degree high power pulse
;p10: 1000u 90 deg soft rectangular water flip-back pulse
;p11: 1000u 90 deg soft rectangular water flip-back pulse
;p16: homospoil/gradient pulse [1 ms]
;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]
;d13: short delay [4 usec]
;d16: delay for homospoil/gradient recovery
;d18: saturation time
;d26 : 1/(4J(YH))
;cnst4: = J(YH)
;inf2: 1/SW(X) = 2 * DW(X)
;in0: 1/(2 * SW(X)) = DW(X)
;nd0: 2
;NS: 2 * n
;DS: 16
;td1: number of experiments
;FnMODE: States-TPPI (or TPPI)
;cpd3: decoupling according to sequence defined by cpdprg3
;pcpd3: f3 channel - 90 degree pulse for decoupling sequence
;for z-only gradients:
;gpz1: 37%
;gpz2: 23%
;gpz3: 49%
;gpz4: 0.5%
;gpz5: 31%
;gpz6: 57%
;gpz7: 11%
;use gradient files:
;gpnam1: SQUA.100
;gpnam2: SQUA.100
;gpnam3: SQUA.100
;gpnam5: SQUA.100
;gpnam6: SQUA.100
;gpnam7: SQUA.100
;preprocessor-flags-start
;preprocessor-flags-end