function [u,v]=wind_components(ws,wd) % wind_components - calculates u(eastward) and v(northward) wind components % from wind speed & direction % NB wind direction is meteorology convention, i.e. direction wind is FROM % % [u,v]=wind_components(ws,wd) % % ws,wd - wind speed and direction (degrees from north) % u,v - east and north components % % Ian Brooks : January 23 1998 u=-1.*ws.*sin(wd.*pi/180); v=-1.*ws.*cos(wd.*pi/180);