% runs coare model for hourly average inputs
% assumes data from HiWinGS_2013_hrMeans_r2.txt is loaded in workspace

sigH = NaN;
zi = NaN;
cp = NaN;
zu = 16;
zu_ship = 15;
zt = 14;
zt_ship = 15;
zq = 14;
zq_ship = 15;

A = coare35vn(WspdR,zu,Tair,zt,RH,zq,Pair,Tsea,Solar,IR,Lat,zi,Pcp,cp,sigH);

% Output:  A=[usr tau hsb hlb hbb hsbb hlwebb tsr qsr zot zoq Cd Ch Ce  L zet dter dqer tkt Urf Trf Qrf RHrf UrfN Rnl Le rhoa UN U10 U10N Cdn_10 Chn_10 Cen_10 RF Qs Evap T10 Q10 RH10 us wbar zo];
%              1   2   3   4   5   6    7      8   9  10  11  12 13 14  15 16   17  18   19  20  21  22  23   24   25 26  27  28  29  30   31     32     33    34 35  36  37  38   39  40  41  42

A = [Yday,A];

dfl = '/Users/byronb/Raw_Data_Files.noindex/High_Wind_Flux_2013/Bulk Met Data Final/HiWinGS_2013_hrC35outputs_r2.txt';
save(dfl,'A','-ascii','-tabs');