function geotail_compare(flag)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
load geotail_mag_930911.txt;
temp = geotail_mag_930911;

yrs = temp(:,1)+1900;
mons = temp(:,2);
days= temp(:,3);
hrs= temp(:,4);
mins= temp(:,5);
secs= temp(:,6);
bxgsm= temp(:,7);
bygsm= temp(:,8);
bzgsm= temp(:,9);
xgsm= temp(:,15);
ygsm= temp(:,16);
zgsm= temp(:,17);
btot = sqrt(bxgsm.*bxgsm + bygsm.*bygsm + bzgsm.*bzgsm)

Mtime  = datenum(yrs,mons,days,hrs,mins,secs);

[alldata,info] = cdfread('ge_k0_mgf_19930911_v02.cdf');
ct =1;
len = length(alldata);
for jj = 1:1:len,
	Mtemptime(ct) = todatenum(alldata{jj,1});
	cxgsm(ct) = double(alldata{jj,8}(1));
	cygsm(ct) = double(alldata{jj,8}(2));
	czgsm(ct) = double(alldata{jj,8}(3));
	cbxgsm(ct) = double(alldata{jj,5}(1))/10.0;
	cbygsm(ct) = double(alldata{jj,5}(2))/10.0;
	cbzgsm(ct) = double(alldata{jj,5}(3))/10.0;
	cbtot(ct)=sqrt((double(cbxgsm(ct))*double(cbxgsm(ct)))/100.0 +(double(cbygsm(ct))*double(cbygsm(ct)))/100.0 +(double(cbzgsm(ct))*double(cbzgsm(ct)))/100.0);
	ct = ct+1;
end


plot(Mtime,btot,'b',Mtemptime,cbtot,'r');