void S_position_offset_list(double posarray[], double grad, double nslices,
double resfrq, int device, int listno, codeint apv1)
{
int k;
double freqarray[MAXSLICE], nucgamma;
if (nslices <= 1)
position_offset(posarray[0],grad,resfrq,device);
else {
nucgamma = nucleus_gamma(device);
for (k=0; k<(int)nslices; k++) {
freqarray[k] = resfrq + nucgamma*grad*posarray[k];
}
create_offset_list(freqarray,(int)nslices,device,listno);
voffset(listno,apv1);
}
}
pulsesequence()
{
/* menu variables */
double slice_str,pe_length,read_str,power,angle,TR,TIR,slice_coffset;
double trim_read,trim_slice,freqir,encode_on,read_on,ir_on,sel_crush;
/* internal variables */
double pe_length_eff,GPmin,GPstep;
double aq_time,area_readout,read_rephase_str,GRR,GRead,freqreadout;
double slice_rephase_str,area_slice,GSR,GSS,GCR2,grad_factor;
double time[5],cycle_time,start_cycle[20],relaxation,finish_time;
double powir,powexcite,sel_length,pe_step,pe_min;
double slice_offset[128],GCR1,ir_str;
double B1_90,B1_excite,extra_TE;
double rf_max,ir_len,cur_angle,ramp_lenr,ramp_lens;
int a;
read_str = getval("Read_Str");
slice_str = getval("Slice_Str");
pe_length = getval("PE_length");
ir_str = getval("IR_Grad_Str");
sel_crush = getval("Sel_Crush");
ramp_lenr = getval("ramp_lenr");
ramp_lens = getval("ramp_lens");
extra_TE = getval("extra_TE");
TIR = getval("TIR");
TR = getval("TR");
power=getval("Power");
angle = getval("angle");
trim_slice = getval("Slice_Trim");
trim_read = getval("Read_Trim");
freqreadout = getval("Read_Offset");
freqir = getval("IR_Offset");
encode_on=getval("encode_on");
read_on=getval("read_on");
ir_on=getval("ir_on");
/* rf_max = getval("rf_max"); */
ir_len = getval("ir_len");
cur_angle = getval("cur_angle");
/* rf length calculation and offset lists #0*/
sel_length = 0.003*(1000.0/(slice_str*thk));
for (a=0;a<ns;a++)
{
slice_offset[a] = pss[a]* 10.0 * slice_str;
}
GCR2 = sel_crush * HZ_TO_GAUSS;
GSS = slice_str * HZ_TO_GAUSS;
/* phase encoding */
pe_length_eff = (2.0/3.14159)*pe_length;
pe_step = 1.0/(lpe*10.0*pe_length_eff);
pe_min = pe_step*(nv/2);
GPmin = pe_min * HZ_TO_GAUSS;
GPstep = -pe_step * HZ_TO_GAUSS;
/* rephase for readout */
aq_time = at;
grad_factor = sw/((double)read_str*lro*10.0);
area_readout = (0.001 + (ramp_lenr*0.001/2.0) + (aq_time/2)) * read_str*grad_factor;
read_rephase_str = -1.0*area_readout/pe_length_eff;
GRR = read_rephase_str * HZ_TO_GAUSS * trim_read;
GRead = read_str * HZ_TO_GAUSS*grad_factor;
/* rephase the slice */
area_slice = ((sel_length/2) + (ramp_lens*0.001/2.0))*slice_str;
slice_rephase_str = -1.0*area_slice/pe_length_eff;
GSR = slice_rephase_str * HZ_TO_GAUSS*trim_slice;
/* calculate the rf levels */
B1_90 = (slice_str*thk/1000.0)* B1_SINC3S;
B1_excite = B1_90*(cur_angle/90.0);
powexcite = (B1_excite/1000.0)*4096.0;
powir = 4096.0;
/* calculate the timing */
if (ir_on==1)
time[0] = 0.011 + ir_len*0.001 + 0.0001;
else
time[0] = 0.0001;
/* ir */
time[1] = sel_length + (2*ramp_lens*0.001) + 0.001; /* slice select */
time[2] = extra_TE+ pe_length; /* phase encoding etc. */
time[3] = aq_time + (2*ramp_lenr*0.001)+0.001; /* acquisition window */
time[4] = pe_length; /* rewind */
cycle_time = time[1]+time[2]+time[3]+time[4];
/* printf("\nCycle Time = %f",cycle_time); */
create_offset_list(slice_offset,ns,OBSch,0);
/*
for (a=0;a<ns;a++)
printf("\nOffset %d = %f",a,slice_offset[a]);
*/
GCR1 = ir_str * HZ_TO_GAUSS;
if (ir_on==1)
start_cycle[0] = 0.011 + TIR;
else
start_cycle[0] = 0.0001;
finish_time = start_cycle[0] + (ns*ne*cycle_time);
/* printf("\n Finish Time = %f",finish_time); */
relaxation = TR - finish_time;
/* START THE PULSE SEQUENCE */
assign(v10,zero);
assign(v11,zero);
initval(ne,v7);
status(A);
obspower(power);
/* begin phase encode loop */
peloop(seqcon[2],nv,v1,v2);
delay(relaxation);
/* ir */
obspwrf(powir);
obsoffset(freqir);
if (ir_on==1)
{
shaped_pulse("HYPSEC4",ir_len*0.001,v10,rof1,rof2);
obl_shapedgradient("","CRUSH","",0.011,0,GCR1,0,1,WAIT);
zero_all_gradients();
delay(TIR-0.011);
}
delay(0.0001);
loop(v7,v8); /* loop over T1 times measured */
msloop(seqcon[1],ns,v5,v6); /* loop over slices */
/* slice selection readout */
obspwrf(powexcite);
voffset(0,v6);
obl_shapedgradient("","","RAMPUP",ramp_lens*0.001,0,0,GSS,1,WAIT);
obl_shapedgradient("","","CONST",sel_length+0.001,0,0,GSS,1,NOWAIT);
delay(0.001);
shaped_pulse("SINC3S",sel_length,v11,rof1,rof2);
obl_shapedgradient("","","RAMPDN",ramp_lens*0.001,0,0,GSS,1,WAIT);
zero_all_gradients();
delay(extra_TE);
/* encode and rephasing */
if ((read_on==1)&&(encode_on==1))
pe3_shapedgradient("SINUSOID",pe_length,GRR,GPmin,GSR,0,GPstep,0,zero,v2,zero);
if ((read_on!=1)&&(encode_on==1))
pe3_shapedgradient("SINUSOID",pe_length,0,GPmin,GSR,0,GPstep,0,zero,v2,zero);
if ((read_on!=1)&&(encode_on!=1))
pe3_shapedgradient("SINUSOID",pe_length,0,0,GSR,0,0,0,zero,zero,zero);
if ((read_on==1)&&(encode_on!=1))
pe3_shapedgradient("SINUSOID",pe_length,GRR,0,GSR,0,0,0,zero,zero,zero);
zero_all_gradients();
/* readout */
obsoffset(freqreadout);
if (read_on==1)
{
obl_shapedgradient("RAMPUP","","",ramp_lenr*0.001,GRead,0,0,1,WAIT);
obl_shapedgradient("CONST","","",at+0.001,GRead,0,0,1,NOWAIT);
delay(0.001);
acquire(np,1.0/(sw));
obl_shapedgradient("RAMPDN","","",ramp_lenr*0.001,GRead,0,0,1,WAIT);
zero_all_gradients();
}
else
{
delay(0.001);
acquire(np,1.0/(sw));
}
/* rewind and crush */
if (encode_on==1)
{
pe3_shapedgradient("SINUSOID",pe_length,0,-GPmin,GCR2,0,-GPstep,0,zero,v2,zero);
zero_all_gradients();
}
else delay(pe_length);
endmsloop(seqcon[1],v6);
endloop(v8);
endpeloop(seqcon[2],v2);
}