%flx=load('C:\data\cwf\Dropbox\DATA\ATOMIC\RBH\flux\ATOMIC_2020_10min_flux_met_sea_data_01_23_jd023.txt');
flx=load('C:\data\cwf\Dropbox\DATA\ATOMIC\RBH\flux\ATOMIC_2020_10min_flux_met_sea_data_01_25_jd025.txt');


    Cd=flx(:,1);%: [144×1 double]		drag coefficient = turbulent exchange coefficient for wind stress adjusted for 10-m height and neutral atmospheric stability, unitless
    Ce=flx(:,2);%: [144×1 double]		Dalton number = turbulent exchange coefficient for evaporation adjusted for 10-m height and neutral atmospheric stability, unitless
    Ch=flx(:,3);%: [144×1 double]		Stanton number = turbulent exchange coefficient for sensible heat assuming adjusted for 10-m height and neutral atmospheric stability, unitless
    Evap=flx(:,4);%: [144×1 double]		evaporation rate, mm/hr
    L=flx(:,5);%: [144×1 double] 	Monin Obukhov length scale, m
    SST=flx(:,6);%: [144×1 double]		calculated skin-level sea surface temperature from Tsnk - calculated cool skin (dter), C
       Stsg=flx(:,7);%: [144×1 double]		ship thermosalinograph SBE 45 salinity, psu
        T10=flx(:,8);%: [144×1 double]		air T adjusted for 10-m height and neutral atmospheric stability, C
         Ta=flx(:,9);%: [144×1 double]		air T at height zt, C
       Tsnk=flx(:,10);%: [144×1 double]		sea snake T at depth zsnk, C
      Tstar=flx(:,11);%: [144×1 double]		temperature scaling parameter, K
       Ttsg=flx(:,12);%: [144×1 double]		ship thermosalinograph SBE 38 T at depth ztsg_ship, C
        U10=flx(:,13);%: [144×1 double]		wind speed adjusted to 10 m, m/s
         Ue=flx(:,14);%: [144×1 double]		eastward wind speed, m/s
         Un=flx(:,15);%: [144×1 double]		northward wind speed, m/s
        cog=flx(:,16);%: [144×1 double]		course over ground, deg
    dT_skin=flx(:,17);%: [144×1 double]		cool-skin temperature depression, C
        day=flx(:,18);%: [144×1 double]		day
    dz_skin=flx(:,19);%: [144×1 double]		cool-skin vertical thickness, m
        hed=flx(:,20);%: [144×1 double]		heading, deg
         hl=flx(:,21);%: [144×1 double]		latent heat flux, W/m2
       hnet=flx(:,22);%: [144×1 double]		net heat flux, W/ms (hnet = lwnet + swnet + hl + hs + hr)
       hour=flx(:,23);%: [144×1 double]		hour
         hr=flx(:,24);%: [144×1 double]		rain heat flux, W/m2
         hs=flx(:,25);%: [144×1 double]		sensible heat flux, W/m2
         jd=flx(:,26);%: [144×1 double]		decimal julian date
        lat=flx(:,27);%: [144×1 double]		latitude, deg
        lon=flx(:,28);%: [144×1 double]		longitude, deg
      lw_dn=flx(:,29);%: [144×1 double]		downwelling long wave radiation, W/m2
     lw_net=flx(:,30);%: [144×1 double]		net long wave radiation, W/m2
      lw_up=flx(:,31);%: [144×1 double]		upwelling long wave radiation, W/m2
     minute=flx(:,32);%: [144×1 double]		minute
      month=flx(:,33);%: [144×1 double]		month
      prate=flx(:,34);%: [144×1 double]		instantaneous precipitation rate from optical rain gauge, mm/hour
        q10=flx(:,35);%: [144×1 double]		air specific humidity adjusted for 10-m height and neutral atmospheric stability, g/kg
         qa=flx(:,36);%: [144×1 double]		air specific humidity at height zt, zq, zp, g/kg
         qs=flx(:,37);%: [144×1 double]		saturation specific humidity at surface, g/kg
      qstar=flx(:,38);%: [144×1 double]		specific humidity scaling parameter, g/kg
       rdir=flx(:,39);%: [144×1 double]		relative wind direction, deg
         rh=flx(:,40);%: [144×1 double]		relative humidity at height zq_ship, %
       rh10=flx(:,41);%: [144×1 double]		relative humidity adjusted for 10-m height and neutral atmospheric stability, %
       rhoa=flx(:,42);%: [144×1 double]		density of air at height zt_ship, kg/m3
       rspd=flx(:,43);%: [144×1 double]		relative wind speed, m/s
        slp=flx(:,44);%: [144×1 double]		sea level pressure, mb
        sog=flx(:,45);%: [144×1 double]		speed over ground, m/s
      sw_dn=flx(:,46);%: [144×1 double]		downwelling short wave radiation, W/m2
     sw_net=flx(:,47);%: [144×1 double]		net short wave radiation, W/m2
      sw_up=flx(:,48);%: [144×1 double]		upwelling short wave radiation, W/m2U10
        tau=flx(:,49);%: [144×1 double]		wind stress, N/m2
      ustar=flx(:,50);%: [144×1 double]		friction velocity that includes gustiness, m/s
       wdir=flx(:,51);%: [144×1 double]		true wind direction, deg
       wspd=flx(:,52);%: [144×1 double]		true wind speed, deg
       year=flx(:,53);%: [144×1 double]		year
        zet=flx(:,54);%: [144×1 double]		Monin-Obukhov stability parameter zu/L, unitless
        
        figure;plot(jd,hs,jd,hl,jd,lw_net,jd,sw_net);xlabel('Yearday (2020)');ylabel('Heat flux (W/m^2)');legend('Hs','Hl','Rln','Rsn');grid;title('ATOMIC Ronald H Brown');
        figure;
        subplot(2,1,1);plot(jd,Tsnk,jd,Ttsg,jd,T10);ylabel('Temperature (C)');legend('Tsnk','Tsg','Ta10');grid;title('ATOMIC Ronald H Brown');
        subplot(2,1,2);plot(jd,U10);xlabel('Yearday (2020)');ylabel('U10 (m/s)');grid;title('ATOMIC Ronald H Brown');
        
