本文整理汇总了C++中decshaped_pulse函数的典型用法代码示例。如果您正苦于以下问题:C++ decshaped_pulse函数的具体用法?C++ decshaped_pulse怎么用?C++ decshaped_pulse使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了decshaped_pulse函数的15个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: pulsesequence
void pulsesequence()
{
char pwpat[MAXSTR],p1pat[MAXSTR];
/* equilibrium period */
getstr("pwpat",pwpat);
if (pwpat[0] == '\0')
{
abort_message("no pwpat? ABORT");
}
getstr("p1pat",p1pat);
if (p1pat[0] == '\0')
{
abort_message("no p1pat? ABORT");
}
status(A);
obspower(zero);
decpower(tpwr);
hsdelay(d1);
rcvroff();
status(B);
if (is_y(rfwg[1])) decshaped_pulse(p1pat,p1,zero,rof1,rof2);
else apshaped_decpulse(p1pat,p1,zero,t1,t2,rof1,rof2);
hsdelay(d2);
status(C);
if (is_y(rfwg[1])) decshaped_pulse(pwpat,pw,oph,rof1,rof2);
else apshaped_decpulse(pwpat,pw,oph,t1,t2,rof1,rof2);
}
示例2: dec_c13_shpulse
double dec_c13_shpulse(shape pwxshape, codeint phase)
{
decpower(pwxshape.pwr);
decshaped_pulse(pwxshape.name, pwxshape.pw, phase, 0.0, 0.0);
decpower(getval("pwClvl"));
return(pwxshape.pw + 2.0*POWER_DELAY +WFG_STOP_DELAY + WFG_START_DELAY);
}
示例3: pulsesequence
//.........这里部分代码省略.........
delay(timeTC -tau1 - 202.0e-6 - gt7-2.0*POWER_DELAY-pwC9a-
WFG_START_DELAY-WFG_STOP_DELAY-2.0e-6-SAPS_DELAY);
}
}
else {
if (dm3[B]=='y') {
rgpulse(pwHd,one,0.0,0.0);
txphase(zero);
delay(2.0e-6);
obsprgon("waltz16", pwHd, 90.0);
xmtron();
dec3rgpulse(1/dmf3,one,10.0e-6,2.0e-6);
dec3unblank();
dec3phase(zero);
delay(2.0e-6);
setstatus(DEC3ch, TRUE, 'w', FALSE, dmf3);
delay(timeTC-tau1-202.0e-6-gt7-2.0*WFG_START_DELAY-4.0e-6-1/dmf3-pwHd-
PRG_START_DELAY-2.0*POWER_DELAY-pwC9a-2.0e-6-SAPS_DELAY);
}
else {
delay(2.0e-6);
rgpulse(pwHd,one,0.0,0.0);
txphase(zero);
delay(2.0e-6);
obsprgon("waltz16", pwHd, 90.0);
xmtron();
delay(timeTC-tau1-202.0e-6-gt7-4.0e-6-pwHd-PRG_START_DELAY-
2.0*POWER_DELAY-pwC9a-WFG_START_DELAY-WFG_STOP_DELAY-SAPS_DELAY);
}
}
decpwrf(rf9);
decshaped_pulse("offC9",pwC9a,zero,0.0,0.0);
initval(phshift9, v9);
decstepsize(1.0);
dcplrphase(v9); /* SAPS_DELAY */
} /* %%%%%%%%%%%%%%%%%%ENDING 13Ca Constant Time EVOLUTION %%%%%%%%%%%%%%%%%%*/
decphase(t5);
decpwrf(rf1);
decrgpulse(pwC1, t5, 2.0e-6, 0.0);
dec2phase(t8);
dcplrphase(zero);
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 */
}
if (h1dec[0]=='y')
{
xmtroff();
obsprgoff();
rgpulse(pwHd,three,2.0e-6,0.0);
txphase(one);
}
delay(2.0e-6);
zgradpulse(gzlvl4, gt4);
delay(2.0e-4);
示例4: pulsesequence
//.........这里部分代码省略.........
if (tpwrsf<4095.0)
{
obspower(tpwrs+6.0);
obspwrf(tpwrsf);
shaped_pulse("H2Osinc", pwHs, two, 5.0e-4, 0.0);
obspower(tpwr);
obspwrf(4095.0);
}
else
{
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 5.0e-4, 0.0);
obspower(tpwr);
}
zgradpulse(gzlvl3, gt3);
dec2phase(t3);
decpwrf(rfst);
decoffset(dof100);
delay(2.0e-4);
dec2rgpulse(pwN, t3, 0.0, 0.0);
decphase(zero);
/* xxxxxxxxxxxxxxxxxx OPTIONS FOR N15 EVOLUTION xxxxxxxxxxxxxxxxxxxxx */
txphase(zero);
dec2phase(t9);
if (NH2only[A]=='y')
{
delay(tau2);
/* optional sech/tanh pulse in middle of t2 */
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(tau2);
dec2rgpulse(2.0*pwN, t9, 0.0, 0.0);
if (mag_flg[A] == 'y')
{
magradpulse(gzcal*gzlvl1, gt1);
}
else
{
zgradpulse(gzlvl1, gt1);
}
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 ( (C13refoc[A]=='y') && (tau2 > 0.5e-3 + WFG2_START_DELAY) )
{ delay(tau2 - 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(tau2 - 0.5e-3);
delay(gt1 + 2.0e-4);
示例5: pulsesequence
//.........这里部分代码省略.........
decphase(zero);
txphase(zero);
decpwrf(rf10);
delay(tau2);
/* WFG3_START_DELAY */
sim3shaped_pulse("", "offC10", "", 2.0*pw, pwC10, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
if(pwC10>2.0*pwN) pwZ=0.0; else pwZ=2.0*pwN - pwC10;
delay(tau2);
decpwrf(rf0);
if (mag_flg[A] == 'y')
{
magradpulse(-icosel2*gzcal*gzlvl2, 1.8*gt1);
}
else
{
zgradpulse(-icosel2*gzlvl2, 1.8*gt1);
}
delay(2.02e-4);
decrgpulse(2.0*pwC, zero, 0.0, 0.0);
decpwrf(rf10);
if (mag_flg[A] == 'y')
{
magradpulse(icosel2*gzcal*gzlvl2, 1.8*gt1);
}
else
{
zgradpulse(icosel2*gzlvl2, 1.8*gt1);
}
delay(2.0e-4 + WFG3_START_DELAY + pwZ);
decshaped_pulse("offC10", pwC10, zero, 0.0, 0.0);
decpwrf(rf0);
decrgpulse(pwC, zero, 2.0e-6, 0.0);
zgradpulse(gzlvl5, gt5);
delay(0.5*del1 - gt5);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
txphase(one);
decphase(t10);
delay(0.5*del1 - gt5);
simpulse(pw, pwC, one, t10, 0.0, 0.0);
zgradpulse(gzlvl6, gt5);
txphase(zero);
decphase(zero);
delay(0.5*del2 - gt5);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl6, gt5);
delay(0.5*del2 - gt5);
simpulse(pw, pwC, zero, zero, 0.0, 0.0);
delay(0.5*del - 0.5*pwC);
simpulse(2.0*pw,2.0*pwC, zero, zero, 0.0, 0.0);
if (mag_flg[A] == 'y')
magradpulse(gzcal*gzlvl1, gt1);
示例6: pulsesequence
void pulsesequence()
{
double tpwr180r = getval("tpwr180r"),
pw180r = getval("pw180r"),
pwx180 = getval("dnbippw"),
pwxlvl180 = getval("dnbippwr"),
pp = getval("pp"),
pplvl = getval("pplvl"),
tauC = getval("tauC"),
tau;
int prgcycle=(int)(getval("prgcycle")+0.5);
char pw180ref[MAXSTR],
dnbipshp[MAXSTR],
bipflg[MAXSTR];
getstr("pw180ref", pw180ref);
getstr("dnbipshp", dnbipshp);
getstr("bipflg",bipflg);
tau = 1.0 / (2.0 * (getval("j1xh")));
assign(ct,v17);
assign(zero,v18);
assign(zero,v19);
if (getflag("prgflg") && getflag("satmode") && (prgcycle > 1.5))
{
hlv(ct,v17);
mod2(ct,v18); dbl(v18,v18);
if (prgcycle > 2.5)
{
hlv(v17,v17);
hlv(ct,v19); mod2(v19,v19); dbl(v19,v19);
}
}
settable(t1,4,phs1);
settable(t2,8,phs2);
settable(t3,4,phs3);
getelem(t1,v17,v1);
getelem(t2,v17,v2);
getelem(t3,v17,v3);
assign(v1,oph);
add(oph,v18,oph);
add(oph,v19,oph);
if (getflag("prgflg") && getflag("satmode"))
assign(v3,v2);
/* equilibrium period */
status(A);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
if (getflag("satmode"))
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
if (getflag("slpsat"))
{
shaped_satpulse("relaxD",satdly,v2);
if (getflag("prgflg"))
shaped_purge(v1,v2,v18,v19);
}
else
{
satpulse(satdly,v2,rof1,rof1);
if (getflag("prgflg"))
purge(v1,v2,v18,v19);
}
}
else
delay(d1);
if (getflag("wet"))
wet4(zero,one);
status(B);
obspower(tpwr);
if (bipflg[1] == 'y')
decpower(pwxlvl180);
else
decpower(pplvl);
rgpulse(pw,v1,rof1,rof1);
if (bipflg[1] == 'y')
{
decshaped_pulse(dnbipshp,pwx180,zero,2.0e-6,2.0e-6);
delay(tau+tauC);
}
else
{
decrgpulse(2*pp,zero,2.0e-6,2.0e-6);
delay(tau);
}
obspower(tpwr180r);
//.........这里部分代码省略.........
示例7: pulsesequence
//.........这里部分代码省略.........
/* shaped pulse */
rgpulse(pw,zero,0.0,0.0); /* 90 deg 1H pulse */
delay(0.2e-6);
zgradpulse(gzlvl5,gt5);
delay(2.0e-6);
delay(taua - gt5 - 2.2e-6); /* taua <= 1/4JNH */
sim3pulse(2*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
txphase(three); dec2phase(zero); decphase(zero);
delay(0.2e-6);
zgradpulse(gzlvl5,gt5);
delay(200.0e-6);
delay(taua - gt5 - 200.2e-6 - 2.0e-6);
if (sel_flg[A] == 'n')
{
rgpulse(pw,three,2.0e-6,0.0);
delay(0.2e-6);
zgradpulse(gzlvl3,gt3);
delay(200.0e-6);
dec2rgpulse(pwN,zero,0.0,0.0);
delay( zeta );
dec2rgpulse(2.0*pwN,zero,0.0,0.0);
decshaped_pulse(spchirp,pwchirp,zero,0.0,0.0);
delay(zeta -WFG_START_DELAY -pwchirp -WFG_STOP_DELAY -2.0e-6);
dec2rgpulse(pwN,zero,2.0e-6,0.0);
}
else
{
rgpulse(pw,one,2.0e-6,0.0);
initval(1.0,v3);
dec2stepsize(45.0);
dcplr2phase(v3);
delay(0.2e-6);
zgradpulse(gzlvl3,gt3);
delay(200.0e-6);
dec2rgpulse(pwN,zero,0.0,0.0);
dcplr2phase(zero);
delay(1.34e-3 - SAPS_DELAY - 2.0*pw);
rgpulse(pw,one,0.0,0.0);
rgpulse(2.0*pw,zero,0.0,0.0);
rgpulse(pw,one,0.0,0.0);
delay( zeta - 1.34e-3 - 2.0*pw);
dec2rgpulse(2.0*pwN,zero,0.0,0.0);
decshaped_pulse(spchirp,pwchirp,zero,0.0,0.0);
示例8: pulsesequence
//.........这里部分代码省略.........
/* ===========Beginning of NOE transfer ======= */
delay(mix - gt4 - gt5 - pwN -pwC - 2.0*gstab);
obsoffset(tof);
zgradpulse(gzlvl4,gt4);
delay(gstab);
decrgpulse(pwC, zero, 2.0e-6, 2.0e-6);
dec2rgpulse(pwN, zero, 2.0e-6, 2.0e-6);
zgradpulse(gzlvl5,gt5);
delay(gstab);
/* H2O relaxes back to +z */
/* =========== End of NOE transfer ======= */
rgpulse(pw, zero, 2.0e-6, 2.0e-6);
zgradpulse(gzlvl6,gt6);
delay(tauNH - gt6 - 4.0e-6 - 2.0*GRADIENT_DELAY);
sim3pulse(2.0*pw, 0.0, 2*pwN, zero, zero, zero, 2.0e-6, 2.0e-6);
delay(tauNH - gt6 - gstab -4.0e-6 - 2.0*GRADIENT_DELAY);
zgradpulse(gzlvl6,gt6);
txphase(one);
delay(gstab);
rgpulse(pw,one,2.0e-6,2.0e-6);
txphase(two);
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 2.0e-6, 2.0e-6);
obspower(tpwr);
zgradpulse(gzlvl7,gt7);
dec2phase(t3);
decoffset(dofcaco); /* offset on 120ppm for CaCO decoupling */
decpwrf(rf200); /* fine power for stC200 */
delay(gstab);
dec2rgpulse(pwN,t3,2.0e-6,2.0e-6);
dec2phase(t4);
delay(tau3);
rgpulse(2.0*pw, zero, 2.0e-6, 2.0e-6);
decshaped_pulse("stC200", 1.0e-3, zero, 2.0e-6, 2.0e-6);
delay(tau3);
delay(gt1 +gstab +8.0e-6 - 1.0e-3 - 2.0*pw);
dec2rgpulse(2.0*pwN, t4, 2.0e-6, 2.0e-6);
dec2phase(t5);
zgradpulse(gzlvl1, gt1);
delay(gstab + WFG_START_DELAY + WFG_STOP_DELAY - 2.0*GRADIENT_DELAY);
sim3pulse(pw, 0.0, pwN, zero, zero, t5, 2.0e-6, 2.0e-6);
dec2phase(zero);
zgradpulse(gzlvl8, gt8);
delay(tauNH - gt8 - 2.0*GRADIENT_DELAY -4.0e-6);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6);
delay(tauNH - gt8 - gstab -4.0e-6 - 2.0*GRADIENT_DELAY);
zgradpulse(gzlvl8, gt8);
txphase(one); dec2phase(one);
delay(gstab);
sim3pulse(pw, 0.0, pwN, one, zero, one, 2.0e-6, 2.0e-6);
txphase(zero); dec2phase(zero);
zgradpulse(gzlvl9, gt9);
delay(tauNH - gt9 - 2.0*GRADIENT_DELAY -4.0e-6);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 2.0e-6, 2.0e-6);
delay(tauNH - gt9 - 2.0*GRADIENT_DELAY -gstab -4.0e-6);
zgradpulse(gzlvl9, gt9);
delay(gstab);
rgpulse(pw, zero, 2.0e-6, 2.0e-6);
delay((gt1/10.0) +gstab -2.0e-6 - 0.5*pw + 2.0*GRADIENT_DELAY + POWER_DELAY);
rgpulse(2.0*pw, zero, 2.0e-6, 2.0e-6);
zgradpulse(icosel*gzlvl2, gt1/10.0);
dec2power(dpwr2);
delay(gstab);
status(C);
setreceiver(t11);
}
示例9: pulsesequence
//.........这里部分代码省略.........
decrgpulse(2.0*pwC, one, 0.2e-6, 0.0);
decphase(zero);
decrgpulse(pwC, zero, 0.2e-6, 0.0);
delay(tau1);
rgpulse(pw,zero,0.0,0.0);
delay(mix - gt4 - gt5 - pwN - 2.0e-3);
obsoffset(tof);
zgradpulse(gzlvl4,gt4);
delay(1.0e-3);
sim3pulse((double)0.0,pwC,pwN,zero,zero,zero,0.0,2.0e-6);
zgradpulse(gzlvl5,gt5);
delay(1.0e-3);
rgpulse(pw,zero,0.0,2.0e-6);
zgradpulse(gzlvl6,gt6);
delay(taunh - gt6 - 2.0e-6);
sim3pulse(2.0*pw,(double)0.0,2*pwN,zero,zero,zero,0.0,0.0);
delay(taunh - gt6 - 500.0e-6);
zgradpulse(gzlvl6,gt6);
txphase(one);
delay(500.0e-6);
rgpulse(pw,one,0.0,0.0);
txphase(two);
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 2.0e-6, 2.0e-6);
obspower(tpwr);
zgradpulse(gzlvl7,gt7);
dec2phase(t3);
decoffset(dofcaco);
decpwrf(rfst);
delay(200.0e-6);
dec2rgpulse(pwN,t3,0.0,0.0);
dec2phase(t4);
delay(tau3);
rgpulse(2.0*pw, zero, 0.0, 0.0);
decshaped_pulse("stC200", 1.0e-3, zero, 0.0, 0.0);
delay(tau3);
delay(gt1 + 202.0e-6 - 1.0e-3 - 2.0*pw);
dec2rgpulse(2.0*pwN, t4, 0.0, 2.0e-6);
dec2phase(t5);
if(mag_flg[A] == 'y')
{
magradpulse(gzcal*gzlvl1, gt1);
}
else
{
zgradpulse(gzlvl1, gt1);
delay(4.0*GRADIENT_DELAY);
}
delay(200.0e-6 + WFG_START_DELAY + WFG_STOP_DELAY - 6.0*GRADIENT_DELAY);
sim3pulse(pw, 0.0, pwN, zero, zero, t5, 0.0, 2.0e-6);
dec2phase(zero);
zgradpulse(gzlvl8, gt8);
delay(taunh - gt8);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(taunh - gt8 - 400.0e-6);
zgradpulse(gzlvl8, gt8);
txphase(one);
dec2phase(one);
delay(400.0e-6);
sim3pulse(pw, 0.0, pwN, one, zero, one, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
zgradpulse(gzlvl9, gt9);
delay(taunh - gt9);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl9, gt9);
delay(taunh - gt9);
rgpulse(pw, zero, 0.0, 0.0);
delay((gt1/10.0) + gstab - 0.5*pw + 6.0*GRADIENT_DELAY + POWER_DELAY);
rgpulse(2.0*pw, zero, 0.0, 0.0);
if(mag_flg[A] == 'y')
{
magradpulse(icosel*gzcal*gzlvl2, gt1/10.0);
}
else
{
zgradpulse(icosel*gzlvl2, gt1/10.0);
delay(4.0*GRADIENT_DELAY);
}
dec2power(dpwr2);
delay(gstab);
status(C);
setreceiver(t11);
}
示例10: pulsesequence
//.........这里部分代码省略.........
if (ix < (int)(2.0*ni)) lk_hold(); /*force lock sampling at end of experiment*/
obspower(shlvl1);
shaped_pulse(shname1,shpw1,zero,2.0e-4,2.0e-6);
zgradpulse(gzlvl1, gt1);
delay(1.0e-4);
if (ni == 0)
{
delay(tauch-gt1-1.0e-4-WFG_START_DELAY+pwC*4.0-adjust);
obspower(shlvl2);
shaped_pulse(shname2,shpw2,zero,2.0e-6,2.0e-6);
obspower(shlvl1);
decrgpulse(pwC,t1,0.0,0.0);
decrgpulse(pwC*2.0,zero,0.0,0.0);
decrgpulse(pwC,zero,0.0,0.0);
}
else
{
delay(tauch-gt1-1.0e-4-adjust);
obspower(shlvl2);
status(B);
if(CT_flg[0] == 'y')
{
/*************************************************/
/**** CT EVOLUTION **********/
/*************************************************/
decrgpulse(pwC,t1,0.0,0.0);
delay((CT_delay-tau1)*0.25-pwC*2.0/3.0);
decpower(shlvl3);
decshaped_pulse(shname3,shpw3,zero,0.0,0.0);
delay((CT_delay+tau1)*0.25-shpw2*0.5);
shapedpulse(shname2,shpw2,zero,0.0,0.0);
delay((CT_delay+tau1)*0.25-shpw2*0.5);
decshaped_pulse(shname3,shpw3,zero,0.0,0.0);
decpower(pwClvl);
delay((CT_delay-tau1)*0.25-pwC*2.0/3.0);
decrgpulse(pwC,zero,0.0,0.0);
}
else
{
/*************************************************/
/**** REAL-TIME EVOLUTION **********/
/*************************************************/
if ((tau1) > shpw2)
{
decrgpulse(pwC,t1,0.0,0.0);
if(Cdecflg[0] == 'y')
{
decpower(getval("Cdecpwr"));
decprgon(Cdecseq, 1.0/getval("Cdecdmf"), getval("Cdecres"));
decon();
}
delay((tau1-shpw2)*0.5);
xmtrphase(zero);
xmtron();
obsunblank();
obsprgon(shname2,1/shdmf2,9.0);
delay(shpw2);
obsprgoff();
obsblank();
示例11: pulsesequence
//.........这里部分代码省略.........
delay(TC - taub
- gt4 - 102e-6
- PRG_STOP_DELAY - POWER_DELAY - pwd1
- 4.0e-6 - POWER_DELAY - 4.0e-6 - WFG_START_DELAY);
/* 2D decoupling off */
dec3off();
dec3prgoff();
dec3blank();
dec3phase(three);
dec3power(dpwr3);
dec3rgpulse(pwd1,three,4.0e-6,0.0);
/* 2D decoupling off */
ifzero(v1);
delay(2.0e-6);
zgradpulse(gzlvl4,gt4);
delay(gstab);
elsenz(v1);
delay(2.0e-6);
zgradpulse(-1.0*gzlvl4,gt4);
delay(gstab);
endif(v1);
initval(1.0,v3);
decstepsize(353.0);
dcplrphase(v3);
decpower(d_ca180);
decshaped_pulse(shca180,pwca180,zero,4.0e-6,0.0);
dcplrphase(zero);
ifzero(v1);
delay(2.0e-6);
zgradpulse(gzlvl4,gt4);
delay(gstab);
elsenz(v1);
delay(2.0e-6);
zgradpulse(-1.0*gzlvl4,gt4);
delay(gstab);
endif(v1);
/* 2D decoupling on */
dec3phase(one);
dec3power(dpwr3);
dec3rgpulse(pwd1,one,4.0e-6,0.0);
dec3phase(zero);
dec3unblank();
dec3power(dpwr3_D); /* keep power down */
dec3prgon(ddseq,pwd,dresD);
dec3on();
/* 2D decoupling on */
delay(TC - taud - WFG_STOP_DELAY - gt4 - 102e-6
- POWER_DELAY - 4.0e-6 - pwd1 - POWER_DELAY
- PRG_START_DELAY);
/* 2D decoupling off */
示例12: pulsesequence
pulsesequence()
{
double tau, rg1 = 2.0e-6,
jXH = getval("jXH"),
mult = getval("mult"),
gt0 = getval("gt0"),
gzlvl0 = getval("gzlvl0"),
gt1 = getval("gt1"),
gzlvl1 = getval("gzlvl1"),
gt2 = getval("gt2"),
gzlvl2 = getval("gzlvl2"),
gt3 = getval("gt3"),
gzlvl3 = getval("gzlvl3"),
Gratio = getval("Gratio"),
gstab = getval("gstab"),
compH = getval("compH"),
compX = getval("compX"),
Hpwr = getval("Hpwr"),
jevol = 0.25 / 140.0; /* default for C-13 */
char pshape[MAXSTR], zz[MAXSTR], ad180[MAXSTR], sspul[MAXSTR];
shape hdx; /* HADAMARD stuff */
void setphases();
getstr("pshape", pshape); /* pulse shape as in wavelib */
getstr("sspul", sspul);
getstr("zz", zz);
getstr("ad180", ad180);
if (jXH > 0.0) jevol = 0.25 / jXH;
tau = gt3 + gstab;
if (mult > 0.5)
tau = 2*jevol;
setphases();
/* MAKE PBOX SHAPES */
if (FIRST_FID)
{
ad180sh = pbox_ad180("ad180", pwx, pwxlvl); /* make adiabatic 180 */
if (strcmp(pwpat,"Hdec") == 0)
{
/* make H decoupling shape */
Hdec = pboxHT_dec(pwpat, pshape, jXH, 20.0, pw, compH, tpwr);
Hpwr = Hdec.pwr;
}
}
hdx = pboxHT_F1i(pshape, pwx*compX, pwxlvl); /* HADAMARD stuff */
if (getval("htcal1") > 0.5) /* enable manual power calibration */
hdx.pwr += getval("htpwr1");
/* THE PULSE PROGRAMM STARTS HERE */
status(A);
delay(5.0e-5);
zgradpulse(gzlvl0,gt0);
if (sspul[0] == 'y')
{
rgpulse(pw,zero,rof1,rof1);
zgradpulse(gzlvl0,gt0);
}
pre_sat();
if (getflag("wet"))
wet4(zero,one);
decoffset(dof);
decpower(pwxlvl);
obspower(tpwr);
obsoffset(tof);
delay(2.0e-5);
status(B);
rgpulse(pw, v2, rg1, rg1);
delay(jevol - ad180sh.pw/2.0 - 3.0*rg1); /* +rg on both sides */
decshaped_pulse(ad180sh.name, ad180sh.pw, zero, rg1, rg1);
rgpulse(2.0*pw, zero, rg1, rg1);
delay(jevol + ad180sh.pw/2.0 - rg1);
rgpulse(pw, one, rg1, rg1);
zgradpulse(-gzlvl2, gt2);
delay(gstab);
decpower(hdx.pwr);
if (strcmp(pwpat,"") == 0)
decshaped_pulse(hdx.name, hdx.pw, zero, rg1, rg1);
else
{
obspower(Hpwr);
simshaped_pulse(pwpat, hdx.name, hdx.pw, hdx.pw, zero, zero, rg1, rg1);
}
zgradpulse(gzlvl2*1.4, gt2);
decpower(pwxlvl); obspower(tpwr);
if((zz[A] == 'y') && (mult == 0))
{
//.........这里部分代码省略.........
示例13: pulsesequence
//.........这里部分代码省略.........
zgradpulse(gzlvl5*1.2,gt5);
dec2phase(t4);
delay(gstab);
sim3pulse(0.0,0.0,2.0*pwN,zero,zero,t4,rof1,rof1); /* 180 15N */
zgradpulse(gzlvl5*1.2,gt5);
delay(gstab);
delay(0.5*taunco - gt5 - gstab - 0.5*pw180onco); /* 1/8J(NCO) */
sim3pulse(0.0*pw,pw180onco,0.0,zero,zero,zero,rof1,rof1); /* 180 CO */
delay(0.5*taunco - gt5 - gstab - 0.5*pw180onco - POWER_DELAY); /* 1/8J(NCO) */
decpwrf(rf180offca);
zgradpulse(gzlvl5*1.2,gt5);
delay(gstab);
decphase(t5);
dec2phase(zero);
dec2rgpulse(pwN,zero,0.0,0.0); /* purge for 15N */
}
zgradpulse(gzlvl3*1.3,gt3);
delay(gstab);
dec2rgpulse(pwN,t1,0.0,0.0); /* read for 15N */
delay(tau1);
decshaped_pulse("offC3",pw180offca,t5,0.0,0.0);
delay((1-kappa)*tau1);
zgradpulse(gzlvl1,gt1);
delay(gstab);
dec2phase(zero);
dec2rgpulse(2.0*pwN,zero,0.0,0.0);
delay(gt1 + gstab + (4.0/PI)*pwN + WFG_START_DELAY + pw180offca - kappa*tau1);
sim3pulse(pw,0.0,0.0,t2,zero,zero,rof1,rof1); /* Pulse for 1H */
zgradpulse(gzlvl5,gt5);
delay(gstab);
delay(tauhn - gt5 - gstab); /* delay=1/4J (NH) */
sim3pulse(2.0*pw,0.0,2.0*pwN,zero,zero,zero,rof1,rof1);
txphase(zero);
dec2phase(zero);
delay(tauhn - gt5 - gstab); /* 1/4J (NH) */
zgradpulse(gzlvl5,gt5);
delay(gstab);
txphase(one);
sim3pulse(pw,0.0,pwN,one,zero,zero,rof1,rof1); /* 90 for 1H and 15N */
示例14: pulsesequence
//.........这里部分代码省略.........
/*** purge */
dec2power(pwNlvl);
dec2pwrf(4095.0);
txphase(zero);
decphase(t2);
delay(gstab);
zgradpulse(gzlvl2, gt2);
delay(3.0*gstab);
/* evolution on t1 */
/* real time here */
if(CT[A]=='n')
{
decrgpulse(pwC, t2, 0.0, 0.0);
decphase(zero);
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay(tau1);
obsunblank();
rgpulse(pw,one,rof1,0.0);
rgpulse(2.0*pw,zero,0.0,0.0);
rgpulse(pw,one,0.0,rof1);
obsblank();
if(refocN15[A]=='y') dec2rgpulse(2.0*pwN,zero,0.0,0.0); /*n15 refocusing */
if(refocCO[A]=='y')
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
delay(tau1);
/*************** CODING with CLUB sandwitches and BIPS */
zgradpulse(gzlvl3*icosel, gt3);
delay(gstab1);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
zgradpulse(-gzlvl3*icosel, gt3);
delay(gstab1 +rof1*2.0+pw*4.0 +pwC*4.0/M_PI +2.0*POWER_DELAY+2.0*PWRF_DELAY);
if(refocN15[A]=='y') delay(2.0*pwN); /*n15 refocusing */
if(refocCO[A]=='y') /* ghost CO pulse */
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
zgradpulse(-gzlvl3*icosel, gt3);
delay(gstab1);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
zgradpulse( gzlvl3*icosel, gt3);
decpower(pwClvl);
decpwrf(4095.0);
示例15: pulsesequence
//.........这里部分代码省略.........
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
if (getflag("slpsat"))
{
shaped_satpulse("relaxD",satdly,zero);
if (getflag("prgflg"))
shaped_purge(v1,zero,v18,v19);
}
else
{
satpulse(satdly,zero,rof1,rof1);
if (getflag("prgflg"))
purge(v1,zero,v18,v19);
}
}
else
delay(d1);
if (getflag("wet"))
wet4(zero,one);
decpower(pwxlvl);
status(B);
if (getflag("nullflg"))
{
rgpulse(0.5 * pw, zero, rof1, rof1);
delay(2 * bird);
decpower(pwxlvl180);
decshaped_pulse(pwx180ad, pwx180, zero, rof1, rof1);
rgpulse(2.0 * pw, zero, rof1, rof1);
delay(2 * bird + 2 * POWER_DELAY);
decshaped_pulse(pwx180adR, pwx180, zero, rof1, rof1);
decpower(pwxlvl);
rgpulse(1.5 * pw, zero, rof1, rof1);
zgradpulse(hsglvl, hsgt);
delay(1e-3);
}
if (getflag("cpmgflg"))
{
rgpulse(pw, v1, rof1, 0.0);
cpmg(v1, v15);
}
else
rgpulse(pw, v1, rof1, rof1);
delay(tau/2 - gtA - gstab);
delay(jdly/4);
zgradpulse(gzlvlA,gtA);
delay(gstab);
obspower(selpwrA);
decpower(pwxlvl180);
decshaped_pulse(pwx180ad, pwx180, zero, rof1, rof1);
shaped_pulse(selshapeA,selpwA,zero,rof1,rof1);
zgradpulse(gzlvlA,gtA);
delay(gstab);
delay(tau/2 - gtA - gstab);
delay(jdly/2);
delay(tau/2 - gtB - gstab);