void pulsesequence()
{
double j1xh,
pwxlvl,
pwx,
tau;
pwxlvl = getval("pwxlvl");
pwx = getval("pwx");
j1xh = getval("j1xh");
tau = 1.0 / (2.0*j1xh);
settable(t1, 4, ph1);
settable(t2, 4, ph2);
getelem(t1, ct, v1);
getelem(t2, ct, oph);
/* BEGIN ACTUAL PULSE SEQUENCE CODE */
status(A);
decpower(pwxlvl);
delay(d1);
status(B);
rgpulse(pw, zero, rof1, rof1);
delay(tau - (2*pw/PI) - 2*rof1);
simpulse(2*pw,pwx,v1,zero,rof1,rof2);
decpower(dpwr);
delay(tau - POWER_DELAY);
status(C);
}
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();
}
pulsesequence()
{
double cycles,
bigtau = getval("bigtau"),
tau = getval("tau"),
satdly = getval("satdly");
char satmode[MAXSTR],
sspul[MAXSTR];
getstr("satmode",satmode);
getstr("sspul",sspul);
settable(t1,4,phs1);
settable(t2,8,phs2);
settable(t3,4,phs3);
getelem(t1,ct,v1);
getelem(t2,ct,v2);
getelem(t3,ct,v4);
assign(v1,oph);
/* calculate 'big tau' values */
cycles = bigtau/(2.0*tau);
cycles = (double)((int)((cycles/2.0) + 0.5)) * 2.0;
initval(cycles,v3);
/* equilibration period */
status(A);
delay(5.0e-5);
if (sspul[0] == 'y')
steadystate();
if (satmode[0] == 'y')
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
satpulse(satdly,v2,rof1,rof1);
}
else
delay(d1);
/* calculate exact delay and phases */
mod2(oph,v5);
incr(v5);
/* spin-echo loop */
status(B);
rgpulse(pw,v1,rof1,0.0);
starthardloop(v3);
delay(tau - p1/2.0 - rof2);
rgpulse(p1,v5,rof2,rof2);
delay(tau - p1/2.0 - rof2);
endhardloop();
/* observation period */
status(C);
}
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();
}
pulsesequence()
{
double gzlvl1,
gt1,
gstab,
hsglvl,
hsgt;
char sspul[MAXSTR];
gzlvl1 = getval("gzlvl1");
gt1 = getval("gt1");
gstab = getval("gstab");
hsglvl = getval("hsglvl");
hsgt = getval("hsgt");
getstr("sspul",sspul);
settable(t1,4,ph1);
settable(t2,4,ph2);
settable(t3,4,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);
if (sspul[0] == 'y')
{
zgradpulse(hsglvl,hsgt);
rgpulse(pw,zero,rof1,rof1);
zgradpulse(hsglvl,hsgt);
}
delay(d1);
if (getflag("wet"))
wet4(zero,one);
status(B);
rgpulse(pw, v1, rof1, rof1);
delay(d2);
zgradpulse(gzlvl1,gt1);
delay(gstab);
rgpulse(pw, v2, rof1, rof2);
zgradpulse(gzlvl1,gt1);
delay(gstab);
status(C);
}
void pulsesequence()
{
double pwx;
char rcvrsflag[MAXSTR];
pwx = getval("pwx");
getstr("rcvrs",rcvrsflag);
/* check decoupling modes */
if ( (dm[C] == 'y') || (dm[D] == 'y') || (h**o[0] == 'y') )
{
printf("dm[C], dm[D] should be set to 'n' and/or h**o should set to 'n'");
psg_abort(1);
}
if (strcmp(rcvrsflag,"yy"))
printf("rcvrs parameter should be set to 'yy'\n");
settable(t1,4,ph1);
getelem(t1,ct,v1);
assign(v1,oph);
settable(t2,4,ph2);
status(A);
obspower(tpwr);
decpower(dpwr);
delay(d1);
status(B);
delay(d2);
rgpulse(pw, t2, rof1, rof2);
status(C);
setactivercvrs("yn");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
status(B);
delay(d2);
decrgpulse(pwx, t2, rof1, rof2);
status(D);
setactivercvrs("ny");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
}
pulsesequence()
{
double cmult = getval("cmult");
settable(t1,8,ph1);
settable(t2,8,ph2);
settable(t3,8,ph3);
getelem(t1,ct,v1);
getelem(t2,ct,v2);
getelem(t3,ct,oph);
/*
mod2(id2,v10);
dbl(v10,v10);
*/
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v10);
add(v1,v10,v1);
add(oph,v10,oph);
assign(v1,v4);
status(A);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
if (satmode[0] == 'y')
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
satpulse(satdly,v4,rof1,rof1);
}
else
delay(d1);
if (getflag("wet"))
wet4(zero,one);
status(B);
rgpulse(pw, v1, rof1, rof1);
delay(d2);
rgpulse(cmult*pw, v2, rof1, rof2);
status(C);
}
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();
}
void pulsesequence()
{
char compshape[MAXSTR],
sspul[MAXSTR],
wet[MAXSTR],
composit[MAXSTR];
getstr("compshape",compshape);
getstr("composit",composit);
getstr("sspul",sspul);
getstr("wet",wet);
settable(t1,4,phs1);
settable(t2,8,phs2);
getelem(t1,ct,oph);
getelem(t2,ct,v2);
assign(oph,v1);
/* equilibrium period */
status(A);
delay(5.0e-5);
if (sspul[0] == 'y')
steadystate();
delay(d1);
if (wet[0] == 'y')
wet4(zero,one);
status(B);
pulse(p1,zero);
hsdelay(d2);
if (composit[0] == 'y')
{
if (rfwg[OBSch-1] == 'y')
shaped_pulse(compshape,4.0*pw+0.8e-6,v1,rof1,rof2);
else
comp90pulse(pw,v1,rof1,rof2);
}
else
rgpulse(pw,v1,rof1,rof2);
status(C);
}
char* intToRoman(int num) {
int ary[4] = {-1};
char *ans;
if(num == 0) {
return NULL;
}
Roman *table = (Roman*)malloc(7*sizeof(Roman));
if(table == 0) {
return NULL;
}
settable(table);
int dig = getdigitnumber(num);
ary[0] = num % 10;
ary[1] = (num /10) % 10;
ary[2] = (num /100) % 10;
ary[3] = (num /1000) % 10;
for(int i = dig+1 ; i < 4; i++) {
ary[i] = -1;
}
ans = malloc(150*sizeof(char));
memset(ans, 0, 150);
for(int i = 3; i >=0; i--) {
int j = 2*i;
if(ary[i]>=0 && ary[i] <=3) {
for(int k = 0; k < ary[i]; k++) {
strcat(ans, table[j].letter);
}
}
else if(ary[i] == 4) {
strcat(ans, table[j].letter);
strcat(ans, table[j+1].letter);
}
else if(ary[i] == 5) {
strcat(ans, table[j+1].letter);
}
else if(ary[i] > 5 && ary[i] <= 8) {
strcat(ans, table[j+1].letter);
for(int k = 0; k < ary[i] - 5; k++) {
strcat(ans, table[j].letter);
}
}
else if(ary[i] == 9) {
strcat(ans, table[j].letter);
strcat(ans, table[j+2].letter);
}
}
return ans;
}
void breakpoint_Relocation::set_target(address x) {
assert(settable(), "must be settable");
jint target_bits = (jint)internal() ? scaled_offset (x)
: runtime_address_to_index(x);
short* p = &live_bits() + 1;
add_long(p, target_bits);
assert(p == instrs(), "new target must fit");
_target = x;
}
pulsesequence()
{
settable(t1,4,phasecycle);
dec3power(dpwr3);
diplexer_override(0);
delay(d1);
dec3rgpulse(pw, t1, rof1, rof2);
setreceiver(t1);
}
pulsesequence()
{
double pd, seqtime;
double minte,ted1,ted2;
double restol, resto_local;
int vph180 = v2; /* Phase of 180 pulse */
init_mri(); /****needed ****/
restol=getval("restol"); //local frequency offset
roff=getval("roff"); //receiver offset
init_rf(&p1_rf,p1pat,p1,flip1,rof1,rof2); /* hard pulse */
calc_rf(&p1_rf,"tpwr1","tpwr1f");
init_rf(&p2_rf,p2pat,p2,flip2,rof1,rof2); /* hard pulse */
calc_rf(&p2_rf,"tpwr2","tpwr2f");
seqtime = at+(p1/2.0)+rof1+d2;
pd = tr - seqtime; /* predelay based on tr */
if (pd <= 0.0) {
abort_message("%s: Requested tr too short. Min tr = %f ms",seqfil,seqtime*1e3);
}
minte = p1/2.0 + p2 + 2*rof2 + rof1;
if(d2 > 0) {
if(d2 < minte+4e-6)
abort_message("%s: TE too short. Min te = %f ms",seqfil,minte*1e3);
}
ted1 = d2/2 - p1/2 - p2/2 + rof2 + rof1;
ted2 = d2/2 - p2/2 + rof2;
resto_local=resto-restol;
status(A);
xgate(ticks);
delay(pd);
/* --- observe period --- */
obsoffset(resto_local);
obspower(p1_rf.powerCoarse);
obspwrf(p1_rf.powerFine);
shapedpulse(p1pat,p1,oph,rof1,rof2);
/* if d2=0 no 180 pulse applied */
if (d2 > 0) {
obspower(p2_rf.powerCoarse);
obspwrf(p2_rf.powerFine);
settable(t2,2,ph180); /* initialize phase tables and variables */
getelem(t2,ct,v6); /* 180 deg pulse phase alternates +/- 90 off the rcvr */
add(oph,v6,vph180); /* oph=zero */
delay(ted1);
shapedpulse(p2pat,p2,vph180,rof1,rof2);
delay(ted2);
}
startacq(alfa);
acquire(np,1.0/sw);
endacq();
}
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();
}
void test(void)
{
int i; unsigned long k;
settable(12345,65435,34221,12345,9983651,95746118);
for(i=1;i<1000001;i++){k=LFIB4;} printf("%lu\n", k-1064612766U);
for(i=1;i<1000001;i++){k=SWB ;} printf("%lu\n", k- 627749721U);
for(i=1;i<1000001;i++){k=KISS ;} printf("%lu\n", k-1372460312U);
for(i=1;i<1000001;i++){k=CONG ;} printf("%lu\n", k-1529210297U);
for(i=1;i<1000001;i++){k=SHR3 ;} printf("%lu\n", k-2642725982U);
for(i=1;i<1000001;i++){k=MWC ;} printf("%lu\n", k- 904977562U);
for(i=1;i<1000001;i++){k=FIB ;} printf("%lu\n", k-3519793928U);
}
pulsesequence()
{
double gzlvl1,pwClvl,
gt1,
gstab;
gzlvl1 = getval("gzlvl1");
gt1 = getval("gt1");
gstab = getval("gstab");
pwClvl = getval("pwClvl");
settable(t1,4,ph1);
settable(t2,4,ph2);
settable(t3,4,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);
obspower(pwClvl);
delay(d1);
status(B);
rgpulse(pw, v1, rof1, rof1);
delay(d2);
zgradpulse(gzlvl1,gt1);
delay(gstab);
rgpulse(pw, v2, rof1, rof2);
zgradpulse(gzlvl1,gt1);
delay(gstab);
status(C);
}
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);
}
void
settables(Page *p)
{
Table *t;
Item *i;
if(p->doc==nil)
return;
for(i=p->items; i!=nil; i=i->next)
if(i->tag == Itabletag)
((Itable *)i)->table->flags |= Ttoplevel;
for(t=p->doc->tables; t!=nil; t=t->next)
settable(t);
}
int breakpoint_Relocation::pack_data() {
short* p = data();
assert(p == &live_bits(), "initializing live_bits");
*p++ = _bits;
jint target_bits = (jint) internal() ? scaled_offset (_target)
: runtime_address_to_index(_target);
if (settable()) {
// save space for set_target later
add_long(p, target_bits);
} else {
add_int (p, target_bits);
}
for (int i = 0; i < instrlen(); i++) {
add_short(p, (short)0x7777); // placeholder value until bytes can be saved
}
return p - data();
}
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);
}
void kv2t(lua_State *L, FieldValueMap &m, const CFieldDict *dict) {
CFieldDict::fconst_iterator fit;
assert(dict);
for (fit = dict->fbegin();fit != dict->fend();fit++)
{
const FIELD_ATTR &fa = fit->second;
if (dict->IsRepeat(fa.m_fid))
{
const CFieldDict *dict2;
size_t rn = 0;
m.GetRepeatCount(fa.m_fid, rn);
lua_pushnumber(L, fa.m_fid);
lua_pushnumber(L, rn);
lua_settable(L, -3);
if (dict->GetRepeat(fa.m_fid, dict2))
{
// cout<<fa.m_fid<<":"<<rn<<endl;
for(unsigned int ri = 0; ri < rn; ri++)
{
char idx[128];
sprintf(idx, "%d,%d", fa.m_fid, ri);
lua_pushstring(L, idx);
lua_newtable(L);
kv2t(L, m(fa.m_fid, ri), dict2);
lua_settable(L, -3);
}
}
}
else
{
settable(L, m, fa.m_fid);
}
}
}//kv2t;
pulsesequence()
{
double hsglvl = getval("hsglvl"),
hsgt = getval("hsgt"),
tau,
evolcorr,
taug,
mult = getval("mult"),
null = getval("null");
int phase1 = (int)(getval("phase")+0.5),
prgcycle = (int)(getval("prgcycle")+0.5),
ZZgsign;
tau = 1/(4*(getval("j1xh")));
if (mult > 0.5)
taug = 2*tau;
else
taug = 0.0;
evolcorr=2*pw+4.0e-6;
ZZgsign=-1;
if (mult == 2) ZZgsign=1;
assign(ct,v17);
assign(zero,v18);
assign(zero,v19);
if (getflag("prgflg") && (satmode[0] == 'y') && (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,ph1);
settable(t2,2,ph2);
settable(t3,8,ph3);
settable(t4,16,ph4);
settable(t5,16,ph5);
getelem(t1, v17, v1);
getelem(t3, v17, v3);
getelem(t4, v17, v4);
getelem(t2, v17, v2);
getelem(t5, v17, oph);
assign(zero,v6);
add(oph,v18,oph);
add(oph,v19,oph);
/*
mod2(id2,v14);
dbl(v14,v14);
*/
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v14);
if (phase1 == 2)
incr(v2);
add(v2,v14,v2);
add(oph,v14,oph);
if (mult > 0.5)
add(oph,two,oph);
status(A);
obspower(tpwr);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
if (satmode[0] == 'y')
{
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(v6,zero,v18,v19);
}
else
{
satpulse(satdly,zero,rof1,rof1);
if (getflag("prgflg"))
purge(v6,zero,v18,v19);
}
}
else
delay(d1);
if (getflag("wet"))
wet4(zero,one);
decpower(pwxlvl);
status(B);
if ((getflag("PFGflg")) && (getflag("nullflg")))
{
rgpulse(0.5*pw,zero,rof1,rof1);
delay(2*tau);
simpulse(2.0*pw,2.0*pwx,zero,zero,rof1,rof1);
delay(2*tau);
rgpulse(1.5*pw,two,rof1,rof1);
zgradpulse(hsglvl,hsgt);
delay(1e-3);
}
else if (null != 0.0)
{
rgpulse(pw,zero,rof1,rof1);
delay(2*tau);
simpulse(2*pw,2*pwx,zero,zero,rof1,rof1);
delay(2*tau);
rgpulse(pw,two,rof1,rof1);
if (satmode[1] == 'y')
{
if (getflag("slpsat"))
shaped_satpulse("BIRDnull",null,zero);
else
satpulse(null,zero,rof1,rof1);
}
else
delay(null);
}
rgpulse(pw,v6,rof1,rof1);
delay(tau);
simpulse(2*pw,2*pwx,zero,zero,rof1,rof1);
delay(tau);
rgpulse(pw,v1,rof1,rof1);
if (getflag("PFGflg"))
{
zgradpulse(hsglvl,2*hsgt);
delay(1e-3);
}
decrgpulse(pwx,v2,rof1,2.0e-6);
if (mult > 0.5)
{
delay(d2/2);
rgpulse(2*pw,zero,2.0e-6,2.0e-6);
delay(d2/2);
delay(taug);
simpulse(mult*pw,2*pwx,zero,zero,rof1,rof1);
delay(taug + evolcorr);
}
else
{
if (d2/2 > 0.0)
delay(d2/2 - (2*pwx/PI) - pw - 4.0e-6);
else
delay(d2/2);
rgpulse(2*pw,zero,2.0e-6,2.0e-6);
if (d2/2 > 0.0)
delay(d2/2 - (2*pwx/PI) - pw - 4.0e-6);
else
delay(d2/2);
}
decrgpulse(pwx,v4,2.0e-6,rof1);
if (getflag("PFGflg"))
{
zgradpulse(ZZgsign*0.6*hsglvl,1.2*hsgt);
delay(1e-3);
}
rgpulse(pw,v3,rof1,rof1);
delay(tau - (2*pw/PI) - 2*rof1);
simpulse(2*pw,2*pwx,zero,zero,rof1, rof2);
decpower(dpwr);
delay(tau - POWER_DELAY);
status(C);
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hy = getcp("HY",0.0,0.0,0,1);
strncpy(hy.fr,"dec",3);
strncpy(hy.to,"dec2",4);
putCmd("frHY='dec'\n");
putCmd("toHY='dec2'\n");
GP inept = getinept("ineptYX");
strncpy(inept.ch1,"dec2",4);
strncpy(inept.ch2,"obs",3);
putCmd("ch1YXinept='dec2'\n");
putCmd("ch2YXinept='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ mix = getdseq("Hmix");
strncpy(mix.t.ch,"dec",3);
putCmd("chHmixtppm='dec'\n");
strncpy(mix.s.ch,"dec",3);
putCmd("chHmixspinal='dec'\n");
// Dutycycle Protection
double simpw1 = inept.pw1;
if (inept.pw2 > inept.pw1) simpw1 = inept.pw2;
double simpw2 = inept.pw3;
if (inept.pw4 > inept.pw3) simpw2 = inept.pw4;
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHY") + 2.0*simpw1 + 2.0*simpw2;
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 = inept.t1 + inept.t2 + inept.t3 + inept.t4 +
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 = inept.t1 + inept.t2 + inept.t3 + inept.t4 +
getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,16,table1);
settable(phHhy,4,table2);
settable(phYhy,4,table3);
settable(ph1Yyxinept,4,table4);
settable(ph1Xyxinept,4,table5);
settable(ph2Yyxinept,4,table6);
settable(ph2Xyxinept,16,table7);
settable(ph3Yyxinept,8,table8);
settable(ph3Xyxinept,4,table9);
settable(phRec,8,table10);
setreceiver(phRec);
// Begin Sequence
txphase(ph1Xyxinept); decphase(phH90); dec2phase(phYhy);
obspwrf(getval("aXyxinept")); decpwrf(getval("aH90")); dec2pwrf(getval("aYhy"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to Y Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhy);
_cp_(hy,phHhy,phYhy);
decphase(zero);
// INEPT Transfer from Y to X
_dseqon(mix);
_ineptref(inept,ph1Yyxinept,ph1Xyxinept,ph2Yyxinept,ph2Xyxinept,ph3Yyxinept,ph3Xyxinept);
_dseqoff(mix);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence()
{
double j1min = getval("j1min"),
j1max = getval("j1max"),
gzlvl0 = getval("gzlvl0"),
gt0 = getval("gt0"),
tau,
tauA,
tauB,
taumb;
int phase1 = (int)(getval("phase")+0.5);
tauA = 1/(2*(j1min + 0.146*(j1max - j1min)));
tauB = 1/(2*(j1max - 0.146*(j1max - j1min)));
taumb = 1/(2*(getval("jnxh")));
tau=1/(j1min+j1max);
settable(t1,1,ph1);
settable(t2,2,ph2);
settable(t3,4,ph3);
settable(t4,1,ph4);
settable(t5,8,ph5);
settable(t6,8,ph6);
settable(t7,4,ph7);
if (getflag("jfilter"))
{
if (getflag("PFGflg"))
{
getelem(t2,ct,v3);
getelem(t7,ct,oph);
getelem(t3,ct,v4);
}
else
{
getelem(t3,ct,v3);
getelem(t6,ct,oph);
getelem(t5,ct,v4);
}
}
else
{
getelem(t2,ct,v3);
getelem(t7,ct,oph);
getelem(t3,ct,v4);
}
if (phase1 == 2)
incr(v3);
/*
mod2(id2,v10);
dbl(v10,v10);
*/
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v10);
add(v3,v10,v3);
add(oph,v10,oph);
status(A);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
if (satmode[0] == 'y')
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
satpulse(satdly,zero,rof1,rof1);
}
else
delay(d1);
if (getflag("wet"))
wet4(zero,one);
decpower(pwxlvl);
status(B);
rgpulse(pw,t1,rof1,rof2);
/* Start of J filter */
if (getflag("jfilter"))
{
if (getflag("PFGflg"))
{
zgradpulse(gzlvl0/2,gt0);
delay(tauA - gt0);
decrgpulse(pwx, zero, rof1, rof1);
zgradpulse(-gzlvl0/3,gt0);
delay(tauB - gt0);
decrgpulse(pwx, zero, rof1, rof1);
zgradpulse(-gzlvl0/6,gt0);
delay(gt0/3);
}
else
{
delay(tau - rof2 - rof1);
decrgpulse(pwx,t2,rof1,rof1);
}
}
/* End of J filter */
delay(taumb);
decrgpulse(pwx,v3,rof1,2.0e-6);
if (d2/2 > 0.0)
delay(d2/2 - 2*pwx/PI - pw - 4.0e-6);
else
delay(d2/2);
rgpulse(pw*2.0,t4,2.0e-6,2.0e-6);
if (d2/2 > 0.0)
delay(d2/2 - 2*pwx/PI - pw - 4.0e-6);
else
delay(d2/2);
decrgpulse(pwx,v4,2.0e-6,rof2);
decpower(dpwr);
status(C);
}
pulsesequence()
{
double gzlvl1 = getval("gzlvl1"),
gt1 = getval("gt1"),
zqfpw1 = getval("zqfpw1"),
zqfpwr1 = getval("zqfpwr1"),
gzlvlzq1 = getval("gzlvlzq1"),
gstab = getval("gstab"),
h1freq_local = getval("h1freq_local"),
flip1 = getval("flip1"),
flip2 = getval("flip2"),
swfactor = 9.0, /* do the adiabatic sweep over 9.0*sw */
gzlvlzq,invsw;
int iphase = (int) (getval("phase") + 0.5),
prgcycle = (int)(getval("prgcycle")+0.5);
char satmode[MAXSTR],
zqfpat1[MAXSTR],
wet[MAXSTR],
antiz_flg[MAXSTR],
alt_grd[MAXSTR];
getstr("satmode",satmode);
getstr("wet",wet);
getstr("zqfpat1",zqfpat1);
getstr("antiz_flg", antiz_flg);
getstr("alt_grd",alt_grd);
invsw = sw*swfactor;
if (invsw > 60000.0) invsw = 60000.0; /* do not exceed 60 kHz */
invsw = invsw/0.97; /* correct for end effects of the cawurst-20 shape */
assign(ct,v17);
assign(zero,v18);
assign(zero,v19);
if (getflag("prgflg") && (satmode[0] == 'y') && (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);
}
}
// sub(ct,ssctr,v12);
settable(t1,8,ph1); getelem(t1,v17,v1);
settable(t2,8,ph2); getelem(t2,v17,v2);
settable(t3,8,ph3); getelem(t3,v17,v3);
settable(t4,8,ph4); getelem(t4,v17,oph);
add(oph,v18,oph);
add(oph,v19,oph);
if (alt_grd[0] == 'y') mod2(ct,v8); /* alternate gradient sign */
if (getflag("Gzqfilt")) add(oph,two,oph);
if (iphase == 2) incr(v1);
/* HYPERCOMPLEX MODE USES REDFIELD TRICK TO MOVE AXIAL PEAKS TO EDGE */
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v6);
if ((iphase==1)||(iphase==2))
{add(v1,v6,v1); add(oph,v6,oph);}
/* BEGIN THE ACTUAL PULSE SEQUENCE */
status(A);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
if (satmode[0] == 'y')
{
if ((d1-satdly) > 0.02)
delay(d1-satdly);
else
delay(0.02);
if (getflag("slpsat"))
{
shaped_satpulse("relaxD",satdly,v6);
if (getflag("prgflg"))
shaped_purge(v1,v6,v18,v19);
}
else
{
satpulse(satdly,v6,rof1,rof1);
if (getflag("prgflg"))
purge(v1,v6,v18,v19);
}
}
else
delay(d1);
if (wet[0] == 'y')
wet4(zero,one);
obsstepsize(45.0);
initval(7.0,v7);
xmtrphase(v7);
status(B);
if (antiz_flg[0] == 'n') rgpulse(flip1*pw/90.0,v1,rof1,1.0e-6);
else rgpulse(flip1*pw/90.0+2.0*pw,v1,rof1,1.0e-6);
xmtrphase(zero);
if (d2 > 0.0)
{
if (antiz_flg[0] == 'n')
delay(d2-1.0e-6-rof1-SAPS_DELAY-(2.0*pw/3.14159)*(flip1+flip2)/90.0);
else
delay(d2-1.0e-6-rof1-SAPS_DELAY-(2.0*pw/3.14159)*(flip1+flip2)/90.0+4.0*pw);
}
else delay(0.0);
if (antiz_flg[0] == 'n') rgpulse(flip2*pw/90.0,v2,rof1,1.0e-6);
else rgpulse(flip2*pw/90.0+2.0*pw,v2,rof1,1.0e-6);
status(C);
if (getflag("Gzqfilt"))
{
obspower(zqfpwr1);
rgradient('z',gzlvlzq1);
delay(100.0e-6);
shaped_pulse(zqfpat1,zqfpw1,zero,rof1,rof1);
delay(100.0e-6);
rgradient('z',0.0);
delay(gstab);
obspower(tpwr);
}
ifzero(v8); zgradpulse(gzlvl1,gt1);
elsenz(v8); zgradpulse(-1.0*gzlvl1,gt1); endif(v8);
delay(gstab);
status(D);
rgpulse(flip2*pw/90.0,v3,rof1,rof2);
}
void pulsesequence()
{
/* DECLARE AND LOAD VARIABLES */
char
ch90shape[MAXSTR],
ch180shape[MAXSTR],
exp_mode[MAXSTR], /* flag to run 3D, or 2D time-shared 15N TROSY /13C HSQC-SE*/
decCACO[MAXSTR],
caco180shape[MAXSTR],
f1180[MAXSTR], /* Flag to start t1 @ halfdwell */
f2180[MAXSTR],
f3180[MAXSTR],
f4180[MAXSTR]; /* do TROSY on N15 and H1 */
int icosel, max_pcyc; /* used to get n and p type */
double
tpwrs,
ni2=getval("ni2"),
ni3=getval("ni3"),
tau1, tau1p,tau2,tau3,tau3p, /*evolution times in indirect dimensions */
tauNH=getval("tauNH"), /* 1/(4Jhn)*/
tauCH=getval("tauCH"), /* 1/(4Jch) */
tauCH1= getval("tauCH1"), /* tauCH/2.0+tauNH/2.0,*/ /* 1/(8Jch) +1/(8Jnh) */
tauCH2= getval("tauCH2"),
swC = getval("swC"), /* spectral widths in 13C methyls */
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
swN = getval("swN"), /* spectral widths in 15N */
pwNlvl = getval("pwNlvl"), /* power for N15 pulses */
pwN = getval("pwN"), /* N15 90 degree pulse length at pwNlvl */
ch90pwr=getval("ch90pwr"),
ch90pw=getval("ch90pw"),
ch90corr=getval("ch90corr"),
ch90dres=getval("ch90dres"),
ch90dmf=getval("ch90dmf"),
ch180pw=getval("ch180pw"),
ch180pwr=getval("ch180pwr"),
caco180pw=getval("caco180pw"),
caco180pwr=getval("caco180pwr"),
mix=getval("mix"),
tpwrsf_d = getval("tpwrsf_d"), /* fine power adustment for first soft pulse(down)*/
tpwrsf_u = getval("tpwrsf_u"), /* fine power adustment for second soft pulse(up) */
pwHs = getval("pwHs"), /* H1 90 degree pulse length at tpwrs */
compH =getval("compH"),
gstab = getval("gstab"),
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gt6 = getval("gt6"),
gt7 = getval("gt7"),
gt8 = getval("gt8"),
gt9 = getval("gt9"),
gt10 = getval("gt10"),
gzlvl0 = getval("gzlvl0"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5"),
gzlvl6 = getval("gzlvl6"),
gzlvl7 = getval("gzlvl7"),
gzlvl8 = getval("gzlvl8"),
gzlvl9 = getval("gzlvl9"),
gzlvl10 = getval("gzlvl10"),
gzlvl11 = getval("gzlvl11");
getstr("f1180",f1180);
getstr("f2180",f2180);
getstr("ch180shape",ch180shape);
getstr("ch90shape",ch90shape);
getstr("decCACO",decCACO);
getstr("caco180shape",caco180shape);
getstr("exp_mode",exp_mode);
tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69)); /*needs 1.69 times more*/
tpwrs = (int) (tpwrs); /*power than a square pulse */
if (tpwrsf_d<4095.0)
tpwrs=tpwrs+6.0; /* add 6dB to let tpwrsf_d control fine power ~2048*/
if( (exp_mode[A]!='2') && (exp_mode[A]!='3') && (exp_mode[A]!='4') )
{text_error("invalid exp_mode, Should be either 2D or 3D or 4D\n "); psg_abort(1); }
/* LOAD PHASE TABLE */
settable(t1,1,phi1);
settable(t2,4,phi2);
settable(t12,4,phi2); {tsadd(t12,2,4);}
settable(t3,1,phi3);
settable(t4,2,phi4);
settable(t5,2,phi5);
settable(t6,4,phi6);
settable(t7,8,phi7);
settable(t8,8,phi8);
/* changing sign */
if( (exp_mode[A]=='4') && (exp_mode[C]=='a') )
{tsadd(t7,2,4); tsadd(t5,2,4); }
settable(t21,1,psi1); /*trosy and SE hsqc in reverse INPET */
settable(t22,1,psi2);
settable(t23,1,psi2c);
if(exp_mode[A]=='2') {settable(t31,2,rec);}
if(exp_mode[A]=='3') {settable(t31,4,rec);}
if(exp_mode[A]=='4') {settable(t31,8,rec);}
if((dm2[A] == 'y') || (dm2[B] == 'y') || (dm2[C] == 'y') || (dm2[D] == 'y'))
{ text_error("incorrect dec2 decoupler flags! Should be 'nnnn' "); psg_abort(1); }
/* special case for swapping t2 and t3 for test purposes */
if( (exp_mode[A]=='4') && (exp_mode[B]=='x') && (ni3=1) )
{
text_error("Acquiring t3 axis in ni2 dimension (instead of t2), set nt to 8! ");
tau3 = 0.5*(d3_index/swC+0.5/swC)-pw-rof1 -pwC*2.0/M_PI ; /* increment corresponds to 13C increment */
tau3p = 0.5*(d3_index/swN+0.5/swN) -pw-rof1 -pwC -pwN*2.0/M_PI -tau3;
if(d3_index % 2) { tsadd(t7,2,4); tsadd(t8,2,4); tsadd(t31,2,4); }
if (phase2 == 2) {tsadd(t7 ,1,4); tsadd(t8 ,1,4);}
tau2=0.0;
}
else {
if (phase2 == 2) {tsadd(t2 ,1,4); tsadd(t12,1,4);}
if (phase3 == 2) {tsadd(t7 ,1,4); tsadd(t8 ,1,4);}
if(d3_index % 2) { tsadd(t2,2,4); tsadd(t12,2,4); tsadd(t31,2,4); }
tau3 = 0.5*(d4_index/swC+0.5/swC)-pw -rof1 -pwC*2.0/M_PI ; /* increment corresponds to 13C increment */
tau3p = 0.5*(d4_index/swN+0.5/swN) -pw -rof1 -pwC -pwN*2.0/M_PI -tau3;
if(d4_index % 2) { tsadd(t7,2,4); tsadd(t8,2,4); tsadd(t31,2,4); }
tau2 = d3;
tau2 += 0.0*(-pw*4.0/M_PI-rof1*2.0);
if((f2180[A] == 'y') && (ni2 > 0.0)) {tau2 += ( 1.0 / (2.0*sw2) ); }
if(tau2 < 0.2e-6) {tau2 = 0.0;}
tau2 = tau2/2.0;
}
/* simultaneous Ntrosy-ChsqcSE, last part */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase1 == 1) {icosel = 1; }
else { tsadd(t21,2,4); tsadd(t22,2,4); tsadd(t23,2,4); icosel = -1; }
if(d2_index % 2) { tsadd(t4,2,4); tsadd(t5,2,4); tsadd(t31,2,4); }
/* ECHO-ANTIECHO + STATES-TPPI t1, t1' in TROSY/HSQC last step */
tau1 = 1.0*d2_index/swC; /* increment corresponds to 13C increment */
tau1p = 1.0*d2_index*(1.0/swN-1.0/swC);
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
/* Destroy 13C magnetization*/
decrgpulse(pwC*1.0, zero, 0.0, 0.0);
zgradpulse(-gzlvl0, gt0);
delay(gstab);
/* NOESY */
if(exp_mode[A]!='2')
{ /* 3-4D */
if(exp_mode[A]=='4')
{ /* full 4D */
/* t3 evolution, the very first HSQC */
/* Hz -> HzXz INEPT */
rgpulse(pw,two,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl7, gt0); delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0); delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab); zgradpulse(gzlvl7, gt0); delay(gstab);
rgpulse(pw, one, rof1, rof1);
/* water defoc-refoc */
delay(gstab); zgradpulse(gzlvl8, gt8); delay(gstab);
/* t3 time */
if((ni3==0))
{
dec2rgpulse(pwN,t7,0.0,0.0);
dec2rgpulse(pwN,two,0.0,0.0);
decrgpulse(pwC, t8, 0.0, 0.0);
decrgpulse(pwC, two, 0.0, 0.0);
rgpulse(pw*2.0, zero, rof1, rof1);
delay(pwN*2.0+pwC*2.0);
}
else
{
dec2rgpulse(pwN,t7,0.0,0.0);
delay(tau3p);
decrgpulse(pwC, t8, 0.0, 0.0);
delay(tau3);
rgpulse(pw*2.0, zero, rof1, rof1);
delay(tau3);
decrgpulse(pwC, two, 0.0, 0.0);
delay(tau3p);
dec2rgpulse(pwN,two,0.0,0.0);
}
/* water defoc-refoc */
delay(gstab); zgradpulse(gzlvl8, gt8); delay(gstab);
/* back inept, water to +Z */
rgpulse(pw,one,rof1,rof1);
zgradpulse(gzlvl9, gt0); delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0); delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab); zgradpulse(gzlvl9, gt0); delay(gstab);
rgpulse(pw, two, rof1, rof1);
/* purge */
zgradpulse(gzlvl10, gt10); delay(2.0*gstab);
} /*end of full 4D */
/************************* t2 evolution, 1H and NOE****************************** */
/* for the case of no flipbacks in NOE part of the experiment, shift first pulse
in t2 time by 45 deg and let water bring itself back at the end
of mixing time by radiation dumping */
if( (exp_mode[D]=='t') )
{ initval(1.0, v10);
obsstepsize(45.0);
xmtrphase(v10);
}
else
{
xmtrphase(zero);
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc",pwHs,t12,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
}
rgpulse(pw,t2,rof1,rof1);
xmtrphase(zero); /* SAPS_DELAY */
delay(tau2);
decrgpulse(2.0*pwC,zero,0.0,0.0); dec2rgpulse(2.0*pwN,zero,0.0,0.0);
delay(tau2);
rgpulse(pw*2.0,zero,rof1,rof1);
delay(pwN*2.0+pwC*2.0 + SAPS_DELAY);
rgpulse(pw,zero,rof1,rof1);
if( (exp_mode[D]!='t') )
{
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc",pwHs,zero,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
}
/* NOESY period */
delay(mix-gt2-4.0*gstab );
zgradpulse(gzlvl2, gt2);
delay(4.0*gstab);
} /* end 3-4 D acquisition */
/* N-TROSY/C-HSQCse */
/* Hz -> HzXz INEPT */
rgpulse(pw,two,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl0, gt0);
delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0);
delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab);
zgradpulse(gzlvl0, gt0);
delay(gstab);
rgpulse(pw, one, rof1, rof1);
/* on HzXz now */
/* water flipback*/
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc",pwHs,two,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
/* purge */
zgradpulse(gzlvl3, gt3);
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);
}
pulsesequence()
{
double slpwrT = getval("slpwrT"),
slpwT = getval("slpwT"),
mixT = getval("mixT"),
trim = getval("trim"),
tauz1 = getval("tauz1"),
tauz2 = getval("tauz2"),
tauz3 = getval("tauz3"),
tauz4 = getval("tauz4"),
selpwrA = getval("selpwrA"),
selpwA = getval("selpwA"),
gzlvlA = getval("gzlvlA"),
gtA = getval("gtA"),
selpwrB = getval("selpwrB"),
selpwB = getval("selpwB"),
gzlvlB = getval("gzlvlB"),
gtB = getval("gtB"),
gstab = getval("gstab"),
selfrq = getval("selfrq");
char selshapeA[MAXSTR], selshapeB[MAXSTR], slpatT[MAXSTR],
alt_grd[MAXSTR];
//synchronize gradients to srate for probetype='nano'
// Preserve gradient "area"
gtA = syncGradTime("gtA","gzlvlA",1.0);
gzlvlA = syncGradLvl("gtA","gzlvlA",1.0);
gtB = syncGradTime("gtB","gzlvlB",1.0);
gzlvlB = syncGradLvl("gtB","gzlvlB",1.0);
getstr("slpatT",slpatT);
getstr("selshapeA",selshapeA);
getstr("selshapeB",selshapeB);
getstr("alt_grd",alt_grd);
/* alternate gradient sign on every 2nd transient */
if (strcmp(slpatT,"mlev17c") &&
strcmp(slpatT,"dipsi2") &&
strcmp(slpatT,"dipsi3") &&
strcmp(slpatT,"mlev17") &&
strcmp(slpatT,"mlev16"))
abort_message("SpinLock pattern %s not supported!.\n", slpatT);
assign(ct,v17);
assign(v17,v6);
if (getflag("zqfilt"))
{ hlv(v6,v6); hlv(v6,v6); }
settable(t1,4,ph1); getelem(t1,v6,v1);
settable(t3,8,ph3); getelem(t3,v6,v11);
settable(t4,8,ph4);
settable(t5,4,ph5); getelem(t5,v6,v5);
settable(t2,4,ph2); getelem(t2,v6,v2);
settable(t7,8,ph7); getelem(t7,v6,v7);
settable(t8,4,ph8); getelem(t8,v6,v8);
settable(t11,16,phs8); getelem(t11,v6,v3); /* 1st echo in ES */
settable(t12,16,phs9); getelem(t12,v6,v4); /* 2nd exho in ES */
if (getflag("zqfilt"))
getelem(t4,v6,oph);
else
assign(v1,oph);
add(oph,v5,oph); mod4(oph,oph);
sub(v2,one,v21);
add(v21,two,v23);
mod4(ct,v10);
if (alt_grd[0] == 'y') mod2(ct,v12); /* alternate gradient sign on every 2nd transient */
/* BEGIN THE ACTUAL PULSE SEQUENCE */
status(A);
if (getflag("lkgate_flg")) lk_sample(); /* turn lock sampling on */
obspower(tpwr);
delay(5.0e-5);
if (getflag("sspul"))
steadystate();
delay(d1);
if (getflag("lkgate_flg")) lk_hold(); /* turn lock sampling off */
status(B);
rgpulse(pw, v1, rof1, rof1);
if (selfrq != tof)
obsoffset(selfrq);
ifzero(v12); zgradpulse(gzlvlA,gtA);
elsenz(v12); zgradpulse(-gzlvlA,gtA); endif(v12);
delay(gstab);
obspower(selpwrA);
shaped_pulse(selshapeA,selpwA,v1,rof1,rof1);
obspower(tpwr);
ifzero(v12); zgradpulse(gzlvlA,gtA);
elsenz(v12); zgradpulse(-gzlvlA,gtA); endif(v12);
delay(gstab);
if (selfrq != tof)
delay(2*OFFSET_DELAY);
ifzero(v12); zgradpulse(gzlvlB,gtB);
elsenz(v12); zgradpulse(-gzlvlB,gtB); endif(v12);
delay(gstab);
obspower(selpwrB);
shaped_pulse(selshapeB,selpwB,v2,rof1,rof1);
obspower(slpwrT);
ifzero(v12); zgradpulse(gzlvlB,gtB);
elsenz(v12); zgradpulse(-gzlvlB,gtB); endif(v12);
delay(gstab);
if (selfrq != tof)
obsoffset(tof);
if (mixT > 0.0)
{
rgpulse(trim,v11,0.0,0.0);
if (dps_flag)
rgpulse(mixT,v21,0.0,0.0);
else
SpinLock(slpatT,mixT,slpwT,v21);
}
if (getflag("zqfilt"))
{
obspower(tpwr);
rgpulse(pw,v7,1.0e-6,rof1);
ifzero(v10); delay(tauz1); endif(v10);
decr(v10);
ifzero(v10); delay(tauz2); endif(v10);
decr(v10);
ifzero(v10); delay(tauz3); endif(v10);
decr(v10);
ifzero(v10); delay(tauz4); endif(v10);
rgpulse(pw,v8,rof1,2.0e-6);
}
ExcitationSculpting(v3,v4,v12);
delay(rof2);
status(C);
}
long read_ALOS_data_SS (FILE *imagefile, FILE *outfile, struct PRM *prm, long *byte_offset,
int *nsub, int *burst_skip, int *num_burst) {
char *data, *shift_data, *gap_data;
int header_size; /* file header size 720 bytes */
int line_prefix_size; /* line header size 412 bytes*/
int record_length0; /* data record size start 10788 bytes */
int totrecl; /* total record length 11200 bytes = line_prefix_size + record_length */
int record_length1; /* data record size curr. 10788 bytes */
int line_suffix_size; /* number of bytes after data 36 bytes */
int kswath0=0; /* subswath at start of file */
int skip_swath; /* number of subswaths to skip to get to desired subswath */
int skiprecl; /* number of lines to skip to get to desired subswath */
int data_length; /* bytes of data */
int n = 1, ishift, shift, shift0;
int k, kburst = 0, totburst = 0;
int j, ngap, nlines = 0, ntot;
int nprfchange, nburstchange;
double tbias = 0.0, get_clock();
double ttot=0.,dt=0.,tgap=0.; /* total time, burst interval, burst gap, fractional gap */
settable(12345);
if (debug) fprintf(stderr,".... reading header \n");
/* read header information */
read_sardata_info(imagefile, prm, &header_size, &line_prefix_size);
assign_sardata_params(prm, line_prefix_size, &line_suffix_size, &record_length0);
totrecl = line_prefix_size + record_length0;
set_file_position(imagefile, byte_offset, header_size);
/* get the sarting burst number and the PRF information for all the bursts and rewind the file */
for (k=0;k<5;k++) totburst = totburst + n_burst[k];
for(j = 0; j < totburst; j++) {
if ( fread((void *) &sdr,sizeof(struct sardata_info), 1, imagefile) !=1 ) break;
fseek(imagefile, record_length0, SEEK_CUR);
if (swap) swap_ALOS_data_info(&sdr);
for (k=0;k<5;k++) {
if(sdr.n_data_pixels == n_data_burst[k]) {
PRF[k]=sdr.PRF;
if(j == 0) kswath0 = k;
}
}
}
/* get the time interval of a burst and rewind the file */
dt = 0.;
for (k=0;k<5;k++) if(PRF[k] > 0.) dt = dt + (double)n_burst[k]/(0.001*PRF[k]);
tgap = dt - (double)n_burst[*nsub]/(0.001*PRF[*nsub]);
ngap = (int)(tgap*(0.001*PRF[*nsub])+0.5);
prm->num_valid_az = 6*(n_burst[*nsub] + ngap); /* use 6 bursts per aperture */
rewind(imagefile);
fseek(imagefile, header_size, SEEK_SET);
if (verbose) fprintf(stderr," totburst start_busrt# burstcycle_time %d %d %f \n", totburst, kswath0, dt);
/* set the total number of lines output based on the number of patches requested or default 1000 */
ntot = *num_burst * (n_burst[*nsub] + ngap);
if (verbose) fprintf(stderr," dt, tgap, n_burst, ngap ,ntot %f %f %d %d %d %d \n", dt, tgap, n_burst[*nsub], ngap, prm->num_valid_az,ntot);
/* allocate the memory for data */
if ((data = (char *) malloc(record_length0)) == NULL) die("couldn't allocate memory for input indata.\n","");
if ((shift_data = (char *) malloc(record_length0)) == NULL) die("couldn't allocate memory for input indata.\n","");
if ((gap_data = (char *) malloc(record_length0)) == NULL) die("couldn't allocate memory for input indata.\n","");
/* seek to beginning of next burst of sub swath nsub, read the line header, reset the parameters, and set the counters */
skip_swath = ((*nsub +5) - kswath0)%5;
skiprecl = 0;
for( k = kswath0; k < kswath0 + skip_swath; k++) skiprecl = skiprecl + n_burst[k%5];
skiprecl = skiprecl + *burst_skip * totburst;
if (verbose) fprintf(stderr, " skip_swath skiprecl %d %d \n",skip_swath, skiprecl);
fseek(imagefile,skiprecl*totrecl, SEEK_CUR);
/* recalculate the PRF at the new skip location */
for(j = 0; j < totburst; j++) {
if ( fread((void *) &sdr,sizeof(struct sardata_info), 1, imagefile) !=1 ) break;
fseek(imagefile, record_length0, SEEK_CUR);
if (swap) swap_ALOS_data_info(&sdr);
for (k=0;k<5;k++) {
if(sdr.n_data_pixels == n_data_burst[k]) {
PRF[k]=sdr.PRF;
}
}
}
/* get the time interval of a burst at the new skip location and rewind the file */
dt = 0.;
for (k=0;k<5;k++) if(PRF[k] > 0.) dt = dt + (double)n_burst[k]/(0.001*PRF[k]);
tgap = dt - (double)n_burst[*nsub]/(0.001*PRF[*nsub]);
ngap = (int)(tgap*(0.001*PRF[*nsub])+0.5);
prm->num_valid_az = 6*(n_burst[*nsub] + ngap); /* use 6 bursts per aperture */
if (verbose) fprintf(stderr," totburst start_busrt# burstcycle_time %d %d %f \n", totburst, kswath0, dt);
/* rewind and seek again to start in the correct location */
rewind(imagefile);
fseek(imagefile, header_size, SEEK_SET);
fseek(imagefile,skiprecl*totrecl, SEEK_CUR);
/* read the line header and set the parameters for this subswath */
fread((void *) &sdr,line_prefix_size, 1, imagefile);
if (swap) swap_ALOS_data_info(&sdr);
prm->clock_start = get_clock(sdr, tbias);
prm->SC_clock_start = ((double) sdr.sensor_acquisition_year)*1000 + prm->clock_start;
prm->prf = sdr.PRF;
if(prm->near_range < 0) prm->near_range = sdr.slant_range;
fseek(imagefile, -1*line_prefix_size, SEEK_CUR);
n = sdr.sequence_number -1;
//m = sdr.sequence_number;
*byte_offset = ftell(imagefile);
if (verbose) fprintf(stderr," n skiprecl byte_offset %d %d %ld \n", n, skiprecl, *byte_offset);
/* now start at the beginning but with the intervals known */
if (verbose) fprintf(stderr,".... reading data (byte %ld) \n",ftell(imagefile));
shift0 = 0;
//m = 0;
/* read the rest of the file */
while ( (fread((void *) &sdr,sizeof(struct sardata_info), 1, imagefile)) == 1 ) {
n++;
if (swap) swap_ALOS_data_info(&sdr);
/* check to make sure there is no prf-change in any of the subswaths */
for (k=0;k<5;k++) {
if(sdr.n_data_pixels == n_data_burst[k]) {
if ((sdr.PRF) != PRF[k]) {
PRF[*nsub] = 0.;
}
}
}
/* detect a different subswath and skip intil the next subswath */
if(sdr.n_data_pixels == n_data_burst[*nsub]) {
kburst++;
if (sdr.sequence_number != n) printf(" missing line: n, seq# %d %d \n", n, sdr.sequence_number);
/* check for changes in record_length and PRF */
record_length1 = sdr.record_length - line_prefix_size;
if (record_length0 != record_length1) die("record_length changed","");
/* if prf changes exit */
if ((sdr.PRF) != PRF[*nsub]) {
fprintf(stderr," ERROR PRF changed, oldPRF, newPRF %f %f \n",PRF[*nsub]*.001,sdr.PRF*.001);
*byte_offset=ftell(imagefile);
nprfchange = (*byte_offset - header_size)/totrecl;
nburstchange = nprfchange/totburst;
fprintf(stderr," rec# burst# %d %d \n",nprfchange,nburstchange);
break;
}
/* check shift to see if it varies from beginning or from command line value */
check_shift(prm, &shift, &ishift, &shift0, record_length1);
if ((verbose) && (n/2000.0 == n/2000)) {
fprintf(stderr," Working on line %d prf %f record length %d slant_range %d \n"
,sdr.sequence_number, 0.001*sdr.PRF, record_length1, sdr.slant_range);
}
/* read data (and trailing bytes) */
if ( fread ((char *) data, record_length1, (size_t) 1, imagefile) != 1 ) break;
data_length = 2*n_data_burst[*nsub];
/* write line header to output data */
fwrite((void *) &sdr, line_prefix_size, 1, outfile);
nlines++;
/* check to see if this is enough data */
if(nlines >= ntot) break;
/* Shimada says the first 13 lines are bad. I only saw 11 bad. shift these lines outside the window */
if(kburst < 12) ishift = data_length;
/* ishift and write the data */
fill_shift_data(shift, ishift, data_length, line_suffix_size, record_length1, data, shift_data, outfile);
/* if kburst it equal to the length of this burst then write the appropriate number of zero lines */
if(kburst == n_burst[*nsub]) {
ttot = ttot + dt;
if(verbose) fprintf(stderr," %d %d %f %f %d %f \n",*nsub+1,kburst,dt,tgap,ngap,ttot);
kburst = 0;
/* write the appropriate number of zero lines */
for(j=0;j<ngap;j++) {
fwrite((void *) &sdr, line_prefix_size, 1, outfile);
nlines++;
/* set the gap data to 15.5 */
for (k=0;k<record_length0;k++) gap_data[k]=NULL_DATA+znew%2;
fwrite((char *) gap_data, record_length1, 1, outfile);
}
}
}
else {
record_length1 = sdr.record_length - line_prefix_size;
if ( fread ((char *) data, record_length1, (size_t) 1, imagefile) != 1 ) break;
}
}
/* calculate end time and fix prf */
prm->prf = 0.001*PRF[*nsub];
prm->clock_stop = get_clock(sdr, tbias);
prm->SC_clock_stop = ((double) sdr.sensor_acquisition_year)*1000 + prm->clock_stop;
/* m is non-zero only in the event of a prf change */
prm->num_lines = nlines-1;
prm->num_patches = (int)((1.0*prm->num_lines)/(1.0*prm->num_valid_az));
if (prm->num_lines == 0) prm->num_lines = 1;
if (verbose) print_params(prm);
free(data);
free(shift_data);
fclose (outfile);
return(*byte_offset);
}
pulsesequence()
{
/* DECLARE AND LOAD VARIABLES */
char f1180[MAXSTR], /* Flag to start t1 @ halfdwell */
mag_flg[MAXSTR], /*magic-angle coherence transfer gradients */
f2180[MAXSTR], /* Flag to start t2 @ halfdwell */
C13refoc[MAXSTR], /* C13 sech/tanh pulse in middle of t1*/
NH2only[MAXSTR]; /* spectrum of only NH2 groups */
int icosel, /* used to get n and p type */
t1_counter, /* used for states tppi in t1 */
t2_counter, /* used for states tppi in t2 */
PRexp, /* projection-reconstruction flag */
ni2 = getval("ni2");
double tau1, /* t1 delay */
mix = getval("mix"), /* NOESY mix time */
tau2, /* t2 delay */
lambda = 0.94/(4.0*getval("JNH")), /* 1/4J H1 evolution delay */
tNH = 1.0/(4.0*getval("JNH")), /* 1/4J N15 evolution delay */
csa, sna,
pra = M_PI*getval("pra")/180.0,
/* temporary Pbox parameters */
bw, pws, ofs, ppm, nst, /* bandwidth, pulsewidth, offset, ppm, # steps */
/* the sech/tanh pulse is automatically calculated by the macro "biocal", */
/* 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 */
compC = getval("compC"), /* adjustment for C13 amplifier compression */
dof100, /* C13 frequency at 100ppm for both aliphatic & aromatic*/
tpwrsf = getval("tpwrsf"), /* fine power adjustment for flipback pulse*/
pwHs = getval("pwHs"), /* H1 90 degree pulse length at tpwrs */
tpwrs, /* power for the pwHs ("H2Osinc") pulse */
pwNlvl = getval("pwNlvl"), /* power for N15 pulses */
pwN = getval("pwN"), /* N15 90 degree pulse length at pwNlvl */
sw1 = getval("sw1"),
sw2 = getval("sw2"),
gzcal=getval("gzcal"),
gt1 = getval("gt1"), /* coherence pathway gradients */
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gt0 = getval("gt0"), /* other gradients */
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gstab = getval("gstab"),
gzlvl0 = getval("gzlvl0"),
gzlvl3 = getval("gzlvl3"),
gzlvl6 = getval("gzlvl6"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5");
getstr("f1180",f1180);
getstr("mag_flg",mag_flg);
getstr("f2180",f2180);
getstr("C13refoc",C13refoc);
getstr("NH2only",NH2only);
csa = cos(pra);
sna = sin(pra);
/* LOAD PHASE TABLE */
settable(t1,2,phi1);
settable(t3,4,phi3);
settable(t9,16,phi9);
settable(t10,1,phi10);
settable(t11,8,rec);
/* INITIALIZE VARIABLES */
/* maximum fine power for pwC pulses (and initialize rfst) */
rf0 = 4095.0;
rfst=0.0;
setautocal(); /* activate auto-calibration flags */
if (autocal[0] == 'n')
{
/* 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);
}
}
}
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);
}
ofs_check(H1ofs, C13ofs, N15ofs, H2ofs);
}
if (C13refoc[A]=='y') rfst = stC200.pwrf;
}
/* 180 degree adiabatic C13 pulse from 0 to 200 ppm */
dof100 = dof + 65.0*dfrq;
/* selective H20 one-lobe sinc pulse */
tpwrs = tpwr - 20.0*log10(pwHs/((compH*pw)*1.69)); /* needs 1.69 times more */
tpwrs = (int) (tpwrs); /* power than a square pulse */
/* CHECK VALIDITY OF PARAMETER RANGES */
if ( (mix - gt4 - gt5) < 0.0 )
{ text_error("mix is too small. Make mix equal to %f or more.\n",(gt4 + gt5));
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 > 20.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);
}
/* PHASES AND INCREMENTED TIMES */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase1 == 2)
tsadd(t1,1,4);
if (phase2 == 1)
{
tsadd(t10,2,4);
icosel = 1;
}
else icosel = -1;
/* Set up f1180 */
PRexp = 0;
if((pra > 0.0) && (pra < 90.0)) PRexp = 1;
if(PRexp) /* set up Projection-Reconstruction experiment */
tau1 = d2*csa;
else
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;
/* Set up f2180 */
if(PRexp)
tau2 = d2*sna;
else
{
tau2 = d3;
if((f2180[A] == 'y') && (ni2 > 1.0))
{
tau2 += ( 1.0 / (2.0*sw2) );
if(tau2 < 0.2e-6) tau2 = 0.0;
}
}
tau2 = tau2/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(t1,2,4);
tsadd(t11,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(t11,2,4);
}
/* Correct inverted signals for NH2 only spectra */
if (NH2only[A]=='y')
{
tsadd(t3,2,4);
}
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
decoffset(dof);
decpwrf(rf0);
txphase(zero);
dec2phase(zero);
delay(d1);
dec2rgpulse(pwN, zero, 0.0, 0.0); /*destroy N15 and C13 magnetization*/
decrgpulse(pwC, zero, 0.0, 0.0);
zgradpulse(gzlvl0, 0.5e-3);
delay(1.0e-4);
dec2rgpulse(pwN, one, 0.0, 0.0);
decrgpulse(pwC, one, 0.0, 0.0);
zgradpulse(0.7*gzlvl0, 0.5e-3);
txphase(t1);
decphase(zero);
dec2phase(zero);
initval(135.0,v1);
obsstepsize(1.0);
xmtrphase(v1);
delay(5.0e-4);
rcvroff();
rgpulse(pw, t1, 50.0e-6, 0.0); /* 1H pulse excitation */
xmtrphase(zero); /* SAPS_DELAY */
txphase(zero);
if (tau1 > (2.0*GRADIENT_DELAY + pwN + 0.64*pw + 5.0*SAPS_DELAY))
{
if (tau1>0.002)
{
zgradpulse(gzlvl6, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) - SAPS_DELAY);
}
else
{
delay(tau1-pwN-0.64*pw);
}
if (C13refoc[A]=='y')
sim3pulse(0.0, 2.0*pwC, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
else
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
if (tau1>0.002)
{
zgradpulse(-1.0*gzlvl6, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) - SAPS_DELAY);
}
else
{
delay(tau1-pwN-0.64*pw);
}
}
else if (tau1 > (0.64*pw + 0.5*SAPS_DELAY))
delay(2.0*tau1 - 2.0*0.64*pw - SAPS_DELAY );
rgpulse(pw, zero, 0.0, 0.0);
delay(mix - gt4 - gt5 -gstab -200.0e-6);
dec2rgpulse(pwN, zero, 0.0, 0.0);
zgradpulse(gzlvl4, gt4);
delay(gstab);
rgpulse(pw, zero, 200.0e-6,0.0); /* HSQC begins */
dec2phase(zero);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
txphase(one);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0);
rgpulse(pw, one, 0.0, 0.0);
txphase(two);
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);
}
else
{ delay(tau2);
rgpulse(2.0*pw, zero, 0.0, 0.0);
delay(gt1 + 2.0e-4 - 2.0*pw);
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);
delay(2.0e-4 - 2.0*GRADIENT_DELAY);
}
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
sim3pulse(pw, 0.0, pwN, zero, zero, t10, 0.0, 0.0);
dec2phase(zero);
zgradpulse(gzlvl5, gt5);
delay(lambda - 1.5*pwN - gt5);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
txphase(one);
dec2phase(one);
delay(lambda - 1.5*pwN - gt5);
sim3pulse(pw, 0.0, pwN, one, zero, one, 0.0, 0.0);
txphase(zero);
dec2phase(zero);
zgradpulse(1.5*gzlvl5, gt5);
delay(lambda - 1.5*pwN - gt5);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(1.5*gzlvl5, gt5);
delay(lambda - pwN - 0.5*pw - gt5);
rgpulse(pw, zero, 0.0, 0.0);
delay((gt1/10.0) + 1.0e-4+ gstab - 0.5*pw + 2.0*GRADIENT_DELAY + POWER_DELAY);
rgpulse(2.0*pw, zero, 0.0, rof1);
dec2power(dpwr2); /* POWER_DELAY */
if (mag_flg[A] == 'y')
{
magradpulse(icosel*gzcal*gzlvl2, gt1/10.0);
}
else
{
zgradpulse(icosel*gzlvl2, gt1/10.0);
}
delay(gstab);
rcvron();
statusdelay(C,1.0e-4-rof1);
setreceiver(t11);
}
pulsesequence()
{
char
f1180[MAXSTR],
f2180[MAXSTR],
mag_flg[MAXSTR]; /* y for magic angle, n for z-gradient only */
int
icosel,
t1_counter,
t2_counter;
double
ni2,
ratio, /* used to adjust t1 semi-constant time increment */
tau1,
tau2,
taua, /* ~ 1/4JCH = 1.5 ms - 1.7 ms] */
taub, /* ~ 3.3 ms */
bigTC, /* ~ 8 ms */
bigTCO, /* ~ 6 ms */
bigTN, /* ~ 12 ms */
tauc, /* ~ 5.4 ms */
taud, /* ~ 2.3 ms */
gstab, /* ~0.2 ms, gradient recovery time */
pwClvl, /* High power level for carbon on channel 2 */
pwC, /* C13 90 degree pulse length at pwClvl */
compH, /* Compression factor for H1 on channel 1 */
compC, /* Compression factor for C13 on channel 2 */
pwNlvl, /* Power level for Nitrogen on channel 3 */
pwN, /* N15 90 degree pulse lenght at pwNlvl */
maxpwCN,
bw, ofs, ppm, /* bandwidth, offset, ppm - temporary Pbox parameters */
pwCa90, /*90 "offC13" pulse at Ca(56ppm) xmtr at CO(174ppm) */
pwCa180, /*180 "offC17" pulse at Ca(56ppm) xmtr at CO(174ppm) */
pwCO90, /* 90 "offC6" pulse at CO(174ppm) xmtr at CO(174ppm)*/
pwCO180, /* 180 "offC8" pulse at CO(174ppm) xmtr at CO(174ppm)*/
pwCab180, /* 180 "offC27" pulse at Cab(46ppm) xmtr at CO(174ppm)*/
tpwrHd, /* Power level for proton decoupling on channel 1 */
pwHd, /* H1 90 degree pulse lenth at tpwrHd. */
waltzB1 = getval("waltzB1"), /* waltz16 field strength (in Hz) */
phi_CO, /* phase correction for Bloch-Siegert effect on CO */
phi_Ca, /* phase correction for Bloch-Siegert effect on Ca */
gt1,
gt2,
gt3,
gt4,
gt5,
gt6,
gt7,
gt0,
gzlvl1,
gzlvl2,
gzlvl3,
gzlvl4,
gzlvl5,
gzlvl6, /* N15 selection gradient level in DAC units */
gzlvl7,
gzlvl0, /* H1 gradient level in DAC units */
gzcal, /* gradient calibration (gcal) */
dfCa180,
dfCab180,
dfC90,
dfCa90,
dfCO180;
/* LOAD VARIABLES */
getstr("f1180",f1180);
getstr("f2180",f2180);
getstr("mag_flg", mag_flg);
gzcal = getval("gzcal");
ni2 = getval("ni2");
taua = getval("taua");
taub = getval("taub");
tauc = getval("tauc");
bigTC = getval("bigTC");
bigTCO = getval("bigTCO");
bigTN = getval("bigTN");
taud = getval("taud");
gstab = getval("gstab");
pwClvl = getval("pwClvl");
pwC = getval("pwC");
compH = getval("compH");
compC = getval("compC");
pwNlvl = getval("pwNlvl");
pwN = getval("pwN");
phi_CO = getval("phi_CO");
phi_Ca = getval("phi_Ca");
gt1 = getval("gt1");
gt2 = getval("gt2");
gt3 = getval("gt3");
gt4 = getval("gt4");
gt5 = getval("gt5");
gt6 = getval("gt6");
gt7 = getval("gt7");
gt0 = getval("gt0");
gzlvl1 = getval("gzlvl1");
gzlvl2 = getval("gzlvl2");
gzlvl3 = getval("gzlvl3");
gzlvl4 = getval("gzlvl4");
gzlvl5 = getval("gzlvl5");
gzlvl6 = getval("gzlvl6");
gzlvl7 = getval("gzlvl7");
gzlvl0 = getval("gzlvl0");
setautocal(); /* activate auto-calibration flags */
if (autocal[0] == 'n')
{
pwCa90 = getval("pwCa90");
pwCa180 = getval("pwCa180");
pwCab180 = getval("pwCab180");
pwCO90 = getval("pwCO90");
pwCO180 = getval("pwCO180");
dfCa180 = (compC*4095.0*pwC*2.0*1.69)/pwCa180; /*power for "offC17" pulse*/
dfCab180 = (compC*4095.0*pwC*2.0*1.69)/pwCab180; /*power for "offC27" pulse*/
dfC90 = (compC*4095.0*pwC*1.69)/pwCO90; /*power for "offC6" pulse */
dfCa90 = (compC*4095.0*pwC)/pwCa90; /*power for "offC13" pulse*/
dfCO180 = (compC*4095.0*pwC*2.0*1.65)/pwCO180; /*power for "offC8" pulse */
dfCa90 = (int) (dfCa90 + 0.5);
dfCa180 = (int) (dfCa180 + 0.5);
dfC90 = (int) (dfC90 + 0.5);
dfCO180 = (int) (dfCO180 + 0.5);
dfCab180 = (int) (dfCab180 +0.5);
/* power level and pulse time for WALTZ 1H decoupling */
pwHd = 1/(4.0 * waltzB1) ;
tpwrHd = tpwr - 20.0*log10(pwHd/(compH*pw));
tpwrHd = (int) (tpwrHd + 0.5);
}
else
{
if(FIRST_FID) /* call Pbox */
{
ppm = getval("dfrq"); bw = 118.0*ppm; ofs = -118.0*ppm;
offC6 = pbox_make("offC6", "sinc90n", bw, 0.0, compC*pwC, pwClvl);
offC8 = pbox_make("offC8", "sinc180n", bw, 0.0, compC*pwC, pwClvl);
offC17 = pbox_make("offC17", "sinc180n", bw, ofs, compC*pwC, pwClvl);
bw = 128.0*ppm;
offC13 = pbox_make("offC13", "square90n", bw, ofs, compC*pwC, pwClvl);
ofs = -128.0*ppm;
offC27 = pbox_make("offC27", "sinc180n", bw, ofs, compC*pwC, pwClvl);
bw = 2.8*7500.0;
wz16 = pbox_Dcal("WALTZ16", 2.8*waltzB1, 0.0, compH*pw, tpwr);
ofs_check(H1ofs, C13ofs, N15ofs, H2ofs);
}
dfC90 = offC6.pwrf; pwCO90 = offC6.pw;
dfCO180 = offC8.pwrf; pwCO180 = offC8.pw;
dfCa90 = offC13.pwrf; pwCa90 = offC13.pw;
dfCa180 = offC17.pwrf; pwCa180 = offC17.pw;
dfCab180 = offC27.pwrf; pwCab180 = offC27.pw;
tpwrHd = wz16.pwr; pwHd = 1.0/wz16.dmf;
}
maxpwCN = 2.0*pwN;
if (pwCab180 > pwN) maxpwCN = pwCab180;
/* LOAD PHASE TABLE */
settable(t1,4,phi1);
settable(t2,2,phi2);
settable(t3,8,phi3);
settable(t4,16,phi4);
settable(t5,1,phi5);
settable(t16,8,rec);
/* CHECK VALIDITY OF PARAMETER RANGES */
if(ni > 64)
{
printf("ni is out of range. Should be: 14 to 64 ! \n");
psg_abort(1);
}
/*
if(ni/sw1 > 2.0*(bigTCO))
{
printf("ni is too big, should be < %f\n", sw1*2.0*(bigTCO));
psg_abort(1);
}
*/
if(ni2/sw2 > 2.0*(bigTN - pwCO180))
{
printf("ni2 is too big, should be < %f\n",2.0*sw2*(bigTN-pwCO180));
psg_abort(1);
}
if((dm[A] == 'y' || dm[B] == 'y' ))
{
printf("incorrect dec1 decoupler flags! Should be 'nnn' ");
psg_abort(1);
}
if((dm2[A] == 'y' || dm2[B] == 'y'))
{
printf("incorrect dec2 decoupler flags! Should be 'nny' ");
psg_abort(1);
}
if( dpwr > 50 )
{
printf("don't fry the probe, DPWR too large! ");
psg_abort(1);
}
/* Phase incrementation for hypercomplex 2D data */
if (phase1 == 1)
{
tsadd(t1, 1, 4);
}
if (phase2 == 2)
{
tsadd(t5,2,4);
icosel = 1;
}
else icosel = -1;
/* Set up f1180 tau1 = t1 */
tau1 = d2;
if ((f1180[A] == 'y') && (ni > 1))
{ tau1 += (1.0/(2.0*sw1)); }
if(tau1 < 0.2e-6) tau1 = 0.0;
tau1 = tau1/4.0;
ratio = 2.0*bigTCO*sw1/((double) ni);
ratio = (double)((int)(ratio*100.0))/100.0;
if (ratio > 1.0) ratio = 1.0;
if((dps_flag) && (ni > 1))
printf("ratio = %f => %f\n",2.0*bigTCO*sw1/((double) ni), ratio);
/* Set up f2180 tau2 = t2 */
tau2 = d3;
if ((f2180[A] == 'y') && (ni2 > 1))
{ tau2 += (1.0/(2.0*sw2)); }
if(tau2 < 0.2e-6) tau2 = 0.0;
tau2 = tau2/4.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(t1,2,4);
tsadd(t16,2,4);
}
if( ix == 1) d3_init = d3 ;
t2_counter = (int)((d3-d3_init)*sw2 + 0.5);
if((t2_counter % 2))
{
tsadd(t2,2,4);
tsadd(t16,2,4);
}
decstepsize(1.0);
initval(phi_CO, v1);
initval(phi_Ca, v2);
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
delay(d1-1.0e-3);
obsoffset(tof);
decoffset(dof);
obspower(tpwr);
decpower(pwClvl);
decpwrf(4095.0);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
rcvroff();
if(gt6 > 0.2e-6)
{
delay(10.0e-6);
decrgpulse(pwC, zero, 1.0e-6, 1.0e-6);
delay(0.2e-6);
zgradpulse(gzlvl6, gt6);
}
decpwrf(dfCa180);
delay(1.0e-3);
rgpulse(pw,zero,1.0e-6,1.0e-6);
delay(2.0e-6);
zgradpulse(gzlvl0,gt0);
delay(taua - gt0 - 2.0e-6 - WFG_START_DELAY);
simshaped_pulse("","offC17",2.0*pw,pwCa180,zero,zero,1.0e-6,1.0e-6);
/* c13 offset on CO, slp 180 on Ca */
delay(taua - gt0 - 500.0e-6 - WFG_STOP_DELAY);
zgradpulse(gzlvl0,gt0);
txphase(one);
delay(500.0e-6);
rgpulse(pw, one, 1.0e-6, 1.0e-6);
decphase(zero);
delay(2.0e-6);
zgradpulse(gzlvl3,gt3);
obspower(tpwrHd);
decpwrf(dfCa90);
delay(200.0e-6);
/* c13 offset on CO, slp 90 on Ca */
decshaped_pulse("offC13", pwCa90, zero, 0.0, 0.0);
delay(taub -PRG_START_DELAY);
obsprgon("waltz16", pwHd, 180.0);
xmtron();
decpwrf(dfC90);
decphase(t1);
delay(bigTC -taub -SAPS_DELAY -PWRF_DELAY);
/* c13 offset on CO, on-res 90 on CO */
decshaped_pulse("offC6", pwCO90, t1, 0.0, 0.0);
/* CO EVOLUTION BEGINS */
decpwrf(dfCO180);
decphase(zero);
delay(bigTCO/2.0 +maxpwCN/2.0 +WFG_STOP_DELAY -2.0*pwCO90/PI -ratio*tau1);
/* c13 offset on CO, on-res 180 on CO */
decshaped_pulse("offC8", pwCO180, zero, 0.0, 0.0);
decpwrf(dfCab180);
delay(bigTCO/2.0 +(2.0 -ratio)*tau1 -PRG_STOP_DELAY);
xmtroff();
obsprgoff();
/* c13 offset on CO, slp 180 at Cab */
sim3shaped_pulse("","offC27","",0.0,pwCab180,2.0*pwN,zero,zero,zero,0.0,0.0);
obsprgon("waltz16", pwHd, 180.0);
xmtron();
decpwrf(dfCO180);
delay(bigTCO/2.0 +(2.0 -ratio)*tau1 -PRG_START_DELAY);
/* c13 offset on CO, on-res 180 on CO */
decshaped_pulse("offC8", pwCO180, zero, 0.0, 0.0);
decpwrf(dfC90);
dcplrphase(v1);
delay(bigTCO/2.0 +maxpwCN/2.0 +WFG_STOP_DELAY -2.0*pwCO90/PI -ratio*tau1 -SAPS_DELAY);
/* CO EVOLUTION ENDS */
decshaped_pulse("offC6", pwCO90, zero, 0.0, 0.0);
/* c13 offset on CO, on-res 90 on CO */
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);
}
