function [flux]=load_flux_file(infile)
lcmd=sprintf('x=load(''%s'');', infile)
eval(lcmd);

flux.jdy=x(:,1); 	%Day-Of-Year at beginning of time average
flux.s1t1=x(:,2);      %psd bow true wind speed, m/s
flux.dir1t1=x(:,3);    %psd bow true wind direction, deg
flux.ts=x(:,4);        %psd seasnake T, C
flux.ts_tsg=x(:,5);    %ship thermosalinograph T, C  bow (NaN is not available)
flux.sal_tsg=x(:,6);   %ship thermosalinograph salinity, psu  bow (NaN is not available)
flux.ta=x(:,7);        %psd air T, C
flux.qs=x(:,8);        %psd air specific humidity at sea surface, g/kg
flux.qa=x(:,9);        %psd air specific humidity, g/kg
flux.psp=x(:,10);      %psd solar flux, w/m^2
flux.pir=x(:,11);      %psd IR flux, w/m^2
flux.org=x(:,12);      %psd optical raingage precip rate, mm/hr
flux.sogm5=x(:,13);    %SOG from gps, m/s
flux.cogm5=x(:,14);    %COG from gps, deg
flux.relsp1=x(:,15);   %relative bow wind speed, m/s
flux.reldir1=x(:,16);  %relative bow wind direction, deg
flux.lat=x(:,17);      %decimal latitude, deg
flux.lon=x(:,18);      %decimal longtude, deg
flux.zt=x(:,19);       %Depth of SST sensor used in heat flux calc, m
flux.sig_sp=x(:,20);   %standard deviation of ship speed, m/s
flux.taub=x(:,21);     %wind stress, coare 3.0, N/m^2
flux.hsb=x(:,22);      %sensible heat flux, coare 3.0, w/m^2
flux.hlb=x(:,23);      %latent heat flux, coare 3.0, w/m^2
flux.rf=x(:,24);       %rain heat flux, w/m^2
flux.ta_im=x(:,25);    %IMET air t, C  from RTD (NaN if not available)
flux.qa_im=x(:,26);    %IMET air specific humidity from RTD, g/kg (NaN if not available)
flux.s_shp1=x(:,27);   %ship port true wind speed, m/s (NaN if not available)
flux.dir_shp1=x(:,28); %IMET port true wind direction, deg (NaN if not available)
flux.psp_im=x(:,29);   %IMET solar flux, w/m^2 (NaN if not available)
flux.pir_im=x(:,30);   %IMET IR flux, w/m^2 (NaN if not available)
flux.pressm=x(:,31);   %Barometric pressure, mb
flux.rh_psd=x(:,32);   %RH (%)
