pro phil_restore,outfile,x0,x1,y0,y1 close,/all ;restore,'20110924.115237.1.hssd' restore,'20110922.103517.1.hssdrf' I=hssd2.spec.i Q=hssd2.spec.q U=hssd2.spec.u V=hssd2.spec.v nwvl=n_elements(reform(I(*,0,0))) xsize=n_elements(reform(I(0,*,0))) ysize=n_elements(reform(I(0,0,*))) ;print,max(I(*,0,y0:y1)),min(I(*,0,y0:y1)) ; ;tvframe,I(*,0,y0:y1),-1000.,200. ;stop openw,1,outfile printf,1,xsize,y1-y0,x1-x0+1 count=0 wvl=findgen(x1-x0+1) for k=y0,y1 do begin for j=0,xsize-1 do begin I(x0:x1,j,k)=I(x0:x1,j,k)-min(I(x0:x1,j,k)) cont=max(I(x0:x1,j,k)) count=count+1 profile=[transpose(wvl),$ transpose(I(x0:x1,j,k))/cont,$ transpose(Q(x0:x1,j,k))/cont,$ transpose(U(x0:x1,j,k))/cont,$ transpose(V(x0:x1,j,k))/cont] printf,1,count printf,1,profile endfor endfor close,1 end