function [wd]=wind_direction_deg(u,v, m,n)
wd=NaN*ones(m,n);

u_wd=find(v>0 & u >0);
alpha=atan(u(u_wd)./v(u_wd));
alpha=pi+alpha;
wd(u_wd)=alpha;

u_wd=find(v>0 & u <0);
alpha=-atan(v(u_wd)./u(u_wd));
alpha=pi/2+alpha;
wd(u_wd)=alpha;

u_wd=find(v<0 & u >0);
alpha=-atan(v(u_wd)./u(u_wd));
alpha=3*pi/2+alpha;
wd(u_wd)=alpha;

u_wd=find(v<0 & u <0);
alpha=atan(u(u_wd)./v(u_wd));
wd(u_wd)=alpha;

u_wd=find(v>0 & u ==0);
alpha=pi;
wd(u_wd)=alpha;

u_wd=find(v<0 & u ==0);
alpha=0;
wd(u_wd)=alpha;

u_wd=find(v==0 & u >0);
alpha=3*pi/2;
wd(u_wd)=alpha;

u_wd=find(v==0 & u <0);
alpha=pi/2;
wd(u_wd)=alpha;

wd=rad2deg(wd);
%wd=wd';