function readfilesreduc,file,h

im=readadfits(file,head)
sz=size(im)

split=strsplit(file,'.',/extract)
ext=split[n_elements(split)-1]

if ext eq 'fits' then begin
    idx = where(im le !values.f_infinity, complement=idx_nan)
    if idx_nan[0] ne -1 then im[idx_nan] = -3.1e+38
    if n_params() eq 2 then begin
        case sz[0] of
            1:h=head
            2:begin
                h=FitsHeaderToAd2Trailer(head)
            end
            3:begin
                ;im=transpose(im,[1,2,0])
                ;h.nbdim=+-3
                h=FitsHeaderToAd3Trailer(head)
                ;h.nbdim=+-3
            end
        endcase
    endif
endif else if sz[0] gt 1 then h=head

return,im

end