PRO read_slice,file,nph,nth,time,ph,lat,$ br,bth,bph,vr,vth,vph,endian=endian if n_elements(endian) le 0 then endian='little' if endian eq 'little' then begin openr,/get_lun,unit,file endif else begin openr,/get_lun,unit,file,/swap_if_little_endian endelse print,'read ',file nph=long(0) nth=long(0) readu,unit,nph,nth time=0.D0 ph=fltarr(nph) lat=fltarr(nth) br=fltarr(nph,nth) bth=fltarr(nph,nth) bph=fltarr(nph,nth) vr=fltarr(nph,nth) vth=fltarr(nph,nth) vph=fltarr(nph,nth) readu,unit,time readu,unit,ph readu,unit,lat readu,unit,br readu,unit,bth readu,unit,bph readu,unit,vr readu,unit,vth readu,unit,vph free_lun,unit return end