pulsesequence()
{
char
shname1[MAXSTR],
f1180[MAXSTR],
SE_flg[MAXSTR];
int icosel = 0,
t1_counter,
phase;
double d2_init=0.0,
pwS4,pwS5,pwS6,pwS7,
lambda = getval("lambda"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl5 = getval("gzlvl5"),
gt1 = getval("gt1"),
gt3 = getval("gt3"),
gt5 = getval("gt5"),
gstab = getval("gstab"),
shpw1,shlvl1=getval("shlvl1"),
pwClvl = getval("pwClvl"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
ni = getval("ni"),
d2 = getval("d2"),
tau1 = getval("tau1");
getstr("shname1", shname1);
getstr("SE_flg",SE_flg);
getstr("f1180",f1180);
phase = (int) (getval("phase") + 0.5);
settable(t3,2,phi3);
settable(t10,1,phi10);
settable(t12,2,phi12);
settable(t13,2,phi13);
/* INITIALIZE VARIABLES */
shpw1 = getval("shpw1");
pwS4 = h_shapedpw("eburp2",4.0,3.5,zero, 0.0, 0.0);
pwS5 = h_shapedpw("pc9f",4.0,3.5,zero, 2.0e-6, 0.0);
pwS6 = h_shapedpw("reburp",4.0,3.5,zero, 0.0, 0.0);
pwS7 = c_shapedpw2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0);
if (SE_flg[0] =='y')
{
if (phase == 2) {tsadd(t10,2,4); icosel = +1;}
else icosel = -1;
}
else
{
if (phase == 2) {tsadd(t3,1,4); icosel=1;}
}
tau1 = d2;
if((f1180[A] == 'y') )
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1;
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); tsadd(t13,2,4); }
status(A);
decpower(pwClvl);
dec2power(pwNlvl);
set_c13offset("co");
zgradpulse(gzlvl3, gt3);
delay(d1-gt3);
lk_hold();
rcvroff();
h_shapedpulse("pc9f",4.0,3.5,zero, 2.0e-6, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_shapedpulse("pc9f_",4.0,3.5,one, 0.0, 0.0);
obspower(shlvl1);
/**************************************************************************/
/** Sensitivity enhanced version **********************/
/**************************************************************************/
if (SE_flg[0] == 'y')
{
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
delay(tau1*0.5);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
c_shapedpulse2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0);
delay(tau1*0.5);
zgradpulse(-gzlvl1, gt1);
delay(pwS4-shpw1-pwS7-gt1);
dec2rgpulse(pwN*2.0,zero,0.0,0.0);
h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0);
dec2rgpulse(pwN, t10, 0.0, 0.0);
}
/**************************************************************************/
/** standard INEPT-based version **********************/
/**************************************************************************/
else
{
if (ni < 1.0)
{
dec2rgpulse(pwN,t3,0.0,0.0);
dec2rgpulse(pwN*2.0, zero, 0.0, 0.0);
dec2rgpulse(pwN, zero, 0.0, 0.0);
}
else
{
if (tau1 < shpw1)
{
shaped_pulse(shname1,shpw1,two,0.0,0.0);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
delay(tau1);
dec2rgpulse(pwN*3.0, zero, 0.0, 0.0);
}
else
{
if (tau1*0.5 < (pwS7+shpw1*0.5))
{
shaped_pulse(shname1,shpw1,two,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
delay(tau1*0.5-shpw1*0.5);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
delay(tau1*0.5-shpw1*0.5);
dec2rgpulse(pwN*3.0, zero, 0.0, 0.0);
}
else
{
shaped_pulse(shname1,shpw1,two,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
delay(tau1*0.5-shpw1*0.5);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
c_shapedpulse2("isnob5",40.0,-125.0,"isnob5",40.0,0.0 , two, 0.0, 0.0);
delay(tau1*0.5-shpw1*0.5-pwS7);
dec2rgpulse(pwN*3.0, zero, 0.0, 0.0);
}
}
}
}
/**************************************************************************/
if (SE_flg[0] == 'y')
{
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
dec2rgpulse(pwN, one, 0.0, 0.0);
h_shapedpulse("eburp2_",4.0,3.5,one, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t10);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_shapedpulse("eburp2",4.0,3.5,zero, 0.0, 0.0);
delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY);
h_shapedpulse("reburp",4.0,3.5,zero, 0.0, 0.0);
zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */
delay(gstab);
}
else
{
h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4);
}
dec2power(dpwr2); /* POWER_DELAY */
lk_sample();
if (SE_flg[0] == 'y')
setreceiver(t13);
else
setreceiver(t12);
statusdelay(C,1.0e-4 );
}
pulsesequence()
{
char
shname1[MAXSTR],
f1180[MAXSTR],
f2180[MAXSTR],
SE_flg[MAXSTR],
dec_flg[MAXSTR],
CT_flg[MAXSTR];
int icosel=1,
t1_counter,
t2_counter,
ni2 = getval("ni2"),
phase;
double d2_init=0.0,
d3_init=0.0,
pwS1,pwS2,pwS3,pwS4,pwS5,pwS6,
lambda = getval("lambda"),
CTdelay = getval("CTdelay"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5"),
gzlvl6 = getval("gzlvl6"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gt6 = getval("gt6"),
gstab = getval("gstab"),
tpwrsf = getval("tpwrsf"),
shlvl1,
shpw1,
pwClvl = getval("pwClvl"),
pwC = getval("pwC"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
dpwr2 = getval("dpwr2"),
d2 = getval("d2"),
shbw = getval("shbw"),
shofs = getval("shofs")-4.77,
timeTN = getval("timeTN"),
tauC = getval("tauC"),
tau1 = getval("tau1"),
tau2 = getval("tau2"),
taunh = getval("taunh");
getstr("shname1", shname1);
getstr("SE_flg",SE_flg);
getstr("dec_flg",dec_flg);
getstr("CT_flg",CT_flg);
getstr("f1180",f1180);
getstr("f2180",f2180);
phase = (int) (getval("phase") + 0.5);
settable(t1,2,phi1);
settable(t3,4,phi3);
settable(t10,1,phi10);
settable(t12,4,phi12);
settable(t13,4,phi13);
/* INITIALIZE VARIABLES */
shpw1 = pw*8.0;
shlvl1 = tpwr;
pwS1 = c13pulsepw("co", "ca", "sinc", 90.0);
pwS2 = c13pulsepw("co", "ca", "sinc", 180.0);
pwS3 = c13pulsepw("ca", "co", "square", 180.0);
pwS4 = h_shapedpw("eburp2",shbw,shofs,zero, 0.0, 0.0);
pwS6 = h_shapedpw("reburp",shbw,shofs,zero, 0.0, 0.0);
pwS5 = h_shapedpw("pc9f",shbw,shofs,zero, 2.0e-6, 0.0);
if (SE_flg[0] == 'y')
{
if ( ni2*1/(sw2)/2.0 > (timeTN-pwS2*0.5-pwS4))
{ printf(" ni2 is too big. Make ni2 equal to %d or less.\n",
((int)((timeTN-pwS2*0.5-pwS4)*2.0*sw2))); psg_abort(1);}
}
else
{
if ( ni2*1/(sw2)/2.0 > (timeTN-pwS2*0.5))
{ printf(" ni2 is too big. Make ni2 equal to %d or less.\n",
((int)((timeTN-pwS2*0.5)*2.0*sw2))); psg_abort(1);}
}
if (CT_flg[0] == 'y')
{
if ( ni*1/(sw1) > (CTdelay))
{ printf(" ni is too big. Make ni equal to %d or less.\n",
((int)(CTdelay*sw1))); psg_abort(1);}
}
if (phase == 1) ;
if (phase == 2) tsadd(t1,1,4);
if (SE_flg[0] =='y')
{
if (phase2 == 2) {tsadd(t10,2,4); icosel = +1;}
else icosel = -1;
}
else
{
if (phase2 == 2) tsadd(t3,1,4);
}
tau1 = d2;
if((f1180[A] == 'y') )
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1;
tau2 = d3;
if((f2180[A] == 'y') )
{ tau2 += ( 1.0 / (2.0*sw2) ); if(tau2 < 0.2e-6) tau2 = 0.0; }
tau2 = tau2;
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t12,2,4); tsadd(t13,2,4); }
if( ix == 1) d3_init = d3;
t2_counter = (int) ( (d3-d3_init)*sw2 + 0.5 );
if(t2_counter % 2)
{ tsadd(t3,2,4); tsadd(t12,2,4); tsadd(t13,2,4); }
status(A);
rcvroff();
decpower(pwClvl);
decoffset(dof);
dec2power(pwNlvl);
dec2offset(dof2);
obspwrf(tpwrsf);
decpwrf(4095.0);
obsoffset(tof);
set_c13offset("co");
zgradpulse(gzlvl2, gt2);
delay(1.0e-4);
delay(d1-gt2);
lk_hold();
h_shapedpulse("pc9f",shbw,shofs,zero, 2.0e-6, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_shapedpulse("pc9f_",shbw,shofs,one, 0.0, 0.0);
obspower(shlvl1);
/**************************************************************************/
dec2rgpulse(pwN,zero,0.0,0.0);
zgradpulse(gzlvl6,gt6);
delay(timeTN-pwS2*0.5-gt6);
sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
delay(timeTN-pwS2*0.5-taunh*0.5-shpw1);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
zgradpulse(gzlvl6, gt6);
delay(taunh*0.5-gt6);
dec2rgpulse(pwN,zero,0.0,0.0);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
/**************************************************************************/
/*** CO -> CA transfer *********************************/
/**************************************************************************/
c13pulse("co", "ca", "sinc", 90.0, zero, 2.0e-6, 0.0);
zgradpulse(-gzlvl4, gt4);
decphase(zero);
delay(tauC - gt4 - 0.5*10.933*pwC);
decrgpulse(pwC*158.0/90.0, zero, 0.0, 0.0);
decrgpulse(pwC*171.2/90.0, two, 0.0, 0.0);
decrgpulse(pwC*342.8/90.0, zero, 0.0, 0.0); /* Shaka 6 composite */
decrgpulse(pwC*145.5/90.0, two, 0.0, 0.0);
decrgpulse(pwC*81.2/90.0, zero, 0.0, 0.0);
decrgpulse(pwC*85.3/90.0, two, 0.0, 0.0);
zgradpulse(-gzlvl4, gt4);
decphase(one);
delay(tauC - gt4 - 0.5*10.933*pwC - WFG_START_DELAY - 2.0e-6);
c13pulse("co", "ca", "sinc", 90.0, one, 2.0e-6, 0.0);
/**************************************************************************/
/* xxxxxxxxxxxxxxxxxxxx 13Ca EVOLUTION xxxxxxxxxxxxxxxxxx */
/**************************************************************************/
set_c13offset("ca");
if (CT_flg[0] == 'y')
{
c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 0.0);
delay(tau1*0.5);
sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 0.0, 0.0);
delay(CTdelay*0.5);
c13pulse("cab", "co", "square", 180.0, zero, 0.0, 0.0);
delay((CTdelay-tau1)*0.5);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
c13pulse("ca", "co", "square", 90.0, zero, 0.0, 0.0);
}
else
{
if (dec_flg[0] == 'y')
{
c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 2.0e-6);
delay(tau1*0.5);
c_shapedpulse2("isnob5",30.0,-31.0,"isnob5",30.0,120.0 , zero, 0.0, 0.0);
dec2rgpulse(pwN*2.0,zero,0.0,0.0);
delay(tau1*0.5);
c13pulse("ca", "co", "square", 180.0, zero, 2.0e-6, 2.0e-6);
dec2rgpulse(pwN*2.0,two,0.0,0.0);
c_shapedpulse2("isnob5",30.0,-31.0,"isnob5",30.0,120.0 , two, 0.0, 0.0);
c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 2.0e-6);
}
else
{
c13pulse("ca", "co", "square", 90.0, t1, 2.0e-6, 2.0e-6);
delay(tau1*0.5);
sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
delay(tau1*0.5);
c13pulse("ca", "co", "square", 180.0, zero, 2.0e-6, 2.0e-6);
c13pulse("co", "ca", "sinc", 180.0, zero, 2.0e-6, 2.0e-6);
if (pwN*2.0 > pwS2) delay(pwN*2.0 - pwS2);
c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 2.0e-6);
}
}
set_c13offset("co");
zgradpulse(gzlvl3, gt3);
delay(1.0e-4);
/**************************************************************************/
/*** CO -> CA transfer *********************************/
/**************************************************************************/
c13pulse("co", "ca", "sinc", 90.0, one, 2.0e-6, 0.0);
zgradpulse(gzlvl4, gt4*0.7);
delay(tauC - gt4*0.7 - 0.5*10.933*pwC);
decrgpulse(pwC*158.0/90.0, zero, 0.0, 0.0);
decrgpulse(pwC*171.2/90.0, two, 0.0, 0.0);
decrgpulse(pwC*342.8/90.0, zero, 0.0, 0.0); /* Shaka 6 composite */
decrgpulse(pwC*145.5/90.0, two, 0.0, 0.0);
decrgpulse(pwC*81.2/90.0, zero, 0.0, 0.0);
decrgpulse(pwC*85.3/90.0, two, 0.0, 0.0);
zgradpulse(gzlvl4, gt4*0.7);
delay(tauC - gt4*0.7 - 0.5*10.933*pwC - WFG_START_DELAY - 2.0e-6);
/* WFG_START_DELAY */
c13pulse("co", "ca", "sinc", 90.0, zero, 2.0e-6, 0.0);
/**************************************************************************/
obspower(shlvl1);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
if (SE_flg[0] == 'y')
{
delay(tau2*0.5);
c13pulse("ca", "co", "square", 180.0, zero, 0.0, 0.0);
delay(taunh*0.5-pwS3);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
delay(timeTN-pwS2*0.5-shpw1-taunh*0.5-gt1-1.0e-4);
zgradpulse(-gzlvl1, gt1);
delay(1.0e-4);
sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
delay(timeTN-pwS2*0.5-tau2*0.5-pwS4);
h_shapedpulse("eburp2",shbw,shofs,zero, 2.0e-6, 0.0);
dec2rgpulse(pwN, t10, 0.0, 0.0);
}
else
{
delay(tau2*0.5);
c13pulse("ca", "co", "square", 180.0, zero, 0.0, 0.0);
delay(taunh*0.5-pwS3);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
delay(timeTN-pwS2*0.5-taunh*0.5-shpw1);
sim3_c13pulse("", "co", "ca", "sinc", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
delay(timeTN-pwS2*0.5-tau2*0.5);
dec2rgpulse(pwN, zero, 0.0, 0.0);
}
/**************************************************************************/
if (SE_flg[0] == 'y')
{
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
dec2rgpulse(pwN, one, 0.0, 0.0);
h_shapedpulse("eburp2_",shbw,shofs,one, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t10);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_shapedpulse("eburp2",shbw,shofs,zero, 0.0, 0.0);
delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY);
h_shapedpulse("reburp",shbw,shofs,zero, 0.0, rof2);
zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */
delay(gstab);
}
else
{
h_shapedpulse("eburp2",shbw,shofs,zero, 2.0e-6, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",shbw,shofs,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4);
}
dec2power(dpwr2); /* POWER_DELAY */
lk_sample();
if (SE_flg[0] == 'y')
setreceiver(t13);
else
setreceiver(t12);
rcvron();
statusdelay(C,1.0e-4 );
}
pulsesequence()
{
char
shname1[MAXSTR],
f1180[MAXSTR],
f2180[MAXSTR],
SE_flg[MAXSTR];
int icosel = 0,
t1_counter,
t2_counter,
ni2 = getval("ni2"),
phase;
double d2_init=0.0,
d3_init=0.0,
pwS1,pwS2,pwS3,pwS4,pwS5,pwS6,pwS7,pwS8,
lambda = getval("lambda"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5"),
gzlvl6 = getval("gzlvl6"),
gt1 = getval("gt1"),
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gt6 = getval("gt6"),
gstab = getval("gstab"),
shlvl1 = getval("shlvl1"),
shpw1 = getval("shpw1"),
pwClvl = getval("pwClvl"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
d2 = getval("d2"),
timeTN = getval("timeTN"),
Delta,
tauNCO = getval("tauNCO"),
tauC = getval("tauC"),
tau1 = getval("tau1"),
tau2 = getval("tau2"),
taunh = getval("taunh");
getstr("shname1", shname1);
getstr("SE_flg",SE_flg);
getstr("f1180",f1180);
getstr("f2180",f2180);
phase = (int) (getval("phase") + 0.5);
settable(t1,2,phi1);
settable(t2,1,phi2);
settable(t3,4,phi3);
settable(t4,8,phi4);
settable(t10,1,phi10);
settable(t12,8,phi12);
settable(t13,8,phi13);
/* INITIALIZE VARIABLES */
Delta = timeTN-tauNCO;
pwS1 = c13pulsepw("ca", "co", "square", 90.0);
pwS2 = c13pulsepw("ca", "co", "square", 180.0);
pwS3 = c13pulsepw("co", "ca", "sinc", 180.0);
pwS7 = c13pulsepw("co", "ca", "sinc", 90.0);
pwS8 = c_shapedpw("reburp",60.0 ,-135.0,zero, 0.0, 0.0);
pwS4 = h_shapedpw("eburp2",4.0,3.5,zero, 0.0, 0.0);
pwS6 = h_shapedpw("reburp",4.0,3.5,zero, 0.0, 0.0);
pwS5 = h_shapedpw("pc9f",4.0,3.5,zero, 2.0e-6, 0.0);
if (SE_flg[0] == 'y')
{
if ( ni2*1/(sw2)/2.0 > (timeTN-Delta-pwS3-pwS4))
{ printf(" ni2 is too big. Make ni2 equal to %d or less.\n",
((int)((timeTN-Delta-pwS3-pwS4)*2.0*sw2))); psg_abort(1);}
}
else
{
if ( ni2*1/(sw2)/2.0 > (timeTN-Delta-pwS3))
{ printf(" ni2 is too big. Make ni2 equal to %d or less.\n",
((int)((timeTN-Delta-pwS3)*2.0*sw2))); psg_abort(1);}
}
if (phase == 1) ;
if (phase == 2) {tsadd(t1,1,4);}
if (SE_flg[0] =='y')
{
if (phase2 == 2) {tsadd(t10,2,4); icosel = +1;}
else icosel = -1;
}
else
{
if (phase2 == 2) {tsadd(t3,1,4); icosel = 1;}
}
tau1 = d2;
if((f1180[A] == 'y') )
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1;
tau2 = d3;
if((f2180[A] == 'y') )
{ tau2 += ( 1.0 / (2.0*sw2) ); if(tau2 < 0.2e-6) tau2 = 0.0; }
tau2 = tau2;
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t12,2,4); tsadd(t13,2,4); }
if( ix == 1) d3_init = d3;
t2_counter = (int) ( (d3-d3_init)*sw2 + 0.5 );
if(t2_counter % 2)
{ tsadd(t3,2,4); tsadd(t12,2,4); tsadd(t13,2,4); }
status(A);
decpower(pwClvl);
dec2power(pwNlvl);
set_c13offset("co");
zgradpulse(gzlvl6, gt6);
delay(1.0e-4);
delay(d1-gt6);
lk_hold();
rcvroff();
h_shapedpulse("pc9f",4.0,3.5,zero, 2.0e-6, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0);
delay(lambda-pwS5*0.5-pwS6*0.4);
h_shapedpulse("pc9f_",4.0,3.5,one, 2.0e-6, 0.0);
obspower(shlvl1);
/**************************************************************************/
/* xxxxxxxxxxxxxxxxxxxxxx N-> CA transfer xxxxxxxxxxxxxxxxxx */
/**************************************************************************/
set_c13offset("ca");
dec2rgpulse(pwN,zero,0.0,0.0);
delay(timeTN);
sim3_c13pulse("", "ca", "co", "square", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
delay(Delta);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
delay(timeTN-Delta-pwS3-taunh*0.5-shpw1);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
delay(taunh*0.5);
dec2rgpulse(pwN,zero,0.0,0.0);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
/**************************************************************************/
/* xxxxxxxxxxxxxxxxxxxxxx CA->CO TRANSFER xxxxxxxxxxxxxxxxxx */
/**************************************************************************/
set_c13offset("ca");
c13pulse("ca", "co", "square", 90.0, t2, 2.0e-6, 0.0);
/*
initval(0.0, v2);
decstepsize(1.0);
dcplrphase(v2);
*/
zgradpulse(gzlvl4, gt4);
delay(tauC*0.5-gt4);
c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,zero,0.0,0.0);
zgradpulse(gzlvl4, gt4);
delay(tauC-gt4);
c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,two,0.0,0.0);
delay(tauC*0.5);
c13pulse("ca", "co", "square", 90.0, one, 0.0, 0.0);
zgradpulse(gzlvl3, gt3*3.5);
delay(1.0e-4);
/**************************************************************************/
/* xxxxxxxxxxxxxxxxxxxxxx 13CO EVOLUTION xxxxxxxxxxxxxxxxxx */
/**************************************************************************/
set_c13offset("co");
c13pulse("co", "ca", "sinc", 90.0, t1, 2.0e-6, 0.0);
delay(tau1*0.5);
sim3_c13pulse(shname1, "ca", "co", "square", "",shpw1, 180.0, 2.0*pwN,
zero, zero, zero, 0.0, 0.0);
delay(tau1*0.5);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
sim3_c13pulse(shname1, "ca", "co", "square", "",shpw1, 180.0, 0.0,
two, zero, zero, 0.0, 0.0);
if (pwN*2.0 > pwS2) delay(pwN*2.0-pwS2);
c13pulse("co", "ca", "sinc", 90.0, t4, 0.0, 0.0);
/***************************************************************************/
/* CA->CO transfer */
/***************************************************************************/
set_c13offset("ca");
c13pulse("ca", "co", "square", 90.0, zero, 2.0e-6, 0.0);
initval(0.0, v2);
decstepsize(1.0);
dcplrphase(v2);
zgradpulse(gzlvl3, gt3*2.0);
delay(tauC*0.5-gt3*2.0);
c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,two,0.0,0.0);
zgradpulse(gzlvl3, gt3*2.0);
delay(tauC-gt3*2.0);
c_shapedpulse2("isnob5",20.0,0.0,"isnob5",20.0,119.0,zero,0.0,0.0);
delay(tauC*0.5);
c13pulse("ca", "co", "square", 90.0, one, 0.0, 0.0);
/* dcplrphase(v2); */
/**************************************************************************/
obspower(shlvl1);
shaped_pulse(shname1,shpw1,zero,0.0,0.0);
dec2rgpulse(pwN,t3,0.0,0.0);
delay(tau2*0.5+taunh*0.5);
shaped_pulse(shname1,shpw1,two,0.0,0.0);
delay(timeTN-shpw1-taunh*0.5-gt1-1.0e-4);
zgradpulse(-gzlvl1, gt1);
delay(1.0e-4);
sim3_c13pulse("", "ca", "co", "square", "", 0.0, 180.0, 2.0*pwN,
zero, zero, zero, 2.0e-6, 2.0e-6);
if (SE_flg[0] == 'y')
{
delay(Delta);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
delay(timeTN-tau2*0.5-pwS4-Delta-pwS3);
h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0);
dec2rgpulse(pwN, t10, 0.0, 0.0);
}
else
{
delay(Delta);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
delay(timeTN-tau2*0.5-Delta-pwS3);
dec2rgpulse(pwN, zero, 0.0, 0.0);
}
/**************************************************************************/
if (SE_flg[0] == 'y')
{
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
dec2rgpulse(pwN, one, 0.0, 0.0);
h_shapedpulse("eburp2_",4.0,3.5,one, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t10);
delay(lambda-pwS4*0.5-pwS6*0.4);
h_shapedpulse("eburp2",4.0,3.5,zero, 0.0, 0.0);
delay((gt1/10.0) + 1.0e-4 +gstab + 2.0*GRADIENT_DELAY + POWER_DELAY);
h_shapedpulse("reburp",4.0,3.5,zero, 0.0, 0.0);
zgradpulse(icosel*gzlvl2, gt1/10.0); /* 2.0*GRADIENT_DELAY */
delay(gstab);
}
else
{
h_shapedpulse("eburp2",4.0,3.5,zero, 2.0e-6, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5);
h_sim3shapedpulse("reburp",4.0,3.5,0.0,2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
delay(lambda-pwS6*0.4 - gt5-POWER_DELAY-1.0e-4);
}
dec2power(dpwr2); /* POWER_DELAY */
lk_sample();
if (SE_flg[0] == 'y')
setreceiver(t13);
else
setreceiver(t12);
rcvron();
statusdelay(C,1.0e-4 );
}
