本文整理汇总了C++中setreceiver函数的典型用法代码示例。如果您正苦于以下问题:C++ setreceiver函数的具体用法?C++ setreceiver怎么用?C++ setreceiver使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了setreceiver函数的15个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: pulsesequence
void pulsesequence()
{
/* DECLARE & READ IN NEW PARAMETERS */
char compshape[MAXSTR];
getstr("compshape",compshape); /* Composit pulse shape */
loadtable("lc1d"); /* Phase table */
/* PULSE SEQUENCE */
status(A);
hsdelay(d1);
status(B);
if (getflag("wet")) wet4(t1,t2);
status(C);
if (getflag("composit"))
{
if (rfwg[OBSch-1] == 'y')
shaped_pulse(compshape,4.0*pw+0.8e-6,t3,rof1,rof2);
else
composite_pulse(pw,t3,rof1,rof2,v1);
}
else
rgpulse(pw,t3,rof1,rof2);
setreceiver(t4);
}
示例2: pulsesequence
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
MPSEQ dec = getblew("blewH",0,0.0,0.0,0,1);
strncpy(dec.ch,"dec",3);
putCmd("chHblew='dec'\n");
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + getval("rd") + getval("ad") + at;
d.dutyoff = d1 + 4.0e-6;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(phRec,4,table4);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
_mpseqon(dec, phHhx);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_mpseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
示例3: pulsesequence
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(phRec,4,table2);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
示例4: pulsesequence
pulsesequence()
{
settable(t1,4,phasecycle);
dec3power(dpwr3);
diplexer_override(0);
delay(d1);
dec3rgpulse(pw, t1, rof1, rof2);
setreceiver(t1);
}
示例5: pulsesequence
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
initval(getval("periods"),v2);
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Set Phase Tables
settable(phX90,4,table1);
settable(phRec,4,table2);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
xgate(1.0);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Apply a Rotorsync Delay
rgpulse(getval("pwX90"),phX90,0.0,0.0);
rotorsync(v2);
rgpulse(getval("pwX90"),phX90,0.0,0.0);
xgate(getval("xperiods"));
rgpulse(getval("pwX90"),phX90,0.0,0.0);
delay(10.0e-6);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
obsunblank(); decunblank(); _unblank34();
}
示例6: pulsesequence
void pulsesequence()
{
int rxgate;
double pp,
pplvl;
pp = getval("pp");
pplvl = getval("pplvl");
rxgate = (rof1 == 0.0);
if (rxgate)
rof1 = 1.0e-6; /* phase switching time */
if (newdecamp)
{
if (rxgate)
rof1 = 40.0e-6;
}
status(A);
hsdelay(d1);
status(B);
settable(t1,4,phasecycle);
pulse(pw, t1);
delay(d2);
if (newdecamp)
{
pplvl = getval("pplvl");
decpower(pplvl); /* sets DEC atten = pplvl */
}
else
{
declvlon(); /* sets dhp = 255 level */
}
simpulse(p1, pp, t1, t1, rof1, rof1);
if (newdecamp)
{
decpower(dpwr); /* sets DEC atten = dpwr */
}
else
{
declvloff();
}
delay(d2);
status(C);
setreceiver(t1);
}
示例7: pulsesequence
void pulsesequence(){
//Define Variables and Get Parameter Values
double pwTune = getval("pwTune");
pwTune = pwTune*6.0;
at = pwTune*2.0;
char atval[MAXSTR];
sprintf(atval,"at = %f\n", at);
putCmd(atval);
int chTune = (int) getval("chTune");
if ((chTune < 1) || (chTune > 4)) {
abort_message("chTune(%d) must be between 1 and 4\n", chTune);
}
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Set Phase Tables
settable(phTune,4,table1);
settable(phRec,4,table2);
setreceiver(t2);
//Begin Sequence
obspwrf(getval("aTune"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2e-6); sp1off(); delay(2.0e-6);
//Begin Phase Detected Pulse
set4Tune(chTune,getval("gain"));
delay(1.0e-4);
ShapedXmtNAcquire("phtran",pwTune,phTune,6.0e-6,chTune);
obsunblank(); decunblank(); _unblank34();
}
示例8: pulsesequence
void pulsesequence()
{
/* equilibrium period */
status(A);
hsdelay(d1);
/* tau delay */
status(B);
if (newdecamp)
{
decpower(tpwr);
decrgpulse(p1, zero, rof1, rof2);
decpower(dpwr);
}
else
{
declvlon();
decrgpulse(p1, zero, rof1, rof2);
declvloff();
}
hsdelay(d2);
/* observe period */
status(C);
settable(t1,4,phasecycle);
if (newdecamp)
{
decpower(tpwr);
decrgpulse(pw, t1, rof1, rof2);
decpower(dpwr);
}
else
{
declvlon();
decrgpulse(pw, t1, rof1, rof2);
declvloff();
}
setreceiver(t1);
}
示例9: pulsesequence
//.........这里部分代码省略.........
tsadd(t10,2,4);
}
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
if (satmode[A] == 'y')
{
if(satmove)
obsoffset(satfrq);
obspower(satpwr);
rgpulse(d1,zero,rof1,rof1);
if(satmove)
obsoffset(tof);
}
else
{
delay(d1);
}
rcvroff();
decphase(t13);
obspower(tpwr); /* Set transmitter power for hard 1H pulses */
decpower(pwcalvl);
status(B);
rgpulse(pw,t21,1.0e-5,0.0); /* First 1H 90 degree pulse */
txphase(t11);
if(tau1 < (2*pwca + pw)) /* first t1-value */
{
delay(tauhc - 2*pwca - 1.0e-6); /* assuming pwca > pw */
decrgpulse(pwca,t13,0.0,0.0); /* 13C composite 180 deg pulse */
simpulse(2*pw, 2*pwca,t11, t14, 1.0e-6, 0.0);
decrgpulse(pwca,t13,1.0e-6,0.0); /* 13C composite 180 deg pulse */
txphase(one);
decphase(t12);
delay(tauhc - 2*pwca - 1.0e-6); /* assuming pwca > pw */
}
else
{
delay(tauhc + tau1/2 - 2*pwca -1.0e-6); /* t1-evol. plus pol. trans. */
decrgpulse(pwca,t13,0.0,0.0); /* 13C composite 180 deg pulse */
decrgpulse(2*pwca,t14,1.0e-6,0.0); /* 13C composite 180 deg pulse */
decrgpulse(pwca,t13,0.0,0.0); /* 13C composite 180 deg pulse */
delay(tau1/2 - 2*pwca - pw); /* continued t1-evol. */
rgpulse(2*pw,t11,0.0,0.0); /* proton echo pulse */
txphase(one);
decphase(t12);
delay(tauhc - pw);
}
simpulse(pw,pwca,one,t12,1.0e-6,0.2e-6);
decphase(zero);
decpower(pwcolvl);
if(tau2/2 > POWER_DELAY)
delay(tau2/2 - POWER_DELAY); /* t2 evolution */
else
delay(tau2/2);
decshaped_pulse(coshape,pwco,zero,0.0,0.0);
decpower(pwcalvl);
delay(delta1 - pwco - POWER_DELAY - WFG_START_DELAY - WFG_STOP_DELAY - pw);
rgpulse(2*pw,zero,0.0,0.0);
if (tau2/2 > (pw + pwca))
delay(tau2/2 - pw -pwca); /* t2 evolution */
else
delay(tau2/2);
decrgpulse(2*pwca,t15,0.0,0.0);
decphase(t16);
decpower(dipsipwr);
delay(delta1 - pwca - POWER_DELAY);
if (ncyc>0.0)
{
decrgpulse(trim,t16,0.0,0.0);
starthardloop(v7);
dipsi3a(); dipsi3b(); dipsi3b(); dipsi3a();
endhardloop();
}
txphase(zero); decphase(zero);
decpower(pwcalvl);
delay(delta2 - POWER_DELAY - pwca);
simpulse(2*pw,2*pwca,zero,zero,0.0,0.0);
txphase(t13);
delay(delta2);
/* start reversed INEPT */
simpulse(pw,pwca,t13,zero,0.0,0.0);
decphase(zero);
txphase(zero);
delay(tauhc - 2*pwca - 1.0e-6); /* delay = 1/4J */
decrgpulse(pwca,zero,0.0,0.0);
simpulse(2*pw,2*pwca,zero,one,1.0e-6,0.0);
decrgpulse(pwca,zero,1.0e-6,rof2);
delay(tauhc - 2*pwca - 1.0e-6 - rof2); /* delay = 1/4J */
decrgpulse(pwca,zero,0.0,0.0); /* Filter out IySz terms */
decrgpulse(pwca,t16,1.0e-6,0.0);
decpower(dpwr); /* Set power for decoupling */
/* BEGIN ACQUISITION */
setreceiver(t10);
status(C);
}
示例10: pulsesequence
//.........这里部分代码省略.........
setstatus(DEC3ch, FALSE, 'c', FALSE, dmf3);
dec3blank();
dec3phase(three);
dec3power(pwDlvl);
dec3rgpulse(pwD,three,4.0e-6,0.0);
/* Turn off D decoupling */
decoffset(dof); /* set carrier back to Co */
decpower(d_chirp);
delay(0.2e-6);
zgradpulse(gzlvl9,gt9);
delay(150.0e-6);
/* t3 period begins */
dec2rgpulse(pwN,t2,2.0e-6,0.0);
dec2phase(t3);
delay(bigTN - tau3);
dec2rgpulse(2.0*pwN,t3,0.0,0.0);
decshaped_pulse(spchirp,pwchirp,zero,0.0,0.0);
txphase(zero);
dec2phase(t4);
delay(0.2e-6);
zgradpulse(gzlvl1,gt1);
delay(500.0e-6);
delay(bigTN - WFG_START_DELAY - pwchirp - WFG_STOP_DELAY
-gt1 -500.2e-6 -2.0*GRADIENT_DELAY);
delay(tau3);
sim3pulse(pw,0.0e-6,pwN,zero,zero,t4,0.0,0.0);
/* t3 period ends */
decpower(d_c90);
decrgpulse(pwc90,zero,4.0e-6,0.0);
decoffset(dof-(174-56)*dfrq);
decrgpulse(pwc90,zero,20.0e-6,0.0);
delay(0.2e-6);
zgradpulse(gzlvl6,gt6);
delay(2.0e-6);
dec2phase(zero);
delay(taub - POWER_DELAY - 4.0e-6 - pwc90 - 20.0e-6 - pwc90 - gt6 - 2.2e-6);
sim3pulse(2*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
decoffset(dof);
delay(0.2e-6);
zgradpulse(gzlvl6,gt6);
delay(200.0e-6);
txphase(one);
dec2phase(one);
delay(taub - gt6 - 200.2e-6);
sim3pulse(pw,0.0e-6,pwN,one,zero,one,0.0,0.0);
delay(0.2e-6);
zgradpulse(gzlvl7,gt7);
delay(2.0e-6);
txphase(zero);
dec2phase(zero);
delay(taub - gt7 - 2.2e-6);
sim3pulse(2*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
delay(0.2e-6);
zgradpulse(gzlvl7,gt7);
delay(200.0e-6);
delay(taub - gt7 - 200.2e-6);
sim3pulse(pw,0.0e-6,pwN,zero,zero,zero,0.0,0.0);
delay(gt2 +gstab -0.5*(pwN -pw) -2.0*pw/PI);
rgpulse(2*pw,zero,0.0,0.0);
delay(2.0e-6);
zgradpulse(icosel*gzlvl2,gt2);
decpower(dpwr);
dec2power(dpwr2);
delay(gstab -2.0e-6 -2.0*GRADIENT_DELAY -2.0*POWER_DELAY);
lk_sample();
/* BEGIN ACQUISITION */
status(C);
setreceiver(t6);
}
示例11: pulsesequence
//.........这里部分代码省略.........
dec2phase(t4);
delay(gstab*2.0);
/* t1 (C) and t1+t1'(N) evolution */
dec2rgpulse(pwN, t4, 0.0, 0.0);
delay(gt4+gstab + gt4+gstab + pwC*3.0
+2.0*(pwHs +2.0*rof1));
if(decCACO[A]=='y'){ delay(2.0*caco180pw);}
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(tau1p); /* t1 */
decrgpulse(pwC,t5,0.0,0.0);
delay(tau1*0.5);
if(decCACO[A]=='y')
{
decpower(caco180pwr);
decshaped_pulse(caco180shape,caco180pw,zero, 0.0, 0.0);
decpower(pwClvl);
}
obspower(ch180pwr); /*180 on methyls*/
shaped_pulse(ch180shape,ch180pw,zero,rof1,rof1);
obspower(tpwr);
delay(tau1*0.5);
zgradpulse(gzlvl4, gt4); /*coding */
delay(gstab + pwHs -ch180pw -2.0*GRADIENT_DELAY -2.0*POWER_DELAY -WFG_START_DELAY- WFG_STOP_DELAY);
decrgpulse(2.0*pwC,zero,0.0,0.0);
if(decCACO[A]=='y')
{
decpower(caco180pwr);
decshaped_pulse(caco180shape,caco180pw,zero, 0.0, 0.0);
decpower(pwClvl);
}
/* delay(ch180pw+2.0*rof1);*/
zgradpulse(gzlvl5, gt4);
delay(gstab - rof1 -2.0*GRADIENT_DELAY -2.0*POWER_DELAY -WFG_START_DELAY- WFG_STOP_DELAY);
/*Water flipback (flipdown actually ) */
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc",pwHs,three,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
/* reverse double INEPT */
sim3pulse(pw, pwC, 0.0, t21, t23, zero, rof1, rof1);
/* rgpulse(pw, t21, rof1, rof1); */
zgradpulse(gzlvl11, gt1);
delay(gstab);
delay(tauCH1 -gt1 -gstab -2.0*pwC
+ (-2.0/M_PI*pwC-0.5*(pwN-pwC) +pwN));
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -tauCH1 - 0.65*(pw + pwN)-rof1 -(pwC-pw)
-(-2.0/M_PI*pwC-0.5*(pwN-pwC) +pwN) );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl11, gt1);
delay(gstab);
delay(tauCH1-gt1 -gstab-2.0*pwC);
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -1.3*pwN -tauCH1);
sim3pulse(pw, pwC, pwN, one, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(tauCH2-2.0*pwC-gt1-gstab);
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -1.3*pwN-tauCH2);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(tauNH-1.6*pwN -POWER_DELAY -ch90pw*ch90corr -gt1-gstab +pwN*0.5 -PRG_START_DELAY);
/* delay(ch90pw*0.5);
sim3pulse(0.0,0.0, pwN, one, zero, t22, 0.0, 0.0);
delay(ch90pw*0.5);*/
/* sim3shaped_pulse(ch90shape,"hard","hard",ch90pw,0.0, pwN, zero,zero, t22,0.0,0.0);*/
/* ch90corr is a fraction of ch90pw to correct for a 1H phase rollcaused by shaped 90 on
CH3 protons for a sinc pulse ch90corr=0.41 seems to be good. */
obspower(ch90pwr);
txphase(one);
obsunblank(); xmtron();
obsprgon(ch90shape,1.0/ch90dmf,ch90dres); /*PRG_START_DELAY */
delay(ch90pw*ch90corr-pwN*0.5);
dec2rgpulse(pwN, t22, 0.0, 0.0);
delay(ch90pw*(1.0-ch90corr)-pwN*0.5);
obsprgoff(); xmtroff(); obsblank(); /*PRG_STOP_DELAY */
obspower(tpwr); /*POWER_DEALY */
delay( gstab +gt6 +2.0*GRADIENT_DELAY
+2.0*POWER_DELAY -0.65*pw -POWER_DELAY
+pwN*0.5 -ch90pw*(1.0-ch90corr)
-PRG_STOP_DELAY);
rgpulse(2.0*pw, zero, rof1, rof1);
dec2power(dpwr2); decpower(dpwr); /* 2.0*POWER_DELAY */
zgradpulse(gzlvl6*icosel, gt6); /* 2.0*GRADIENT_DELAY */
delay(gstab);
status(C);
setreceiver(t31);
}
示例12: pulsesequence
//.........这里部分代码省略.........
decpwrf(dfCa90);
decphase(t3); dcplrphase(v2);
delay(bigTC -2.0*SAPS_DELAY -PWRF_DELAY);
/* c13 offset on CO, slp 90 at Ca */
decshaped_pulse("offC13", pwCa90, t3, 0.0, 0.0);
xmtroff();
decpwrf(dfCO180);
decphase(zero); dcplrphase(zero);
dec2phase(t2);
delay(2.0e-5);
zgradpulse(gzlvl4,gt4);
delay(2.0e-6);
obsprgon("waltz16", pwHd, 180.0);
xmtron();
txphase(zero);
delay(150.0e-6);
dec2rgpulse(pwN, t2, 0.0, 0.0);
/* N15 EVOLUTION BEGINS HERE */
delay(bigTN/2.0 -tau2);
decshaped_pulse("offC8", pwCO180, zero, 0.0, 0.0); /* c13 offset on CO, on-res 180 on CO */
decpwrf(dfCa180);
dec2phase(t4);
delay(bigTN/2.0 -tau2);
dec2rgpulse(2.0*pwN, t4, 0.0, 0.0);
decshaped_pulse("offC17", pwCa180, zero, 0.0, 0.0); /* c13 offset on CO, slp 180 at Ca */
decpwrf(dfCO180);
delay(bigTN/2.0 +tau2 -pwCa180 -WFG_START_DELAY -WFG_STOP_DELAY);
decshaped_pulse("offC8", pwCO180, zero, 0.0, 0.0);
/* c13 offset on CO, on-res 180 on CO */
delay(bigTN/2.0 +tau2 -tauc -PRG_STOP_DELAY);
dec2phase(t5);
xmtroff();
obsprgoff();
obspower(tpwr);
if (mag_flg[A] == 'y') magradpulse(gzcal*gzlvl1, gt1);
else zgradpulse(gzlvl1, gt1);
delay(tauc -gt1 -2.0*GRADIENT_DELAY);
/* N15 EVOLUTION ENDS HERE */
sim3pulse(pw,0.0, pwN, zero,zero, t5, 0.0, 0.0);
dec2phase(zero);
delay(2.0e-6);
zgradpulse(0.8*gzlvl5, gt5);
delay(taud - gt5 - 2.0e-6);
sim3pulse(2.0*pw,0.0, 2.0*pwN, zero,zero, zero, 0.0, 0.0);
delay(taud - gt5 - 500.0e-6);
zgradpulse(0.8*gzlvl5, gt5);
txphase(one);
decphase(one);
delay(500.0e-6);
sim3pulse(pw,0.0, pwN, one,zero, one, 0.0, 0.0);
delay(2.0e-6);
txphase(zero);
decphase(zero);
zgradpulse(gzlvl5, gt5);
delay(taud - gt5 - 2.0e-6);
sim3pulse(2.0*pw,0.0, 2.0*pwN, zero,zero, zero, 0.0, 0.0);
delay(taud - gt5 - 2.0*POWER_DELAY - 500.0e-6);
zgradpulse(gzlvl5, gt5);
decpower(dpwr);
dec2power(dpwr2);
delay(500.0e-6);
rgpulse(pw, zero, 0.0, 0.0);
delay(gstab +gt2 +2.0*GRADIENT_DELAY);
rgpulse(2.0*pw, zero, 0.0, 0.0);
if (mag_flg[A] == 'y') magradpulse(icosel*gzcal*gzlvl2, gt2);
else zgradpulse(icosel*gzlvl2, gt2);
delay(0.5*gstab);
rcvron();
statusdelay(C, 0.5*gstab);
setreceiver(t16);
}
示例13: pulsesequence
//.........这里部分代码省略.........
{
dutycycle = dutycycle + d2 + p180;
timeoff = timeoff + d2;
}
if (toss[0] == 'y')
{
dutycycle = dutycycle + 2.142/srate;
if (dm[2] != 'y') timeoff = timeoff + 2.142/srate - 8*pwx;
}
if (pdp[0] == 'y')
{
dutycycle = dutycycle + 2.0/srate;
timeoff = timeoff + pdpd2;
}
dutycycle = timeoff/dutycycle;
if ((dutycycle < 0.8) || (dm[0] == 'y'))
{
fprintf(stdout, "Duty cycle is %5.2f%%.\n", (1.0 - dutycycle) * 100.0);
fprintf(stdout, "ABORT! The duty cycle must be less than 20%%.\n");
psg_abort(1);
}
/*begin pulse sequence*/
if (xpol[0] == 'n')
{
settable(t1,4,table1);
settable(t3,4,table3);
settable(t4,4,table4);
status(A);
setreceiver(t1);
obs_pw_ovr(TRUE);
dec_pw_ovr(TRUE);
declvloff();
decpwrf(dipolr);
obspwrf(tpwrm);
delay(d1);
if(p180 > 0.0)
{
decoff();
if (dm[1] == 'y') decon();
rgpulse(p180, zero, rof1, 0.0);
delay(d2);
}
rcvroff();
rgpulse(pw, t1,rof1, 0.0);
decoff();
if (dm[2] == 'y') decon();
}
else
{
settable(t1,4,table1);
settable(t2,4,table2);
settable(t3,4,table3);
settable(t4,4,table4);
status(A);
setreceiver(t1);
obs_pw_ovr(TRUE);
dec_pw_ovr(TRUE);
示例14: pulsesequence
//.........这里部分代码省略.........
decpwrf(rf180onco);
delay(tau1);
if (lambda>0.0)
{
delay(lambda*tau1);
sim3pulse(0.0,pw180onco,2.0*pwN,zero,t5,zero,0.0,0.0);
delay(lambda*tau1);
}
else decrgpulse(pw180onco,t5,0.0,0.0);
decpwrf(rf180offca); /* Set decoupler power to rf180offca */
simshaped_pulse("","offC3",0.0,pw180offca,zero,zero,0.0,0.0); /* 180 CA */
decphase(zero);
decpwrf(rf90onco);
/* start reverse transfer from CO to N by INEPT */
decrgpulse(pw90onco,zero,0.0,0.0); /* 90 for 13C' */
zgradpulse(gzlvl3*0.7,gt3);
delay(gstab);
dec2rgpulse(pwN,t1,0.0,0.0); /* 90 for 15N */
delay((taunco - tauhn) - kappa*tau2); /* 1/4J(NCO) - 1/4J(NH) - kt2/2 */
dec2rgpulse(2.0*pwN,zero,0.0,0.0); /* 180 for 15N */
delay((1-kappa)*tau2); /* (1-k)t2/2 */
decpwrf(rf180onco);
decrgpulse(pw180onco,zero,0.0,0.0); /* 180 for 13C' */
decpwrf(rf180offca);
delay(taunco - tauhn - gt1 - gstab - pw180onco - pw180offca - 3.0*POWER_DELAY);
zgradpulse(gzlvl1,gt1);
delay(gstab);
decshaped_pulse("offC3",pw180offca,zero,0.0,0.0);
delay(tau2); /* t2/2 */
/* start TROSY transfer from N to HN */
rgpulse(pw,t2,rof1,0.0); /* 180 for 1H */
zgradpulse(gzlvl5,gt5);
delay(gstab);
delay(tauhn - gt5 - gstab );
decpwrf(rf180onco);
sim3pulse(2.0*pw,pw180onco,2.0*pwN,zero,zero,zero,rof1,rof1);
delay(tauhn - gt5 - gstab );
zgradpulse(gzlvl5,gt5);
delay(gstab);
decpwrf(rf90onco);
sim3pulse(pw,pw90onco,pwN,one,t6,zero,rof1,0.0);
/* shaped pulse for water flip-back */
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc_u",pwHs,t2,rof1,0.0);
obspwrf(4095.0);
obspower(tpwr);
/* shaped pulse */
zgradpulse(gzlvl5*0.9,gt5);
delay(gstab);
delay(tauhn - gt5 - gstab - POWER_DELAY - pwHs);
sim3pulse(2.0*pw,0.0,2.0*pwN,zero,zero,zero,rof1,rof1);
delay(tauhn - gt5 - gstab );
decpower(dpwr);
zgradpulse(gzlvl5*0.9,gt5);
delay(gstab);
dec2rgpulse(pwN,t3,0.0,0.0); /* 90 for 15N */
dec2power(dpwr2);
delay((gt1/10.0) - pwN + gstab -POWER_DELAY);
rgpulse(2.0*pw, zero, rof1, rof1);
zgradpulse(gzlvl2*icosel,gt1/10.0);
delay(gstab);
/* acquire data */
status(C);
setreceiver(t7);
}
示例15: pulsesequence
//.........这里部分代码省略.........
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);
else
zgradpulse(gzlvl1, gt1);
rcvron();
if ((STUD[A]=='n') && (dm[C] == 'y'))
decpower(dpwr);
if ( dm3[B] == 'y' ) /* turns off 2H decoupling */
{
delay(0.5*del-40.0e-6 -gt1 -1/dmf3);
setstatus(DEC3ch, FALSE, 'c', FALSE, dmf3);
dec3rgpulse(1/dmf3,three,2.0e-6,2.0e-6);
dec3blank();
lk_autotrig(); /* resumes lock pulsing */
lk_sample();
if (mag_flg[A] == 'y')
statusdelay(C,40.0e-6 - 2.0*VAGRADIENT_DELAY - POWER_DELAY);
else
statusdelay(C,40.0e-6 - 2.0*GRADIENT_DELAY - POWER_DELAY);
}
else
{
delay(0.5*del-40.0e-6 -gt1);
if (mag_flg[A] == 'y')
statusdelay(C,40.0e-6 - 2.0*VAGRADIENT_DELAY - POWER_DELAY);
else
statusdelay(C,40.0e-6 - 2.0*GRADIENT_DELAY - POWER_DELAY);
}
if ((STUD[A]=='y') && (dm[C] == 'y'))
{decpower(studlvl);
decunblank();
decon();
decprgon(stCdec,1/stdmf, 1.0);
startacq(alfa);
acquire(np, 1.0/sw);
decprgoff();
decoff();
decblank();
}
setreceiver(t11);
}