% To read in 4 second averages of IDM files to get B field and V field from text file

filename2 = '/hao/aim3/dmsp/dmsp_idm_4s/2005/l1505013a.txt_w_Bin'

% Read IDM data

IDM=fopen(filname2,'r');

%Need to read in columns somehow

lindat2 = fgetl(IDM);
 idmdat = sscanf(lindat2,'%f');
  nn = nn + 1;
  uts(nn)  = idmdat(1);
  babs(nn) = idmdat(16); 
  bn(nn)   = idmdat(17);
  be(nn)   = idmdat(18);
  bd(nn)   = idmdat(19);
  vy(nn)   = idmdat(2);
  vz(nn)   = idmdat(3);