function [flux]=load_etl_flux_file(infile)
% this is the processed file created daily
% by the etl system
x = ones(1,1): NaN;
lcmd=sprintf('x=load(''%s'');', infile);
eval(lcmd);

flux.yday=x(:,1); 	   %Day-Of-Year at beginning of time average
flux.wspd_true=x(:,2);      %psd bow true wind speed, m/s
flux.wdir_true=x(:,3);    %psd bow true wind direction, deg
flux.sst=x(:,4);        %psd seasnake T, C
%these are nans stratus 2011
% 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.atmp=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.swr=x(:,10);      %psd solar flux, w/m^2
flux.lwr=x(:,11);      %psd IR flux, w/m^2
flux.rainrate=x(:,12);      %psd optical raingage precip rate, mm/hr
% % these are NaNs, sadly
% flux.sogm5=x(:,13);    %SOG from gps, m/s
flux.cog=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.bpr=x(:,31);   %Barometric pressure, mb  (lloks like db??)
flux.hrh=x(:,32);   %RH (%)
