disp('eval_nau')

U_ship=[];DIR_ship=[];GYRO_ship=[];TSEA_ship=[];RH_ship=[];TAIR_ship=[];
TAIR_mean=[];TSEA_mean=[];RH_mean=[];T_sonic=[];U_sonic=[];V_sonic=[];W_sonic=[];
%Q_ophir=[];

for hhh=0:23;
   disp(['  '])
   disp(['hour = ',num2str(hhh)])
   disp(['  '])
   
   eval_ship;
   eval_mean;
   eval_son;
   
   if hhh==0 | hhh==12;
      disp('2')
      eval_ophir;
      disp('3')
   end
      
   U_ship=[U_ship imum];
   DIR_ship=[DIR_ship imdm];
   GYRO_ship=[GYRO_ship lrgm];
   TSEA_ship=[TSEA_ship tsg];
   RH_ship=[RH_ship imrh];
   TAIR_ship=[TAIR_ship imta];
   
   TAIR_mean=[TAIR_mean Tvais];
   TSEA_mean=[TSEA_mean Tsea];
   RH_mean=[RH_mean Rhvais];
   
   T_sonic=[T_sonic Tslow];
   U_sonic=[U_sonic U];
   V_sonic=[V_sonic V];
   W_sonic=[W_sonic w];
   
   clear imum imdm lrgm tsg imrh imta Tvais Tsea Rhvais Tslow U V w 
   
   %   Q_ophir=[Q_ophir abshum];
   
   dndx=ddd+(0:length(U_ship)-1)./length(U_ship);
   
end;

figure;plot(dndx,sqrt(U_sonic.^2+V_sonic.^2),dndx,U_ship);title('Wind Speed:  Blue=ETL, Green=Ship')
figure;plot(dndx,RH_mean,dndx,RH_ship);title('Relative Humidity:  Blue=ETL, Green=Ship')
figure;plot(dndx,TSEA_mean,dndx,TSEA_ship);title('Sea Temperature:  Blue=ETL, Green=Ship')
figure;plot(dndx,TAIR_mean,dndx,TAIR_ship);title('Air Temperature:  Blue=ETL, Green=Ship')

