% 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

lindat = fgetl(SSJ);
 weidat = sscanf(lindat,'%f');
 if abs(idmdat(8)) > 45
  nn = nn + 1;
  uts(nn)  = idmdat(1);
  babs(nn) = idmdat(6); 
  bn(nn)   = idmdat();
  be(nn)   = idmdat();
  bd(nn)   = idmdat();
  vy(nn)   = idmdat();
  vz(nn)   = idmdat();
 end
end
