disp('rdson1_08')   %('eval_epi_rhb_07')

%  DEW  5/15/2007 just do sonics
%changed by cwf to inlcde imet PIR  12/2004
U_ship=[];DIR_ship=[];GYRO_ship=[];TSEA_ship=[];TAIR_ship=[];Xaxis_ship=[];Xaxis_mean=[];RH_ship=[];
LW_ship=[];SOL_ship=[];LAT_ship=[];LON_ship=[];IR_ship=[];
Xaxis_sonic=[];T_sonic=[];U_sonic=[];V_sonic=[];
W_sonic=[];DIR_sonic=[];TAIR_mean=[];TSEA_mean=[];RH_mean=[];LW_mean=[];SOL_mean=[];QA_mean=[];

for hhh=6:6  %23;
   disp(['hour = ',num2str(hhh)])
   
   
   %%%%%%%%%%%%%%%%%%%%%%%%%
   %rdson2_08    % does spectra, but not working yet    eval_son_rhb_07
   
   
   imum=[]; imdm=[]; lrgm=[]; tsg=[]; imrh=[]; imta=[]; Tvais=[]; Tsea=[]; Rhvais=[]; Tslow=[]; U=[]; V=[]; w =[];irm=[];
   clear eval_ship_error eval_mean_error eval_son_error
end;
%tdk=273.15;
%jjx=min(length(Xaxis_ship),length(QA_mean));
%rlclr=(.52+.13/60*abs(LAT_ship(1:jjx)')+(.082-.03/60*abs(LAT_ship(1:jjx)')).*sqrt(QA_mean(1:jjx))).*(5.67e-8*(TAIR_mean(1:jjx)+tdk).^4); 

clear jdy;
