function [l, u, ln] = assign_long_info_met_PISTON(the_fields, nfields)

%%% function to assign long info, units, and long variable names to the variables of cruise data
l = cell(nfields,1);     
u = cell(nfields,1);     
ln = cell(nfields,1);    
        

%%% save calibration info from data loggers, serial numbers
%%% we can make sure Sergio's aware we need that info
%%% calibrate more consistently, when cals are done and by whom
%%% each cruise, get list of each serial number and cal... used to be in
%%% log book manually. Sergio might keep a separate log of serials. To
%%% submit variations in whatever sensors that we put on. 

%% long info
l(ismember(the_fields,'ID_Ftsg_s')) = {'seawater flow rate ID for TSG: 1 = good, 0 = bad (F < 1.75 L/min or F > 3.5) main lab TSG'};
l(ismember(the_fields,'ID_Ftsg1_s'))= {'seawater flow rate ID for bow TSG: 1 = good, 0 = bad (F < 1.75 L/min or F > 3.5) main lab TSG'};
l(ismember(the_fields,'ID_Ftsg2_s'))= {'seawater flow rate ID for seachest TSG: 1 = good, 0 = bad (F < 1.75 L/min or F > 3.5) main lab TSG'};
l(ismember(the_fields,'ID_tsg_s'))  = {'0 = do not use tsg data during switch over, 1 = before switch, 2 = after switch'};
l(ismember(the_fields,'Ctsg_s'))    = {'ship thermosalinograph SBE 45 conductivity at depth ztsg_ship, S/m'};
l(ismember(the_fields,'Ctsg1_s'))   = {'ship bow thermosalinograph SBE 45 conductivity at depth ztsg1_ship, S/m'};
l(ismember(the_fields,'Ctsg2_s'))   = {'ship seachest thermosalinograph SBE 45 conductivity at depth ztsg2_ship, S/m'};
l(ismember(the_fields,'Ftsg_s'))    = {'ship thermpsalinograph flow rate at depth ztsg_ship, L/min'};
l(ismember(the_fields,'Ftsg1_s'))   = {'ship bow thermpsalinograph flow rate at depth ztsg1_ship, L/min'};
l(ismember(the_fields,'Ftsg2_s'))   = {'ship seachest thermpsalinograph flow rate at depth ztsg2_ship, L/min'};
l(ismember(the_fields,'FLtsg_s'))   = {'ship thermpsalinograph fluorometer dry value at depth ztsg_ship, micrograms/L'};
l(ismember(the_fields,'FLtsg1_s'))  = {'ship bow thermpsalinograph fluorometer dry value at depth ztsg1_ship, micrograms/L'};
l(ismember(the_fields,'FLtsg2_s'))  = {'ship seachest thermpsalinograph fluorometer dry value at depth ztsg2_ship, micrograms/L'};
l(ismember(the_fields,'Stsg_s'))    = {'ship thermosalinograph SBE 45 salinity at depth ztsg_ship, psu'};
l(ismember(the_fields,'Stsg1_s'))   = {'ship bow thermosalinograph SBE 45 salinity at depth ztsg1_ship, psu'};
l(ismember(the_fields,'Stsg2_s'))   = {'ship seachest thermosalinograph SBE 45 salinity at depth ztsg2_ship, psu'};
l(ismember(the_fields,'Ta'))        = {'PSL air T at height zt, C'};
l(ismember(the_fields,'Ta_s'))      = {'ship air T at height zt_ship, C'};
l(ismember(the_fields,'Td_s'))      = {'ship air dewpoint T at height zt_ship, C'};
l(ismember(the_fields,'Tskin'))     = {'skin T of seawater from ROSR infrared radiometer, C (true SST felt by atmosphere)'};
l(ismember(the_fields,'Tsnk'))      = {'PSL sea snake T at depth zsnk, C'};
l(ismember(the_fields,'Tsonic'))    = {'PSL air Tv virtual temp of 3D sonic anemometer at height zu, C'};
l(ismember(the_fields,'Ttsg_int_s'))= {'ship thermosalinograph SBE 45 internal T, C (use for instrument calibration only)'};
l(ismember(the_fields,'Ttsg_s'))    = {'ship thermosalinograph SBE 38 T at depth ztsg_ship, C'};
l(ismember(the_fields,'Ttsg1_s'))   = {'ship bow thermosalinograph SBE 38 T at depth ztsg1_ship, C'};
l(ismember(the_fields,'Ttsg2_s'))   = {'ship seachest thermosalinograph SBE 38 T at depth ztsg2_ship, C'};
l(ismember(the_fields,'aspir_TRH')) = {'PSL fan speed of aspirator on T/RH Vaisala sensor, rpm'};
l(ismember(the_fields,'cog'))       = {'PSL course over ground, deg'};
l(ismember(the_fields,'cog_s'))     = {'ship course over ground, deg'};
l(ismember(the_fields,'day'))       = {'day'};
l(ismember(the_fields,'heave'))     = {'PSL heave, m'};
l(ismember(the_fields,'heave_s'))   = {'ship heave, m'};
l(ismember(the_fields,'hed'))       = {'PSL heading, deg'};
l(ismember(the_fields,'hed_s'))     = {'ship heading, deg'};
l(ismember(the_fields,'hedE'))      = {'PSL heading eastward, deg'};
l(ismember(the_fields,'hedN'))      = {'PSL heading northward, deg'};
l(ismember(the_fields,'hedE_s'))    = {'ship heading northward, deg'};
l(ismember(the_fields,'hedN_s'))    = {'ship heading northward, deg'};
l(ismember(the_fields,'hour'))      = {'hour'};
l(ismember(the_fields,'jd'))        = {'decimal julian date'};
l(ismember(the_fields,'lat'))       = {'PSL lat, deg'};
l(ismember(the_fields,'lat_s'))     = {'ship lat, deg'};
l(ismember(the_fields,'lon'))       = {'PSL lon, deg'};
l(ismember(the_fields,'lon_s'))     = {'ship lon, deg'};
l(ismember(the_fields,'lw_dn_1'))   = {'PSL downwelling longwave radiation unit #1, W/m2'};
l(ismember(the_fields,'lw_dn_2'))   = {'PSL downwelling longwave radiation unit #2, W/m2'};
l(ismember(the_fields,'lw_dn_clr')) = {'calculated clear-sky downwelling longwave radiation, W/m2'};
l(ismember(the_fields,'lw_dn_s'))   = {'ship downwelling longwave radiation, W/m2'};
l(ismember(the_fields,'lw_T_case_1')) = {'PSL T of lw unit #1 case, C'};
l(ismember(the_fields,'lw_T_dome_1')) = {'PSL T of lw unit #1 dome, C'};
l(ismember(the_fields,'lw_therm_1')) = {'PSL thermopile output of lw unit #1, mV'};
l(ismember(the_fields,'lw_T_case_2')) = {'PSL T of lw unit #2 case, C'};
l(ismember(the_fields,'lw_T_dome_2')) = {'PSL T of lw unit #2 dome, C'};
l(ismember(the_fields,'lw_therm_2')) = {'PSL thermopile output of lw unit #2, mV'};
l(ismember(the_fields,'lw_T_case_s')) = {'ship T of lw unit case, C'};
l(ismember(the_fields,'lw_T_dome_s')) = {'ship T of lw unit dome, C'};
l(ismember(the_fields,'lw_therm_s')) = {'ship thermopile output of lw unit, mV'};
l(ismember(the_fields,'minute'))    = {'minute'};
l(ismember(the_fields,'month'))     = {'month'};
l(ismember(the_fields,'pitch'))     = {'PSL pitch, deg'};
l(ismember(the_fields,'pitch_s'))   = {'ship pitch, deg'};
l(ismember(the_fields,'prate'))     = {'PSL instantaneous precipitation rate from optical rain gauge, mm/hour'};
l(ismember(the_fields,'prate_s'))   = {'ship precipitation rate from ?, mm/hour'};
l(ismember(the_fields,'paccum'))    = {'PSL precipitation accumulation from optical rain gauge, mm'};
l(ismember(the_fields,'paccum_s'))  = {'ship precipitation accumulation from ?, mm'};
l(ismember(the_fields,'P'))         = {'PSL barometric pressure at height zp, mb'};
l(ismember(the_fields,'P_s'))       = {'ship barometric pressure at height zp_ship, mb'};
l(ismember(the_fields,'qa'))        = {'PSL air specific humidity at height zt, zq, zp, g/kg'};
l(ismember(the_fields,'qa_s'))      = {'ship air specific humidity at height zt, zq, zp, g/kg'};
l(ismember(the_fields,'qsnk'))      = {'saturation specific humidity calculated with sea snake value Tsnk, g/kg'};
l(ismember(the_fields,'raw_p'))     = {'PSL optical rain gauge raw uncorrected measurement'};
l(ismember(the_fields,'raw_p_despike')) = {'PSL optical rain gauge raw uncorrected measurement despiked'};
l(ismember(the_fields,'rdir'))      = {'PSL relative wind direction, deg, 3D sonic anemometer'};
l(ismember(the_fields,'rdir_raw'))  = {'PSL relative wind direction, deg, 3D sonic anemometer'};
l(ismember(the_fields,'rdir_f_s'))  = {'ship relative wind direction, deg, foremast rotating propellor'};
l(ismember(the_fields,'rdir_p_s'))  = {'ship relative wind direction, deg, port 2D sonic anemometer'};
l(ismember(the_fields,'rdir_s_s'))  = {'ship relative wind direction, deg, starboard 2D sonic anemometer'};
l(ismember(the_fields,'rdir_s'))    = {'ship relative wind direction, deg, 2D sonic anemometer'};
l(ismember(the_fields,'rh'))        = {'PSL relative humidity at height zq, %'};
l(ismember(the_fields,'rh_s'))      = {'ship relative humidity at height zq_ship, %'};
l(ismember(the_fields,'roll'))      = {'PSL roll, deg'};
l(ismember(the_fields,'roll_s'))    = {'ship roll, deg'};
l(ismember(the_fields,'rspd'))      = {'PSL relative wind speed, m/s 3D sonic anemometer'};
l(ismember(the_fields,'rspd_raw'))  = {'PSL raw relative wind speed, m/s 3D sonic anemometer'};
l(ismember(the_fields,'rspd_dir'))  = {'PSL relative wind speed, m/s 3D sonic anemometer'};
l(ismember(the_fields,'rspd_f_s'))  = {'ship relative wind speed, m/s, foremast rotating propellor'};
l(ismember(the_fields,'rspd_p_s'))  = {'ship relative wind speed, m/s, port 2D sonic anemometer'};
l(ismember(the_fields,'rspd_s_s'))  = {'ship relative wind speed, m/s, starboard 2D sonic anemometer'};
l(ismember(the_fields,'rspd_s'))    = {'ship relative wind speed, m/s, 2D sonic anemometer'};
l(ismember(the_fields,'rspd_std_s'))= {'ship standard deviation relative wind speed, m/s, 2D sonic anemometer'};
l(ismember(the_fields,'slp'))       = {'PSL barometric pressure for sea level altitude, mb'};
l(ismember(the_fields,'slp_reported_s')) = {'ship reported barometric pressure for sea level altitude, mb'};
l(ismember(the_fields,'slp_s'))     = {'ship corrected barometric pressure for sea level altitude, mb'};
l(ismember(the_fields,'sog'))       = {'PSL speed over ground, m/s'};
l(ismember(the_fields,'sogE'))      = {'PSL eastward component of speed over ground, m/s'};
l(ismember(the_fields,'sogE_s'))    = {'ship eastward component of speed over ground, m/s'};
l(ismember(the_fields,'sogN'))      = {'PSL northward component of speed over ground, m/s'};
l(ismember(the_fields,'sogN_s'))    = {'ship northward component of speed over ground, m/s'};
l(ismember(the_fields,'sog_s'))     = {'ship speed over ground, m/s'};
l(ismember(the_fields,'spd_log_u_s')) = {'ship eastward actual speed log, m/s'};
l(ismember(the_fields,'spd_log_v_s')) = {'ship northward actual speed log, m/s'};
l(ismember(the_fields,'spd_log_s'))  = {'ship actual speed log, m/s'};
l(ismember(the_fields,'spd_log_std_s'))  = {'standard deviation of ship actual speed log, m/s'};
l(ismember(the_fields,'sw_dn_2'))   = {'PSL downwelling shortwave radiation from unit #1, W/m2'};
l(ismember(the_fields,'sw_dn_1'))   = {'PSL downwelling shortwave radiation from unit #2, W/m2'};
l(ismember(the_fields,'sw_dn_clr')) = {'calculated clear-sky downwelling shortwave radiation, W/m2'};
l(ismember(the_fields,'sw_dn_s'))   = {'ship downwelling shortwave radiation, W/m2'};
l(ismember(the_fields,'t'))         = {'matlab datetime'};
l(ismember(the_fields,'wdir'))      = {'PSL true wind direction, deg, 3D sonic anemometer'};
l(ismember(the_fields,'wdir_raw'))  = {'PSL true wind direction, deg, 3D sonic anemometer'};
l(ismember(the_fields,'wdir_f_s'))  = {'ship true wind direction, deg, foremast rotating propellor'};
l(ismember(the_fields,'wdir_p_s'))  = {'ship true wind direction, deg, port 2D sonic anemometer'};
l(ismember(the_fields,'wdir_s_s'))  = {'ship true wind direction, deg, starboard 2D sonic anemometer'};
l(ismember(the_fields,'wdir_s'))    = {'ship true wind direction, deg, 2D sonic anemometer'};
l(ismember(the_fields,'wspd'))      = {'PSL wind speed, m/s, 3D sonic anemometer'};
l(ismember(the_fields,'wspd_raw'))  = {'PSL wind speed, m/s, 3D sonic anemometer'};
l(ismember(the_fields,'wspd_f_s'))  = {'ship true wind speed, m/s, foremast rotating propellor'};
l(ismember(the_fields,'wspd_p_s'))  = {'ship true wind speed, m/s, port 2D sonic anemometer'};
l(ismember(the_fields,'wspd_s_s'))  = {'ship true wind speed, m/s, starboard 2D sonic anemometer'};
l(ismember(the_fields,'wspd_s'))    = {'ship true wind speed, m/s, 2D sonic anemometer'};
l(ismember(the_fields,'Tsnk_best'))    = {'Tsnk filled with either ROSR (plus cool skin) or TSG, C'};
l(ismember(the_fields,'Ta_best'))      = {'Ta from PSL, with holes filled by ship, C'};
l(ismember(the_fields,'qa_best'))      = {'qa from PSL, with holes filled by ship, g/kg'};
l(ismember(the_fields,'rh_best'))      = {'rh from PSL, with holes filled by ship, %'};
l(ismember(the_fields,'lw_dn_best'))   = {'lw_dn from PSL, W/m2'};
l(ismember(the_fields,'sw_dn_best'))   = {'sw_dn from PSL, with holes filled by ship, C'};
l(ismember(the_fields,'slp_best'))     = {'slp from PSL, with holes filled by ship, mb'};
l(ismember(the_fields,'wspd_best'))    = {'wspd from PSL, with holes filled by ship, m/s'};
l(ismember(the_fields,'rspd_best'))    = {'rspd from PSL, with holes filled by ship, m/s'};
l(ismember(the_fields,'sog_std_s'))    = {'standard deviation of ship SOG, m/s'};
l(ismember(the_fields,'sog_min_s'))    = {'min of ship SOG, m/s'};
l(ismember(the_fields,'sog_max_s'))    = {'max of ship SOG, m/s'};
l(ismember(the_fields,'cog_std_s'))    = {'standard deviation of ship COG, deg'};
l(ismember(the_fields,'cog_min_s'))    = {'min of ship COG, deg'};
l(ismember(the_fields,'cog_max_s'))    = {'max of ship COG, deg'};
l(ismember(the_fields,'hed_std_s'))    = {'standard deviation of ship heading, deg'};
l(ismember(the_fields,'hed_min_s'))    = {'min of ship heading, deg'};
l(ismember(the_fields,'hed_max_s'))    = {'max of ship heading, deg'};
l(ismember(the_fields,'cog_std'))      = {'standard deviation of PSL COG, deg'};
l(ismember(the_fields,'sog_std'))      = {'standard deviation of PSL SOG, m/s'};
l(ismember(the_fields,'hed_std'))      = {'standard deviation of PSL heading, deg'};
l(ismember(the_fields,'cog_max'))      = {'max of PSL COG, deg'};
l(ismember(the_fields,'sog_max'))      = {'max of PSL SOG, m/s'};
l(ismember(the_fields,'hed_max'))      = {'max of PSL heading, deg'};
l(ismember(the_fields,'cog_min'))      = {'min of PSL COG, deg'};
l(ismember(the_fields,'sog_min'))      = {'min of PSL SOG, m/s'};
l(ismember(the_fields,'hed_min'))      = {'min of PSL heading, deg'};
l(ismember(the_fields,'shed_std'))      = {'standard deviation of sine of PSL heading, radians'};
l(ismember(the_fields,'ched_std'))      = {'standard deviation of cosine of PSL heading, radians'};
l(ismember(the_fields,'hnet'))         = {'net heat flux, positive warming ocean, W/m2'};
l(ismember(the_fields,'lw_net'))       = {'net long wave flux, positive warming ocean, W/m2'};
l(ismember(the_fields,'sw_net'))       = {'net short wave flux, positive warming ocean, W/m2'};
l(ismember(the_fields,'dT'))           = {'SST - T10 = air sea T difference PSL, C'};
l(ismember(the_fields,'dq'))           = {'qs - q10 = air sea q difference PSL, g/kg'};
l(ismember(the_fields,'qs'))          = {'sea surface specific humidity PSL, g/kg'};
l(ismember(the_fields,'qs_s'))          = {'sea surface specific humidity ship, g/kg'};
l(ismember(the_fields,'rhoa'))         = {'air density by PSL instruments, kg/m3'};
l(ismember(the_fields,'rhoa_s'))       = {'air density by ship, kg/m3'};
l(ismember(the_fields,'H2O'))          = {'H2O number concentration PSL Vaisala, mmol/m3'};
l(ismember(the_fields,'H2O_s'))        = {'H2O number concentration ship met, mmol/m3'};
l(ismember(the_fields,'licor_AGC'))    = {'Licor box T, aka AGC value used for quality control. Good data have AGC < 60.'};
l(ismember(the_fields,'licor_rhoa'))   = {'Licor air density, kg/m3'};
l(ismember(the_fields,'licor_qa'))     = {'Licor specific humidity, g/kg'};
l(ismember(the_fields,'licor_H2O'))    = {'Licor H2O number concentration, mmol/m3'};
l(ismember(the_fields,'licor_H2O_MR')) = {'Licor H2O mixing ratio, mol/mol'};
l(ismember(the_fields,'licor_Tbox'))   = {'Licor box T, C'};
l(ismember(the_fields,'licor_Pbox'))   = {'Licor box P, kPa'};
l(ismember(the_fields,'licor_H2O_std'))= {'Licor H2O standard deviation, mmol/m3'};
l(ismember(the_fields,'licor_Pbox_std'))= {'Licor Pbox stanard deviation, kPa'};
l(ismember(the_fields,'rUn_wxt'))      = {'WXT met station northward positive wind speed, m/s'};
l(ismember(the_fields,'rUw_wxt'))      = {'WXT met station westward positive wind speed, m/s'};
l(ismember(the_fields,'wspd_wxt'))      = {'WXT met station true wind speed, m/s'};
l(ismember(the_fields,'wdir_wxt'))      = {'WXT met station true wind direction, deg'};
l(ismember(the_fields,'rspd_wxt'))      = {'WXT met station relatve wind speed, m/s'};
l(ismember(the_fields,'rdir_wxt'))      = {'WXT met station relative wind direction, deg'};
l(ismember(the_fields,'Ta_wxt'))       = {'WXT met station air T, C'};
l(ismember(the_fields,'qa_wxt'))       = {'WXT met station air specific humidity, g/kg'};
l(ismember(the_fields,'rh_wxt'))       = {'WXT met station relative humidity, m/s'};
l(ismember(the_fields,'rhoa_wxt'))     = {'WXT met station air density, m/s'};
l(ismember(the_fields,'H2O_wxt'))      = {'WXT met station H2O concentration, mmol/m3'};
l(ismember(the_fields,'prate_wxt'))    = {'WXT met station precip rate, mm/hr'};
l(ismember(the_fields,'paccum_wxt'))   = {'WXT met station precip accumlulation, mm'};
l(ismember(the_fields,'slp_wxt'))      = {'WXT met station sea level pressure, mb'};
l(ismember(the_fields,'P_wxt'))        = {'WXT met station measured pressure, mb'};
l(ismember(the_fields,'S_0p85'))       = {'salinity measured by OSU Otter at 0.85 m depth, psu'};
l(ismember(the_fields,'ID_wh_otter'))  = {'OSU otter present = 1'};
l(ismember(the_fields,'ID_modflow'))   = {'1 = moderate to poor wind directions: absolute value of rdir > 90 deg'};
l(ismember(the_fields,'ID_badflow'))   = {'1 = bad wind directions: absolute value of rdir > 120 deg'};
l(ismember(the_fields,'Ptsg_s'))       = {'seawater pressure at ship TSG, dbar'};
l(ismember(the_fields,'Sigtsg_s'))     = {'seawater density at ship TSG, kg/m3'};
l(ismember(the_fields,'wspd_min_s'))   = {'min ship true wind speed, m/s'};
l(ismember(the_fields,'wspd_max_s'))   = {'max ship true wind speed, m/s'};
l(ismember(the_fields,'wspd_std_s'))   = {'standard deviation ship true wind speed, m/s'};
l(ismember(the_fields,'wdir_min_s'))   = {'min ship true wind direction, deg'};
l(ismember(the_fields,'wdir_max_s'))   = {'max ship true wind direction, deg'};
l(ismember(the_fields,'wdir_std_s'))   = {'standard deviation ship true wind direction, deg'};
l(ismember(the_fields,'lw_dn_std'))    = {'standard deviation of PSL lw_dn_best, W/m2'};
l(ismember(the_fields,'sw_dn_std'))    = {'standard deviation of PSL lw_dn_best, W/m2'};
l(ismember(the_fields,'usr'))          = {'u_star PSL, m/s'};
l(ismember(the_fields,'usr_s'))        = {'u_star ship, m/s'};
l(ismember(the_fields,'tau'))          = {'bulk wind stress PSL, N/m2'};
l(ismember(the_fields,'tau_s'))        = {'bulk wind stress ship, N/m2'};
l(ismember(the_fields,'hs'))          = {'bulk sensible heat flux PSL, W/m2'};
l(ismember(the_fields,'hs_s'))        = {'bulk sensible heat flux ship, W/m2'};
l(ismember(the_fields,'hl'))          = {'bulk latent heat flux PSL, W/m2'};
l(ismember(the_fields,'hl_s'))        = {'bulk latent heat flux ship, W/m2'};
l(ismember(the_fields,'hlwebb'))       = {'Webb correction PSL W/m2 for latent heat flux'};
l(ismember(the_fields,'hlwebb_s'))     = {'Webb correction ship W/m2 for latent heat flux'};
l(ismember(the_fields,'tsr'))          = {'T_star temperature scaling parameter PSL, K'};
l(ismember(the_fields,'tsr_s'))        = {'T_star temperature scaling parameter ship, K'};
l(ismember(the_fields,'qsr'))          = {'q_star specific humidity scaling parameter PSL, g/kg'};
l(ismember(the_fields,'qsr_s'))        = {'q_star specific humidity scaling parameter ship, g/kg'};
l(ismember(the_fields,'zot'))          = {'z over t roughness length with respect to T PSL, m'};
l(ismember(the_fields,'zot_s'))        = {'z over t roughness length with respect to T ship, m'};
l(ismember(the_fields,'zoq'))          = {'z over t roughness length with respect to q PSL, m'};
l(ismember(the_fields,'zoq_s'))        = {'z over t roughness length with respect to q ship, m'};
l(ismember(the_fields,'Cd'))           = {'drag coefficient, momentum transfer coefficient PSL'};
l(ismember(the_fields,'Cd_s'))         = {'drag coefficient, momentum transfer coefficient ship'};
l(ismember(the_fields,'Ch'))           = {'Stanton number, sensible heat transfer coefficient PSL'};
l(ismember(the_fields,'Ch_s'))         = {'Stanton number, sensible heat transfer coefficient PSL'};
l(ismember(the_fields,'Ce'))           = {'Dalton number, latent heat transfer coefficient PSL'};
l(ismember(the_fields,'Ce_s'))         = {'Dalton number, latent heat transfer coefficient PSL'};
l(ismember(the_fields,'L'))            = {'Monin-Obukhov length scale ship, m'};
l(ismember(the_fields,'L_s'))          = {'Monin-Obukhov length scale ship, m'};
l(ismember(the_fields,'zet'))          = {'Monin-Obukhov stability parameter zu/L PSL, unitless'};
l(ismember(the_fields,'zet_s'))        = {'Monin-Obukhov stability parameter zu/L ship, unitless'};
l(ismember(the_fields,'dT_skin'))      = {'T difference skin and snake PSL, C'};
l(ismember(the_fields,'dT_skin_s'))    = {'T difference skin and ship TSG, C'};
l(ismember(the_fields,'dq_skin'))      = {'q difference skin and snake PSL, g/kg'};
l(ismember(the_fields,'dq_skin_s'))    = {'q difference skin and ship TSG, g/kg'};
l(ismember(the_fields,'dz_skin'))          = {'cool skin thickness to PSL snake, m'};
l(ismember(the_fields,'dz_skin_s'))        = {'cool skin thickness to ship TSG, m'};
l(ismember(the_fields,'hrain'))           = {'rain heat flux PSL, W/m2'};
l(ismember(the_fields,'hrain_s'))         = {'rain heat flux ship, W/m2'};
l(ismember(the_fields,'CdN10'))       = {'drag coefficient neutral stability 10 m height, momentum transfer coefficient PSL'};
l(ismember(the_fields,'CdN10_s'))     = {'drag coefficient neutral stability 10 m height, momentum transfer coefficient ship'};
l(ismember(the_fields,'ChN10'))       = {'Stanton number neutral stability 10 m height, sensible heat transfer coefficient PSL'};
l(ismember(the_fields,'ChN10_s'))     = {'Stanton number neutral stability 10 m height, sensible heat transfer coefficient PSL'};
l(ismember(the_fields,'CeN10'))       = {'Dalton number neutral stability 10 m height, latent heat transfer coefficient PSL'};
l(ismember(the_fields,'CeN10_s'))     = {'Dalton number neutral stability 10 m height, latent heat transfer coefficient PSL'};
l(ismember(the_fields,'U10'))          = {'wind speed adjusted to 10 m PSL, m/s'};
l(ismember(the_fields,'U10_s'))        = {'wind speed adjusted to 10 m ship, m/s'};
l(ismember(the_fields,'erate'))         = {'Evaporation rate PSL, mm/hr'};
l(ismember(the_fields,'erate_s'))       = {'Evaporation rate ship mm/hr'};
l(ismember(the_fields,'T10'))          = {'Ta adjusted to 10 m PSL, m/s'};
l(ismember(the_fields,'T10_s'))        = {'Ta adjusted to 10 m ship, m/s'};
l(ismember(the_fields,'q10'))          = {'qa adjusted to 10 m PSL, g/kg'};
l(ismember(the_fields,'q10_s'))        = {'qa adjusted to 10 m ship, g/kg'};
l(ismember(the_fields,'rh10'))         = {'rh adjusted to 10 m PSL, %'};
l(ismember(the_fields,'rh10_s'))       = {'rh adjusted to 10 m ship, %'};
l(ismember(the_fields,'dz_warm'))      = {'warm layer thickness, m'};
l(ismember(the_fields,'dT_warm'))      = {'warming across entire warm layer, C'};
l(ismember(the_fields,'dsea_warm'))    = {'heating due to warm layer at selected depth of TSG, C'};
l(ismember(the_fields,'usid_a'))       = {'ustar from inertial dissipation method a, m/s'};
l(ismember(the_fields,'usid_b'))       = {'ustar from inertial dissipation method a, m/s'};
l(ismember(the_fields,'qsid_a'))       = {'qstar from inertial dissipation method a, g/kg'};
l(ismember(the_fields,'qsid_b'))       = {'qstar from inertial dissipation method a, g/kg'};
l(ismember(the_fields,'tsid_a'))       = {'Tstar from inertial dissipation method a, K'};
l(ismember(the_fields,'tsid_b'))       = {'Tstar from inertial dissipation method a, K'};
l(ismember(the_fields,'Lid_a'))        = {'L output zu over zet inertial dissipation method a, m/s'};
l(ismember(the_fields,'Lid_b'))        = {'L output zu over zet inertial dissipation method a, m/s'};
l(ismember(the_fields,'hs_cov'))       = {'direct covariance sensible heat flux, W/m2'};
l(ismember(the_fields,'hs_cov_sds'))   = {'direct covariance sensible heat flux (Simon de Szoeke method), W/m2'};
l(ismember(the_fields,'hl_cov'))       = {'direct covariance latent heat flux, W/m2'};
l(ismember(the_fields,'hl_covW'))      = {'direct covariance latent heat flux (Webb method), W/m2 '};
l(ismember(the_fields,'hl_cov_sds'))   = {'direct covariance latent heat flux (Simon de Szoeke method), W/m2'};
l(ismember(the_fields,'wu_cov'))       = {'w prime u prime bar covariance, m2/s2'};
l(ismember(the_fields,'wv_cov'))       = {'w prime v prime bar covariance, m2/s2'};
l(ismember(the_fields,'wq_cov'))       = {'w prime q prime bar covariance , m/s'};
l(ismember(the_fields,'wT_cov'))       = {'w prime T prime bar covariance (Simon de Szoeke method), m K/s'};
l(ismember(the_fields,'wq_sds'))       = {'w prime q prime bar covariance (Simon de Szoeke method), m/s'};
l(ismember(the_fields,'wT_sds'))       = {'w prime T prime bar covariance (Simon de Szoeke method), m K/s'};
l(ismember(the_fields,'upat_std'))     = {'platform velocity std deviations stern-to-bow direction for u, m/s'};
l(ismember(the_fields,'vplat_std'))    = {'platform velocity std deviations port-to-starboard for v, m/s'};
l(ismember(the_fields,'wplat_std'))    = {'platform velocity std deviations z-direction for w, m/s'};
l(ismember(the_fields,'tilt'))         = {'tilt angle from first rotation, radians'};
l(ismember(the_fields,'azm'))          = {'azimuth angle from first rotation, radians'};
l(ismember(the_fields,'missing'))      = {'sonic anemometer missing points'};
l(ismember(the_fields,'badPnts'))      = {'sonic anemometer bad points'};
l(ismember(the_fields,'uvar'))         = {'variance of u (stern to bow) from sonic anemometer, m2/s2'};
l(ismember(the_fields,'vvar'))         = {'variance of v (port to starboard) from sonic anemometer, m2/s2'};
l(ismember(the_fields,'wvar'))         = {'variance of w (stern to bow) from sonic anemometer, m2/s2'};
l(ismember(the_fields,'Tvar'))         = {'variance of Tv from sonic anemometer, C2'};
l(ismember(the_fields,'Ubar'))         = {'streamwise mean flow, rotated into flow so no cross flow, m/s'};
l(ismember(the_fields,'Tson_noise'))   = {'noise of sonic anemometer for Tv, C'};
l(ismember(the_fields,'shed_std_s'))   = {'standard deviation of sine of heading from ship, radians'};
l(ismember(the_fields,'ched_std_s'))   = {'standard deviation of cosine of heading from ship, radians'};
l(ismember(the_fields,'Cu2a'))         = {'Cu method a structure function parameter'};
l(ismember(the_fields,'Cu2b'))         = {'Cu method b structure function parameter'};
l(ismember(the_fields,'Cw2a'))         = {'Cw method a structure function parameter'};
l(ismember(the_fields,'Cw2b'))         = {'Cw method b structure function parameter'};
l(ismember(the_fields,'Ct2a'))         = {'Ct method a structure function parameter'};
l(ismember(the_fields,'Ct2b'))         = {'Ct method b structure function parameter'};
l(ismember(the_fields,'Cq2a'))         = {'Cq method a structure function parameter'};
l(ismember(the_fields,'Cq2b'))         = {'Cq method b structure function parameter'};
l(ismember(the_fields,'cspd'))         = {'current speed at near-surface ADCP ship, m/s'};
l(ismember(the_fields,'cdir'))         = {'current direction (direction to) at near-surface ADCP ship, deg'};
l(ismember(the_fields,'wspd_sfc'))     = {'wind speed with respect to moving surface, m/s'};
l(ismember(the_fields,'wdir_sfc'))     = {'wind direction with espect to moving surface, deg'};
l(ismember(the_fields,'SST'))          = {'COARE calculated SST equal to Tsnk-dT_skin, C'};
l(ismember(the_fields,'year'))         = {'year'};
l(ismember(the_fields,'rdir_std'))     = {'PSL standard deviation of relative wind speed, deg'};
l(ismember(the_fields,'rdir_std_s'))     = {'ship standard deviation of relative wind speed, deg'};
l(ismember(the_fields,'Pu'))           = {'power spectrum u, s^-1'};
l(ismember(the_fields,'Pv'))           = {'power spectrum v, s^-1'};
l(ismember(the_fields,'Pw'))           = {'power spectrum w, s^-1'};
l(ismember(the_fields,'Pt'))           = {'power spectrum T, s^-1'};
l(ismember(the_fields,'Pq'))           = {'power spectrum q, s^-1'};
l(ismember(the_fields,'Cuw'))          = {'cospectrum uw, s^-1'};
l(ismember(the_fields,'Cvw'))          = {'cospectrum vw, s^-1'};
l(ismember(the_fields,'Ctw'))          = {'cospectrum Tw, s^-1'};
l(ismember(the_fields,'Cqw'))          = {'cospectrum qw, s^-1'};
l(ismember(the_fields,'mu'))           = {'mu decorr q, s'};
l(ismember(the_fields,'mu_q'))         = {'mu decorr, s'};
l(ismember(the_fields,'lag'))          = {'lag times sonic, s'};
l(ismember(the_fields,'lag'))          = {'lag times sonic, s'};
l(ismember(the_fields,'rUn'))          = {'PSL sonic relative streamwise wind, m/s'};
l(ismember(the_fields,'rUw'))          = {'PSL sonic relative cross stream wind, m/s'};
l(ismember(the_fields,'Ue'))          = {'PSL eastward positive true wind, m/s'};
l(ismember(the_fields,'Un'))          = {'PSL northward positive true wind, m/s'};
l(ismember(the_fields,'good_motion'))= {'flag 1 when motion conditions good for covariance'};
l(ismember(the_fields,'good_motion_licor'))= {'flag 1 when motion and licor conditions good for covariance'};
l(ismember(the_fields,'good_motion_id'))= {'flag 1 when motion conditions good for inertial dissipation'};
l(ismember(the_fields,'good_motion_id_licor'))= {'flag 1 when motion and licor conditions good for inertial dissipation'};
l(ismember(the_fields,'tau_cov'))      = {'covariance stress, N/m2'};
l(ismember(the_fields,'tau_ida'))      = {'inertial dissipation tau method a, N/m2'};
l(ismember(the_fields,'tau_idb'))      = {'inertial dissipation tau method b, N/m2'};
l(ismember(the_fields,'hs_ida'))       = {'inertial dissipation hs method a, W/m2'};
l(ismember(the_fields,'hb'))           = {'buoyancy flux, W/m2'};
l(ismember(the_fields,'hb_son'))       = {'buoyancy flux sonic anemometer, W/m2'};
l(ismember(the_fields,'lw_net'))       = {'net longwave flux, W/m2'};
l(ismember(the_fields,'lw_up'))        = {'upward longwave flux, W/m2'};
l(ismember(the_fields,'rhoa10'))       = {'air density adjusted to 10 m, kg/m3'};
l(ismember(the_fields,'sw_net'))       = {'net shortwave flux, W/m2'};
l(ismember(the_fields,'sw_up'))        = {'upward shortwave flux, W/m2'};
l(ismember(the_fields,'zo'))           = {'roughness length, m'};
l(ismember(the_fields,'dT_warm'))      = {'warm layer T difference top to base, deg C'};
l(ismember(the_fields,'dz_warm'))      = {'warm layer thickness or depth, m'};
l(ismember(the_fields,'dT_warm_at_input'))= {'warm layer T difference at input level, deg C'};
l(ismember(the_fields,'Le'))           = {'latent heat of vaporization, J/K'};
l(ismember(the_fields,'Tref'))         = {'Tair adjusted to reference height, deg C'};
l(ismember(the_fields,'U10N'))         = {'wind speed adjusted to neutral at 10 m, m/s'};
l(ismember(the_fields,'UN'))           = {'wind speed adjusted to neutral at measurement height, m/s'};
l(ismember(the_fields,'UrefN'))        = {'wind speed adjusted to neutral at reference height, m/s'};
l(ismember(the_fields,'Uref'))         = {'wind speed adjusted to reference height, m/s'};
l(ismember(the_fields,'rhref'))        = {'rh adjusted to reference height, %'};
l(ismember(the_fields,'qref'))         = {'qair adjusted to reference height, g/kg'};
l(ismember(the_fields,'P10'))          = {'air pressure adjusted to 10 m, mb'};
l(ismember(the_fields,'gust'))         = {'gustiness factor, m/s'};
l(ismember(the_fields,'hb_s'))           = {'ship buoyancy flux, W/m2'};
l(ismember(the_fields,'hb_son_s'))       = {'ship buoyancy flux sonic anemometer, W/m2'};
l(ismember(the_fields,'lw_net_s'))       = {'ship net longwave flux, W/m2'};
l(ismember(the_fields,'lw_up_s'))        = {'ship upward longwave flux, W/m2'};
l(ismember(the_fields,'rhoa10_s'))       = {'ship air density adjusted to 10 m, kg/m3'};
l(ismember(the_fields,'sw_net_s'))       = {'ship net shortwave flux, W/m2'};
l(ismember(the_fields,'sw_up_s'))        = {'ship upward shortwave flux, W/m2'};
l(ismember(the_fields,'zo_s'))           = {'ship roughness length, m'};
l(ismember(the_fields,'dT_warm_s'))      = {'ship warm layer T difference top to base, deg C'};
l(ismember(the_fields,'dz_warm_s'))      = {'ship warm layer thickness or depth, m'};
l(ismember(the_fields,'dT_warm_at_input_s'))= {'ship warm layer T difference at input level, deg C'};
l(ismember(the_fields,'Le_s'))           = {'ship latent heat of vaporization, J/K'};
l(ismember(the_fields,'Tref_s'))         = {'ship Tair adjusted to reference height, deg C'};
l(ismember(the_fields,'U10N_s'))         = {'ship wind speed adjusted to neutral at 10 m, m/s'};
l(ismember(the_fields,'UN_s'))           = {'ship wind speed adjusted to neutral at measurement height, m/s'};
l(ismember(the_fields,'UrefN_s'))        = {'ship wind speed adjusted to neutral at reference height, m/s'};
l(ismember(the_fields,'Uref_s'))        = {'ship wind speed adjusted to reference height, m/s'};
l(ismember(the_fields,'rhref_s'))        = {'ship rh adjusted to reference height, %'};
l(ismember(the_fields,'qref_s'))         = {'ship qair adjusted to reference height, g/kg'};
l(ismember(the_fields,'P10_s'))          = {'ship air pressure adjusted to 10 m, mb'};
l(ismember(the_fields,'gust_s'))         = {'ship gustiness factor, m/s'};
%% units
u(ismember(the_fields,'ID_Ftsg_s')) = {'flag'};
u(ismember(the_fields,'ID_Ftsg1_s'))= {'flag'};
u(ismember(the_fields,'ID_Ftsg2_s'))= {'flag'};
u(ismember(the_fields,'ID_tsg_s'))  = {'flag'};
u(ismember(the_fields,'Ctsg_s'))    = {'S/m'};
u(ismember(the_fields,'Ctsg1_s'))   = {'S/m'};
u(ismember(the_fields,'Ctsg2_s'))   = {'S/m'};
u(ismember(the_fields,'Ftsg_s'))    = {'L/min'};
u(ismember(the_fields,'Ftsg1_s'))   = {'L/min'};
u(ismember(the_fields,'Ftsg2_s'))   = {'L/min'};
u(ismember(the_fields,'FLtsg_s'))   = {'micrograms/L'};
u(ismember(the_fields,'FLtsg1_s'))  = {'micrograms/L'};
u(ismember(the_fields,'FLtsg2_s'))  = {'micrograms/L'};
u(ismember(the_fields,'Stsg_s'))    = {'psu'};
u(ismember(the_fields,'Stsg1_s'))   = {'psu'};
u(ismember(the_fields,'Stsg2_s'))   = {'psu'};
u(ismember(the_fields,'Ta'))        = {'deg_C'};
u(ismember(the_fields,'Ta_s'))      = {'deg_C'};
u(ismember(the_fields,'Td'))        = {'deg_C'};
u(ismember(the_fields,'Td_s'))      = {'deg_C'};
u(ismember(the_fields,'Tskin'))     = {'deg_C'};
u(ismember(the_fields,'Tsnk'))      = {'deg_C'};
u(ismember(the_fields,'Tsonic'))    = {'deg_C'};
u(ismember(the_fields,'Ttsg_int_s'))= {'deg_C'};
u(ismember(the_fields,'Ttsg_s'))    = {'deg_C'};
u(ismember(the_fields,'Ttsg1_s'))   = {'deg_C'};
u(ismember(the_fields,'Ttsg2_s'))   = {'deg_C'};
u(ismember(the_fields,'aspir_TRH')) = {'rpm'};
u(ismember(the_fields,'cog'))       = {'deg_azimuth'};
u(ismember(the_fields,'cog_s'))     = {'deg_azimuth'};
u(ismember(the_fields,'day'))       = {'day'};
u(ismember(the_fields,'heave'))     = {'m'};
u(ismember(the_fields,'heave_s'))   = {'m'};
u(ismember(the_fields,'hed'))       = {'deg_azimuth'};
u(ismember(the_fields,'hed_s'))     = {'deg_azimuth'};
u(ismember(the_fields,'hedE'))      = {'deg_azimuth'};
u(ismember(the_fields,'hedN'))      = {'deg_azimuth'};
u(ismember(the_fields,'hedE_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'hedN_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'hour'))      = {'hour'};
u(ismember(the_fields,'jd'))        = {'decimal_date'};
u(ismember(the_fields,'lat'))       = {'deg_lat'};
u(ismember(the_fields,'lat_s'))     = {'deg_lat'};
u(ismember(the_fields,'lon'))       = {'deg_lon'};
u(ismember(the_fields,'lon_s'))     = {'deg_lon'};
u(ismember(the_fields,'lw_dn_1'))   = {'W/m2'};
u(ismember(the_fields,'lw_dn_2'))   = {'W/m2'};
u(ismember(the_fields,'lw_dn_clr')) = {'W/m2'};
u(ismember(the_fields,'lw_dn_s'))   = {' W/m2'};
u(ismember(the_fields,'lw_T_case_1')) = {'deg_C'};
u(ismember(the_fields,'lw_T_dome_1')) = {'deg_C'};
u(ismember(the_fields,'lw_therm_1')) = {'W/m2'};
u(ismember(the_fields,'lw_T_case_2')) = {'deg_C'};
u(ismember(the_fields,'lw_T_dome_2')) = {'deg_C'};
u(ismember(the_fields,'lw_therm_2')) = {'W/m2'};
u(ismember(the_fields,'lw_T_case_s')) = {'deg_C'};
u(ismember(the_fields,'lw_T_dome_s')) = {'deg_C'};
u(ismember(the_fields,'lw_therm_s')) = {'W/m2'};
u(ismember(the_fields,'minute'))    = {'min'};
u(ismember(the_fields,'month'))     = {'month'};
u(ismember(the_fields,'pitch'))     = {'deg_from_level'};
u(ismember(the_fields,'pitch_s'))   = {'deg_from_level'};
u(ismember(the_fields,'prate'))     = {'mm/hour'};
u(ismember(the_fields,'prate_s'))   = {'mm/hour'};
u(ismember(the_fields,'paccum'))    = {'mm'};
u(ismember(the_fields,'paccum_s'))  = {'mm'};
u(ismember(the_fields,'P'))         = {'mb'};
u(ismember(the_fields,'P_s'))       = {'mb'};
u(ismember(the_fields,'qa'))        = {'g/kg'};
u(ismember(the_fields,'qa_s'))      = {'g/kg'};
u(ismember(the_fields,'qsnk'))      = {'g/kg'};
u(ismember(the_fields,'raw_p'))     = {'V'};
u(ismember(the_fields,'raw_p_despike')) = {'V'};
u(ismember(the_fields,'rdir'))      = {'deg_azimuth'};
u(ismember(the_fields,'rdir'))      = {'deg_azimuth'};
u(ismember(the_fields,'rdir_f_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'rdir_p_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'rdir_s_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'rdir_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'rh'))        = {'%'};
u(ismember(the_fields,'rh_s'))      = {'%'};
u(ismember(the_fields,'roll'))      = {'deg_from_level'};
u(ismember(the_fields,'roll_s'))    = {'deg_from_level'};
u(ismember(the_fields,'rspd'))      = {'m/s'};
u(ismember(the_fields,'rspd_raw'))  = {'m/s'};
u(ismember(the_fields,'rspd_f_s'))  = {'m/s'};
u(ismember(the_fields,'rspd_p_s'))  = {'m/s'};
u(ismember(the_fields,'rspd_s_s'))  = {'m/s'};
u(ismember(the_fields,'rspd_s'))    = {'m/s'};
u(ismember(the_fields,'rspd_std_s'))= {'m/s'};
u(ismember(the_fields,'slp'))       = {'mb'};
u(ismember(the_fields,'slp_reported_s')) = {'mb'};
u(ismember(the_fields,'slp_s'))     = {'mb'};
u(ismember(the_fields,'sog'))       = {'m/s'};
u(ismember(the_fields,'sogE'))      = {'m/s'};
u(ismember(the_fields,'sogE_s'))    = {'m/s'};
u(ismember(the_fields,'sogN'))      = {'m/s'};
u(ismember(the_fields,'sogN_s'))    = {'m/s'};
u(ismember(the_fields,'sog_s'))     = {'m/s'};
u(ismember(the_fields,'spd_log_u_s')) = {'m/s'};
u(ismember(the_fields,'spd_log_v_s')) = {'m/s'};
u(ismember(the_fields,'spd_log_s'))  = {'m/s'};
u(ismember(the_fields,'spd_log_std_s'))  = {'m/s'};
u(ismember(the_fields,'sw_dn_2'))   = {'W/m2'};
u(ismember(the_fields,'sw_dn_1'))   = {'W/m2'};
u(ismember(the_fields,'sw_dn_clr')) = {'W/m2'};
u(ismember(the_fields,'sw_dn_s'))   = {'W/m2'};
u(ismember(the_fields,'t'))         = {'matlab'};
u(ismember(the_fields,'wdir'))      = {'deg_azimuth'};
u(ismember(the_fields,'wdir_raw'))  = {'deg_azimuth'};
u(ismember(the_fields,'wdir_f_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'wdir_p_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'wdir_s_s'))  = {'deg_azimuth'};
u(ismember(the_fields,'wdir_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'wspd'))      = {'m/s'};
u(ismember(the_fields,'wspd_raw'))  = {'m/s'};
u(ismember(the_fields,'wspd_f_s'))  = {'m/s'};
u(ismember(the_fields,'wspd_p_s'))  = {'m/s'};
u(ismember(the_fields,'wspd_s_s'))  = {'m/s'};
u(ismember(the_fields,'wspd_s'))    = {'m/s'};
u(ismember(the_fields,'year'))      = {'year'};
u(ismember(the_fields,'Tsnk_best'))    = {'deg_C'};
u(ismember(the_fields,'Ta_best'))      = {'deg_C'};
u(ismember(the_fields,'qa_best'))      = {'g/kg'};
u(ismember(the_fields,'rh_best'))      = {'%'};
u(ismember(the_fields,'lw_dn_best'))   = {'W/m2'};
u(ismember(the_fields,'sw_dn_best'))   = {'deg_C'};
u(ismember(the_fields,'slp_best'))     = {'mb'};
u(ismember(the_fields,'wspd_best'))    = {'m/s'};
u(ismember(the_fields,'rspd_best'))    = {'m/s'};
u(ismember(the_fields,'sog_std_s'))    = {'m/s'};
u(ismember(the_fields,'sog_min_s'))    = {'m/s'};
u(ismember(the_fields,'sog_max_s'))    = {'m/s'};
u(ismember(the_fields,'cog_std_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'cog_min_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'cog_max_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'hed_std_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'hed_min_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'hed_max_s'))    = {'deg_azimuth'};
u(ismember(the_fields,'cog_std'))      = {'deg_azimuth'};
u(ismember(the_fields,'sog_std'))      = {'m/s'};
u(ismember(the_fields,'hed_std'))      = {'deg_azimuth'};
u(ismember(the_fields,'cog_min'))      = {'deg_azimuth'};
u(ismember(the_fields,'sog_min'))      = {'m/s'};
u(ismember(the_fields,'hed_min'))      = {'deg_azimuth'};
u(ismember(the_fields,'cog_max'))      = {'deg_azimuth'};
u(ismember(the_fields,'sog_max'))      = {'m/s'};
u(ismember(the_fields,'hed_max'))      = {'deg_azimuth'};
u(ismember(the_fields,'shed_std'))      = {'radians'};
u(ismember(the_fields,'ched_std'))      = {'radians'};
u(ismember(the_fields,'hnet'))         = {'W/m2'};
u(ismember(the_fields,'lw_net'))       = {'W/m2'};
u(ismember(the_fields,'sw_net'))       = {'W/m2'};
u(ismember(the_fields,'dT'))           = {'deg_C'};
u(ismember(the_fields,'dq'))           = {'g/kg'};
u(ismember(the_fields,'qs'))          = {'g/kg'};
u(ismember(the_fields,'qs_s'))          = {'g/kg'};
u(ismember(the_fields,'rhoa'))         = {'kg/m3'};
u(ismember(the_fields,'rhoa_s'))       = {'kg/m3'};
u(ismember(the_fields,'H2O'))          = {'mmol/m3'};
u(ismember(the_fields,'H2O_s'))        = {'mmol/m3'};
u(ismember(the_fields,'licor_AGC'))    = {'deg_C'};
u(ismember(the_fields,'licor_rhoa'))   = {'kg/m3'};
u(ismember(the_fields,'licor_qa'))     = {'g/kg'};
u(ismember(the_fields,'licor_H2O'))    = {'mmol/m3'};
u(ismember(the_fields,'licor_H2O_MR')) = {'mol/mol'};
u(ismember(the_fields,'licor_Tbox'))   = {'deg_C'};
u(ismember(the_fields,'licor_Pbox'))   = {'kPa'};
u(ismember(the_fields,'licor_H2O_std'))= {'mmol/m3'};
u(ismember(the_fields,'licor_Pbox_std'))= {'kPa'};
u(ismember(the_fields,'rUn_wxt'))      = {'m/s'};
u(ismember(the_fields,'rUw_wxt'))      = {'m/s'};
u(ismember(the_fields,'wspd_wxt'))      = {'m/s'};
u(ismember(the_fields,'wdir_wxt'))      = {'deg_azimuth'};
u(ismember(the_fields,'rspd_wxt'))      = {'m/s'};
u(ismember(the_fields,'rdir_wxt'))      = {'deg_azimuth'};
u(ismember(the_fields,'Ta_wxt'))       = {'deg_C'};
u(ismember(the_fields,'qa_wxt'))       = {'g/kg'};
u(ismember(the_fields,'rh_wxt'))       = {'m/s'};
u(ismember(the_fields,'rhoa_wxt'))     = {'m/s'};
u(ismember(the_fields,'H2O_wxt'))      = {'mmol/m3'};
u(ismember(the_fields,'prate_wxt'))    = {'mm/hr'};
u(ismember(the_fields,'paccum_wxt'))   = {'mm'};
u(ismember(the_fields,'slp_wxt'))      = {'mb'};
u(ismember(the_fields,'P_wxt'))        = {'mb'};
u(ismember(the_fields,'S_0p85'))       = {'psu'};
u(ismember(the_fields,'ID_wh_otter'))  = {'flag'};
u(ismember(the_fields,'ID_modflow'))   = {'flag'};
u(ismember(the_fields,'ID_badflow'))   = {'flag'};
u(ismember(the_fields,'Ptsg_s'))       = {'dbar'};
u(ismember(the_fields,'Sigtsg_s'))     = {'kg/m3'};
u(ismember(the_fields,'wspd_min_s'))   = {'m/s'};
u(ismember(the_fields,'wspd_max_s'))   = {'m/s'};
u(ismember(the_fields,'wspd_std_s'))   = {'m/s'};
u(ismember(the_fields,'wdir_min_s'))   = {'deg_azimuth'};
u(ismember(the_fields,'wdir_max_s'))   = {'deg_azimuth'};
u(ismember(the_fields,'wdir_std_s'))   = {'deg_azimuth'};
u(ismember(the_fields,'lw_dn_std'))    = {'W/m2'};
u(ismember(the_fields,'sw_dn_std'))    = {'W/m2'};
u(ismember(the_fields,'usr'))          = {'m/s'};
u(ismember(the_fields,'usr_s'))        = {'m/s'};
u(ismember(the_fields,'tau'))          = {'N/m2'};
u(ismember(the_fields,'tau_s'))        = {'N/m2'};
u(ismember(the_fields,'hs'))          = {'W/m2'};
u(ismember(the_fields,'hs_s'))        = {'W/m2'};
u(ismember(the_fields,'hl'))          = {'W/m2'};
u(ismember(the_fields,'hl_s'))        = {'W/m2'};
u(ismember(the_fields,'hlwebb'))       = {'W/m2'};
u(ismember(the_fields,'hlwebb_s'))     = {'W/m2'};
u(ismember(the_fields,'tsr'))          = {'K'};
u(ismember(the_fields,'tsr_s'))        = {'K'};
u(ismember(the_fields,'qsr'))          = {'g/kg'};
u(ismember(the_fields,'qsr_s'))        = {'g/kg'};
u(ismember(the_fields,'zot'))          = {'m'};
u(ismember(the_fields,'zot_s'))        = {'m'};
u(ismember(the_fields,'zoq'))          = {'m'};
u(ismember(the_fields,'zoq_s'))        = {'m'};
u(ismember(the_fields,'Cd'))           = {'unitless'};
u(ismember(the_fields,'Cd_s'))         = {'unitless'};
u(ismember(the_fields,'Ch'))           = {'unitless'};
u(ismember(the_fields,'Ch_s'))         = {'unitless'};
u(ismember(the_fields,'Ce'))           = {'unitless'};
u(ismember(the_fields,'Ce_s'))         = {'unitless'};
u(ismember(the_fields,'L'))            = {'m'};
u(ismember(the_fields,'L_s'))          = {'m'};
u(ismember(the_fields,'zet'))          = {'unitless'};
u(ismember(the_fields,'zet_s'))        = {'unitless'};
u(ismember(the_fields,'dT_skin'))      = {'deg_C'};
u(ismember(the_fields,'dT_skin_s'))    = {'deg_C'};
u(ismember(the_fields,'dq_skin'))      = {'g/kg'};
u(ismember(the_fields,'dq_skin_s'))    = {'g/kg'};
u(ismember(the_fields,'dz_skin'))          = {'m'};
u(ismember(the_fields,'dz_skin_s'))        = {'m'};
u(ismember(the_fields,'hrain'))           = {'W/m2'};
u(ismember(the_fields,'hrain_s'))         = {'W/m2'};
u(ismember(the_fields,'CdN10'))       = {'unitless'};
u(ismember(the_fields,'CdN10_s'))     = {'unitless'};
u(ismember(the_fields,'ChN10'))       = {'unitless'};
u(ismember(the_fields,'ChN10_s'))     = {'unitless'};
u(ismember(the_fields,'CeN10'))       = {'unitless'};
u(ismember(the_fields,'CeN10_s'))     = {'unitless'};
u(ismember(the_fields,'U10'))          = {'m/s'};
u(ismember(the_fields,'U10_s'))        = {'m/s'};
u(ismember(the_fields,'erate'))         = {'mm/hr'};
u(ismember(the_fields,'erate_s'))       = {'mm/hr'};
u(ismember(the_fields,'T10'))          = {'m/s'};
u(ismember(the_fields,'T10_s'))        = {'m/s'};
u(ismember(the_fields,'q10'))          = {'g/kg'};
u(ismember(the_fields,'q10_s'))        = {'g/kg'};
u(ismember(the_fields,'rh10'))         = {'%'};
u(ismember(the_fields,'rh10_s'))       = {'%'};
u(ismember(the_fields,'dz_warm'))      = {'m'};
u(ismember(the_fields,'dT_warm'))      = {'deg_C'};
u(ismember(the_fields,'dsea_warm'))    = {'deg_C'};
u(ismember(the_fields,'usid_a'))       = {'m/s'};
u(ismember(the_fields,'usid_b'))       = {'m/s'};
u(ismember(the_fields,'qsid_a'))       = {'g/kg'};
u(ismember(the_fields,'qsid_b'))       = {'g/kg'};
u(ismember(the_fields,'tsid_a'))       = {'K'};
u(ismember(the_fields,'tsid_b'))       = {'K'};
u(ismember(the_fields,'Lid_a'))        = {'m/s'};
u(ismember(the_fields,'Lid_b'))        = {'m/s'};
u(ismember(the_fields,'hs_cov'))       = {'W/m2'};
u(ismember(the_fields,'hs_cov_sds'))   = {'W/m2'};
u(ismember(the_fields,'hl_cov'))       = {'W/m2'};
u(ismember(the_fields,'hl_covW'))      = {'W/m2 '};
u(ismember(the_fields,'hl_cov_sds'))   = {'W/m2'};
u(ismember(the_fields,'wu_cov'))       = {'m2/s2'};
u(ismember(the_fields,'wv_cov'))       = {'m2/s2'};
u(ismember(the_fields,'wq_cov'))       = {'m/s'};
u(ismember(the_fields,'wT_cov'))       = {'K/s'};
u(ismember(the_fields,'wq_sds'))       = {'m/s'};
u(ismember(the_fields,'wT_sds'))       = {'K/s'};
u(ismember(the_fields,'upat_std'))     = {'m/s'};
u(ismember(the_fields,'vplat_std'))    = {'m/s'};
u(ismember(the_fields,'wplat_std'))    = {'m/s'};
u(ismember(the_fields,'tilt'))         = {'radians'};
u(ismember(the_fields,'azm'))          = {'radians'};
u(ismember(the_fields,'missing'))      = {'unitless'};
u(ismember(the_fields,'badPnts'))      = {'unitless'};
u(ismember(the_fields,'uvar'))         = {'m2/s2'};
u(ismember(the_fields,'vvar'))         = {'m2/s2'};
u(ismember(the_fields,'wvar'))         = {'m2/s2'};
u(ismember(the_fields,'Tvar'))         = {'C^2'};
u(ismember(the_fields,'Ubar'))         = {'m/s'};
u(ismember(the_fields,'Tson_noise'))   = {'deg_C'};
u(ismember(the_fields,'shed_std_s'))   = {'radians'};
u(ismember(the_fields,'ched_std_s'))   = {'radians'};
u(ismember(the_fields,'Cu2a'))         = {'unitless'};
u(ismember(the_fields,'Cu2b'))         = {'unitless'};
u(ismember(the_fields,'Cw2a'))         = {'unitless'};
u(ismember(the_fields,'Cw2b'))         = {'unitless'};
u(ismember(the_fields,'Ct2a'))         = {'unitless'};
u(ismember(the_fields,'Ct2b'))         = {'unitless'};
u(ismember(the_fields,'Cq2a'))         = {'unitless'};
u(ismember(the_fields,'Cq2b'))         = {'unitless'};
u(ismember(the_fields,'cspd'))         = {'m/s'};
u(ismember(the_fields,'cdir'))         = {'deg_azimuth'};
u(ismember(the_fields,'wspd_sfc'))     = {'m/s'};
u(ismember(the_fields,'wdir_sfc'))     = {'deg_azimuth'};
u(ismember(the_fields,'SST'))          = {'deg_C'};
u(ismember(the_fields,'rdir_std'))     = {'deg_azimuth'};
u(ismember(the_fields,'rdir_std_s'))     = {'deg_azimuth'};
u(ismember(the_fields,'Pu'))           = {'seconds^-1'};
u(ismember(the_fields,'Pv'))           = {'seconds^-1'};
u(ismember(the_fields,'Pw'))           = {'seconds^-1'};
u(ismember(the_fields,'Pt'))           = {'seconds^-1'};
u(ismember(the_fields,'Pq'))           = {'seconds^-1'};
u(ismember(the_fields,'Cuw'))          = {'seconds^-1'};
u(ismember(the_fields,'Cvw'))          = {'seconds^-1'};
u(ismember(the_fields,'Ctw'))          = {'seconds^-1'};
u(ismember(the_fields,'Cqw'))          = {'seconds^-1'};
u(ismember(the_fields,'mu'))           = {'seconds'};
u(ismember(the_fields,'mu_q'))         = {'seconds'};
u(ismember(the_fields,'lag'))          = {'seconds'};
u(ismember(the_fields,'rUn'))          = {'m/s'};
u(ismember(the_fields,'rUw'))          = {'m/s'};
u(ismember(the_fields,'Ue'))          = {'m/s'};
u(ismember(the_fields,'Un'))          = {'m/s'};
u(ismember(the_fields,'good_motion'))= {'flag'};
u(ismember(the_fields,'good_motion_licor'))= {'flag'};
u(ismember(the_fields,'good_motion_id'))= {'flag'};
u(ismember(the_fields,'good_motion_id_licor'))= {'flag'};
u(ismember(the_fields,'kk'))= {'flag'};
u(ismember(the_fields,'good_motion_licor'))= {'flag'};
u(ismember(the_fields,'good_motion_id'))= {'flag'};
u(ismember(the_fields,'good_motion_id_licor'))= {'flag'};
u(ismember(the_fields,'good_motion'))= {'flag'};
u(ismember(the_fields,'good_motion_licor'))= {'flag'};
u(ismember(the_fields,'good_motion_id'))= {'flag'};
u(ismember(the_fields,'good_motion_id_licor'))= {'flag'};
u(ismember(the_fields,'tau_cov'))      = {'N/m2'};
u(ismember(the_fields,'tau_ida'))      = {'N/m2'};
u(ismember(the_fields,'tau_idb'))      = {'N/m2'};
u(ismember(the_fields,'hs_ida'))       = {'W/m2'};
u(ismember(the_fields,'hs_idb'))       = {'W/m2'};
u(ismember(the_fields,'hl_ida'))       = {'W/m2'};
u(ismember(the_fields,'hl_idb'))       = {'W/m2'};
u(ismember(the_fields,'hb'))           = {'W/m2'};
u(ismember(the_fields,'hb_son'))       = {'W/m2'};
u(ismember(the_fields,'lw_net'))       = {'W/m2'};
u(ismember(the_fields,'lw_up'))        = {'W/m2'};
u(ismember(the_fields,'rhoa10'))       = {'kg/m3'};
u(ismember(the_fields,'sw_net'))       = {'W/m2'};
u(ismember(the_fields,'sw_up'))        = {'W/m2'};
u(ismember(the_fields,'zo'))           = {'m'};
u(ismember(the_fields,'dT_warm'))      = {'deg C'};
u(ismember(the_fields,'dz_warm'))      = {'m'};
u(ismember(the_fields,'dT_warm_at_input'))= {'deg C'};
u(ismember(the_fields,'Le'))           = {'J/K'};
u(ismember(the_fields,'Tref'))         = {'deg C'};
u(ismember(the_fields,'U10N'))         = {'m/s'};
u(ismember(the_fields,'UN'))           = {'m/s'};
u(ismember(the_fields,'UrefN'))        = {'m/s'};
u(ismember(the_fields,'Uref'))        = {'m/s'};
u(ismember(the_fields,'rhref'))        = {'%'};
u(ismember(the_fields,'qref'))         = {'g/kg'};
u(ismember(the_fields,'P10'))          = {'mb'};
u(ismember(the_fields,'gust'))         = {'m/s'};
u(ismember(the_fields,'hb_s'))           = {'W/m2'};
u(ismember(the_fields,'hb_son_s'))       = {'W/m2'};
u(ismember(the_fields,'lw_net_s'))       = {'W/m2'};
u(ismember(the_fields,'lw_up_s'))        = {'W/m2'};
u(ismember(the_fields,'rhoa10_s'))       = {'kg/m3'};
u(ismember(the_fields,'sw_net_s'))       = {'W/m2'};
u(ismember(the_fields,'sw_up_s'))        = {'W/m2'};
u(ismember(the_fields,'zo_s'))           = {'m'};
u(ismember(the_fields,'dT_warm_s'))      = {'deg C'};
u(ismember(the_fields,'dz_warm_s'))      = {'m'};
u(ismember(the_fields,'dT_warm_at_input_s'))= {'deg C'};
u(ismember(the_fields,'Le_s'))           = {'J/K'};
u(ismember(the_fields,'Tref_s'))         = {'deg C'};
u(ismember(the_fields,'U10N_s'))         = {'m/s'};
u(ismember(the_fields,'UN_s'))           = {'m/s'};
u(ismember(the_fields,'UrefN_s'))        = {'m/s'};
u(ismember(the_fields,'Uref_s'))        = {'m/s'};
u(ismember(the_fields,'rhref_s'))        = {'%'};
u(ismember(the_fields,'qref_s'))         = {'g/kg'};
u(ismember(the_fields,'P10_s'))          = {'mb'};
u(ismember(the_fields,'gust_s'))         = {'m/s'};

%% long_name

ln(ismember(the_fields,'SST'))       = {'SST_equal_to_Tsnk-dT_skin'};
ln(ismember(the_fields,'ID_Ftsg_s')) = {'seawater_flow_rate_ID_for_TSG'};
ln(ismember(the_fields,'ID_Ftsg1_s'))= {'seawater_flow_rate_ID_for_bow_TSG'};
ln(ismember(the_fields,'ID_Ftsg2_s'))= {'seawater_flow_rate_ID_for_seachest_TSG'};
ln(ismember(the_fields,'ID_tsg_s'))  = {'see_info_for_explaination'};
ln(ismember(the_fields,'Ctsg_s'))    = {'ship_thermosalinograph_conductivity'};
ln(ismember(the_fields,'Ctsg1_s'))   = {'ship_bow_thermosalinograph_conductivity'};
ln(ismember(the_fields,'Ctsg2_s'))   = {'ship_seachest_thermosalinograph_conductivity'};
ln(ismember(the_fields,'Ftsg_s'))    = {'ship_thermosalinograph_flow_rate'};
ln(ismember(the_fields,'Ftsg1_s'))   = {'ship_bow_thermosalinograph_flow_rate'};
ln(ismember(the_fields,'Ftsg2_s'))   = {'ship_seachest_thermosalinograph_flow_rate'};
ln(ismember(the_fields,'FLtsg_s'))   = {'ship_thermosalinograph_fluorometer_dry_value'};
ln(ismember(the_fields,'FLtsg1_s'))  = {'ship_bow_thermosalinograph_fluorometer_dry_value'};
ln(ismember(the_fields,'FLtsg2_s'))  = {'ship_seachest_thermosalinograph_fluorometer_dry_value'};
ln(ismember(the_fields,'Stsg_s'))    = {'ship_thermosalinograph_salinity'};
ln(ismember(the_fields,'Stsg1_s'))   = {'ship_bow_thermosalinograph_salinity'};
ln(ismember(the_fields,'Stsg2_s'))   = {'ship_seachest_thermosalinograph_salinity'};
ln(ismember(the_fields,'Ta'))        = {'PSL_air_T'};
ln(ismember(the_fields,'Ta_s'))      = {'ship_air_T'};
ln(ismember(the_fields,'Td'))        = {'PSL_air_dewpoint_T'};
ln(ismember(the_fields,'Td_s'))      = {'ship_air_dewpoint_T'};
ln(ismember(the_fields,'Tskin'))     = {'skin_T_of_seawater_from_ROSR_infrared_radiometer'};
ln(ismember(the_fields,'Tsnk'))      = {'PSL_sea_snake_T'};
ln(ismember(the_fields,'Tsonic'))    = {'PSL_air_Tv_virtual_temp_3D_sonic_anemometer'};
ln(ismember(the_fields,'Ttsg_int_s'))= {'ship_thermosalinograph_internal_T'};
ln(ismember(the_fields,'Ttsg_s'))    = {'ship_thermosalinograph_T'};
ln(ismember(the_fields,'Ttsg1_s'))   = {'ship_bow_thermosalinograph_T'};
ln(ismember(the_fields,'Ttsg2_s'))   = {'ship_seachest_thermosalinograph_T'};
ln(ismember(the_fields,'aspir_TRH')) = {'PSL_fan_speed_aspirator_on_T_RH_Vaisala_sensor'};
ln(ismember(the_fields,'cog'))       = {'PSL_course_over_ground'};
ln(ismember(the_fields,'cog_s'))     = {'ship_course_over_ground'};
ln(ismember(the_fields,'day'))       = {'day'};
ln(ismember(the_fields,'heave'))     = {'PSL_heave'};
ln(ismember(the_fields,'heave_s'))   = {'ship_heave'};
ln(ismember(the_fields,'hed'))       = {'PSL_heading'};
ln(ismember(the_fields,'hed_s'))     = {'ship_heading'};
ln(ismember(the_fields,'hedE'))       = {'PSL_heading_eastward'};
ln(ismember(the_fields,'hedN'))       = {'PSL_heading_northward'};
ln(ismember(the_fields,'hedE_s'))     = {'ship_heading_eastward'};
ln(ismember(the_fields,'hedN_s'))     = {'ship_heading_northward'};
ln(ismember(the_fields,'hour'))      = {'hour'};
ln(ismember(the_fields,'jd'))        = {'decimal_julian_date'};
ln(ismember(the_fields,'lat'))       = {'PSL_lat'};
ln(ismember(the_fields,'lat_s'))     = {'ship_lat'};
ln(ismember(the_fields,'lon'))       = {'PSL_lon'};
ln(ismember(the_fields,'lon_s'))     = {'ship_lon'};
ln(ismember(the_fields,'lw_dn_1'))   = {'PSL_downwelling_longwave_radiation_unit_1'};
ln(ismember(the_fields,'lw_dn_2'))   = {'PSL_downwelling_longwave_radiation_unit_2'};
ln(ismember(the_fields,'lw_dn_clr')) = {'calculated_clear_sky_downwelling_longwave_radiation'};
ln(ismember(the_fields,'lw_dn_s'))   = {'ship_downwelling_longwave_radiation'};
ln(ismember(the_fields,'lw_T_case_1')) = {'PSL_T_of_lw_unit_1_case'};
ln(ismember(the_fields,'lw_T_dome_1')) = {'PSL_T_of_lw_unit_1_dome'};
ln(ismember(the_fields,'lw_therm_1')) = {'PSL_thermopile_of_lw_unit_1'};
ln(ismember(the_fields,'lw_T_case_2')) = {'PSL_T_of_lw_unit_2_case'};
ln(ismember(the_fields,'lw_T_dome_2')) = {'PSL_T_of_lw_unit_2_dome'};
ln(ismember(the_fields,'lw_therm_2')) = {'PSL_thermopile_of_lw_unit_2'};
ln(ismember(the_fields,'lw_T_case_s')) = {'ship_T_of_lw_unit_case'};
ln(ismember(the_fields,'lw_T_dome_s')) = {'ship_T_of_lw_unit_dome'};
ln(ismember(the_fields,'lw_therm_s')) = {'ship_thermopile_of_lw_unit'};
ln(ismember(the_fields,'minute'))    = {'minute'};
ln(ismember(the_fields,'month'))     = {'month'};
ln(ismember(the_fields,'pitch'))     = {'PSL_pitch'};
ln(ismember(the_fields,'pitch_s'))   = {'ship_pitch'};
ln(ismember(the_fields,'prate'))     = {'PSL_precipitation_rate'};
ln(ismember(the_fields,'prate_s'))   = {'ship_precipitation_rate'};
ln(ismember(the_fields,'paccum'))    = {'PSL_precipitation_accumulation'};
ln(ismember(the_fields,'paccum_s'))  = {'ship_precipitation_accumulation'};
ln(ismember(the_fields,'P'))         = {'PSL_barometric_pressure'};
ln(ismember(the_fields,'P_s'))       = {'ship_barometric_pressure'};
ln(ismember(the_fields,'qa'))        = {'PSL_air_specific_humidity'};
ln(ismember(the_fields,'qa_s'))      = {'ship_air_specific_humidity'};
ln(ismember(the_fields,'qsnk'))      = {'saturation_specific_humidity_calculated_with_sea_snake_value_Tsnk'};
ln(ismember(the_fields,'raw_p'))     = {'PSL_optical_rain_gauge_raw_uncorrected_measurement'};
ln(ismember(the_fields,'raw_p_despike')) = {'PSL_optical_rain_gauge_raw_uncorrected_measurement_despiked'};
ln(ismember(the_fields,'rdir'))      = {'PSL_relative_wind_direction_3D_sonic_anemometer'};
ln(ismember(the_fields,'rdir_raw'))  = {'PSL_relative_wind_direction_3D_sonic_anemometer'};
ln(ismember(the_fields,'rdir_f_s'))  = {'ship_relative_wind_direction_foremast_rotating_propellor'};
ln(ismember(the_fields,'rdir_p_s'))  = {'ship_relative_wind_direction_port_2D_sonic_anemometer'};
ln(ismember(the_fields,'rdir_s_s'))  = {'ship_relative_wind_direction_starboard_2D_sonic_anemometer'};
ln(ismember(the_fields,'rdir_s'))    = {'ship relative_wind_direction_2D sonic_anemometer'};
ln(ismember(the_fields,'rh'))        = {'PSL_relative_humidity'};
ln(ismember(the_fields,'rh_s'))      = {'ship_relative_humidity'};
ln(ismember(the_fields,'roll'))      = {'PSL_roll'};
ln(ismember(the_fields,'roll_s'))    = {'ship_roll'};
ln(ismember(the_fields,'rspd'))      = {'PSL_relative_wind_speed_3D_sonic_anemometer'};
ln(ismember(the_fields,'rspd_raw'))  = {'PSL_raw_relative_wind_speed_3D_sonic_anemometer'};
ln(ismember(the_fields,'rspd_f_s'))  = {'ship_relative_wind_speed_foremast_rotating_propellor'};
ln(ismember(the_fields,'rspd_p_s'))  = {'ship_relative_wind_speed_port_2D_sonic_anemometer'};
ln(ismember(the_fields,'rspd_s_s'))  = {'ship_relative_wind_speed_starboard_2D_sonic_anemometer'};
ln(ismember(the_fields,'rspd_s'))    = {'ship_relative_wind_speed_2D_sonic_anemometer'};
ln(ismember(the_fields,'rspd_std_s'))    = {'ship_standard_deviation_relative_wind_speed_2D_sonic_anemometer'};
ln(ismember(the_fields,'slp'))       = {'PSL_barometric_pressure_for_sea_level_altitude'};
ln(ismember(the_fields,'slp_reported_s')) = {'ship_reported_barometric_pressure_for_sea_level_altitude'};
ln(ismember(the_fields,'slp_s'))     = {'ship_corrected_barometric_pressure_for_sea_level_altitude'};
ln(ismember(the_fields,'sog'))       = {'PSL_speed_over_ground'};
ln(ismember(the_fields,'sogE'))      = {'PSL_eastward_component_of_speed_over_ground'};
ln(ismember(the_fields,'sogE_s'))    = {'ship_eastward_component_of_speed_over_ground'};
ln(ismember(the_fields,'sogN'))      = {'PSL_northward_component_of_speed_over_ground'};
ln(ismember(the_fields,'sogN_s'))    = {'ship_northward_component_of_speed_over_ground'};
ln(ismember(the_fields,'sog_s'))     = {'ship_speed_over_ground'};
ln(ismember(the_fields,'spd_log_u_s')) = {'ship_eastward_actual_speed_log'};
ln(ismember(the_fields,'spd_log_v_s')) = {'ship_northward_actual_speed_log'};
ln(ismember(the_fields,'spd_log_s')) = {'ship_actual_speed_log'};
ln(ismember(the_fields,'spd_log_std_s')) = {'ship_standard_deviation_actual_speed_log'};
ln(ismember(the_fields,'sw_dn_2'))   = {'PSL_downwelling_shortwave_radiation_from_unit_1'};
ln(ismember(the_fields,'sw_dn_1'))   = {'PSL_downwelling_shortwave_radiation_from_unit_2'};
ln(ismember(the_fields,'sw_dn_clr')) = {'calculated_clear_sky_downwelling_shortwave_radiation'};
ln(ismember(the_fields,'sw_dn_s'))   = {'ship_downwelling_shortwave_radiation'};
ln(ismember(the_fields,'t'))         = {'matlab_datetime'};
ln(ismember(the_fields,'wdir'))      = {'PSL_true_wind_direction_3D_sonic_anemometer'};
ln(ismember(the_fields,'wdir_raw'))  = {'PSL_true_wind_direction_3D_sonic_anemometer'};
ln(ismember(the_fields,'wdir_f_s'))  = {'ship_true_wind_direction_foremast_rotating_propellor'};
ln(ismember(the_fields,'wdir_p_s'))  = {'ship_true_wind_direction_port_2D_sonic_anemometer'};
ln(ismember(the_fields,'wdir_s_s'))  = {'ship_true_wind_direction_starboard_2D_sonic_anemometer'};
ln(ismember(the_fields,'wdir_s'))    = {'ship_true_wind_direction_2D_sonic_anemometer'};
ln(ismember(the_fields,'wspd'))      = {'PSL_wind_speed_3D_sonic_anemometer'};
ln(ismember(the_fields,'wspd_raw'))  = {'PSL_wind_speed_3D_sonic_anemometer'};
ln(ismember(the_fields,'wspd_f_s'))  = {'ship_true_wind_speed_foremast_rotating_propellor'};
ln(ismember(the_fields,'wspd_p_s'))  = {'ship_true_wind_speed_port_2D_sonic_anemometer'};
ln(ismember(the_fields,'wspd_s_s'))  = {'ship_true_wind_speed_starboard_2D_sonic_anemometer'};
ln(ismember(the_fields,'wspd_s'))    = {'ship_true_wind_speed_2D_sonic_anemometer'};
ln(ismember(the_fields,'year'))      = {'year'};
ln(ismember(the_fields,'Tsnk_best'))    = {'PSL_Tsnk_filled_ROSR_TSG'};
ln(ismember(the_fields,'Ta_best'))      = {'PSL_Ta_filled_ship'};
ln(ismember(the_fields,'qa_best'))      = {'PSL_qa_filled_ship'};
ln(ismember(the_fields,'rh_best'))      = {'PSL_rh_filled_ship'};
ln(ismember(the_fields,'lw_dn_best'))   = {'PSL_lw_dn_best'};
ln(ismember(the_fields,'sw_dn_best'))   = {'PSL_sw_dn_best'};
ln(ismember(the_fields,'slp_best'))     = {'PSL_slp_filled_ship'};
ln(ismember(the_fields,'wspd_best'))    = {'PSL_wspd_filled_ship'};
ln(ismember(the_fields,'rspd_best'))    = {'PSL_rspd_filled_ship'};
ln(ismember(the_fields,'sog_std_s'))    = {'ship_standard deviation_SOG'};
ln(ismember(the_fields,'sog_min_s'))    = {'ship_min_SOG'};
ln(ismember(the_fields,'sog_max_s'))    = {'ship_max_SOG'};
ln(ismember(the_fields,'cog_std_s'))    = {'ship_standard deviation_COG'};
ln(ismember(the_fields,'cog_min_s'))    = {'ship_min_COG'};
ln(ismember(the_fields,'cog_max_s'))    = {'ship_max_COG'};
ln(ismember(the_fields,'hed_std_s'))    = {'ship_standard_deviation_heading'};
ln(ismember(the_fields,'hed_min_s'))    = {'ship_min_heading'};
ln(ismember(the_fields,'hed_max_s'))    = {'ship_max_heading'};
ln(ismember(the_fields,'cog_std'))      = {'PSL_standard_deviation_of_COG'};
ln(ismember(the_fields,'sog_std'))      = {'PSL_standard_deviation_of_SOG'};
ln(ismember(the_fields,'hed_std'))      = {'PSL_standard_deviation_of_heading'};
ln(ismember(the_fields,'cog_min'))      = {'PSL_min_of_COG'};
ln(ismember(the_fields,'sog_min'))      = {'PSL_min_of_SOG'};
ln(ismember(the_fields,'hed_min'))      = {'PSL_min_of_heading'};
ln(ismember(the_fields,'cog_max'))      = {'PSL_max_of_COG'};
ln(ismember(the_fields,'sog_max'))      = {'PSL_max_of_SOG'};
ln(ismember(the_fields,'hed_max'))      = {'PSL_max_of_heading'};
ln(ismember(the_fields,'ched_std'))      = {'PSL_cosine_standard_deviation_of_heading'};
ln(ismember(the_fields,'shed_std'))      = {'PSL_sine_standard_deviation_of_heading'};
ln(ismember(the_fields,'hnet'))         = {'PSL_net_heat_flux_positive_warming_ocean'};
ln(ismember(the_fields,'lw_net'))       = {'PSL_net_long_wave_flux_positive_warming_ocean'};
ln(ismember(the_fields,'sw_net'))       = {'PSL_net_short_wave_flux_positive_warming_ocean'};
ln(ismember(the_fields,'dT'))           = {'PSL_SST-T10_air_sea_T_difference'};
ln(ismember(the_fields,'dq'))           = {'PSL_qs-q10_air_sea_q_difference'};
ln(ismember(the_fields,'qs'))          = {'PSL_sea surface specific humidity'};
ln(ismember(the_fields,'qs_s'))          = {'ship_sea surface specific humidity'};
ln(ismember(the_fields,'rhoa'))         = {'PSL_air_density'};
ln(ismember(the_fields,'rhoa_s'))       = {'ship_air_density'};
ln(ismember(the_fields,'H2O'))          = {'PSL_H2O_number_concentration_Vaisala'};
ln(ismember(the_fields,'H2O_s'))        = {'ship_H2O_number_concentration_met'};
ln(ismember(the_fields,'licor_AGC'))   = {'Licor_box_T_AGC_value'};
ln(ismember(the_fields,'licor_rhoa'))   = {'Licor_air_density'};
ln(ismember(the_fields,'licor_qa'))     = {'Licor_specific_humidity'};
ln(ismember(the_fields,'licor_H2O'))    = {'Licor_H2O_number_concentration'};
ln(ismember(the_fields,'licor_H2O_MR')) = {'Licor_H2O_mixing_ratio'};
ln(ismember(the_fields,'licor_Tbox'))   = {'Licor_box_T'};
ln(ismember(the_fields,'licor_Pbox'))   = {'Licor_box_P'};
ln(ismember(the_fields,'licor_H2O_std'))= {'Licor_H2O_standard_deviation'};
ln(ismember(the_fields,'licor_Pbox_std'))= {'Licor_Pbox_stanard_deviation'};
ln(ismember(the_fields,'rUn_wxt'))      = {'WXT_northward_positive_wind_speed'};
ln(ismember(the_fields,'rUw_wxt'))      = {'WXT_westward_positive_wind_speed'};
ln(ismember(the_fields,'wspd_wxt'))      = {'WXT_true_wind_speed'};
ln(ismember(the_fields,'wdir_wxt'))      = {'WXT_true_wind_direction'};
ln(ismember(the_fields,'rspd_wxt'))      = {'WXT_relatve_wind_speed'};
ln(ismember(the_fields,'rdir_wxt'))      = {'WXT_relative_wind_direction'};
ln(ismember(the_fields,'Ta_wxt'))       = {'WXT_air_T'};
ln(ismember(the_fields,'qa_wxt'))       = {'WXT_air_specific_humidity'};
ln(ismember(the_fields,'rh_wxt'))       = {'WXT_relative_humidity'};
ln(ismember(the_fields,'rhoa_wxt'))     = {'WXT_air_density'};
ln(ismember(the_fields,'H2O_wxt'))      = {'WXT_H2O_concentration'};
ln(ismember(the_fields,'prate_wxt'))    = {'WXT_precip_rate'};
ln(ismember(the_fields,'paccum_wxt'))   = {'WXT_precip_accumlulation'};
ln(ismember(the_fields,'slp_wxt'))      = {'WXT_sea_level_pressure'};
ln(ismember(the_fields,'P_wxt'))        = {'WXT_measured_pressure'};
ln(ismember(the_fields,'S_0p85'))       = {'salinity_OSU_Otter_0p85_m_depth'};
ln(ismember(the_fields,'ID_wh_otter'))  = {'OSU_otter_present'};
ln(ismember(the_fields,'ID_modflow'))   = {'poor_wind_direction_over_90_deg'};
ln(ismember(the_fields,'ID_badflow'))   = {'bad_wind_directions_over_120_deg'};
ln(ismember(the_fields,'Ptsg_s'))       = {'ship_seawater_pressure_at_TSG'};
ln(ismember(the_fields,'Sigtsg_s'))     = {'ship_seawater_density_at_TSG'};
ln(ismember(the_fields,'wspd_min_s'))   = {'ship_min_true_wind_speed'};
ln(ismember(the_fields,'wspd_max_s'))   = {'ship_max_true_wind_speed'};
ln(ismember(the_fields,'wspd_std_s'))   = {'ship_standard_deviation_true_wind_speed'};
ln(ismember(the_fields,'wdir_min_s'))   = {'ship_min_true_wind_direction'};
ln(ismember(the_fields,'wdir_max_s'))   = {'ship_max_true_wind_direction'};
ln(ismember(the_fields,'wdir_std_s'))   = {'ship_standard_deviation_true_wind_direction'};
ln(ismember(the_fields,'lw_dn_std'))    = {'PSL_standard_deviation_of_lw_dn_best'};
ln(ismember(the_fields,'sw_dn_std'))    = {'PSL_standard_deviation_of_sw_dn_best'};
ln(ismember(the_fields,'usr'))          = {'PSL_u_star'};
ln(ismember(the_fields,'usr_s'))        = {'ship_u_star'};
ln(ismember(the_fields,'tau'))          = {'PSL_bulk_wind_stress'};
ln(ismember(the_fields,'tau_s'))        = {'ship_bulk_wind_stress'};
ln(ismember(the_fields,'hs'))          = {'PSL_bulk_sensible_heat_flux'};
ln(ismember(the_fields,'hs_s'))        = {'ship_bulk_sensible_heat_flux'};
ln(ismember(the_fields,'hl'))          = {'PSL_bulk_latent_heat_flux'};
ln(ismember(the_fields,'hl_s'))        = {'ship_bulk_latent_heat_flux'};
ln(ismember(the_fields,'hlwebb'))       = {'PSL_Webb_correction_for_latent_heat_flux'};
ln(ismember(the_fields,'hlwebb_s'))     = {'ship_Webb_correction_for_latent_heat_flux'};
ln(ismember(the_fields,'tsr'))          = {'PSL_T_star_temperature_scaling_parameter'};
ln(ismember(the_fields,'tsr_s'))        = {'ship_T_star_temperature_scaling_parameter'};
ln(ismember(the_fields,'qsr'))          = {'PSL_q_star_specific_humidity_scaling_parameter'};
ln(ismember(the_fields,'qsr_s'))        = {'ship_q_star_specific_humidity_scaling_parameter'};
ln(ismember(the_fields,'zot'))          = {'PSL_z_over_T_roughness_length_with_respect_to_T'};
ln(ismember(the_fields,'zot_s'))        = {'ship_z_over_T_roughness_length_with_respect_to_T'};
ln(ismember(the_fields,'zoq'))          = {'PSL_z_over_q_roughness_length_with_respect_to_q'};
ln(ismember(the_fields,'zoq_s'))        = {'ship_z_over_q_roughness_length_with_respect_to_q'};
ln(ismember(the_fields,'Cd'))           = {'PSL_drag_coefficient_momentum_transfer_coefficient'};
ln(ismember(the_fields,'Cd_s'))         = {'ship_drag_coefficient_momentum_transfer_coefficient'};
ln(ismember(the_fields,'Ch'))           = {'PSL_Stanton_number_sensible_heat_transfer_coefficient'};
ln(ismember(the_fields,'Ch_s'))         = {'ship_Stanton_number_sensible_heat_transfer_coefficient'};
ln(ismember(the_fields,'Ce'))           = {'PSL_Dalton_number_latent_heat_transfer_coefficient'};
ln(ismember(the_fields,'Ce_s'))         = {'ship_Dalton_number_latent_heat_transfer_coefficient'};
ln(ismember(the_fields,'L'))            = {'PSL_Monin-Obukhov_length_scale'};
ln(ismember(the_fields,'L_s'))          = {'ship_Monin-Obukhov_length_scale'};
ln(ismember(the_fields,'zet'))          = {'PSL_Monin-Obukhov_stability_parameter_zu_over_L'};
ln(ismember(the_fields,'zet_s'))        = {'ship_Monin-Obukhov_stability_parameter_zu_over_L'};
ln(ismember(the_fields,'dT_skin'))      = {'PSL_T_difference_skin_to_snake'};
ln(ismember(the_fields,'dT_skin_s'))    = {'ship_T_difference_skin_to_TSG'};
ln(ismember(the_fields,'dq_skin'))      = {'PSL_q_difference_skin_to_snake'};
ln(ismember(the_fields,'dq_skin_s'))    = {'ship_q_difference_skin_to_TSG'};
ln(ismember(the_fields,'dz_skin'))          = {'PSL_cool_skin_thickness_to_snake'};
ln(ismember(the_fields,'dz_skin_s'))        = {'ship_cool_skin_thickness_to_TSG'};
ln(ismember(the_fields,'hrain'))           = {'PSL_rain_heat_flux'};
ln(ismember(the_fields,'hrain_s'))         = {'ship_rain_heat_flux'};
ln(ismember(the_fields,'CdN10'))       = {'PSL_drag_coefficient_neutral_stability_10m_height'};
ln(ismember(the_fields,'CdN10_s'))     = {'ship_drag_coefficient_neutral_stability_10m_height'};
ln(ismember(the_fields,'ChN10'))       = {'PSL_Stanton_number_neutral_stability_10m_height'};
ln(ismember(the_fields,'ChN10_s'))     = {'ship_Stanton_number_neutral_stability_10m_height'};
ln(ismember(the_fields,'CeN10'))       = {'PSL_Dalton_number_neutral_stability_10m_height'};
ln(ismember(the_fields,'CeN10_s'))     = {'ship_Dalton_number_neutral_stability_10m_height'};
ln(ismember(the_fields,'U10'))          = {'PSL_wind_speed_adjusted_10m'};
ln(ismember(the_fields,'U10_s'))        = {'ship_wind_speed_adjusted_10m'};
ln(ismember(the_fields,'erate'))         = {'PSL_evaporation_rate'};
ln(ismember(the_fields,'erate_s'))       = {'ship_evaporation_rate'};
ln(ismember(the_fields,'T10'))          = {'PSL_Ta_adjusted_10m'};
ln(ismember(the_fields,'T10_s'))        = {'ship_Ta_adjusted_10m'};
ln(ismember(the_fields,'q10'))          = {'PSL_qa_adjusted_10m'};
ln(ismember(the_fields,'q10_s'))        = {'ship_qa_adjusted_10m'};
ln(ismember(the_fields,'rh10'))         = {'PSL_rh_adjusted_10m'};
ln(ismember(the_fields,'rh10_s'))       = {'ship_rh_adjusted_10m'};
ln(ismember(the_fields,'dz_warm'))      = {'warm_layer_thickness'};
ln(ismember(the_fields,'dT_warm'))      = {'warming_across_entire_warm_layer'};
ln(ismember(the_fields,'dsea_warm'))    = {'heating_due_to_warm_layer_at_selected_depth_of_TSG'};
ln(ismember(the_fields,'usid_a'))       = {'ustar_inertial_dissipation_method_a'};
ln(ismember(the_fields,'usid_b'))       = {'ustar_inertial_dissipation_method_b'};
ln(ismember(the_fields,'qsid_a'))       = {'qstar_inertial_dissipation_method_a'};
ln(ismember(the_fields,'qsid_b'))       = {'qstar_inertial_dissipation_method_b'};
ln(ismember(the_fields,'tsid_a'))       = {'Tstar_inertial_dissipation_method_a'};
ln(ismember(the_fields,'tsid_b'))       = {'Tstar_inertial_dissipation_method_b'};
ln(ismember(the_fields,'Lid_a'))        = {'L_output_inertial_dissipation_method_a'};
ln(ismember(the_fields,'Lid_b'))        = {'L_output_inertial_dissipation_method_b'};
ln(ismember(the_fields,'hs_cov'))       = {'direct_covariance_sensible_heat_flux'};
ln(ismember(the_fields,'hs_cov_sds'))   = {'direct_covariance sensible_heat_flux_SimondeSzoeke_method'};
ln(ismember(the_fields,'hl_cov'))       = {'direct_covariance_latent_heat flux'};
ln(ismember(the_fields,'hl_covW'))      = {'direct_covariance_latent_heat_flux_Webb_method'};
ln(ismember(the_fields,'hl_cov_sds'))   = {'direct_covariance_latent_heat_flux_SimondeSzoeke_method'};
ln(ismember(the_fields,'wu_cov'))       = {'wu_covariance'};
ln(ismember(the_fields,'wv_cov'))       = {'wv_ovariance'};
ln(ismember(the_fields,'wq_cov'))       = {'wq_covariance'};
ln(ismember(the_fields,'wT_cov'))       = {'wT_covariance_SimondeSzoeke_method'};
ln(ismember(the_fields,'wq_sds'))       = {'wq_covariance_SimondeSzoeke_method'};
ln(ismember(the_fields,'wT_sds'))       = {'wT_covariance_SimondeSzoeke_method'};
ln(ismember(the_fields,'upat_std'))     = {'platform_velocity_std_deviations_stern_to_bow_direction_u'};
ln(ismember(the_fields,'vplat_std'))    = {'platform_velocity_std_deviations_port_to_starboard_direction_v'};
ln(ismember(the_fields,'wplat_std'))    = {'platform_velocity_std_deviations_vertical_direction_w'};
ln(ismember(the_fields,'tilt'))         = {'tilt_angle_from_first_rotation'};
ln(ismember(the_fields,'azm'))          = {'azimuth_angle_from_first_rotation'};
ln(ismember(the_fields,'missing'))      = {'sonic_anemometer_missing points'};
ln(ismember(the_fields,'badPnts'))      = {'sonic_anemometer_bad_points'};
ln(ismember(the_fields,'uvar'))         = {'variance_stern_to_bow_from_sonic_anemometer_u'};
ln(ismember(the_fields,'vvar'))         = {'variance_port_to_starboard_from_sonic anemometer_v'};
ln(ismember(the_fields,'wvar'))         = {'variance_vertical_from_sonic_anemometer_w'};
ln(ismember(the_fields,'Tvar'))         = {'variance_Tv_from_sonic_anemometer'};
ln(ismember(the_fields,'Ubar'))         = {'streamwise_mean_flow'};
ln(ismember(the_fields,'Tson_noise'))   = {'noise_sonic_anemometer_Tv'};
ln(ismember(the_fields,'shed_std_s'))   = {'ship_standard_deviation_sine_heading'};
ln(ismember(the_fields,'ched_std_s'))   = {'ship_standard_deviation_cosine_heading'};
ln(ismember(the_fields,'Cu2a'))         = {'Cu_method a_structure_function_parameter'};
ln(ismember(the_fields,'Cu2b'))         = {'Cu_method_b_structure_function_parameter'};
ln(ismember(the_fields,'Cw2a'))         = {'Cw_method a_structure_function_parameter'};
ln(ismember(the_fields,'Cw2b'))         = {'Cw_method_b_structure_function_parameter'};
ln(ismember(the_fields,'Ct2a'))         = {'Ct_method a_structure_function_parameter'};
ln(ismember(the_fields,'Ct2b'))         = {'Ct_method_b_structure_function_parameter'};
ln(ismember(the_fields,'Cq2a'))         = {'Cq_method a_structure_function_parameter'};
ln(ismember(the_fields,'Cq2b'))         = {'Cq_method_b_structure_function_parameter'};
ln(ismember(the_fields,'cspd'))         = {'current_speed_ADCP'};
ln(ismember(the_fields,'cdir'))         = {'current_direction_ADCP'};
ln(ismember(the_fields,'wspd_sfc'))     = {'wind_speed_with_respect_to_moving_surface'};
ln(ismember(the_fields,'wdir_sfc'))     = {'wind_direction_with_respect_to_moving_surface'};
ln(ismember(the_fields,'rdir_std'))     = {'PSL_standard_deviation_relative_wind_direction'};
ln(ismember(the_fields,'rdir_std_s'))   = {'ship_standard_deviation_relative_wind_direction'};
ln(ismember(the_fields,'Pu'))           = {'power_spectrum_u'};
ln(ismember(the_fields,'Pv'))           = {'power_spectrum_v'};
ln(ismember(the_fields,'Pw'))           = {'power_spectrum_w'};
ln(ismember(the_fields,'Pt'))           = {'power_spectrum_T'};
ln(ismember(the_fields,'Pq'))           = {'power_spectrum_q'};
ln(ismember(the_fields,'Cuw'))          = {'cospectrum_uw'};
ln(ismember(the_fields,'Cvw'))          = {'cospectrum_vw'};
ln(ismember(the_fields,'Ctw'))          = {'cospectrum_Tw'};
ln(ismember(the_fields,'Cqw'))          = {'cospectrum_qw'};
ln(ismember(the_fields,'mu'))           = {'mu_decorr_q'};
ln(ismember(the_fields,'mu_q'))         = {'mu_decorr'};
ln(ismember(the_fields,'lag'))          = {'lag_times_sonic'};
ln(ismember(the_fields,'rUn'))          = {'PSL_sonic_relative_streamwise_wind'};
ln(ismember(the_fields,'rUw'))          = {'PSL_sonic_relative_crossstream_wind'};
ln(ismember(the_fields,'Ue'))           = {'PSL_eastward_true_wind'};
ln(ismember(the_fields,'Un'))           = {'PSL_northward_true_wind'};
ln(ismember(the_fields,'good_motion'))= {'flag_1_when_motion_conditions_good_for_covariance'};
ln(ismember(the_fields,'good_motion_licor'))= {'flag_1_when_motion_and_licor_conditions_good_for_covariance'};
ln(ismember(the_fields,'good_motion_id'))= {'flag_1_when_motion_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'good_motion_id_licor'))= {'flag_1_when_motion_and_licor_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'kk'))= {'when_motion_conditions_good_for_covariance'};
ln(ismember(the_fields,'hh'))= {'when_motion_and_licor_conditions_good_for_covariance'};
ln(ismember(the_fields,'jj'))= {'when_motion_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'mm'))= {'when_motion_and_licor_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'kk_1'))= {'when_motion_conditions_good_for_covariance'};
ln(ismember(the_fields,'hh_1'))= {'when_motion_and_licor_conditions_good_for_covariance'};
ln(ismember(the_fields,'jj_1'))= {'when_motion_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'mm1'))= {'when_motion_and_licor_conditions_good_for_inertial_dissipation'};
ln(ismember(the_fields,'tau_cov'))     = {'covariance_stress'};
ln(ismember(the_fields,'tau_ida'))     = {'inertial_dissipation_tau_method_a'};
ln(ismember(the_fields,'tau_idb'))     = {'inertial_dissipation_tau_method_b'};
ln(ismember(the_fields,'hs_ida'))      = {'inertial_dissipation_hs_method_a'};
ln(ismember(the_fields,'hs_idb'))      = {'inertial_dissipation_hs_method_b'};
ln(ismember(the_fields,'hl_ida'))      = {'inertial_dissipation_hl_method_a'};
ln(ismember(the_fields,'hl_idb'))      = {'inertial_dissipation_hl_method_b'};
ln(ismember(the_fields,'hb'))           = {'buoyancy_flux'};
ln(ismember(the_fields,'hb_son'))       = {'buoyancy_flux_sonic_anemometer'};
ln(ismember(the_fields,'lw_net'))       = {'net_longwave_flux'};
ln(ismember(the_fields,'lw_up'))        = {'upward_longwave_flux'};
ln(ismember(the_fields,'rhoa10'))       = {'air_density_adjusted_to_10 m'};
ln(ismember(the_fields,'sw_net'))       = {'net_shortwave_flux'};
ln(ismember(the_fields,'sw_up'))        = {'upward_shortwave_flux'};
ln(ismember(the_fields,'zo'))           = {'roughness_length'};
ln(ismember(the_fields,'dT_warm'))      = {'warm_layer_T_difference_top_to_base'};
ln(ismember(the_fields,'dz_warm'))      = {'warm_layer_thickness_or_depth'};
ln(ismember(the_fields,'dT_warm_at_input'))= {'warm_layer_T_difference_at_input_level'};
ln(ismember(the_fields,'Le'))           = {'latent_heat_of_vaporization'};
ln(ismember(the_fields,'Tref'))         = {'Tair_adjusted_to_reference_height'};
ln(ismember(the_fields,'U10N'))         = {'wind_speed_adjusted_to_neutral_at_10_m'};
ln(ismember(the_fields,'UN'))           = {'wind_speed_adjusted_to_neutral_at_measurement_height'};
ln(ismember(the_fields,'UrefN'))        = {'wind_speed_adjusted_to_neutral_at_reference_height'};
ln(ismember(the_fields,'Uref'))         = {'wind_speed_adjusted_to_reference_height'};
ln(ismember(the_fields,'rhref'))        = {'rh_adjusted_to_reference_height'};
ln(ismember(the_fields,'qref'))         = {'qair_adjusted_to_reference_height'};
ln(ismember(the_fields,'P10'))          = {'air_pressure_adjusted_to_10_m'};
ln(ismember(the_fields,'gust'))         = {'gustiness_factor'};
ln(ismember(the_fields,'hb_s'))           = {'ship_buoyancy_flux'};
ln(ismember(the_fields,'hb_son_s'))       = {'ship_buoyancy_flux_sonic_anemometer'};
ln(ismember(the_fields,'lw_net_s'))       = {'ship_net_longwave_flux'};
ln(ismember(the_fields,'lw_up_s'))        = {'ship_upward_longwave_flux'};
ln(ismember(the_fields,'rhoa10_s'))       = {'ship_air_density_adjusted_to_10 m'};
ln(ismember(the_fields,'sw_net_s'))       = {'ship_net_shortwave_flux'};
ln(ismember(the_fields,'sw_up_s'))        = {'ship_upward_shortwave_flux'};
ln(ismember(the_fields,'zo_s'))           = {'ship_roughness_length'};
ln(ismember(the_fields,'dT_warm_s'))      = {'ship_warm_layer_T_difference_top_to_base'};
ln(ismember(the_fields,'dz_warm_s'))      = {'ship_warm_layer_thickness_or_depth'};
ln(ismember(the_fields,'dT_warm_at_input_s'))= {'ship_warm_layer_T_difference_at_input_level'};
ln(ismember(the_fields,'Le_s'))           = {'ship_latent_heat_of_vaporization'};
ln(ismember(the_fields,'Tref_s'))         = {'ship_Tair_adjusted_to_reference_height'};
ln(ismember(the_fields,'U10N_s'))         = {'ship_wind_speed_adjusted_to_neutral_at_10_m'};
ln(ismember(the_fields,'UN_s'))           = {'ship_wind_speed_adjusted_to_neutral_at_measurement_height'};
ln(ismember(the_fields,'UrefN_s'))        = {'ship_wind_speed_adjusted_to_neutral_at_reference_height'};
ln(ismember(the_fields,'Uref_s'))         = {'ship_wind_speed_adjusted_to_reference_height'};
ln(ismember(the_fields,'rhref_s'))        = {'ship_rh_adjusted_to_reference_height'};
ln(ismember(the_fields,'qref_s'))         = {'ship_qair_adjusted_to_reference_height'};
ln(ismember(the_fields,'P10_s'))          = {'ship_air_pressure_adjusted_to_10_m'};
ln(ismember(the_fields,'gust_s'))         = {'ship_gustiness_factor'};
