window, 0, xs=1600, ys=800

fnum = 1
offset = 0

for fnum = 1,30 do begin
for offset = 0,16 do begin
print, fnum + offset/16.
x = proc2(fnum,800,800*offset)
tvscl, 20*alog10(abs(x[2001:3600,*]-mean(x[2001:3600,*]))+5)
pjpg,STRING(FORMAT='("noise_rec",I02,"_",I02,".jpg")',fnum,offset)

tmp =  20*alog10(total(abs(fft(x[2001:3600,*],dim=1)),2))
plot, tmp[800:1599], yrange=[20,60], xstyle = 1
pjpg,STRING(FORMAT='("spec_rec",I02,"_",I02,".jpg")',fnum,offset)

endfor
endfor
END