void pulsesequence()
{
double gzlvl0 = getval("gzlvl0"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gstab = getval("gstab"),
compH = getval("compH");
char sspul[MAXSTR], pshape[MAXSTR];
int iphase;
shape hdx;
getstr("sspul",sspul);
iphase = (int)(getval("phase")+0.5);
/* Make HADAMARD encoding waveforms */
getstr("pshape", pshape);
hdx = pboxHT_F1e(pshape, pw*compH, tpwr);
/* Setup Phase Cycling */
settable(t1,8,ph1);
settable(t2,8,ph2);
settable(t3,8,ph3);
getelem(t1,ct,v1);
getelem(t2,ct,v2);
getelem(t3,ct,oph);
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v10);
add(v1,v10,v1);
add(oph,v10,oph);
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);
obspower(tpwr);
status(B);
pbox_pulse(&hdx, v1, rof1, 2.0e-6);
obspower(tpwr);
rgpulse(pw, v2, 2.0e-6, rof1);
delay(gt1 + gstab + 2*GRADIENT_DELAY);
rgpulse(2*pw, v2, rof1, rof1);
zgradpulse(-gzlvl1,gt1);
delay(gstab);
rgpulse(pw, v2, rof1, rof2);
zgradpulse(gzlvl2,gt1);
delay(gstab - 2*GRADIENT_DELAY);
status(C);
}
pulsesequence()
{
char mixpat[MAXSTR], pshape[MAXSTR], httype[MAXSTR], sspul[MAXSTR];
double rg1 = 2.0e-6,
mix = getval("mix"), /* mixing time */
mixpwr = getval("mixpwr"), /* mixing pwr */
compH = getval("compH"),
gt0 = getval("gt0"), /* gradient pulse in sspul */
gt2 = getval("gt2"), /* gradient pulse preceding mixing */
gzlvl0 = getval("gzlvl0"), /* gradient level for gt0 */
gzlvl2 = getval("gzlvl2"), /* gradient level for gt2 */
gstab = getval("gstab"); /* delay for gradient recovery */
shape hdx;
void setphases();
getstr("sspul", sspul);
getstr("pshape", pshape);
getstr("httype", httype);
getstr("mixpat", mixpat);
setlimit("mixpwr", mixpwr, 48.0);
(void) setphases();
if (httype[0] == 'i')
assign(zero,v2);
/* MAKE PBOX SHAPES */
if (FIRST_FID)
hhmix = pbox_mix("HHmix", mixpat, mixpwr, pw*compH, tpwr);
/* HADAMARD stuff */
if(httype[0] == 'e') /* excitation */
hdx = pboxHT_F1e(pshape, pw*compH, tpwr);
else if(httype[0] == 'r') /* refocusing */
hdx = pboxHT_F1r(pshape, pw*compH, tpwr);
else if(httype[0] == 'd') /* DPFGSE */
{
hdx = pboxHT_F1r(pshape, pw*compH, tpwr);
if (FIRST_FID)
ref180 = hdx;
}
else /* httype[0] == 'i' */ /* inversion */
hdx = pboxHT_F1i(pshape, pw*compH, tpwr);
if (getval("htcal1") > 0.5) /* Optional fine power calibration */
hdx.pwr += getval("htpwr1");
/* THE PULSE PROGRAMM STARTS HERE */
status(A);
delay(5.0e-5);
zgradpulse(gzlvl0,gt0);
if (sspul[A] == 'y')
{
rgpulse(pw,zero,rof1,rof1);
zgradpulse(gzlvl0,gt0);
}
pre_sat();
if (getflag("wet"))
wet4(zero,one);
delay(1.0e-5);
obspower(tpwr);
delay(1.0e-5);
status(B);
if (httype[0] == 'i') /* longitudinal encoding */
{
ifzero(v1);
delay(2.0*(pw+rg1));
zgradpulse(gzlvl2,gt2);
delay(gstab);
elsenz(v1);
rgpulse(2.0*pw,v3,rg1,rg1);
zgradpulse(gzlvl2,gt2);
delay(gstab);
endif(v1);
pbox_pulse(&hdx, zero, rg1, rg1);
zgradpulse(gzlvl2,gt2);
delay(gstab);
}
else /* transverse encoding */
{
if (httype[0] == 'e')
pbox_pulse(&hdx, oph, rg1, rg1);
else
{
rgpulse(pw,oph,rg1,rg1);
if (httype[0] == 'd') /* DPFGSE */
{
zgradpulse(2.0*gzlvl2,gt2);
delay(gstab);
pbox_pulse(&ref180, oph, rg1, rg1);
zgradpulse(2.0*gzlvl2,gt2);
delay(gstab);
}
zgradpulse(gzlvl2,gt2);
delay(gstab);
pbox_pulse(&hdx, v2, rg1, rof2);
zgradpulse(gzlvl2,gt2);
delay(gstab - POWER_DELAY);
}
}
if (mix)
pbox_spinlock(&hhmix, mix, v2);
if (httype[0] == 'i')
{
zgradpulse(0.87*gzlvl2,gt2);
delay(gstab);
obspower(tpwr);
txphase(v3);
rgpulse(pw,v3,rg1,rof2);
}
status(C);
}
pulsesequence()
{
double rg1 = 2.0e-6,
compH = getval("compH"),
gt0=getval("gt0"),
gzlvl0=getval("gzlvl0"),
gt1=getval("gt1"),
gzlvl1=getval("gzlvl1"),
qlvl=getval("qlvl"),
grise=getval("grise"),
gstab=getval("gstab"),
taud2=getval("taud2"),
tau1=getval("tau1");
int icosel, iphase;
char pshape[MAXSTR], xptype[MAXSTR];
shape hdx;
getstr("pshape", pshape);
iphase = (int) (0.5 + getval("phase"));
qlvl++ ;
icosel=1; /* Default to N-type experiment */
strcpy(xptype, "N-type");
if (iphase == 2)
{
icosel=-1;
strcpy(xptype, "P-type");
}
if (FIRST_FID) fprintf(stdout,"%s COSY\n", xptype);
/* Make HADAMARD encoding waveforms */
hdx = pboxHT_F1e(pshape, pw*compH, tpwr);
status(A);
delay(1.0e-4);
zgradpulse(gzlvl0,gt0);
pre_sat();
if (getflag("wet"))
wet4(zero,one);
obspower(tpwr);
delay(1.0e-5);
status(B);
pbox_pulse(&hdx, oph, rg1, rg1);
obspower(tpwr);
zgradpulse(gzlvl1,gt1);
delay(grise);
rgpulse(pw, oph, rg1, rg1);
delay(taud2);
zgradpulse(gzlvl1,gt1);
delay(grise+taud2);
rgpulse(pw,oph,rg1,rof2);
delay(tau1);
zgradpulse(gzlvl1*qlvl*(double)icosel,gt1);
delay(grise+gstab);
status(C);
}
