pulsesequence()
{
void makeHHdec(), makeCdec(); /* utility functions */
int ihh=1, /* used in HH decoupling to improve water suppression */
t1_counter;
char C13refoc[MAXSTR], /* C13 sech/tanh pulse in middle of t1 */
Hdecflg[MAXSTR], /* HH-h**o decoupling flag */
Cdecflg[MAXSTR], /* low power C-13 decoupling flag */
TROSY[MAXSTR],
wtg3919[MAXSTR];
double tauxh, tau1,
pwNt = 0.0, /* pulse only active in the TROSY option */
gsign = 1.0,
/* temporary Pbox parameters */
bw, pws, ofs, ppm, nst, /* bandwidth, pulsewidth, offset, ppm, # steps */
gzlvl3=getval("gzlvl3"),
gt3=getval("gt3"),
JNH = getval("JNH"),
pwN = getval("pwN"),
pwNlvl = getval("pwNlvl"),
pwHs, tpwrs=0.0, compH=1.0, /* H1 90 degree pulse length at tpwrs */
sw1 = getval("sw1"),
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
rfst = 4095.0, /* fine power for the stCall pulse */
compC = getval("compC"), /* adjustment for C13 amplifier compr-n */
tpwrsf = getval("tpwrsf"); /* adjustment for soft pulse power*/
/* INITIALIZE VARIABLES */
getstr("C13refoc",C13refoc);
getstr("TROSY",TROSY);
getstr("wtg3919",wtg3919);
getstr("Hdecflg", Hdecflg);
getstr("Cdecflg", Cdecflg);
if(wtg3919[0] != 'y') /* selective H20 one-lobe sinc pulse needs 1.69 */
{ /* times more power than a square pulse */
pwHs = getval("pwHs");
compH = getval("compH");
}
else
pwHs = pw*2.385+7.0*rof1+d3*2.5;
tauxh = ((JNH != 0.0) ? 1/(4*(JNH)) : 2.25e-3);
setautocal(); /* activate auto-calibration flags */
if (autocal[0] == 'n')
{
if (C13refoc[A]=='y')
{
/* 180 degree adiabatic C13 pulse from 0 to 200 ppm */
rfst = (compC*4095.0*pwC*4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35));
rfst = (int) (rfst + 0.5);
if ( 1.0/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35)) < pwC )
{ text_error( " Not enough C13 RF. pwC must be %f usec or less.\n",
(1.0e6/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35))) ); psg_abort(1); }
}
if(wtg3919[0] != 'y') /* selective H20 one-lobe sinc pulse needs 1.69 */
{ /* times more power than a square pulse */
if (pwHs > 1e-6) tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69));
else tpwrs = 0.0;
tpwrs = (int) (tpwrs);
}
}
else /* if autocal = 'y'(yes), 'q'(quiet), r(read), or 's'(semi) */
{
if(FIRST_FID) /* call Pbox */
{
if (C13refoc[A]=='y')
{
ppm = getval("dfrq"); ofs = 0.0; pws = 0.001; /* 1 ms long pulse */
bw = 200.0*ppm; nst = 1000; /* nst - number of steps */
stC200 = pbox_makeA("stC200", "sech", bw, pws, ofs, compC*pwC, pwClvl, nst);
C13ofs = 100.0;
}
if(wtg3919[0] != 'y')
H2Osinc = pbox_Rsh("H2Osinc", "sinc90", pwHs, 0.0, compH*pw, tpwr);
ofs_check(H1ofs, C13ofs, N15ofs, H2ofs);
}
if (C13refoc[A]=='y') rfst = stC200.pwrf;
if (wtg3919[0] != 'y')
{
pwHs = H2Osinc.pw; tpwrs = H2Osinc.pwr-1.0; /* 1dB correction applied */
}
}
if (tpwrsf<4095.0) tpwrs = tpwrs + 6.0;
tauxh = ((JNH != 0.0) ? 1/(4*(JNH)) : 2.25e-3);
if(Cdecflg[0] == 'y') makeCdec(); /* make shapes for HH h**o-decoupling */
if(Hdecflg[0] == 'y') makeHHdec();
if(Hdecflg[0] != 'n') ihh = -3;
/* check validity of parameter range */
if((dm[A] == 'y' || dm[B] == 'y' || dm[C] == 'y' ))
{ text_error("incorrect Dec1 decoupler flags! "); psg_abort(1); }
if((dm2[A] == 'y' || dm2[B] == 'y') )
{ text_error("incorrect Dec2 decoupler flags! "); psg_abort(1); }
if( dpwr > 0 )
{ text_error("don't fry the probe, dpwr too large! "); psg_abort(1); }
if( dpwr2 > 50 )
{ text_error("don't fry the probe, dpwr2 too large! "); psg_abort(1); }
if ((TROSY[A]=='y') && (dm2[C] == 'y'))
{ text_error("Choose either TROSY='n' or dm2='n' ! "); psg_abort(1); }
/* LOAD VARIABLES */
if(ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5);
tau1 = d2/2.0 - pw;
if(tau1 < 0.0) tau1 = 0.0;
/* LOAD PHASE TABLES */
settable(t6, 4, recT);
if (TROSY[A] == 'y')
{ gsign = -1.0;
pwNt = pwN;
assign(zero,v7);
assign(two,v8);
settable(t1, 1, phT1);
settable(t2, 4, phT2);
settable(t3, 1, phT4);
settable(t4, 1, phT4);
settable(t5, 4, recT); }
else
{ assign(one,v7);
assign(three,v8);
settable(t1, 4, phi1);
settable(t2, 2, phi2);
settable(t3, 8, phi3);
settable(t4, 1, phi4);
settable(t5, 8, rec); }
if ( phase1 == 2 ) /* Hypercomplex in t1 */
{ if (TROSY[A] == 'y')
{ tsadd(t3, 2, 4); tsadd(t5, 2, 4); }
else tsadd(t2, 1, 4); }
if(t1_counter %2) /* calculate modification to phases based on */
{ tsadd(t2,2,4); tsadd(t5,2,4); tsadd(t6,2,4); } /* current t1 values */
if(wtg3919[0] != 'y')
{ add(one,v7,v7); add(one,v8,v8); }
/* sequence starts!! */
status(A);
obspower(tpwr);
dec2power(pwNlvl);
decpower(pwClvl);
decpwrf(rfst);
if(Hdecflg[0] != 'n')
{
delay(5.0e-5);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(1.5*gzlvl3, 0.5e-3);
delay(5.0e-4);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(-gzlvl3, 0.5e-3);
}
delay(d1);
rcvroff();
status(B);
rgpulse(pw, zero, rof1, rof1);
zgradpulse(0.3*gzlvl3,gt3);
txphase(zero);
dec2phase(zero);
delay(tauxh-gt3); /* delay=1/4J(XH) */
sim3pulse(2*pw,0.0,2*pwN,t4,zero,zero,rof1,rof1);
zgradpulse(0.3*gzlvl3,gt3);
dec2phase(t2);
delay(tauxh-gt3 ); /* delay=1/4J(XH) */
rgpulse(pw, t1, rof1, rof1);
zgradpulse(0.5*gsign*ihh*gzlvl3,gt3);
delay(200.0e-6);
decphase(zero);
if (TROSY[A] == 'y')
{
txphase(t3);
if ( phase1 == 2 )
dec2rgpulse(pwN, t6, rof1, 0.0);
else
dec2rgpulse(pwN, t2, rof1, 0.0);
if ( (C13refoc[A]=='y') && (d2 > 1.0e-3 + 2.0*WFG2_START_DELAY) )
{
delay(d2/2.0 - 0.5e-3 - WFG2_START_DELAY);
decshaped_pulse("stC200", 1.0e-3, zero, 0.0, 0.0);
delay(d2/2.0 - 0.5e-3 - WFG2_STOP_DELAY);
}
else
delay(d2);
rgpulse(pw, t3, 0.0, rof1);
zgradpulse(0.65*gzlvl3,gt3);
delay(tauxh-gt3 );
sim3pulse(2*pw,0.0,2*pwN,zero,zero,zero,rof1,rof1);
zgradpulse(0.65*gzlvl3,gt3);
delay(tauxh-gt3 );
sim3pulse(pw,0.0,pwN,zero,zero,t3,rof1,rof1);
}
else
{
txphase(t4);
dec2rgpulse(pwN, t2, rof1, 0.0);
if ( (C13refoc[A]=='y') && (tau1 > 0.5e-3 + WFG2_START_DELAY) )
{
delay(tau1 - 0.5e-3 - WFG2_START_DELAY);
simshaped_pulse("", "stC200", 2.0*pw, 1.0e-3, t4, zero, 0.0, 0.0);
dec2phase(t3);
delay(tau1 - 0.5e-3 - WFG2_STOP_DELAY);
}
else
{
delay(tau1);
rgpulse(2.0*pw, t4, 0.0, 0.0);
dec2phase(t3);
delay(tau1);
}
dec2rgpulse(pwN, t3, 0.0, 0.0);
zgradpulse(0.5*gzlvl3,gt3);
delay(200.0e-6);
rgpulse(pw, two, rof1, rof1);
}
zgradpulse(gzlvl3,gt3);
txphase(v7); dec2phase(zero);
delay(tauxh-gt3-pwHs-rof1+5.0e-5);
if(wtg3919[0] == 'y')
{
rgpulse(pw*0.231,v7,rof1,rof1);
delay(d3);
rgpulse(pw*0.692,v7,rof1,rof1);
delay(d3);
rgpulse(pw*1.462,v7,rof1,rof1);
delay(d3/2-pwN);
dec2rgpulse(2*pwN, zero, rof1, rof1);
txphase(v8);
delay(d3/2-pwN);
rgpulse(pw*1.462,v8,rof1,rof1);
delay(d3);
rgpulse(pw*0.692,v8,rof1,rof1);
delay(d3);
rgpulse(pw*0.231,v8,rof1,rof1);
}
else
{
obspower(tpwrs); if (tpwrsf<4095.0) obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHs, v7, rof1, 0.0);
obspower(tpwr); if (tpwrsf<4095.0) obspwrf(4095.0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, v8, zero, zero, 0.0, 0.0);
obspower(tpwrs); if (tpwrsf<4095.0)obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHs, v7, rof1, 0.0);
obspower(tpwr); if (tpwrsf<4095.0)obspwrf(4095.0);
}
zgradpulse(gzlvl3,gt3);
if(Cdecflg[0] == 'y')
{
delay(tauxh-gt3-pwHs-rof1-pwNt-3.0*POWER_DELAY-PRG_START_DELAY);
dec2rgpulse(pwNt, zero, rof1, rof1);
dec2power(dpwr2);
rcvron();
statusdelay(C,5.0e-5);
setreceiver(t5);
pbox_decon(&Cdseq);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
else
{
delay(tauxh-gt3-pwHs-rof1-pwNt-POWER_DELAY);
dec2rgpulse(pwNt, zero, rof1, rof1);
dec2power(dpwr2);
rcvron();
statusdelay(C,5.0e-5);
setreceiver(t5);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
}
pulsesequence()
{
void makeHHdec(), makeCdec(); /* utility functions */
int ihh=1, /* used in HH decouling to improve water suppression */
t1_counter;
char C13refoc[MAXSTR], /* C13 sech/tanh pulse in middle of t1 */
Hdecflg[MAXSTR], /* HH-h**o decoupling flag */
Cdecflg[MAXSTR], /* low power C-13 decoupling flag */
IPAP[MAXSTR], /* Flag for anti-phase spectrum */
wtg3919[MAXSTR];
double tauxh, tau1,
pwNt = 0.0, /* pulse only active in the TROSY option */
gsign = 1.0,
maxHpwr = 51.0, /* maximum allowed H-H decoupling power */
maxCpwr = 45.0, /* maximum allowed C-H decoupling power */
gzlvl0=getval("gzlvl0"),
gzlvl1=getval("gzlvl1"),
gzlvl2=getval("gzlvl2"),
gzlvl3=getval("gzlvl3"),
gt0=getval("gt0"),
gt1=getval("gt1"),
gt2=getval("gt2"),
gt3=getval("gt3"),
JNH = getval("JNH"),
pwN = getval("pwN"),
pwNlvl = getval("pwNlvl"),
pwHs, tpwrs=0.0, compH, /* H1 90 degree pulse length at tpwrs */
sw1 = getval("sw1"),
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
rfst = 4095.0, /* fine power for the stCall pulse */
compC = getval("compC"); /* adjustment for C13 amplifier compr-n */
/* INITIALIZE VARIABLES */
getstr("C13refoc",C13refoc);
getstr("IPAP",IPAP); /* IPAP = 'y' for AP; IPAP = 'n' for IP */
getstr("wtg3919",wtg3919);
getstr("Hdecflg", Hdecflg);
getstr("Cdecflg", Cdecflg);
tauxh = ((JNH != 0.0) ? 1/(4*(JNH)) : 2.25e-3);
if (C13refoc[A]=='y') /* 180 degree adiabatic C13 pulse from 0 to 200 ppm */
{
rfst = (compC*4095.0*pwC*4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35));
rfst = (int) (rfst + 0.5);
if ( 1.0/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35)) < pwC )
{
text_error( " Not enough C13 RF. pwC must be %f usec or less.\n",
(1.0e6/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35))) );
psg_abort(1);
}
}
if(wtg3919[0] != 'y') /* selective H20 one-lobe sinc pulse needs 1.69 */
{ pwHs = getval("pwHs"); /* times more power than a square pulse */
compH = getval("compH");
if (pwHs > 1e-6) tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69));
else tpwrs = 0.0;
tpwrs = (int) (tpwrs); }
else
pwHs = pw*2.385+7.0*rof1+d3*2.5;
if(Cdecflg[0] == 'y') makeCdec(maxCpwr); /* make shapes for HH h**o-decoupling */
if(Hdecflg[0] == 'y') makeHHdec(maxHpwr);
if(Hdecflg[0] != 'n') ihh = -3;
/* check validity of parameter range */
if((dm[A] == 'y' || dm[B] == 'y' || dm[C] == 'y' ))
{ text_error("incorrect Dec1 decoupler flags! "); psg_abort(1); }
if((dm2[A] == 'y' || dm2[B] == 'y') )
{ text_error("incorrect Dec2 decoupler flags! "); psg_abort(1); }
if( dpwr > 0 )
{ text_error("don't fry the probe, dpwr too large! "); psg_abort(1); }
if( dpwr2 > 50 )
{ text_error("don't fry the probe, dpwr2 too large! "); psg_abort(1); }
/* LOAD VARIABLES */
if(ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5);
tau1 = d2/2.0;
if(tau1 < 0.0) tau1 = 0.0;
/* LOAD PHASE TABLES */
settable(t1, 2, phi1);
if (IPAP[A]=='y') settable(t2,4,phi2A);
else settable(t2, 4, phi2);
settable(t3, 16, phi3);
settable(t4, 16, phi4);
if (IPAP[A]=='y') settable(t5,8,recA);
else settable(t5, 4, rec);
assign(one,v7);
assign(three,v8);
if ( phase1 == 2 ) /* Hypercomplex in t1 */
{
if (IPAP[A] == 'y')
tsadd(t3,1,4);
tsadd(t2, 1, 4);
}
if(t1_counter %2) /* calculate modification to phases based on */
{ /* current t1 values */
tsadd(t2,2,4);
tsadd(t5,2,4);
if (IPAP[A] == 'y')
tsadd(t3,2,4);
}
if(wtg3919[0] != 'y')
{ add(one,v7,v7); add(one,v8,v8); }
/* sequence starts!! */
status(A);
obspower(tpwr);
dec2power(pwNlvl);
decpower(pwClvl);
decpwrf(rfst);
if(Hdecflg[0] != 'n')
{
delay(5.0e-5);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(1.5*gzlvl0, 0.5e-3);
delay(5.0e-4);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(-gzlvl0, 0.5e-3);
}
delay(d1);
rcvroff();
status(B);
dec2rgpulse(pwN,zero,rof1,rof1);
zgradpulse(gzlvl0,gt0);
delay(1e-3);
rgpulse(pw, zero, rof1, rof1);
zgradpulse(gzlvl1,gt1);
txphase(zero);
dec2phase(zero);
delay(tauxh-gt1); /* delay=1/4J(XH) */
sim3pulse(2*pw,0.0,2*pwN,zero,zero,zero,rof1,rof1);
zgradpulse(gzlvl1,gt1);
dec2phase(t2);
delay(tauxh-gt1 ); /* delay=1/4J(XH) */
rgpulse(pw, t1, rof1, rof1);
zgradpulse(0.5*gsign*ihh*gzlvl2,gt2);
delay(200.0e-6);
decphase(zero);
txphase(t4);
dec2rgpulse(pwN, t2, rof1, 0.0);
if (IPAP[A] == 'y')
{
delay(tauxh-pwN);
sim3pulse(2*pw,0.0,2*pwN,zero,zero,t3,rof1,rof1);
delay(tauxh-pwN);
rgpulse(pw,t4,rof1,rof1);
}
if ( (C13refoc[A]=='y') && (tau1 > 0.5e-3 + WFG2_START_DELAY) )
{
delay(tau1 - 0.5e-3 - WFG2_START_DELAY);
simshaped_pulse("", "stC200", 0.0, 1.0e-3, zero, zero, 0.0, 0.0);
dec2phase(zero);
delay(tau1 - 0.5e-3 - WFG2_STOP_DELAY);
}
else
delay(2.0*tau1);
dec2rgpulse(pwN, zero, 0.0, 0.0);
zgradpulse(0.5*gzlvl2,gt2);
delay(200.0e-6);
rgpulse(pw, zero, rof1, rof1);
zgradpulse(gzlvl3,gt3);
txphase(v7); dec2phase(zero);
delay(tauxh-gt3-pwHs-rof1+5.0e-5);
if(wtg3919[0] == 'y')
{
rgpulse(pw*0.231,v7,rof1,rof1);
delay(d3);
rgpulse(pw*0.692,v7,rof1,rof1);
delay(d3);
rgpulse(pw*1.462,v7,rof1,rof1);
delay(d3/2-pwN);
dec2rgpulse(2*pwN, zero, rof1, rof1);
txphase(v8);
delay(d3/2-pwN);
rgpulse(pw*1.462,v8,rof1,rof1);
delay(d3);
rgpulse(pw*0.692,v8,rof1,rof1);
delay(d3);
rgpulse(pw*0.231,v8,rof1,rof1);
}
else
{
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, v7, rof1, 0.0);
obspower(tpwr);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, v8, zero, zero, 0.0, 0.0);
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, v7, rof1, 0.0);
obspower(tpwr);
}
zgradpulse(gzlvl3,gt3);
if(Cdecflg[0] == 'y')
{
delay(tauxh-gt3-pwHs-rof1-pwNt-3.0*POWER_DELAY-PRG_START_DELAY);
dec2rgpulse(pwNt, zero, rof1, rof2);
dec2power(dpwr2);
rcvron();
statusdelay(C,5.0e-5);
setreceiver(t5);
pbox_decon(&Cdseq);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
else
{
delay(tauxh-gt3-pwHs-rof1-pwNt-POWER_DELAY);
dec2rgpulse(pwNt, zero, rof1, rof2);
dec2power(dpwr2);
rcvron();
statusdelay(C,5.0e-5);
setreceiver(t5);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
}
void pulsesequence()
{
char c1d[MAXSTR]; /* option to record only 1D C13 spectrum */
int ncyc;
double tau1 = 0.002, /* t1 delay */
post_del = 0.0,
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
compC = getval("compC"),
compH = getval("compH"),
mixpwr = getval("mixpwr"),
jCH = getval("jCH"),
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gzlvl0 = getval("gzlvl0"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
grec = getval("grec"),
phase = getval("phase");
getstr("c1d",c1d);
ncyc=1;
if(jCH > 0.0)
tau1 = 0.25/jCH;
dbl(ct, v1); /* v1 = 0 */
mod4(v1,oph);
hlv(ct,v2);
add(v2,v1,v2);
if (phase > 1.5)
incr(v1); /* hypercomplex phase increment */
initval(2.0*(double)((int)(d2*getval("sw1")+0.5)%2),v10);
add(v1,v10,v1);
add(oph,v10,oph);
mod4(v1,v1); mod4(v2,v2); mod4(oph,oph);
assign(zero,v3);
if(FIRST_FID)
{
HHmix = pbox_mix("HHmix", "DIPSI2", mixpwr, pw*compH, tpwr);
if(c1d[A] == 'n')
{
opx("CHdec"); setwave("WURST2 30k/1.2m"); pbox_par("steps","600"); cpx(pwC*compC, pwClvl);
CHdec = getDsh("CHdec");
}
}
ncyc = (int) (at/HHmix.pw) + 1;
post_del = ncyc*HHmix.pw - at;
/* BEGIN PULSE SEQUENCE */
status(A);
zgradpulse(gzlvl0, gt0);
rgpulse(pw, zero, 0.0, 0.0); /* destroy H-1 magnetization*/
zgradpulse(gzlvl0, gt0);
delay(1.0e-4);
obspower(tpwr);
txphase(v1);
decphase(zero);
dec2phase(zero);
presat();
obspower(tpwr);
delay(1.0e-5);
status(B);
if(c1d[A] == 'y')
{
rgpulse(pw,v1,0.0,0.0); /* 1H pulse excitation */
delay(d2);
rgpulse(pw,two,0.0,0.0); /* 1H pulse excitation */
assign(oph,v3);
}
else
{
decunblank(); pbox_decon(&CHdec);
rgpulse(pw,v1,0.0,0.0); /* 1H pulse excitation */
txphase(zero);
delay(d2);
pbox_decoff(); decblank();
decpower(pwClvl); decpwrf(4095.0);
delay(tau1 - POWER_DELAY);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
txphase(one); decphase(one); dec2phase(one);
delay(tau1);
simpulse(pw, pwC, one, one, 0.0, 0.0);
txphase(zero); decphase(zero); dec2phase(zero);
delay(tau1);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
delay(tau1);
simpulse(0.0, pwC, zero, zero, 0.0, 0.0);
}
zgradpulse(gzlvl1, gt1);
delay(grec);
simpulse(0.0, pwC, zero, v3, 0.0, rof2);
txphase(v2);
obsunblank(); pbox_xmtron(&HHmix);
status(C);
setactivercvrs("ny");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
delay(post_del);
pbox_xmtroff(); obsblank();
zgradpulse(gzlvl2, gt2);
obspower(tpwr);
delay(grec);
rgpulse(pw,zero,0.0,rof2); /* 1H pulse excitation */
status(D);
setactivercvrs("yn");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
}
pulsesequence()
{
/* DECLARE AND LOAD VARIABLES */
void makeHHdec(), makeCdec(); /* utility functions */
char f1180[MAXSTR], /* Flag to start t1 @ halfdwell */
mag_flg[MAXSTR], /* magic-angle coherence transfer gradients */
C13refoc[MAXSTR], /* C13 sech/tanh pulse in middle of t1*/
NH2only[MAXSTR], /* spectrum of only NH2 groups */
T1[MAXSTR], /* insert T1 relaxation delay */
T1rho[MAXSTR], /* insert T1rho relaxation delay */
T2[MAXSTR], /* insert T2 relaxation delay */
TROSY[MAXSTR], /* do TROSY on N15 and H1 */
Hdecflg[MAXSTR], /* HH-h**o decoupling flag */
Cdecflg[MAXSTR]; /* low power C-13 decoupling flag */
int icosel, /* used to get n and p type */
ihh=1, /* used in HH decouling to improve water suppression */
t1_counter, /* used for states tppi in t1 */
rTnum, /* number of relaxation times, relaxT */
rTcounter; /* to obtain maximum relaxT, ie relaxTmax */
double tau1, /* t1 delay */
lambda = 0.91/(4.0*getval("JNH")), /* 1/4J H1 evolution delay */
tNH = 1.0/(4.0*getval("JNH")), /* 1/4J N15 evolution delay */
relaxT = getval("relaxT"), /* total relaxation time */
rTarray[1000], /* to obtain maximum relaxT, ie relaxTmax */
maxrelaxT = getval("maxrelaxT"), /* maximum relaxT in all exps */
ncyc, /* number of pulsed cycles in relaxT */
pwr_dly, /* power delay */
/* the sech/tanh pulse is automatically calculated by the macro "proteincal", */
/* and is called directly from your shapelib. */
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
rf0, /* maximum fine power when using pwC pulses */
rfst, /* fine power for the stCall pulse */
compH = getval("compH"), /* adjustment for H1 amplifier compression */
compN = getval("compN"), /* adjustment for N15 amplifier compression */
compC = getval("compC"), /* adjustment for C13 amplifier compression */
calH = getval("calH"), /* multiplier on a pw pulse for H1 calibration */
tpwrsf = getval("tpwrsf"), /* fine power adustment for soft pulse */
pwHs = getval("pwHs"), /* H1 90 degree pulse length at tpwrs */
pwHH = 0.0, /* pwHH = pwHs for HH h**o-decoupling */
tpwrs, /* power for the pwHs ("H2Osinc") pulse */
pwNlvl = getval("pwNlvl"), /* power for N15 pulses */
pwN = getval("pwN"), /* N15 90 degree pulse length at pwNlvl */
calN = getval("calN"), /* multiplier on a pwN pulse for calibration */
slNlvl, /* power for N15 spin lock */
slNrf = 1500.0, /* RF field in Hz for N15 spin lock at 600 MHz */
sw1 = getval("sw1"),
gt1 = getval("gt1"), /* coherence pathway gradients */
gzcal = getval("gzcal"), /* dac to G/cm conversion */
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
BPpwrlimits, /* =0 for no limit, =1 for limit */
gt0 = getval("gt0"), /* other gradients */
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gstab = getval("gstab"),
gzlvl0 = getval("gzlvl0"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5");
P_getreal(GLOBAL,"BPpwrlimits",&BPpwrlimits,1);
getstr("f1180",f1180);
getstr("mag_flg",mag_flg);
getstr("C13refoc",C13refoc);
getstr("NH2only",NH2only);
getstr("T1",T1);
getstr("T1rho",T1rho);
getstr("T2",T2);
getstr("TROSY",TROSY);
getstr("Hdecflg", Hdecflg);
getstr("Cdecflg", Cdecflg);
/* LOAD PHASE TABLE */
settable(t3,2,phi3);
settable(t4,1,phx);
if (TROSY[A]=='y')
{settable(t1,1,ph_x);
settable(t9,1,phx);
settable(t10,1,phy);
settable(t11,1,phx);
settable(t12,2,recT);}
else
{settable(t1,1,phx);
settable(t9,8,phi9);
settable(t10,1,phx);
settable(t11,1,phy);
settable(t12,4,rec);}
/* INITIALIZE VARIABLES */
/* maximum fine power for pwC pulses (and initialize rfst) */
rf0 = 4095.0; rfst=0.0;
/* 180 degree adiabatic C13 pulse from 0 to 200 ppm */
if (C13refoc[A]=='y')
{rfst = (compC*4095.0*pwC*4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35));
rfst = (int) (rfst + 0.5);
if ( 1.0/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35)) < pwC )
{ text_error( " Not enough C13 RF. pwC must be %f usec or less.\n",
(1.0e6/(4000.0*sqrt((30.0*sfrq/600.0+7.0)/0.35))) ); psg_abort(1); }}
/* selective H20 one-lobe sinc pulse */
if(pwHs > 1e-6)
tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69)); /* needs 1.69 times more */
else /* power than a square pulse */
tpwrs = 0.0;
tpwrs = (int) (tpwrs);
if (tpwrsf<4095.0) tpwrs = tpwrs + 6.0;
if (tpwrsf < 4095.0)
{
tpwrs = tpwrs + 6.0;
pwr_dly = POWER_DELAY + PWRF_DELAY;
}
else pwr_dly = POWER_DELAY;
/* power level for N15 spinlock (90 degree pulse length calculated first) */
slNlvl = 1/(4.0*slNrf*sfrq/600.0) ;
slNlvl = pwNlvl - 20.0*log10(slNlvl/(pwN*compN));
slNlvl = (int) (slNlvl + 0.5);
/* use 1/8J times for relaxation measurements of NH2 groups */
if ( (NH2only[A]=='y') && ((T1[A]=='y') || (T1rho[A]=='y') || (T2[A]=='y')) )
{ tNH = tNH/2.0; }
/* reset calH and calN for 2D if inadvertently left at 2.0 */
if (ni>1.0) {calH=1.0; calN=1.0;}
/* make shapes and set up parameters for HH h**o-decoupling */
if(Cdecflg[0] == 'y') makeCdec();
if(Hdecflg[0] == 'y') makeHHdec();
if(Hdecflg[0] != 'n')
{
pwHH = pwHs;
pwHs = 0.0;
}
/* CHECK VALIDITY OF PARAMETER RANGES */
if ((TROSY[A]=='y') && (gt1 < -2.0e-4 + pwHs + 1.0e-4 + 2.0*POWER_DELAY))
{ text_error( " gt1 is too small. Make gt1 equal to %f or more.\n",
(-2.0e-4 + pwHs + 1.0e-4 + 2.0*POWER_DELAY) ); psg_abort(1); }
if((dm[A] == 'y' || dm[B] == 'y' || dm[C] == 'y' ))
{ text_error("incorrect dec1 decoupler flags! Should be 'nnn' "); psg_abort(1); }
if((dm2[A] == 'y' || dm2[B] == 'y'))
{ text_error("incorrect dec2 decoupler flags! Should be 'nny' "); psg_abort(1); }
if( dpwr2 > 50 )
{ text_error("don't fry the probe, DPWR2 too large! "); psg_abort(1); }
if( pw > 50.0e-6 )
{ text_error("dont fry the probe, pw too high ! "); psg_abort(1); }
if( pwN > 100.0e-6 )
{ text_error("dont fry the probe, pwN too high ! "); psg_abort(1); }
/* RELAXATION TIMES AND FLAGS */
/* evaluate maximum relaxT, relaxTmax chosen by the user */
rTnum = getarray("relaxT", rTarray);
relaxTmax = rTarray[0];
for (rTcounter=1; rTcounter<rTnum; rTcounter++)
if (relaxTmax < rTarray[rTcounter]) relaxTmax = rTarray[rTcounter];
/* compare relaxTmax with maxrelaxT */
if (maxrelaxT > relaxTmax) relaxTmax = maxrelaxT;
if ( ((T1rho[A]=='y') || (T2[A]=='y')) && (relaxTmax > d1) )
{ text_error("Maximum relaxation time, relaxT, is greater than d1 ! "); psg_abort(1);}
if ( ((T1[A]=='y') && (T1rho[A]=='y')) || ((T1[A]=='y') && (T2[A]=='y')) ||
((T1rho[A]=='y') && (T2[A]=='y')) )
{ text_error("Choose only one relaxation measurement ! "); psg_abort(1); }
if ( ((T1[A]=='y') || (T1rho[A]=='y')) &&
((relaxT*100.0 - (int)(relaxT*100.0+1.0e-4)) > 1.0e-6) )
{ text_error("Relaxation time, relaxT, must be zero or multiple of 10msec"); psg_abort(1);}
if ( (T2[A]=='y') &&
(((relaxT+0.01)*50.0 - (int)((relaxT+0.01)*50.0+1.0e-4)) > 1.0e-6) )
{ text_error("Relaxation time, relaxT, must be odd multiple of 10msec"); psg_abort(1);}
if ( ((T1rho[A]=='y') || (T2[A]=='y')) && (relaxTmax > 0.25) && (ix==1) )
{ printf("WARNING, sample heating will result for relaxT>0.25sec"); }
if ( ((T1rho[A]=='y') || (T2[A]=='y')) && (relaxTmax > 0.5) )
{ text_error("relaxT greater than 0.5 seconds will heat sample"); psg_abort(1);}
if ( ((NH2only[A]=='y') || (T1[A]=='y') || (T1rho[A]=='y') || (T2[A]=='y'))
&& (TROSY[A]=='y') )
{ text_error("TROSY not implemented with NH2 spectrum, or relaxation exps."); psg_abort(1);}
if ((TROSY[A]=='y') && (dm2[C] == 'y'))
{ text_error("Choose either TROSY='n' or dm2='n' ! "); psg_abort(1); }
/* PHASES AND INCREMENTED TIMES */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (TROSY[A]=='y')
{ if (phase1 == 2) icosel = -1;
else { tsadd(t4,2,4); tsadd(t10,2,4); icosel = +1; }
}
else { if (phase1 == 2) {tsadd(t10,2,4); icosel = +1;}
else icosel = -1;
}
if(Hdecflg[0] != 'n') ihh = icosel;
/* Set up f1180 */
tau1 = d2;
if((f1180[A] == 'y') && (ni > 1.0))
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1/2.0;
/* Calculate modifications to phases for States-TPPI acquisition */
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t3,2,4); tsadd(t12,2,4); }
/* Correct inverted signals for NH2 only spectra */
if ((NH2only[A]=='y') && (T1[A]=='n') && (T1rho[A]=='n') && (T2[A]=='n'))
{ tsadd(t3,2,4); }
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
decpwrf(rf0);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
if(Hdecflg[0] != 'n')
{
delay(5.0e-5);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(1.5*gzlvl0, 0.5e-3);
delay(5.0e-4);
rgpulse(pw,zero,rof1,0.0);
rgpulse(pw,one,0.0,rof1);
zgradpulse(-gzlvl0, 0.5e-3);
}
delay(d1);
/* xxxxxxxxxxxxxxxxx CONSTANT SAMPLE HEATING FROM N15 RF xxxxxxxxxxxxxxxxx */
if (T1rho[A]=='y')
{dec2power(slNlvl);
dec2rgpulse(relaxTmax-relaxT, zero, 0.0, 0.0);
dec2power(pwNlvl);}
if (T2[A]=='y')
{ncyc = 8.0*100.0*(relaxTmax - relaxT);
if (BPpwrlimits > 0.5)
{
dec2power(pwNlvl-3.0); /* reduce for probe protection */
pwN=pwN*compN*1.4;
}
if (ncyc > 0)
{initval(ncyc,v1);
loop(v1,v2);
delay(0.625e-3 - pwN);
dec2rgpulse(2*pwN, zero, 0.0, 0.0);
delay(0.625e-3 - pwN);
endloop(v2);}
if (BPpwrlimits > 0.5)
{
dec2power(pwNlvl); /* restore normal value */
pwN=getval("pwN");
}
}
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
rcvroff();
if (TROSY[A]=='n')
dec2rgpulse(pwN, zero, 0.0, 0.0); /*destroy N15 magnetization*/
zgradpulse(gzlvl0, 0.5e-3);
delay(1.0e-4);
if (TROSY[A]=='n') dec2rgpulse(pwN, one, 0.0, 0.0);
zgradpulse(0.7*gzlvl0, 0.5e-3);
decpwrf(rfst);
txphase(t1);
delay(5.0e-4);
if ( dm3[B] == 'y' ) /* begins optional 2H decoupling */
{
lk_hold();
dec3rgpulse(1/dmf3,one,10.0e-6,2.0e-6);
dec3unblank();
dec3phase(zero);
delay(2.0e-6);
setstatus(DEC3ch, TRUE, 'w', FALSE, dmf3);
}
rgpulse(calH*pw,t1,0.0,0.0); /* 1H pulse excitation */
txphase(zero);
dec2phase(zero);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0 - pwHH);
if(Hdecflg[0] != 'n')
{
obspower(tpwrs);
if (tpwrsf<4095.0) obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHH, two, 5.0e-5, 0.0);
obspower(tpwr);
if (tpwrsf<4095.0) obspwrf(4095.0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
obspower(tpwrs);
if (tpwrsf<4095.0) obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHH, two, 5.0e-5, 0.0);
obspower(tpwr);
if (tpwrsf<4095.0) obspwrf(4095.0);
}
else
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
txphase(one);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0 - pwHH);
rgpulse(pw, one, 0.0, 0.0);
txphase(two);
obspower(tpwrs);
if (tpwrsf<4095.0) obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHs, two, 5.0e-5, 0.0);
obspower(tpwr);
if (tpwrsf<4095.0) obspwrf(4095.0);
if (TROSY[A]=='y')
zgradpulse(ihh*gzlvl3, gt3);
else
zgradpulse(-ihh*gzlvl3, gt3);
dec2phase(t3);
delay(2.0e-4);
dec2rgpulse(calN*pwN, t3, 0.0, 0.0);
txphase(zero);
decphase(zero);
/* xxxxxxxxxxxxxxxxxx OPTIONS FOR N15 RELAXATION xxxxxxxxxxxxxxxxxxxx */
if ( (T1[A]=='y') || (T1rho[A]=='y') || (T2[A]=='y') )
{
dec2phase(one);
zgradpulse(gzlvl4, gt4); /* 2.0*GRADIENT_DELAY */
delay(tNH - gt4 - 2.0*GRADIENT_DELAY);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, one, 0.0, 0.0);
zgradpulse(gzlvl4, gt4); /* 2.0*GRADIENT_DELAY */
delay(tNH - gt4 - 2.0*GRADIENT_DELAY);
}
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
if (T1[A]=='y')
{
dec2rgpulse(pwN, one, 0.0, 0.0);
dec2phase(three);
zgradpulse(gzlvl0, gt0); /* 2.0*GRADIENT_DELAY */
delay(2.5e-3 - gt0 - 2.0*GRADIENT_DELAY - pw);
rgpulse(2.0*pw, zero, 0.0, 0.0);
delay(2.5e-3 - pw);
ncyc = (100.0*relaxT);
initval(ncyc,v4);
if (ncyc > 0)
{loop(v4,v5);
delay(2.5e-3 - pw);
rgpulse(2.0*pw, two, 0.0, 0.0);
delay(2.5e-3 - pw);
delay(2.5e-3 - pw);
rgpulse(2.0*pw, zero, 0.0, 0.0);
delay(2.5e-3 - pw);
endloop(v5);}
dec2rgpulse(pwN, three, 0.0, 0.0);
}
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
/* Theory suggests 8.0 is better than 2PI as RF */
/* field multiplier and experiment confirms this.*/
if (T1rho[A]=='y') /* Shift evolution of 2.0*pwN/PI for one pulse */
{ /* at end left unrefocused as for normal sequence*/
delay(1.0/(8.0*slNrf) - pwN);
decrgpulse(pwN, zero, 0.0, 0.0);
dec2power(slNlvl);
/* minimum 5ms spinlock to dephase */
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0); /* spins not locked */
sim3pulse(2.0*pw, 0.0, 2.0*pw, zero, zero, zero, 0.0, 0.0);
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0);
ncyc = 100.0*relaxT;
initval(ncyc,v4); if (ncyc > 0)
{loop(v4,v5);
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0);
sim3pulse(2.0*pw, 0.0, 2.0*pw, two, zero, zero, 0.0, 0.0);
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0);
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0);
sim3pulse(2.0*pw, 0.0, 2.0*pw, zero, zero, zero, 0.0, 0.0);
dec2rgpulse((2.5e-3-pw), zero, 0.0, 0.0);
endloop(v5);}
dec2power(pwNlvl);
decrgpulse(pwN, zero, 0.0, 0.0);
delay(1.0/(8.0*slNrf) + 2.0*pwN/PI - pwN);
}
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
if (T2[A]=='y')
{
dec2phase(zero);
initval(0.0,v3); initval(180.0,v4);
if (BPpwrlimits > 0.5)
{
dec2power(pwNlvl-3.0); /* reduce for probe protection */
pwN=pwN*compN*1.4;
}
ncyc = 100.0*relaxT;
initval(ncyc,v5);
loop(v5,v6);
initval(3.0,v7);
loop(v7,v8);
delay(0.625e-3 - pwN);
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(0.625e-3 - pwN);
endloop(v8);
delay(0.625e-3 - pwN - SAPS_DELAY);
add(v4,v3,v3); obsstepsize(1.0); xmtrphase(v3); /* SAPS_DELAY */
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(0.625e-3 - pwN - pw);
rgpulse(2*pw, zero, 0.0, 0.0);
delay(0.625e-3 - pwN - pw );
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
xmtrphase(zero); /* SAPS_DELAY */
delay(0.625e-3 - pwN - SAPS_DELAY);
initval(3.0,v9);
loop(v9,v10);
delay(0.625e-3 - pwN);
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(0.625e-3 - pwN);
endloop(v10);
endloop(v6);
if (BPpwrlimits > 0.5)
{
dec2power(pwNlvl); /* restore normal value */
pwN=getval("pwN");
}
}
/* xxxxxxxxxxxxxxxxxx OPTIONS FOR N15 EVOLUTION xxxxxxxxxxxxxxxxxxxxx */
txphase(zero);
dec2phase(t9);
if ( (NH2only[A]=='y') || (T1[A]=='y') || (T1rho[A]=='y') || (T2[A]=='y') )
{
delay(tau1);
/* optional sech/tanh pulse in middle of t1 */
if (C13refoc[A]=='y') /* WFG_START_DELAY */
{decshaped_pulse("stC200", 1.0e-3, zero, 0.0, 0.0);
delay(tNH - 1.0e-3 - WFG_START_DELAY - 2.0*pw);}
else
{delay(tNH - 2.0*pw);}
rgpulse(2.0*pw, zero, 0.0, 0.0);
if (tNH < gt1 + 1.99e-4) delay(gt1 + 1.99e-4 - tNH);
delay(tau1);
dec2rgpulse(2.0*pwN, t9, 0.0, 0.0);
if (mag_flg[A] == 'y') magradpulse(gzcal*gzlvl1, gt1);
else zgradpulse(gzlvl1, gt1); /* 2.0*GRADIENT_DELAY */
txphase(t4);
dec2phase(t10);
if (tNH > gt1 + 1.99e-4) delay(tNH - gt1 - 2.0*GRADIENT_DELAY);
else delay(1.99e-4 - 2.0*GRADIENT_DELAY);
}
else if (TROSY[A]=='y')
{
if ( (C13refoc[A]=='y') && (tau1 > 0.5e-3 + WFG2_START_DELAY) )
{delay(tau1 - 0.5e-3 - WFG2_START_DELAY); /* WFG2_START_DELAY */
decshaped_pulse("stC200", 1.0e-3, zero, 0.0, 0.0);
delay(tau1 - 0.5e-3);}
else delay(2.0*tau1);
if (mag_flg[A] == 'y') magradpulse(gzcal*gzlvl1, gt1);
else zgradpulse(gzlvl1, gt1); /* 2.0*GRADIENT_DELAY */
delay(2.0e-4 - 2.0*GRADIENT_DELAY);
dec2rgpulse(2.0*pwN, t9, 0.0, 0.0);
txphase(three);
delay(gt1 + 2.0e-4 - pwHs - 1.0e-4 - 2.0*pwr_dly);
obspower(tpwrs);
if (tpwrsf<4095.0) obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHs, three, 5.0e-5, 0.0);
obspower(tpwr);
if (tpwrsf<4095.0) obspwrf(4095.0);
txphase(t4);
delay(5.0e-5);
}
else
{ /* fully-coupled spectrum */
if (dm2[C]=='n') {rgpulse(2.0*pw, zero, 0.0, 0.0); pw=0.0;}
if ( (C13refoc[A]=='y') && (tau1 > 0.5e-3 + WFG2_START_DELAY) )
{delay(tau1 - 0.5e-3 - WFG2_START_DELAY); /* WFG2_START_DELAY */
simshaped_pulse("", "stC200", 2.0*pw, 1.0e-3, zero, zero, 0.0, 0.0);
delay(tau1 - 0.5e-3);
delay(gt1 + 2.0e-4);}
else
{delay(tau1);
rgpulse(2.0*pw, zero, 0.0, 0.0);
delay(gt1 + 2.0e-4 - 2.0*pw);
delay(tau1);}
pw=getval("pw");
dec2rgpulse(2.0*pwN, t9, 0.0, 0.0);
if (mag_flg[A] == 'y') magradpulse(gzcal*gzlvl1, gt1);
else zgradpulse(gzlvl1, gt1); /* 2.0*GRADIENT_DELAY */
txphase(t4);
dec2phase(t10);
delay(2.0e-4 - 2.0*GRADIENT_DELAY);
}
if (T1rho[A]=='y') delay(POWER_DELAY);
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
if (TROSY[A]=='y') rgpulse(pw, t4, 0.0, 0.0);
else sim3pulse(pw, 0.0, pwN, t4, zero, t10, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
zgradpulse(gzlvl5, gt5);
if (TROSY[A]=='y') delay(lambda - 0.65*(pw + pwN) - gt5);
else delay(lambda - 1.3*pwN - gt5);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
txphase(one);
dec2phase(t11);
delay(lambda - 1.3*pwN - gt5);
sim3pulse(pw, 0.0, pwN, one, zero, t11, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
zgradpulse(1.5*gzlvl5, gt5);
delay(lambda - 1.3*pwN - gt5);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t10);
zgradpulse(1.5*gzlvl5, gt5);
if (TROSY[A]=='y') delay(lambda - 1.6*pwN - gt5);
else delay(lambda - 0.65*pwN - gt5);
if (TROSY[A]=='y') dec2rgpulse(pwN, t10, 0.0, 0.0);
else rgpulse(pw, zero, 0.0, 0.0);
delay((gt1/10.0) + 1.0e-4 +gstab - 0.65*pw + 2.0*GRADIENT_DELAY + POWER_DELAY);
if ( dm3[B] == 'y' ) /* turns off 2H decoupling */
{
setstatus(DEC3ch, FALSE, 'c', FALSE, dmf3);
dec3rgpulse(1/dmf3,three,2.0e-6,2.0e-6);
dec3blank();
lk_autotrig(); /* resumes lock pulsing */
}
rgpulse(2.0*pw, zero, 0.0, 0.0);
dec2power(dpwr2); /* POWER_DELAY */
if (mag_flg[A] == 'y') magradpulse(icosel*gzcal*gzlvl2, 0.1*gt1);
else zgradpulse(icosel*gzlvl2, 0.1*gt1); /* 2.0*GRADIENT_DELAY */
if(Cdecflg[0] == 'y')
{
delay(gstab-2.0*POWER_DELAY-PRG_START_DELAY+rof2);
rcvron();
statusdelay(C,1.0e-4);
if (dm3[B] == 'y')
{
delay(1/dmf3);
lk_sample();
}
setreceiver(t12);
pbox_decon(&Cdseq);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
else
{
delay(gstab+rof2);
rcvron();
statusdelay(C,1.0e-4);
if (dm3[B] == 'y')
{
delay(1/dmf3);
lk_sample();
}
setreceiver(t12);
if(Hdecflg[0] == 'y')
homodec(&HHdseq);
}
}
