vector_avg_wind.zip_matlab wind_vector_vector_avg_wind_wind_wind

function [ws_m,wd_m]=vector_avg_wind(a_time,doy,ws,wd,delta_a) % vector average for wind [doy ws wd] into a_time based on the time interval delta_a ; % output: averaged wind spd and wind dir; % east/north wind components of each data point Ve=sin(wd/180*pi).*ws; Vn=cos(wd/180*pi).*ws; aa=[doy Ve Vn]; bb=a_merge(a_time,aa,delta_a); Ve_avg=bb(:,2); Vn_avg=bb(:,3); len=length(a_time); for i=1:len ws_avg(i,1)=sqrt(Ve_avg(i)^2+Vn_avg(i)^2); wd_avg(i,1)=atan(Ve_avg(i)/Vn_avg(i))*180/pi; % if Vn_avg>=0 % wd_avg(i,1)=wd_avg(i,1)+180; % else % if Ve_avg(i)>=0 % wd_avg(i,1)=wd_avg(i,1)+360; % else % wd_avg(i,1)=wd_avg(i,1)+0; % end % end if Ve_avg(i)>=0 & Vn_avg(i)>=0 %0-90 deg wd_avg(i,1)=atan(Ve_avg(i)/Vn_avg(i))*180/pi; end if Ve_avg(i)>=0 & Vn_avg(i)<0 %90-180 deg wd_avg(i,1)=180+atan(Ve_avg(i)/Vn_avg(i))*180/pi; end if Ve_avg(i)==0 & Vn_avg(i)<0 %180 deg wd_avg(i,1)=180; end if Ve_avg(i)<0 & Vn_avg(i)<=0 %180-270 deg wd_avg(i,1)=atan(Ve_avg(i)/Vn_avg(i))*180/pi+180; end if Ve_avg(i)<0 & Vn_avg(i)>0 %270-360 deg wd_avg(i,1)=360+atan(Ve_avg(i)/Vn_avg(i))*180/pi; end end ws_m = ws_avg; wd_m = wd_avg;