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