disp('eval_epi_sjohnson_05')

U_ship=[];DIR_ship=[];GYRO_ship=[];TSEA_ship=[];TAIR_ship=[];Xaxis_ship=[];Xaxis_mean=[];RH_ship=[];
LW_ship=[];SOL_ship=[];LAT_ship=[];LON_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=0:23;
   disp(['hour = ',num2str(hhh)])
   %eval_scs_altair_04;
   %U_ship=[U_ship 
   %    imum];
   %Xaxis_ship=[Xaxis_ship 
   %    gpstime];
   %DIR_ship=[DIR_ship 
   %    imdm];
   %GYRO_ship=[GYRO_ship 
   %    lrgm];
   %TSEA_ship=[TSEA_ship 
   %    tsgm];
   %RH_ship=[RH_ship 
   %    rhm];
   %TAIR_ship=[TAIR_ship 
   %    tam];
   %SOL_ship=[SOL_ship 
   %    solarm];
   LAT_ship=gpslat;
   LON_ship=gpslon;
   
   eval_pc_day_sjohnson_05;
   TAIR_mean=[TAIR_mean Tvais];
   TSEA_mean=[TSEA_mean Tsea];
   RH_mean=[RH_mean Rhvais];
   LW_mean=[LW_mean pir];
   SOL_mean=[SOL_mean psp];
   QA_mean=[QA_mean qvais];
   Xaxis_mean=[Xaxis_mean ddd+(hhh+(0:length(Rhvais)-1)./length(Rhvais))/24.];

   eval_son_sjohnson_05;
   T_sonic=[T_sonic Tslow];
   U_sonic=[U_sonic U];
   Xaxis_sonic=[Xaxis_sonic ddd+(hhh+(0:length(U)-1)./length(U))/24.];
   V_sonic=[V_sonic V];
   W_sonic=[W_sonic w];
   DIR_sonic=[DIR_sonic dir];
      
   if hhh==1 | hhh==12;
      plot_f=1;
      eval_ophir_sjohnson_05;
   end;

   imum=[]; imdm=[]; lrgm=[]; tsg=[]; imrh=[]; imta=[]; Tvais=[]; Tsea=[]; Rhvais=[]; Tslow=[]; U=[]; V=[]; w =[];
   clear eval_ship_error eval_mean_error eval_son_error
end;
tdk=273.15;
jjx=min(length(LAT_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;
jdy=Xaxis_mean';qa=qmean(1)*ones(length(Xaxis_mean),1);
latd=LAT_ship';
lond=LON_ship';
stdt=ddd;endt=stdt+1;
iv=ones(length(jdy),1)*2.5;
rico05_Rsclear_1;

Xaxis_ship=Xaxis_ship./24+ddd;

figure;plot(Xaxis_sonic,sqrt(U_sonic.^2+V_sonic.^2),'b');xlabel('Time of day');ylabel('Relative Wind Speed (m/s)');title('Relative Wind Speed ETL (blue), True Wind Speed SCS (green)');
figure;plot(Xaxis_sonic,DIR_sonic,'.b');xlabel('Time of day');ylabel('Relative Wind Direction (Deg)');title('Relative Wind Dir ETL (blue), True Wind Dir Ship (green)');
figure;plot(Xaxis_mean,RH_mean,'b');ylabel('Relative Humidity in %');xlabel('Time of day');title('Rel Hum (%):  Blue=ETL, Green dots=Ship');
figure;plot(Xaxis_mean,TSEA_mean,'b');xlabel('Time of day');ylabel('Sea Temperature (C)');title('Sea Temperature:  Blue=ETL, Green dots=Ship');
figure;plot(Xaxis_mean,TAIR_mean,'b');ylabel('Air Temperature (C)');xlabel('Time of day');title('Air Temp:  Blue=ETL, Green=Ship');
figure;plot(Xaxis_mean,LW_mean,'b',Xaxis_mean,rlclr,'r');ylabel('IR Flux (W/m^2)');xlabel('Time of day');title('IR Flux:  Blue=ETL, Red=Clearsky');
figure;plot(Xaxis_mean,Rscl,'g',Xaxis_mean,SOL_mean,'b');ylabel('Solar Flux (W/m^2)');xlabel('Time of day');title('Solar Flux:  Blue=ETL, Grn=Clear sky,  Red=Ship');

%figure;plot(Xaxis_sonic,sqrt(U_sonic.^2+V_sonic.^2),'b');xlabel('Time of day');ylabel('Relative Wind Speed (m/s)');title('Relative Wind Speed ETL (blue), True Wind Speed SCS (green)');
%figure;plot(Xaxis_sonic,DIR_sonic,'.b');xlabel('Time of day');ylabel('Relative Wind Direction (Deg)');title('Relative Wind Dir ETL (blue), True Wind Dir Ship (green)');
%figure;plot(Xaxis_mean,RH_mean,'b');ylabel('Relative Humidity in %');xlabel('Time of day');title('Rel Hum (%):  Blue=ETL, Green dots=Ship');
%figure;plot(Xaxis_mean,TSEA_mean,'b');xlabel('Time of day');ylabel('Sea Temperature (C)');title('Sea Temperature:  Blue=ETL, Green dots=Ship');
%figure;plot(Xaxis_mean,TAIR_mean,'b');ylabel('Air Temperature (C)');xlabel('Time of day');title('Air Temp:  Blue=ETL, Green=Ship');
%figure;plot(Xaxis_mean,LW_mean,'b',Xaxis_mean,rlclr,'r');ylabel('IR Flux (W/m^2)');xlabel('Time of day');title('IR Flux:  Blue=ETL, Red=Clearsky');
%figure;plot(Xaxis_mean+11.5/24,Rscl,'g',Xaxis_mean,SOL_mean,'.b');ylabel('Solar Flux (W/m^2)');xlabel('Time of day');title('Solar Flux:  Blue=ETL, Grn=Clear sky,  Red=Ship');
%clear Xaxis_ship Xaxis_mean Xaxis_sonic
