function rfinv_array ; ; generate rf array for inversion and 5 excitation pulses ; nt= 640 t= findgen(nt)/nt-0.5 bbx= sin(t*2*!pi*1.3)/t win_g= (exp(-(2*t)^2)- exp(-1))/(1-exp(-1)) bbx= bbx*win_g bbx[nt/2]= 8.1683 bbx/=total(bbx) bbx*=!pi/2.*250. rfx= dcomplexarr(nt,5) for i=0,4 do rfx[*,i]= bbx*exp(-complex(0,1)*(t+0.5)*(i-2)*5.25*2*!pi) ; (siemens) inversion pulse bb= hypersech_pulse(1.,3.6/5.12,2*!pi,2560,10.24); scaled to amp 1. and freq 1kHz (2pi/ms) b1= bb[*,1]*2*!pi*0.75*exp(complex(0,1)*bb[*,3]*0.75) ; scale to B1 amp 750Hz in rad/us, phase for 750Hz omega rfinv= dcomplexarr(2560,6) ; combined array for sim rfinv[*,0]= b1 rfinv[0:639,1:*]= rfx return,rfinv end