function format_imp_mag_0(flag)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% function format_imp_mag_0
%
% purpose: format the imp 8 magnetometer files
% into matlab time format and Bx, By, Bz, (GSM)
% X, Y, Z (GSM RE)
%
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Given 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ct = 1;
ii = 1;
format long e;
mm =1;
months = strvcat('Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep,','Oct','Nov','Dec')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Load in data
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%for yy=1974:1974,
	yy=1974
	yyyy=num2str(yy)	
	ssone = ['IMP8_',yyyy,'.dat'];
	[fid,message] = fopen(ssone,'r');
	if fid == -1
		disp(message);
	end

	line = fgetl(fid);
	yyymmdd = fscanf(fid,'%s',1);
        hhmmss = fscanf(fid,'%s',1);
	yrstr = yyymmdd(1:4);
	daystr=yyymmdd(10:11);
	monstr = yyymmdd(6:8);
        yyymmddhhmmss = [daystr,'-',monstr,'-',yrstr,' ',hhmmss];

	monthdoing = monstr;
	where = strmatch(monstr,months)
 	ssout = ['imp8',yyyy(3:4),num2str(where,'%2.2d'),'.dat']; 
        [fidout,message] = fopen(ssout,'w');
        if fidout == -1
                disp(message);
        end

        Mtime(ii) = datenum(yyymmddhhmmss);
        bxgsm(ii) = fscanf(fid,'%f',1);
        bygsm(ii) = fscanf(fid,'%f',1);
        bzgsm(ii) = fscanf(fid,'%f',1);
        dummy = fscanf(fid,'%f',1);
        xgsm(ii) = fscanf(fid,'%f',1);
        ygsm(ii) = fscanf(fid,'%f',1);
        zgsm(ii) = fscanf(fid,'%f',1);
        dummy = fscanf(fid,'%f',1);
        dummy = fscanf(fid,'%f',1);
        dummy = fscanf(fid,'%f',1);
        dummy = fscanf(fid,'%f',1);

	startmin = str2num(hhmmss(4:5)) + 1;
	starthhmmss = strcat(hhmmss(1:3),num2str(startmin),':00');
	startyyymmddhhmmss = [daystr,'-',monstr,'-',yrstr,' ',starthhmmss];
	NewMtime(ct) =datenum(startyyymmddhhmmss);
	[newyear(ct),newmonth(ct),newday(ct),newhr(ct),newmin(ct),newsec(ct)] = datevec(NewMtime(ct));
	%%keyboard
   for ii=2:3390,
   %%while( monstr == monthdoing)
	%%disp('In while loop.');
	yyymmdd = fscanf(fid,'%s',1);
	hhmmss = fscanf(fid,'%s',1);
	yrstr = yyymmdd(1:4);
        daystr=yyymmdd(10:11);
        monstr = yyymmdd(6:8);
        yyymmddhhmmss = [daystr,'-',monstr,'-',yrstr,' ',hhmmss];

	Mtime(ii) = datenum(yyymmddhhmmss);
	bxgsm(ii) = fscanf(fid,'%f',1);
	bygsm(ii) = fscanf(fid,'%f',1);
	bzgsm(ii) = fscanf(fid,'%f',1);
	dummy = fscanf(fid,'%f',1);
	xgsm(ii) = fscanf(fid,'%f',1);
	ygsm(ii) = fscanf(fid,'%f',1);
	zgsm(ii) = fscanf(fid,'%f',1);
	dummy = fscanf(fid,'%f',1);
	dummy = fscanf(fid,'%f',1);
	dummy = fscanf(fid,'%f',1);
	dummy = fscanf(fid,'%f',1);
	cow = [ Mtime(ii-1) , (Mtime(ii), NewMtime(ct)]
	%keyboard
	if ( (Mtime(ii-1) < NewMtime(ct)) & (Mtime(ii) >= NewMtime(ct)) )
		(Mtime(ii-1) < NewMtime(ct)) & (Mtime(ii) >= NewMtime(ct))
		newbxgsm(ct)= interp1(Mtime(ii-1:ii),bxgsm(ii-1:ii),NewMtime(ct),'linear');	
		newbygsm(ct)= interp1(Mtime(ii-1:ii),bygsm(ii-1:ii),NewMtime(ct),'linear');	
		newbzgsm(ct)= interp1(Mtime(ii-1:ii),bzgsm(ii-1:ii),NewMtime(ct),'linear');	
		newxgsm(ct)= interp1(Mtime(ii-1:ii),xgsm(ii-1:ii),NewMtime(ct),'linear');	
		newygsm(ct)= interp1(Mtime(ii-1:ii),ygsm(ii-1:ii),NewMtime(ct),'linear');	
		newzgsm(ct)= interp1(Mtime(ii-1:ii),zgsm(ii-1:ii),NewMtime(ct),'linear');	
		count = fprintf(fidout,'%2.2d %2.2d %4.4d %2.2d %2.2d %5.3f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f\n',...
                        newday(ct),newmonth(ct),newyear(ct),newhr(ct),newmin(ct),newsec(ct),...
                        newbxgsm(ct),newbygsm(ct),newbzgsm(ct),newxgsm(ct),newygsm(ct),newzgsm(ct));
		ct = ct+1;
		NewMtime(ct)=NewMtime(ct-1)+datenum('00:01:00');
		bob = num2str(NewMtime(ct));
		[newyear(ct),newmonth(ct),newday(ct),newhr(ct),newmin(ct),newsec(ct)] = datevec(NewMtime(ct));
	end
	ii = ii + 1
   end
	disp('Exited while loop. Closing file.');
	status = fclose (fidout);	

	keyboard
	status = fclose (fid);
	figure;
	plot(Mtime,bxgsm,'b',NewMtime(1:ct-1),newbxgsm,'k*')
%%end